三维图形绘制MATLAB程序

三维图形绘制

clc,clear,close all

a=-0.98;b=0.98;c=-1;d=1;n=10;

x=linspace(a,b,n); y=linspace(c,d,n);

[X,Y]=meshgrid(x,y);

plot(X,Y,'+')

clear,clf,

a=-1;b=1;c=-15;d=15;n=20;eps1=0.01;

x=linspace(a,b,n);y=linspace(c,d,n);

[X,Y]=meshgrid(x,y);

for i=1:n %计算函数值z ，并作定义域裁剪for j=1:n

if (1-X(i,j))

z(i,j)=NaN; %作定义域裁剪，定义域以外的函数值为NaN else

z(i,j)=1000*sqrt(1-X(i,j))^-1.*log(X(i,j)-Y(i,j));

end

end

end

zz=-20*ones(1,n);plot3(x,x,zz),grid off,hold on %画定义域的边界线

mesh(X,Y,z) %绘图，读者可用meshz, surf, meshc在此替换之view([-56.5 38]);

xlabel('x'),ylabel('y'),zlabel('z'), box on %把三维图形封闭在箱体里