MATLAB作业2

1. （1）编写M 函数文件，实现分段绘制曲面，绘制每个曲面为一个子函数，

并为图形窗口添加标题，设置标题的颜色和字号为红色，13号，设置文字颜色为红色。

???????-<+≤+<->+=+-------15457.0117575.015457.0),(5.175.375.065.175.375.022222

2y x e y x e y x e y x z x x y x y x x y

(2)使用句柄函数实现函数的调用

解：

function y = draw( )

%draw 分段绘制曲面

x1=-2:0.1:2;

x2=-2:0.1:2;

[xx1,xx2]=meshgrid(x1,x2); %产生栅格点坐标

[n1,m1]=size(xx1); %求维数

z=zeros(n1,m1); %生成n1*m1的全0矩阵

h_draw1=str2func('draw1') %创建子函数draw1的句柄

h_draw2=str2func('draw2') %创建子函数draw1的句柄

h_draw3=str2func('draw3') %创建子函数draw1的句柄

for n=1:n1

for m=1:m1

if xx1(n,m)+xx2(n,m)>1

z(n,m)=feval(h_draw1,xx1(n,m),xx2(n,m));

elseif (xx1(n,m)+xx2(n,m)>-1)&(xx1(n,m)+xx2(n,m)<=1)

z(n,m)=feval(h_draw2,xx1(n,m),xx2(n,m));

else

z(n,m)=feval(h_draw3,xx1(n,m),xx2(n,m));

end

end

end

surf(xx1,xx2,z)

function z=draw1(x,y) %以下为子函数

z=0.5457*exp(-0.75*y^2-3.75*x^2-1.5*x);

function z=draw2(x,y)

z=0.7575*exp(-y^2-6*x^2);

function z=draw3(x,y)

z=0.5457*exp(-0.75*y^2-3.75*x^2+1.5*x);

结果为：

为sphere

.解：

function yuan

[m,n,p]=sphere(30);

t=abs(p);

surf(m,n,p,t);

hold on;

mesh(2*m,2*n,2*p),colormap(hot);

hold off;

hidden off;

grid off

axis equal