MATLAB课后习题及练习.pdf
7.6习 题1.分别绘制下列函数图形:r=3(l-cos 0)(极坐标)t=0:.01:2*pi;figure polar(t,abs(3*(l-cos(t);运行结果：(2)y(t)=1.25e-0.25t+cos(3t)x=0:pi/180:2*pi;yl=1.25*exp(-0.25*x);y2=cos(3*x);plot(x,yl+y2)运行结果：2.绘制函数y(t)=l-2e tsint(0 WtW8)的图形，且 在 x 轴上标注Time”,y 轴上标注 Amplitude”,图形的标题为 aDecaying Oscillating Exponentialw x=0:0,01:8;y=l-2*exp(-l*x).*sin(x);figure plot(x,y)xlabelCTimeJjylabelCAmplitude);title(Decaying Oscillating Exponential);运行结果：3.在同一图中绘制下列两条曲线(x 0,25内)：(1)y l(t)=2.6e(-0.5x)cos(0.6x)+0.8；(2)y 2(t)=1.6cos(3x)+sin(x)要求用不同的颜色和线型分别表示y i(t)和 y 2(t),并给图形加注解。x=0:0.1:25;yl=2.6*exp(-0.5*x).*cos(0.6*x)+0.8;y2=1.6*cos(3*x)+sin(x);figure plot(x,yLr-X,x,y2,b-)；Iegend(yl=2.6*exp(-0.5*x).*cos(0.6*x)+0.87y2=1.6*cos(3*x)+sin(x)4.在一个图形窗口下绘制两个子图，分别显示下列曲线:(1)y=sin2xcos3x；(2)y=0.4 x要求给x 轴、y 轴加标注，每个子图加标题。x=linspace(0z2*pi,100);yl=sin(2*x).*cos(3*x);y2=0.4*x;subplot。,2,l);plot(x,yl);xlabel(x);ylabel(y)；title(,y=sin(2*x)*cos(3*x);subplot(l,2,2);plot(x,y2);xlabel(,x,);ylabel(y);title(0.4*x,);5.画出下列二元函数z(x,y)的图形：Z(xzy)=l/(x+l)A2+(y+l)A2+l-l/(x-l)A2+(y-l)A2+l(-3=x=3,-3=y=title(恻视图)subplot(2,2,3),plot3(x,y/z),view(0,0)title(正视图)subplot(2,2,4),plot3(x,y,z),view(0,90)title。俯视图)运行结果：三维视图侧视图(2)三维曲面图的绘制：MATLAB绘制网线图和网面图的函数分别是mesh()和 surf(),其具体操作步骤是：用 函 数 meshgrid()生成平面网格点矩阵 X,Y；由 X H 计算函数数值矩阵Z；用 mesh()绘制网线图，用 surf()绘制网面图。绘制椭圆抛物面：clear all,close all;x=-4:0.2:4;y=x;X/Y=meshgrid(x,y);Z=X.A2/9+Y.A2/9;mesh(XzYzZ);title。椭圆抛物面网线图,)figure(2)surf(XXZ);title(椭圆抛物面网面图,)绘制阔边帽面：clear all,close all;x=-7.5:0.5:7.5;y=x;X,Y=meshgrid(xzy);R=sqrt(X.A2+Y.A2)+eps;%避开零点，以免零做除数Z=sin(R)./R;mesh(XXZ);title(悯边帽面网线图)figure(2)surf(XXZ)；title(阔边帽面网面图)练习3:考虑以下问题:设 z=(X.A2).*exp(-l*(X，2+Y.A2);,求定义域 x=-2,2,y=-2,2内的z 值(网格取0.1),并用网面图表示出来。解：x=-2:0.1:2;y=-2:0.1:2;X,Y=meshgrid(x,y);Z=(X.A2).*exp(-l*(X.A2+Y.A2);surf(X,YZ);运行结果：3.图像文件的读取和显示 x=imread(cameraman.tif);imshow(x)y=255-double(x);y=uint8(y);figure imshow(y)imwrite(y,reverse.tif)%首先读取图像文件%对图像进行反色处理%将图像数据保存为文件
