matlab指导书.docx

1、matlab指导书实验一 数值运算实验1 实验目的（1）学习MATLAB语言的基本矩阵运算（2）学习MATLAB语言的点运算（3）学习复杂运算2 实验内容在下面的实验操作中，认真记录每项操作的作用和目的。（1） 基本矩阵运算1） 创建数值矩阵a=1 2 3;4 5 6;7 8 9; 观察aa(3,2)a(:,1)键入t=0:10;u=0:0.1:10;观察矩阵变量t、u的值。键入a(:,3)=2;3;4;a观察矩阵a的变化。键入b=1 1+2i;3+4i 3;观察复数矩阵。2） 创建特殊矩阵键入a=ones(3,3);b=zeros(2,2);c=eye(4);magic(4);观察特殊矩阵。

2、3） 练习矩阵运算a=0 1 0;0 0 1;-6 -11 -6;b=1 2;3 4;5 6;c=1 1 0;0 1 1;作矩阵乘运算v1=c*av2=a*bv3=c*a*bv4=b*cv5=c*b矩阵乘方运算a2a(1/2)矩阵加减运算a1=a+b*ca2=c*b-a(1:2,1:2)a3=a(1:2,2:3)+c*b矩阵右除ar=c/a矩阵左除al=ab4） 练习矩阵特征运算完成以下特征运算inv(a),tril(a),rank(a)（2） MATLAB语言的点运算1） 练习点乘与点除a1=1 2;3 4;a2=0.2*a1;观察a1 a2a1.*a2 a1./a22） 由点运算完成标量函

