matlab典型例子(吐血整理)

一、在同一坐标系内，分别用红色实线和蓝色点划线绘制曲线y1(t)=e0.1sin(0.3t)和y2(t)=2e0.5cos(0.6t)，其中t[0,3]，并用星号标记出两条曲线的交叉点。请写出实现上述功能的M文件程序代码，并给出每条语句注释说明。t=0:0.01:3%t取值[0,3]步长为0.01y1=exp(0.1)*sin(0.3*pi*t)%写出公式y1(t)=e0.1sin(0.3t)plot(t,y1,'r-')%显示该图形holdon%继续显示y2=2*exp(-0.5)*cos(0.6*pi*t)%写出公式y2(t)=2e0.5cos(0.6t)plot(t,y2,'b-.')%显示该图形r0=abs(y2-exp(0.1)*sin(0.3*pi*t))=0.015%设定差值的绝对值=0.015yy=r0.*y1;xx=r0.*t;plot(xx(r0~=0),yy(r0~=0),'r*')%用*标记交叉点结果如下：00.511.522.53-1.5-1-0.500.511.5二、求解多项式3x48x2+x10除以多项式x2+x+3；计算方程3x48x2+x10=0的解。请编写实现上述运算的M程序代码，并给出每条语句注释说明。a=[30-81-10]%3x48x2+x10的系数矩阵b=[113]%x2+x+3的系数矩阵[q,r]=deconv(a,b)%求根和余。q为根，r为余结果如下：q=3-3-14r=0002432a=[30-81-10]%多项式3x4-8x2+x-10=0的矩阵b=roots(a)%求根结果如下：b=-1.93191.85810.0369+0.9629i0.0369-0.9629i三、编写M文件，计算下题：取自变量x的区间为[6,6]，数据间隔为0.01，有函数y=sin(x)+0.5,z=f(x,y)=3sin(x)y/10，绘制三维曲线图z。然后新建一个图形窗口，利用图形窗口分割方法分别画出该三维曲线图在4个不同视角(37.5,30)、(16,60)、(100,0)、(8,12)下的三维图形，并对每个子图形标注坐标轴名称和图标题(标题内容为该图形的视角值)。要求给出程序代码，以及每条语句的注释说明。(其中视角坐标值(az,el)中az--方位角；el--俯视角)x=-6:0.01:6%变量区间和间隔设置为0.01y=sin(x)+0.5z=3*sin(x)-y/10%列出公式[az,el]=viewplot3(x,y,z)xlabel('X');ylabel('Y');zlabel('Z')%显示坐标轴图形如下：-10-50510-0.500.511.5-3-2-10123Xaz=332.5el=30YZx=-6:0.01:6%'变量区间[-6,6]的间隔为0.01y=sin(x)+0.5z=3*sin(x)-y/10%列出公式subplot(2,2,1),plot3(x,y,z),view(-37.5,30);xlabel('X');ylabel('Y');zlabel('Z');title('az=-37.5el=30')%第一个窗口显示的视角为(-37.5,30)subplot(2,2,2),plot3(x,y,z),view(-16,60);xlabel('X');ylabel('Y');zlabel('Z');title('az=-16el=60')%第二个窗口显示的视角为('az=-16el=60')subplot(2,2,3),plot3(x,y,z),view(-100,0);xlabel('X');ylabel('Y');zlabel('Z');title('az=-100el=0')%第三个窗口显示的视角为('az=-100el=0')subplot(2,2,4),plot3(x,y,z),view(-8,-12);xlabel('X');ylabel('Y');zlabel('Z');title('az=-8el=-12')%第四个窗口显示的视角为('az=-8el=-12')结果如下：-10010-202-505Xaz=-37.5el=30YZ-10010-202-505Xaz=-16el=60YZ-10010-1012-4-2024Yaz=-100el=0XZ-10010-202-505YXaz=-8el=-12Z四、A、B、C三个城市，上半年每个月的国民生产总值（单位：亿元）如右表。请编写M文件分别计算：(1)三个城市这6个月的平均国民生产总值。(2)这三个城市的月国民生产总值最大的城市及其最大值。(3)写出生成如右下图的二维累加直方图的M文件程序代码。要求给出每条语句的注释说明。各城市每月国民生产总值原始数据城市1月2月3月4月5月6月A728010590120113B609011410211095C5140708596110a=[728010590120113;609011410211095;5140708596110]b=mean(a,2)%平均每行c=max(b)%求出最大运行结果：b=96.6667%对应A95.1667%对应B75.3333%对应Cc=96.6667%A城最大month=[123456]%月数status=[607251;908040;11410570;1029085;11012096;95113110]%每月的国民生产总值数据bar(month,status,'stack')%画出直方图legend('B市','A市','C市')%增添字符xlabel('x');ylabel('∑y')%显示坐标单击showplottools修改为题目所要求的颜色结果如下：123456050100150200250300350x∑yB市A市C市