实验作业2MATLAB作图

1、实验作业2MATLAB作图1、用 plot,fplot 绘制函数 y=cos（tan（二 x）的图形。解：代码如下：x=li nspace(-0.5,0.5,30);y=cos(ta n( pi*x);plot(x,y)结果：H Ftgure 1I =回 加fplot( 'cos(ta n(pi*x)'，卜0.5,0.5)结果:BQ Figure 1I d | 回S3File Edit View Insert Tools Desktop Window Help生h乱Q酸動運乙TG&日囲回注意事项：fplot('cos(tan(pi*x)'，卜0.5,0.

2、5)不能写成fplot( 'cos(tan(pi*x)', (-0.5Q.5 )否则是错误程序。2. 用ezplot绘制函数 exy-sin(x+y)=0在-3, 3上的图形。解：代码如下：ezplot('exp(x*y)-si n(x+y)=0',-3,3)3. 用 matlab 以动画的方式绘制出摆线 x=a(t-sint) y=a(1-cost) ( a自己赋值)的渐屈线clear;clc;close;a=1;syms tx=a*(t-s in( t);y=a*(1-cos(t);ezplot(x,y,0,2*pi), grid on ;hold on;

3、dy=diff(y)/diff(x);dyy=diff(dy)/diff(x);xx=x-(1+dyA2)*dy/dyy;% 渐屈线的坐标 yy=y+(1+dyA2)/dyy;M=50;t=0;xxx=subs(xx);yyy=subs(yy);H1=plot(xxx,yyy,'r');hold on; grid on; axis(0,7,-2.5,2.5); x1=subs(x);y1=subs(y);H2=plot(x1,xxx,y1,yyy,'k-');H3=plot(x1,y1,'ko');H4=plot(xxx,yyy,'ro&

4、#39;);tt=li nspace(0,2*pi,M);for i=1:Mpause(0.2);t=tt(1:i);xxx=subs(xx);yyy=subs(yy);x1=subs(x);y1=subs(y);set(H1,'xdata',xxx,'ydata',yyy); set(H2,'xdata',x1(i),xxx(i),'ydata',y1(i),yyy(i); set(H3,'xdata',x1(i),'ydata',y1(i);set(H4,'xdata',xxx(

5、i),'ydata',yyy(i);end绘制函数x=a(t-sint )在tqo,2眄上的图形。 y =a(1 _cost)程序：a=in put('please in put a=');t=li nspace(0,2*pi,100);x=a*(sin(t);y=a*(t-cos(t);plot(x,y);please input a=104.5. 在同一平面中的两个窗口分别画出心形线和马鞍面。要求： 1)在图形上加格栅、图例和标注；2)定制坐标；3)以不同角度观察马鞍面。解：心型线代码如下：clear;clc;i=-pi:0.1:pi;x=2.*(sin(i

6、)-sin(2*i)./2);y=2.*(cos(i)-cos(i)42);plot(x,y);axis(-3 3 -4.2 1);%制定图形坐标title(' 心形线 ');% 当前图顶端加图例xlabel('x');%当前图x轴加图例ylabel('y');%当前图y轴加图例grid on;% 加格栅到当前图上心型线运行结果显示如下图：1心形线图3.1 （心型线）马鞍线代码如下:clear;clc;a=10;b=10;x=-1.5:0.1:1.5;y=-1.5:0.1:1.5;X,Y=meshgrid(x,y);%x行 y 列的矩阵Z=X.A

7、2/a-Y.A2/b;figure (2);%新建窗口subplot(221);% 分割当前窗口，作图显示在该分割块上surf (X,Y,Z);title(' 马鞍面');xlabel('x');ylabel('y');subplot(222);surf (X,Y,Z);title(' 马鞍面');xlabel('x');ylabel('y');view(20,30);% 改变视角到(20 , 30)subplot(223);surf (X,Y,Z);title(' 马鞍面');xlabel('x');ylabel('y');view(50,20);% 改变视角到(50,20 )subplot(224);surf (X,Y,Z);title(' 马鞍面');xlabel('x');ylabel('y');view(60,10);% 改变视角到(60,10 )马鞍线运行结果显示如下图：图3.2 （马鞍面）7.以不同的视角观察球面x2+y2+z2=r2和圆柱面x2+y2=rx所围区域.解：代码如下：clear;clc;r=2;x0,y0,z0=sphere(50);