MATLAB编程设计与应用课程设计.docVIP

《MATLAB编程设计与应用》课程设计 所在学院： 理学院 专业名称： 应用物理学 所 在 班 级： 物理09-1 学 生 姓 名： 学 生 学 号： 2011 年 6 月 一，程序的结构和流程控制： 用Matlab程序的一般语句（循环语句for…end,条件循环语句while…end和条件转向语句if…else…end）分别求1～100的奇数和。 程序一 程序二： 程序三 sum=0; for x=1:1:100; if mod(x,2) sum=sum+x; else continue; end end sum clear x=1; sum=0; while x101 sum=sum+x; x=x+2; end sum clear sum=0; for x=1:2:99; sum=sum+x; end sum 二，2D绘制、函数图形的绘制： 1，用plot语句绘制函数图形，熟悉分区绘图命令（subplot），建立图形命令（figure），以及用title命令和一些指定的线性、颜色和标记对图形作注释的功能。 x=[-pi:pi/20:pi]; y1=sin(x); y2=cos(x); figure(1) subplot(2,2,1),plot(x,y1,+) grid on,title(sin(x)) subplot(2,2,2),plot(x,y2,k) grid on,title(cos(x)) subplot(2,2,3),plot(x,y1,b*,x,y2,g-.) grid on,title(sin(x) and cos(x)) subplot(2,2,4),plot(x,y1,-,x,y1,ro) grid on,title(sin(x)) 2，用subplot命令将画图区域分为三份，分别用ploy，comet和fplot命令做y=sin（3x+2）的图像。并用图形标题设置命令（title），X轴标签设置命令（xlabel），y轴标签设置命令（ylabel）对图形作注释。 程序： x=-2*pi:pi/20:2*pi; y=sin(3*x+2); subplot(3,1,1) plot(x,y) title(y=sin(3x+2)) xlabel(x [-2*pi,2*pi]) ylabel(y=sin(3x+2)) subplot(3,1,2) comet(x,y,g) title(y=sin(3x+2)) xlabel(x [-2*pi,2*pi]) ylabel(y=sin(3x+2)) subplot(3,1,3) fplot(sin(3*x+2),[-2*pi,2*pi]) title(y=sin(3x+2)) xlabel(x [-2*pi,2*pi]) ylabel(y=sin(3x+2)) 用行结果： ,3，练习使用ezplot命令绘制参数曲线图形： 参数方程及程序 用行结果 参数方程： 程序： subplot(2,1,1) ezplot((cos(x))^3,(sin(x)^3),[-4*pi,4*pi]) subplot(2,1,2) ezplot((x-sin(x)),(1-cos(x)),[-4*pi,4*pi] 参数方程： 程序： subplot(2,1,1) ezplot(x*(cos(x))^3,x*(sin(x)^3),[-pi,pi]) subplot(2,1,2) ezplot(x*(x-sin(x)),x*(1-cos(x)),[-pi,pi]) 4，用polar命令绘制极坐标图形，并熟悉gtext命令，用鼠标将文本点击到图形相应的位置上。 在同一坐标区域内绘制三叶玫瑰线和四叶玫瑰线。 程序： subplot(2,2,1) ezpolar(sin(2*x),[0,2*pi]) gtext(??ò??μ1???r=sin(2x)) subplot(2,2,2) ezpolar(cos(2*x),[0,2*pi]) gtext(??ò??μ1???r=cos(2x)) subplot(2,2,3) ezpolar(sin(3*x),[0,2*pi]) gtext(èyò??μ1???r=sin(3x)) subplot(2,2,4) ezpolar(cos(3*x),[0,2*pi]) gtext(èyò??μ1???