《MATLAB课后习题及练习.docx》由会员分享,可在线阅读,更多相关《MATLAB课后习题及练习.docx(10页珍藏版)》请在第一文库网上搜索。
1、7.6i.分别绘制下列函数图形:r=3(1cos)(极坐标)t=O:.O1:2*pi;figurepo1ar(t,abs(3*(1-cos(t);运行结果:(2)y(t)=1.25e-O.25t+cos(3t)x=0:pi/180:2*pi;y1=1.25*exp(-0.25*x);y2=cos(3*x);p1ot(x,y1+y2)运行结果:23 .绘制函数y(t)=1-2etsint(0t8)的图形,且在x轴上标注“Time”,y轴上标注Amp1itude”,图形的标题为DecayingOsci11atingExponentia1”0x=0:0.01:8;y=1-2*exp(-1*x).*s
2、in(x);figurePot(,y)x1abe1(Time)jyabe1(Amp1itude,);tit1e(DecayingOsci11atingExponentia1);运行结果:45 .在同一图中绘制下列两条曲线(XW0,25内):(1) y1(t)=2.6e(-0.5x)cos(0.6x)+0.8;y2(t)=16cos(3x)+sin(x)要求用不同的颜色和线型分别表示y1(t)和y2,并给图形加注解。x=0:0.1:25;y1=2.6*exp(-0.5*x).*cos(0.6*x)+0.8;y2=1.6*cos(3*x)+sin(x);figurep1ot(x,y1,r-X,x2
3、b-.);Iegend(y1=2.6*exp(-0.5*x).*cos(0.6*x)+0.87y2=1.6*cos(3*x)+sin(x),)4.在一个图形窗口下绘制两个子图,分别显示下列曲线:(1) y=si2xcos3x;(2)y=0.4X要求给X轴、y轴加标注,每个子图加标题。X=IinSPaCe(C),2*pi,100);y1=sin(2*x).*cos(3*x);y2=0.4*x;SUbP1Ot(1,2,1);PIOt(X,y1);x1abe1(xyabe1(,y)jtit1e(y=sin(2*x)*cos(3*x);SUbP1ot(1,2,2);P1Ot(X,y2);x1abe1(
4、,x,by1abe1(y)jtit1e(,0.4*x,);5.画出下列二元函数z(x,y)的图形:Z(x,y)=1/(x+1)A2+(y+1)A2+1-1(x-1)2+(y-1)2+1(-3=x=3,-3=y=d=1;(2) a=5j,b=4c=3d=1a=-2:0.1:2;x,y=meshgrid(a);a=5;b=4;c=3;d=1;z1=c.*sqrt(d-(x.A2)./(a.A2)-(y.A2)./(b.A2);SUbP1Ot(1,2,1),SUrf(X,y,zI)a=5j;z2=c.*sqrt(d-(x.A2)./(a.A2)-(y.A2)./(b.A2);SUbP1Ot(1,2,
5、2),SUrf(X,y,z2)实验六T1AB的数据可视化一、实验目的掌握MAT1AB二维、三维图形绘制,掌握图形属性的设置和图形修饰;掌握图像文件的读取和显示。二、实验内容(1)二维图形绘制。(2)三维曲线和三维曲面绘制。(3)图像文件的读取和显示。三、实验步骤(3) 维图形绘制(1)二维图形绘制主要使用函数p1ot。c1eara11;x=1inspace(0,2*pi,100);y1=sin(x);p1ot(x,y1)ho1don%保持原有的图形y2=cos(x);p1ot(x,y2)注:ho1don用于保持图形窗口中原有的图形,ho1doff解除保持。(2)函数PIOt的参数也可以是矩阵。
6、c1osea11%关闭所有图形窗口X=IinSPaCe(0,2*pi,100);y1=sin(x);y2=cos(x);A=y1;y2,;%把矩阵转置B=x;x,;p1ot(B,A)(3)选用绘图线形和颜色。c1osea11%关闭所有图形窗口p1ot(xzy1zg+,x,y2,r:)gridon%添加网格线(4)添加文字标注。tit,正弦曲线和余弦曲线,)yabe(幅度)XIabeI(时间,)1egend(,sin(x)z,cos(x),)gtet(,1eftarrowsin)%可用鼠标选择标注的位置,%1eftarrow产生左箭头,为转义符修改坐标轴范围。axisequa1axisnorma
7、1ais(0pi01.5)运行结果:正弦曲线和余弦曲线*n()COft(X)t,1w*+十+4sn+(6)子图和特殊图形绘制。subp1ot(2z2,1)t1=0:0.1:3;y1=exp(-t1);bar(t1zy1);subp1ot(2z2z2)t2=0:0.2:2*pi;y2=sin(t2);stem(t2,y2);subp1ot(2z2,3)t3=0:0.1:3;y3=t3.2+1;stairs(t3zy3);subp1ot(2,2,4)t4=0:.01:2*pi;y4=abs(cos(2*t4);po1ar(t4,y4);运行结果:练习1:写出图A2的绘制方法。提示:按照以下的步骤进
8、行产生曲线的数据(共有3组数据);(2)选择合适的线形、标记、颜色(正弦曲线为红色,余弦曲线为紫色卜(3)添加图例及文字说明信息;添加坐标轴说明与图标题。解:XI=IinSPaCe(0,4*pi,100);x2=0:0.5:16;y1=sin(x);y2=cos(x);y3=0;figureP1Ot(X1yI,r-,x1,y2,m-,x2,y3,k:.)yabe1(,gx=cos(z);y=sin(z);p1ot3(xzyzz);练习2:利用子图函数,绘制以上的空间螺旋线的俯视图、左侧视图和前视图。z=0:0.1:6*pi;x=cos(z);y=sin(z);SUbP1Ot(2,2,1),P1
9、ot3(x,y,z);tit1e(三维视图)subp1ot(2,2,2),p1ot3(x,y,z),view(90,0)tit1e(恻视图)subp1ot(2,2z3)Pot3(x,y,z),view(0,0)tit1e(正视图)subp1ot(2,2,4),p1ot3(x,z),view(0,90)tit1ed俯视图)O1-二-1-05005三维视图005运行结果:确困总物面冏面困(2)三维曲面图的绘制:MAT1AB绘制网线图和网面图的函数分别是mesh()和surf(),其具体操作步骤是:用函数meshgrid()生成平面网格点矩阵X,Y;由XM计算函数数值矩阵Z;用mesh()绘制网线图
10、,用surf()绘制网面图。绘制椭圆抛物面:c1eara11,c1osea11;x=-4:0.2:4;y=x;X,Y=meshgrid(x,y);Z=X.29+Y.29;mesh(XZ);tit1e(椭圆抛物面网线图,)figure(2)surf(X,Z);tit1e(,椭圆抛物面网面图,)椭因财物面网经图绘制阔边帽面:c1eara11,c1osea11;x=-7.5:0.5:7.5;y=x;XzY=meshgrid(x,y);R=sqrt(X.A2+Y.A2)+eps;%避开零点,以免零做除数Z=sin(R).R;mesh(XZ);tit1e(,阔边帽面网线图,)figure(2)surf(
11、X,Z);tit1e(悯边帽面网面图,)练习3:考虑以下问题:设Z=(XC2)*exp(-1*(X.A2+YC2);,求定义域x=-2,2,y=-2,2内的Z值(网格取0.D,并用网面图表示出来。解:x=-2:0.1:2;y=-2:0.1:2;X,Y=meshgrid(x,y);Z=(X.A2).*exp(-1*(X.A2+Y.A2);SUrf(XM;运行结果:3.图像文件的读取和显示x=imread(cameraman.tif);imshow(x)y=255-doub1e(x);y=uint8(y);figureimshow(y)imwrite(y,reverse.tif,)%首先读取图像文件%对图像进行反色处理%将图像数据保存为文件