3、数运算与作图正余弦函数的点运算t=0:2*pi/180:2*pi;y1=sin(t);y2=cos(t);y=y1.*y2;plot(t,y y1 y2);复变函数的点运算。w=0.1:0.1:2;g1=(1+0.5*w*i)/(1-0.5*w*i);g1g2=(1+0.5*w*i)./(1-0.5*w*i);g2plot(g2);xlabel(real g2(w);ylabel(inag g2(w)axis(square);（3） 多项式运算1） 建立多项式向量ap=1 2 2 1;b=-1 -2 -3;bp=poly(b)2） 练习多项式乘与求根p=conv(ap,bp)roots(p)3

4、） 练习多项式运算a=1 2 3 4;b=1 -1;c=a+zeros(1,length(a)-length(b) b;poly2str(c,x)polyvalm(a,3)（4） 代数方程组1）恰定方程组A=1 2;2 3;b=8;13;方法1：逆矩阵求解X=inv(a)*b方法2：矩阵左除求解X=ab3 实验报告要求按上述步骤完成数值运算实验，并按练习实例做好记录。实验二 MATLAB语言绘图1 实验目的（1）学习MATLAB的各种二维绘图（2）学习MATLAB的三维绘图（3）MATLAB的绘图修饰2 实验内容（1）基本二维绘图1）向量绘图x=0:2*pi/100:2*pi;y1=sin(2

5、*x);y2=cos(2*x);plot(x,y1)plot(x,y2)plot(x,y1,x,y2)plot(x,y1);hold on;plot(x,y2);hold off;plot(x,y1 y2)plot(x,y1,c:,x,y2,ro)figure(1);plot(x,y1);figure(2);plot(x,y2);subplot(2,2,1);plot(x,y1);subplot(2,2,2);plot(x,y2);subplot(2,2,3);plot(x,y1,x,y1+y2);subplot(2,2,4);plot(x,y2,x,y1-y2);2）函数绘图fplot(si

6、n,0 4*pi)f=sin(x);fplot(f,0 4*pi)fplot(sin(1/x),0.01 0.1,1e-3)fplot(tan(x),sin(x),cos(x),-2*pi,2*pi,-2*pi,2*pi)3）符号函数快捷绘图syms xf=exp(-0.5*x)*sin(x);ezplot(f,0,10)（2）多种二维绘图1）半对数绘图w=logspace(-1,1);g=20*log10(1./(1+2*w*i);p=angle(1./(1+2*w*i)*180/pi;subplot(2,2,1);semilogx(w,g);grid;subplot(2,1,2);semi

7、logx(w,p);grid;2）极坐标绘图t=0:2*2*pi/180:2*pi;mo=cos(2*t);polar(t,mo);3）直方图t=0:2*pi/8:2*pi;y=sin(t);bar(t,y)4）离散棒图t=0:2*pi/8:2*pi;y=sin(t);stem(t,y)5）阶梯图t=0:2*pi/8:2*pi;y=sin(t);stairs(t,y)（3）图形注释y1=dsolve(D2u+2*Du+10*u=0,Du(0)=1,u(0)=0,x);y2=dsolve(D2u+2*Du+10*u=1,Du(0)=0,u(0)=0,x);ezplot(y1,0,5);hold

8、on;ezplot(y2,0,5);axis(0,5,-0.1,0.2)title(二阶系统时间响应)；xlable(时间t);ylabel(响应幅值y);gtext(零输入响应)；gtext(零状态响应)；grid;hold off;（4）三维绘图1）三维线图t=0:pi/50:10*pi;plot3(sin(t),cos(t),t);comet3(sin(t),cos(t),t);2）单变量高度网线图Z2=1 1;1 -1;Z4=Z2 Z2;Z2 Z2;Z8=Z4 Z4;Z4 Z4;Mesh(Z8);3）3变量马鞍面网线图x-4:0.5:4,y=x;X,Y=meshgrid(x,y);Z=

9、X.2-Y.2;Mesh(X,Y,Z);4）圆锥面网线图t=0:0.1:0.9;t=1:0.1:2;r=t1 t2+2;X,Y,Z=cylinder(r,40);Mesh(X,Y,Z):5）视角修饰t=0:0.1:0.9;t=1:0.1:2;r=t1 t2+2;X,Y,Z=cylinder(r,30);Mesh(X,Y,Z);subplot(2,2,1);mesh(X,Y,Z);view(0,0);subplot(2,2,2);mesh(X,Y,Z);view(-20,20);subplot(2,2,3);mesh(X,Y,Z);view(-30,30);subplot(2,2,4);mesh

10、(X,Y,Z);view(-40,40);3 实验报告要求按照上述步骤进行实验，并按实验记录要求完成报告实验三 控制系统阶跃响应与脉冲响应实验1 实验目的（1）观察学习控制系统的单位阶跃响应和单位脉冲响应（2）记录单位阶跃响应和单位脉冲响应曲线（3）掌握时间响应分析的一般方法2 实验步骤（1）运行MATLAB（2）建立系统模型1）传递函数模型TF2）ZPK模型3）MATLAB的阶跃响应函数3 实验内容已知二阶系统： 1） 建立系统模型，观察阶跃响应曲线和单位脉冲响应，并计算系统的闭环根、阻尼比，无阻尼振荡频率，并作记录。2） 修改参数，分别实验=1，=2的响应曲线，并作记录。4 实验报告要求（

11、1）分析系统的阻尼比和无阻尼振荡频率对系统阶跃响应和脉冲响应的影响；（2）分析响应曲线的零初值、非零初值与系统模型的关系；（3）分析响应曲线的稳态值与系统模型的关系；（4）分析系统零点对阶跃响应曲线和单位脉冲响应曲线的影响；实验要求：记录给定系统的响应曲线，对并上述问题作出分析实验四 控制系统根轨迹分析1实验目的（1） 利用计算机完成控制系统的根轨迹作图。（2） 了解控制系统根轨迹作图的一般规律。（3） 利用根轨迹进行系统分析。2 实验步骤（1） 运行matlab。（2） 练习根轨迹的相关函数。rlocus(sys)rlocus(sys,k)r=rlocus(sys)r,k=rlocus(sy

12、s)（3） 实验内容1） 给定如下系统的开环传递函数，作出它们的根轨迹，并完成给定要求。要求：准确记录根轨迹的起点、终点和根轨迹的条数。 确定根轨迹的分离点与相应的根轨迹增益。 确定临界稳定时根轨迹增益k。 要求：确定根轨迹与虚轴交点并确定系统稳定的根轨迹增益k。 要求：确定系统具有最大超调量时Mp,max时的根轨迹增益，作时域仿真验证。 确定系统阶跃响应无超调量时的根轨迹增益取值范围，并作时域仿真验证。3 实验要求：（1） 记录给定系统与显示的根轨迹图；（2） 完成上述各要求，分析闭环极点在s平面上的位置与系统动态性能的关系。实验五 控制系统Simulink仿真实验1 实验目的1. 学习结构图编程，掌握结构图系统文件的设计方法2. 对于给定的控制系统，通过Simulink仿真观察系统的时域响应。2 实验步骤（1）运行matlab，执行simulink。（2）建立图示仿真结构。（3） 构图单元设置1） 激活信号单元2） 激活加法器单元3） 激活连续单元模块4） 激活示波器单元5） 进行仿真参数设置3 实验要求：（1） 作闭环结构时域仿真