實(shí)驗(yàn)作業(yè)2MATLAB作圖

1、.實(shí)驗(yàn)作業(yè)2 MATLAB作圖1、用plot,fplot 繪制函數(shù)y=cos(tan(x)的圖形。解：代碼如下:x=linspace(-0.5,0.5,30); y=cos(tan(pi*x); plot(x,y)結(jié)果：fplot(cos(tan(pi*x),-0.5,0.5)結(jié)果：注意事項(xiàng)： fplot(cos(tan(pi*x),-0.5,0.5)不能寫(xiě)成fplot(cos(tan(pi*x),（-0.5,0.5）)否則是錯(cuò)誤程序。2. 用ezplot繪制函數(shù) exy-sin(x+y)=0 在-3, 3上的圖形。解：代碼如下：ezplot(exp(x*y)-sin(x+y)=0,-3,3)

2、3. 用matlab以動(dòng)畫(huà)的方式繪制出擺線 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+dy2)*dy/dyy;%漸屈線的坐標(biāo)yy=y+(1+dy2)/dyy;M=50;t=0;xxx=subs(xx);yyy=subs(yy);H1=plot(xxx,yyy,r);hold on;grid on;

3、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,

4、x1(i),ydata,y1(i); set(H4,xdata,xxx(i),ydata,yyy(i);end 繪制函數(shù)在上的圖形。程序：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=104.5. 在同一平面中的兩個(gè)窗口分別畫(huà)出心形線和馬鞍面。要求：1)在圖形上加格柵、圖例和標(biāo)注； 2)定制坐標(biāo)； 3)以不同角度觀察馬鞍面。解：心型線代碼如下：clear;clc;i=-pi:0.1:pi;x=2.*(sin(i)-sin(2*i)./2);y=

5、2.*(cos(i)-cos(i).2);plot(x,y);axis(-3 3 -4.2 1);%制定圖形坐標(biāo)title(心形線);%當(dāng)前圖頂端加圖例xlabel(x);%當(dāng)前圖x軸加圖例ylabel(y);%當(dāng)前圖y軸加圖例grid on;%加格柵到當(dāng)前圖上心型線運(yùn)行結(jié)果顯示如下圖： 圖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);%分割當(dāng)前窗口，作圖顯示在該分割塊上su

6、rf (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）馬鞍線運(yùn)行結(jié)果顯示如下圖：圖3.2（馬鞍面） 7.以不同的視角觀察球面x2+y2+z2=r2 和圓柱面 x2+y2=rx所圍區(qū)域.解：代碼如下：clear;clc;r=2;x0,y0,z0=sphere(50);x=r*x0;y=r*y0;z=r*z0;sur