matlab数学实验课

for j=2:N-1 if imag(z(k,j)~=0) z1(k,j)=0; end
if all(imag(z([k-1,k-1],[j-1,j1])))~=0
z1(k,j)=NaN; end
end end surf(x,y,z1),hold if u==1
z2=-z1; surf(x,y,z2); axis([-abs(a),abs(a),abs(b),abs(b),-abs(c),abs(c)]); end xlabel('x'),ylabel('y'),zlabel('z') title(['[a,b,c,d,N]=[',int2str(a),',',int 2str(b),',',int2str(c),',',int2str(d),',',i nt2str(N),']'])
空间解析几何 MATLAB实验
小组成员： 李姿周斯阳 熊向南刘志宇 化铮文涛榕
MATLAB实验
实验目的 1、学会用Matlab绘制简单的空间曲线. 2、学会用Matlab绘制简单的空间曲面.
实验内容 1、编写MATLAB程序绘制二次曲面图形 2、编写MATLAB程序绘制旋转曲面图形 3、编写MATLAB程序绘制空间曲线图形
x cos(t)
ysin(t)Fra bibliotekz t
编写M文件xuanzhuanqumian.m，如下：
t=0:pi/50:10*pi;
plot3(sin(t),cos(t),t)
xlabel('x'),ylabel('y')
1. *sec(u). 表示正号 2. x=a *cos(v);
y=b *sin(v); z=c *tan(v);
ADD YOUR TITLE HERE
79% 9%
42% 24%
ADD YOUR KEY TEXT HERE ADD YOUR KEY TEXT HERE ADD YOUR KEY TEXT HERE ADD YOUR KEY TEXT HERE ADD YOUR KEY TEXT HERE
单击此处添加文字信息 92% A good book may be among the best of friends. It is the
same today that it always was, and it will never change. A good book may be among the best of friends. It is the same today that it always was, and it will never change. A good book may be among the best of friends. It is the same today that it always was, and it will never change.
z=2*sin(v);
x2 y2 z2
绘制旋转单叶曲面
1
994
即代入： x=3*sec(u).*cos(v); y=3*sec(u).*tan(u); z=2*sin(v);
x2 y2 z2
绘制旋转双叶曲面
1
944
即代入： x=3*sec(u).*cos(v); y=2*tan(u); z=2*sin(v);
二次曲面图形
解：设二次曲面的方程为
x2 a2
y2 b2
z2 c2
d
编写M文件kongjianqumian.m，如下：
1 2 3
椭圆抛物面 椭球面 双叶双曲面
a=input('a=');b=input('b=');c=inpu t('c=');d=input('d=');N=input('N='); xgrid=linspace(abs(a),abs(a),N);ygrid=linspace(abs(b),abs(b),N); [x,y]=meshgrid(xgrid,ygrid); z=c*sqrt(d-y.*y/b/bx.*x/a/a);u=1;z1=real(z); for k=2:N-1
绘制抛物柱面 y2 2x
即代入：
x=3*sec(u).*cos(v);
y2
y=3*sec(u).*tan(u);
z=2*sin(v);
旋转曲面图形
解：设旋转曲面的方程为
x2
y2
z2
a2 b2 c2 1
1 旋转单叶双曲面
编写M文件xuanzhuanqumian.m，如下：
2
旋转双叶双曲面
u=-pi/3:pi/30;pi/3;
v=0:pi/30:2*pi;
[u,v]=meshgrid(u,v);
x=a*cos(v);
y=b*tan(u);
z=c*sin(v); mesh(x,y,z)
1. *sec(u). 表示正号 2. x=a *cos(v);
y=b *sin(v); z=c *tan(v);
空间曲线图形
1 螺旋线
解：设空间曲线的方程为
抛物柱面图形
解：设抛物柱面的方程为
x2 a2
y2 b2
z2 c2
d
编写M文件paowuzhumian.m，如下：
1 抛物柱面
x=-2*pi:0.1:2*pi; y=a*x.^2; n=size(x,2); for i=-2:0.1:2
z=i*ones(1,n); plot3(x,y,z) grid on hold on end xlabel('x'),ylabel('y')
绘制螺旋线
x cos(t)
y
sin(t)
z t
即代入： x=3*sec(u).*cos(v); y=2*tan(u); z=2*sin(v);
小组总结
遇到的困难：
收获的内容：
thanks
ADD YOUR TITLE HERE
绘制椭圆抛物面 y2 2x
即代入：
x=3*sec(u).*cos(v);
y2
y=3*sec(u).*tan(u);
z=2*sin(v);
绘制椭球面 y2 2x
即代入：
x=3*sec(u).*cos(v);
y2
y=3*sec(u).*tan(u);
z=2*sin(v);
绘制双叶双曲面 y2 2x
即代入：
x=3*sec(u).*cos(v);
y2
y=3*sec(u).*tan(u);