matlab实验一

开课学院、实验室：实验时间：年月日

图4.1（上上）图4.2（上）图5.1.1（下）图5.1.2(下下)

5．作出下列曲面的3维图形，

1）)sin(22y x z +π=；

第一次程序:x=-7.5:0.5:7.5;y=x;[X,Y]=meshgrid(x,y);Z=sin(pi*sqrt(X.^2+Y .^2));mesh(X,Y ,Z);图如图5.1.1 第二次程序x=-2:0.01:2;y=x;[X,Y]=meshgrid(x,y);Z=sin(pi*sqrt(X.^2+Y .^2));mesh(X,Y,Z);图如图5.1.2所示 第三次程序x=-7.5:0.01:7.5;y=x;[X,Y]=meshgrid(x,y);Z=sin(pi*sqrt(X.^2+Y.^2));mesh(X,Y,Z); 如图5.1.3所示

图5.1.3（上）图5.2（下）

2）环面：⎪⎩

⎪

⎨⎧=+=+=,

sin ,sin )cos 1(,

cos )cos 1(u z v u y v u x )2,0()2,0(ππ∈∈v u 。

程序：

u=linspace(0,2*pi,30);v=u;[U,V]=meshgrid(u,v);X=(1+cos(U)).*cos(V);Y=(1+cos(U)).*sin(V);Z=sin(U);mesh(X,Y,Z), 图如图5.2所示

3）分别作出单位球面在参数为两种不同取值范围的图形,注意坐标轴的单位长度要相等。提示：附加命

令rotate3d 可实现3维图形旋转。

a) cos sin ,sin sin ,cos ,x u v y u v z v =⎧⎪=⎨⎪=⎩ (0,1.6)(0,)u v ππ∈∈;b) cos sin ,

sin sin ,cos ,

x u v y u v z v =⎧⎪

=⎨⎪=⎩

(0,2)(0.5,)u v πππ∈∈

程序：u1=linspace(0,1.6*pi,30);v1=linspace(0,pi,30);u2=linspace(0,2*pi,30);v2=linspace(0.5*pi,pi,30);...

[U1,V1]=meshgrid(u1,v1);[U2,V2]=meshgrid(u2,v2);X1=cos(U1).*sin(V1);Y1=sin(U1).*sin(V1);Z1=cos(V1);. X2=cos(U2).*sin(V2);Y2=sin(U2).*sin(V2);Z2=cos(V2);

subplot(1,2,1),mesh(X1,Y1,Z1),title('a');subplot(1,2,2),mesh(X2,Y2,Z2),title('b');图如图5.3所示

图5.3（上）

图5.4（上）图5.5（下）

4）z =y 2 绕z 轴的旋转面图形

程序：x=-5:0.01:5;y=x;[x,y]=meshgrid(x,y);z=x.*x+y.*y;mesh(x,y,z),图如图5.4 5) y = -2

z ,0

x=eps:0.01:5-eps;y=-x;[x,y]=meshgrid(x,y);z=sqrt(-y);mesh(x,y,z)图如图5.5 6．建立一个命令M-文件：求所有的“水仙花数”，所谓“水仙花数”是指一个三位数，其各位数字的立

方和等于该数本身。例如，153是一个水仙花数，因为153=13+53+33

。 程序：for a=1:9; for b=0:9; for c=0:9;

if a*100+b*10+c==a^3+b^3+c^3; d=a*100+b*10+c, end end end end

运行结果：sxh d =153 d = 370 d =371 d =407 7．编写函数M-文件sq.m ：用迭代法求a =

x 的值。求平方根的迭代公式为

)a (211n

n n x x x +=

+ 迭代的终止条件为前后两次求出的x 的差的绝对值小于10-5

。

程序: function f=sq(a); x=1;for k=1:10000 x=0.5*(x+a/x); end