程序我是这样写的,不知道怎么设定sqrt(X.^2+Y.^2)的范围在【R1,R2】就直接没设定,所以做出来的图很奇怪,和原图不一样。
syms a b R1 R2;
a=60;
R1=13;
R2=100;
x=linspace(-100,100,100);
y=linspace(-100,100,100);
[X,Y] = meshgrid(x,y);
b=sqrt(X.^2+Y.^2);
Z =(R1-b).*(b-R2)./a;
subplot(1,2,1);
surf(X,Y,real(Z));
subplot(1,2,2);
imagesc(real(Z));
我做出来的图是这样的:
别人的图是这样的: