qq_37486958的博客

画图主要是用极坐标，我们使用画复数的函数来画。 clc;%清空上次的 %draw_Bessel%画这个东西 a=4;%圆形半径% %定义角向和径向的步长和范围 fai=[-4:0.1:4]; r=[0:0.02:a]; %将fai和r组合起来并生成复平面 [R,T] = meshgrid(fai,r); Z = R.*exp(1i*T); %和画具体图有关的 m=1;%阶数% k_11=4.6...

画的是《重力学与固体潮》那本书的斜阶梯的重力异常 function [gravity]=draw_gravity_dixing(h1,h2,alpha,xigma,x0,z,direction) %h1=100;%上表面 %h2=200;%下表面 %x0=0;%移动原点位置，向右为正 %xigma=3*10^3;%剩余密度 pi=3.1415926535 %alpha=alpha*pi;%倾角 ...

h=xlabel('x'); set(h,'Fontsize',14); h=ylabel('df(x)/dx'); set(h,'Fontsize',14); h=title('精确解和二阶差分对比'); set(h,'Fontsize',14); h=legend('精确的一阶偏导','二阶差分得到的偏导'); set(h,'Fontsize',14); 这样就行了 ...

读写： https://seisman.github.io/SAC_Docs_zh/commands/write/ 保存文件： https://seisman.github.io/SAC_Docs_zh/graphics/save-image/