我的课程设计

长 春 大 学

┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ 装 ┊ ┊ ┊ ┊ ┊ 订 ┊ ┊ ┊ ┊ ┊ 线 ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊

学 生 实 习 报 告

[实习名称] 数学软件课程实习

[实习时间] 2011年6月20日—2011年7月1日 [机器型号] 清华同方

[系统环境] Microsoft Windows 2000 [开发工具] matlab [所用机时] 80学时 [实习地点] 计算中心机房 [实习项目][1]绘制下列曲线。 ①

100y?1?x2 ② y?1e2??x22

2?x?t?22③ x?y?1 ④ ?3 ??y?5t代码x=0:pi/100:2*pi; subplot(2,2,1); y=100./(1+x.*x); plot(x,y);

subplot(2,2,2) x=0:pi/100:2*pi;

y=exp(-x.^2/2)*(1/2*pi); plot(x,y);

subplot(2,2,3);

ezplot('x^2+y^2=1',[-1,1]) subplot(2,2,4); t=0:pi/100:2*pi; x=t.*t; y=t.*t*5; plot(x,y); 结果

共 28 页 第 1 页

长 春 大 学

┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ 装 ┊ ┊ ┊ ┊ ┊ 订 ┊ ┊ ┊ ┊ ┊ 线 ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊

学 生 实 习 报 告

[实习项目][2] 绘制下列极坐标图。 ① ??5cos??4 ② ??③ ??12?

5??7 ④ ???2 cos?3代码theta=0:0.01:2*pi;

subplot(2,2,1);

rho=5*(cos(theta))+4; polar(theta,rho,'k'); subplot(2,2,2);

rho=12./(theta.^0.5); polar(theta,rho,'k'); subplot(2,2,3);

rho=5./cos(theta)-7; polar(theta,rho,'k'); subplot(2,2,4);

共 28 页 第 2 页

长 春 大 学

┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ 装 ┊ ┊ ┊ ┊ ┊ 订 ┊ ┊ ┊ ┊ ┊ 线 ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊

学 生 实 习 报 告

rho=pi/3*(theta.*theta); polar(theta,rho,'k'); 结果

[实习项目] 【3】绘制下列三维图形。

?x?cost?x?(1?cosu)cosv??①?y?sint ② ?y?(1?cosu)sinv

?z?t?z?sinu??③ z?5 ④ 半径为10的球面 代码subplot(2,2,1); t=0:pi/100:2*pi; x=cos(t); y=sin(t); z=t;

plot3(x,y,z,'p'); subplot(2,2,2); u=0:pi/100:2*pi; v=0:pi/100:2*pi;

共 28 页 第 3 页

长 春 大 学

┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ 装 ┊ ┊ ┊ ┊ ┊ 订 ┊ ┊ ┊ ┊ ┊ 线 ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊

学 生 实 习 报 告

x=(1+cos(u)).*cos(v); y=(1+cos(u)).*sin(v); z=sin(u);

plot3(x,y,z,'p'); subplot(2,2,3);

[x,y]=meshgrid(-10:0.1:10); z=5*ones(size(x)); plot3(x,y,z); subplot(2,2,4);

[x,y,z]=sphere(50); meshz(x,y,z);

surf(10*x,10*y,10*z); 结果

[实习项目] 【4】在同一图形窗口采用子图形式分别绘制正方形、圆、三角形和六边形。 代码

subplot(2,2,1);

共 28 页 第 4 页