实验作业2----MATLAB作图

实验作业2MATLAB作图1、用plot,fplot绘制函数y=cos(tan(x))的图形。解：代码如下:x=linspace(-0.5,0.5,30);y=cos(tan(pi*x));plot(x,y)结果：fplot('[cos(tan(pi*x))]',[-0.5,0.5])结果：注意事项：fplot('[cos(tan(pi*x))]',[-0.5,0.5])不能写成fplot('[cos(tan(pi*x))]',（-0.5,0.5）)否则是错误程序。2.用ezplot绘制函数exy-sin(x+y)=0在[-3,3]上的图形。解：代码如下：ezplot('exp(x*y)-sin(x+y)=0',[-3,3])3.用matlab以动画的方式绘制出摆线x=a(t-sint)y=a(1-cost)（a自己赋值）的渐屈线clear;clc;close;a=1;symstx=a*(t-sin(t));y=a*(1-cos(t));ezplot(x,y,[0,2*pi]),gridon;holdon;dy=diff(y)/diff(x);dyy=diff(dy)/diff(x);xx=x-(1+dy^2)*dy/dyy;%渐屈线的坐标yy=y+(1+dy^2)/dyy;M=50;t=0;xxx=subs(xx);yyy=subs(yy);H1=plot(xxx,yyy,'r');holdon;gridon;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');tt=linspace(0,2*pi,M);fori=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(i),'ydata',yyy(i));end绘制函数)cos1()sin(tayttax在]2,0[t上的图形。程序：a=input('pleaseinputa=');t=linspace(0,2*pi,100);x=a*(sin(t));y=a*(t-cos(t));plot(x,y);pleaseinputa=104.5.在同一平面中的两个窗口分别画出心形线和马鞍面。要求：1)在图形上加格栅、图例和标注；2)定制坐标；3)以不同角度观察马鞍面。解：心型线代码如下：clear;clc;i=-pi:0.1:pi;x=2.*(sin(i)-sin(2*i)./2);y=2.*(cos(i)-cos(i).^2);plot(x,y);axis([-33-4.21]);%制定图形坐标title('心形线');%当前图顶端加图例xlabel('x');%当前图x轴加图例ylabel('y');%当前图y轴加图例gridon;%加格栅到当前图上心型线运行结果显示如下图：图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.^2/a-Y.^2/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);x=r*x0;y=r*y0;z=r*z0;surf(x,y,z);holdon;f=@(x,y)x.^2+y.^2-r*x;ezsurf(f,[-1.5,2.5],[-2,1.5]);axisequal;view(30,20);运行结果显示如下图：图3.3