Question

用matlab在(x,y)=［-8,8］×［-8,8］区域中绘制曲面图形并比较线性、立方及样条插值的结果。

图形为z=［sin根号(x的平方+y的平方）］÷［根号(x的平方+y的平方）］，可以把步骤详细的写下来吗？谢谢

Answer 1

x=-8:8;
y=x;
[X,Y]=meshgrid(x,y);
Z=sin(sqrt(X.^2+Y.^2))./sqrt(X.^2+Y.^2+(X.^2+Y.^2==0)*eps);
figure(1)
mesh(X,Y,Z)
title('粗糙图')

x1=-8:0.5:8;
y1=x1;
[X1,Y1]=meshgrid(x1,y1);

figure(2)
Z1=interp2(X,Y,Z,X1,Y1,'linear');
mesh(X1,Y1,Z1)
title('线性插值精细图')

figure(3)
Z1=interp2(X,Y,Z,X1,Y1,'cubic');
mesh(X1,Y1,Z1)
title('立方插值精细图')

figure(4)
Z1=interp2(X,Y,Z,X1,Y1,'spline');
mesh(X1,Y1,Z1)
title('样条插值精细图')