实验作业2MATLAB作图.docVIP

实验作业2 MATLAB作图 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;syms tx=a*(t-sin(t));y=a*(1-cos(t));ezplot(x,y,[0,2*pi]),grid on;hold on;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);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);tt=linspace(0,2*pi,M);for i=1:M pause(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 绘制函数在上的图形。 程序：a=input(please input a=); t=linspace(0,2*pi,100); x=a*(sin(t)); y=a*(t-cos(t)); plot(x,y); please input a=10 4. 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([-3 3 -4.2 1]);%制定图形坐标 title(心形线);%当前图顶端加图例 xlabel(x);%当前图x轴加图例 ylabel(y);%当前图y轴加图例 grid on;%加格栅到当前图上 心型线运行结果显示如下图： 图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);