MATLAB 实习报告.doc

1.计算矩阵 与 之和程序：cleara=5 3 5;3 7 4;7 9 8;b=2 4 2;6 7 9;8 3 6;c=b+a运行结果：c =7 7 -79 14 1315 12 142. 求的共轭转置。程序：x=4+8i 3+5i 2-7i 1+4i 7-5i ;3+2i 7-6i 9+4i 3-9i 4+4i;x运行结果：ans = 4.0000- 8.0000i 3.0000 - 2.0000i3.0000-5.0000i 7.0000+6.0000i2.0000+7.0000i 9.0000-4.0000i1.0000-4.0000i 3.0000+9.0000i3. 对于，如果，求解X。程序：clearA=4 9 2;7 6 4;3 5 7;B=37;26;28;X=AB运行结果：X =-0.51184.04271.33184.已知：，分别计算a的数组平方和矩阵平方，并观察其结果。程序：cleara=1 2 3;4 5 6;7 8 9;a.2a2运行结果：ans =1 4 916 25 3649 64 81ans =30 36 4266 81 96102 126 1505. 计算多项式乘法(x2+2x+2)(x2+5x+4)程序：a=1 2 2;b=1 5 4;h=conv(a,b)运行结果：h =1 7 16 18 86. 对下式进行部分分式展开：程序：a=1 3 4 2 7 2;b=3 2 5 4 6;r,s,k=residue(b,a)运行结果：r =1.1274 + 1.1513i1.1274 - 1.1513i-0.0232 - 0.0722i-0.0232 + 0.0722i0.7916 s =-1.7680 + 1.2673i-1.7680 - 1.2673i0.4176 + 1.1130i0.4176 - 1.1130i-0.2991 k = 7.计算多项式的微分和积分。程序：clearf=sym(4*x4-12*x3-14*x2+5*x+9)diff(f)clearf=sym(4*x4-12*x3-14*x2+5*x+9);int(f)运行结果：ans =16*x3-36*x2-28*x+5ans =4/5*x5-3*x4-14/3*x3+5/2*x2+9*x8. 因式分解： 程序：syms x;f=x4-5*x3+5*x2+5*x-6;factor(f)运行结果： ans = (x-1)*(x-2)*(x-3)*(x+1)9. ，用符号微分求df/dx。程序：syms a x;f=a,x2,1/x;exp(a*x),log(x),sin(x);df=diff(f) 运行结果：df = 0, 2*x, -1/x2 a*exp(a*x), 1/x, cos(x)10. 求代数方程组关于x,y的解。程序：S=solve(a*x2+b*y+c=0,b*x+c=0,x,y);disp(S.x=),disp(S.x)disp(S.y=),disp(S.y)运行结果：S.x=-c/bS.y=-c*(a*c+b2)/b311. 有一组测量数据满足，t的变化范围为0-10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线。程序:t=0:0.5:10;y1=exp(-0.1*t);y2=exp(-0.2*t);y3=exp(-0.5*t);plot(t,y1,-.b,t,y2,:*r,t,y3,-g)运行结果：12. 在5.1题结果图中添加标题和图例框。程序：t=0:0.5:10;y1=exp(-0.1*t);y2=exp(-0.2*t);y3=exp(-0.5*t);plot(t,y1,-.b,t,y2,:*r,t,y3,-g)legend(a=0.1,a=0.2,a=0.5)title(y=exp(-a*x)运行结果：13. x= 66 49 71 56 38，绘制饼图，并将第五个切块分离出来。程序：x=66 49 71 56 38;y=0 0 0 0 1;pie(x,y)运行结果：14. ，当x和y的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口中绘制出三维线图、网线图、表面图。程序：x,y=meshgrid(-2:0.2:2);z=x.*exp(-x.2-y.2);mesh(x,y,z)subplot(2,2,1),plot3(x,y,z)title(plot3(x,y,z)subplot(2,2,2),mesh(x,y,z)title(mesh(x,y,z)subplot(2,2,3),surf(x,y,z)title(surf(x,y,z)运行结果：15. 有传递函数如下的控制系统，用Simulink建立系统模型，并对系统的阶跃响应进行仿真。程序： 运行结果：16. 建立一个简单模型，用信号发生器产生一个幅度为2V、频率为0.5Hz的正弦波，并叠加一个0.1V的噪声信号，将叠加后的信号显示在示波器上并传送到工作空间。程序：运行结果：17. 已知五个数据点: (1,5.5), (2,43.1), (3,128), (4,290.7), (5,498.4)(1)用三次曲线拟合上述数据点(2)在同一图形中绘出数据点和拟合曲线程序：clearx=1 2 3 4 5;y=5.5 43.1 128 290.7 498.4;n=2;p=polyfit(x,y,n);x1=linspace(0,6,20);y1=polyval(p,x1);plot(x,y,o-b,x1,y1,r-.,linewidth,2);运行结果：18. 解线性方程组a1x+b1y+c1z=d1a2x+b2y+c2z=d2a3x+b3y+d3z=d3，其中ai, bi, ci, di为常数, x, y, z为变量程序：g1=sym(a1*x+b1*y+c1*z=d1); g2=sym(a2*x+b2*y+c2*z=d2);g3=sym(a3*x+b3*y+c3*z=d3); f=solve(g1,g2,g3)x,y,z=solve(a1*x+b1*y+c1*z=d1,a2*x+b2*y+c2*z=d2,a3*x+b3*y+c3*z=d3)运行结果：x =-(-b1*c2*d3+b1*d2*c3-c1*b3*d2+c1*b2*d3-d1*b2*c3+d1*b3*c2)/(a1*b2*c3-a1*b3*c2-a2*b1*c3+b3*a2*c1-b2*a3*c1+a3*b1*c2)y =(-a1*c2*d3+a1*d2*c3-d2*a3*c1+c2*a3*d1+a2*c1*d3-a2*d1*c3)/(a1*b2*c3-a1*b3*c2-a2*b1*c3+b3*a2*c1-b2*a3*c1+a3*b1*c2)z =(b3*a2*d1-a1*b3*d2+a1*b2*d3-b2*a3*d1-a2*b1*d3+a3*b1*d2)/(a1*b2*c3-a1*b3*c2-a2*b1*c3+b3*a2*c1-b2*a3*c1+a3*b1*c2)19. 绘制下列各种函数图形(1)绘制下列极坐标图形 r=3(1-cosq) r=2(1+cos q) r=2(1+sin q) r=cos3 q(2)求函数z的三维图形.定义区间与z的函数表达式如下:-5 x 5 , -5 y 5程序：(1)x=linspace(0,2*pi);r=3*(1-cos(x);polar(x,r);x=linspace(0,2*pi);r=2*(1+cos(x);polar(x,r);x=linspace(0,2*pi);r=2*(1+sin(x)polar(x,r);x=linspace(0,2*pi);r=cos(3*x);polar(x,r);运行结果：(2)x=linspace(-5,5);y=linspace(-5,5); x,y=meshgrid(x,y); z=1(x+1).2+(y+1).2+1-1(x-1).2+(y-1).2+1; mesh(x,y,z);运行结果：20. 用MATLAB命令绘制下列数学函数的图形(1) 绘制如下函数图形y(t)=1-2exp(-t)sin(t) (0 t 8)且在x轴写上“Time”,y轴写上“Amplitude”标号,图形标题为“Decaying-oscillating Exponential”(2) 绘制如下图形y(t)=5exp(-0.2t)cos(0.9t-30O)+0.8exp(-2t) (0 t 30)(3)在0 t10区间内绘制下面图形y(t)=1.23cos(2.83t+240O)+0.625x(t)=0.625t(4)在0 t20区间内,且在同一图中绘制下述函数图形y1(t)=2.62exp(-0.25t)cos(2.22t+174O)+0.6 y2(t)=2.62exp(-0.25t)+0.6 y3(t)=0.6程序：(1) t=0:0.2:8;y=1-2.*exp(-t).*sin(t);plot(t,y,linewidth,3);title(Decaying-oscillating Exponential);xlabel(Time)ylabel(Amplitude)(2) t=0:0.2:30;y=5*exp(-0.2*t).*cos(0.9*t-300)+0.8.*exp(-2*t);plot(t,y,linewidth,3);(3) t=0:0.2:10;y=1.23*cos(2.83*t+2400)+0.625;x=0.625*t;plot(t,y,-or,t,x,b*);(4) t=