MATLAB基础知识考试复习总结

1、第一章MATLAB基础1系统仿真是根据被研究的真实系统的数学模型研究系统性能的一门学科，现在尤指利用电脑去研究数学模型行为的方法，即数值仿真。MATLAB集计算，可视化及编程于一身。其主要产品模块构成：1MATLAB2MATLABtoolbox3MATLABCompiler4simulink(5)stateflow(6)Real-TimeWorkshop。3MATLAB5言被称为第四代电脑语言。有以下几个主要特点：编程效率高；使用方便；高效方便的科学计算；先进的可视化工具；开放性、可扩展性强；运行时动态连接外部C或FORTRAN用函数；在独立C或FORTRAN序中调用MATLAB函数；输入输出

2、各种MATLAB及其他标准格式的数据文件；创建图文并茂的技术文档；特殊应用工具箱；高效仿真工具Smulink。只能是英文、数字和下连线求复数Z实部求复数Z的模回校函数0imag(Z)求复数Z虚部angle(Z)求复数Z的相角单位是弧度mdata=csvread(engdata.txt)clf活除图形窗clear活除工作空间保存的变量exit、quit关闭、退出MATLAB4变量命名规则：变量名、函数名对字母大小写敏感；变量名第一个字母必须是英文字母5real(z)abs(z)callbackclc活除指令窗cd设置当前工作目录edit打开M文件编辑器6c=3+5ic=3+5*i(a=3;b=5

3、;c=a+b*i)1.3e-4A.转置S.*BS./BB.SA.An7CommanHistory历史指令窗记录着用户在命令指令窗中所输入过的所有指令行，且所有这些被记录的指令行都能被复制，并送到指令窗中再运行。8WorkspaceBrowser(X作空间浏览器也叫内存浏览器，他保存了指令窗所使用过的全部变量除非有意删除),可通过该浏览器对内存变量进行操作。10点击MATLAEg面工具条上的？图标，或选择下拉菜单项【Help】，都能提供帮助；MATLAB还提供现场帮助，用鼠标点亮指令并点击右键，在弹出的菜单中选择【HelpOnSelection】。第二章数据及其运算1简单数组生成方法：逐个元素输

4、入法；冒号生成法x=a:intc:b;x=linspace(a,b,n)=a:(b-a)/(n-1):b;logspace(w1,w2,n)2diag(3,3,3)产生对角形数组ans=300030003eye(3)产生单位数组ans=100010001magic(3)ans=816产生魔方数组2 5792ones产生全1数组ans=1111 11111rand产生0、1间随机数组zeros(3)产生全0数组ans=000000000zeros(3,2)ans=0000003一位数组寻访x=2.00001.04721.73213.0000+5.0000i;x1=x(3)x2=x(124)x3=

5、x(2:end)x4=x(4:-1:1)4二维数组寻访：A：，j)表示A矩阵的第j歹0全部元素；Ai,:表示A矩阵第i行全部元素;A1:3,2:4表示对A矩阵取第13行，第24列中所有元素。例：A=rand(3,5)A1=A(1,:)A2=A(1:2,2:5)A3=A(1,3,2,5)5非数NotaNumber）指的是0/0、00/8、0x8之类的运算，用Nanan表示。非数参与运算所得结果也是非数，具有传递性；非数没有大小，不能比较两个非数的大小。6在MATLAB,空数组除了用表示外，某维或假设干维长度均为0的数组都是空数组|或非弟二草数据和数组的可视化%生成101x1）的时间采样列向量%生

6、成1x7的行向量生成101x7）的矩阵绘制曲线，横坐标为每列元素对应的下标以t为横坐标，X为纵坐标，按Y的歹U数绘制曲线1例t=(0:pi/50:2*pi);k=0.4:0.1:1;X=cos(t)*k;plot(X)；Plot(t,X)多层叠绘axis0,1,-1,17坐标范围holdon,holdoffalpha(0)完全透明，alpha(0.5)半透明，alpha(1)完全不透明plotyy(X1,Y1,X2,Y2)双纵坐标3view(az,el)View(vx,xy,xz)shadingflat,interp,faceted通过方位角、俯视角设置视点，例view-82,58通过直角坐标

7、设置视点7逻辑运算符：线型符号-:-.含义r实线虚线点划线双划线色彩符号bgrcmykw含义蓝绿红青品红黄里八、白%2gridonboxontitle(s)xlable(s不画分格线使坐标开启在x,y）处写字符注释纵坐标轴名gridoffboxofftext(x,y,s)ylable(s)画出分格线使坐标封闭书写图名)横坐标轴名3 M文件包括M脚本文件和M函数文件；M脚本文件的执行过程与在指令窗中直接输入指令的效果一样，但效率更高；M函数文件与脚本文件类似之处在丁他们都是一个扩展名为“.m”的文件；M函数文件通常由以下几个部分组成：函数定义行、H1行、函数帮助文本、函数体、注释。legend(

8、s1,s2,在图的右上角建立图例3 例：t=(0:0.02:2)*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z,b-,x,y,z,bd)boxonlegend(链,宝石)4 x=-4:4;y=x;X,Y=meshgrid(x,y);曲面图网线图；boxon%曲线图Z=X.A2+Y.A2;subplot(1,3,1),surf(X,Y,Z)%subplot(1,3,2),mesh(X,Y,Z)%subplot(1,3,3),plot3(x,y,x.A2+y.A2)5x,y=meshgrid(-3:0.1:3,-2:0.1:2);z=(x.A2+2*x).

9、*exp(-x.A2-y.A2-x.*y);subplot(1,2,1),mesh(x,y,z)axis(-3,3,-2,2,-0.5,1.0)title(toushi)hiddenoff%透视被叠压图形subplot(1,2,2),mesh(x,y,z)title(xiaoyin)hiddenon%消隐被叠压图形axis(-3,3,-2,2,-0.5,1.0)第四章MATLABg程1forx=array说明：循环体执行的次数由数组array的歹U数决定(commands)例：endb=sqrt(1-xA2);end2whileexpression说明：当expression为真时，执行循环体

10、直commands到为倡iend3ifexpressionifexpressioncommandscommandsendelseend4 ezplot(f)其中f=f(x,y)为用符号函数表示的隐函数例：ezplot(xA2+x*y+yA2-10)5x,y=ginput(n)该指令用鼠标从图形上获取n个点的坐标x,y6例：用冒泡法对数组A=1689356进行升序排序解：r=length(A);fori=1:r-1forj=i+1:rifA(i)A(j)%A(i)sys=zpk(-0.6,-3-2-1,5)Zero/pole/gain:5(s+0.6)(s+3)(s+2)(s+1)z,p,k=z

11、pkdata(sys,v)p=-3-23 -1k=5num,den=tfdata(sys,v)num=0053den=16116pzmap(sys)注：该指令将在S平面上用符号“。表示零点，符号“X”表示极点。例：延时模型。G(s)=eA-0.5t(5s+3/sA3+6sA2+11s+6)num=53;sys=tf(num,den,inputdelay,tao)den=16116;sys=zpk(z,p,k,inputdelay,tao)sys=tf(num,den,inputdelay,0.5)Transferfunction:5 s+3exp(-0.5*s)*sA3+6sA2+11s+6申

12、联模型：sys=series(sys1,sys2)=sys1*sys2并联模型：sys=parallel(sys1,sys2)=sys1+sys2反馈连接：sys=feedback(sys1,sys2,sign)第六章系统时间响应及其仿真1对连续时间系统进行数字动态仿真，主要有两种方法：基于数值积分的仿真方法；基于离散相似法的仿真方法。2数值积分方法所得到的离散数值解只是精确解的近似，其误差来自于两个方面，即舍入误差和局部截断误差。舍入误差与计算步长h成反比，局部截断误差与h成正比，显然选择一个合适的积分步长可使总误差到达最小。3数值积分常用方法：欧拉法、梯度法、龙格-库塔法、Gear法ste

13、p(sys)step(sys,Tfinal)step(sys,T)阶跃响应仿真说明：sys是系统模型；Tfinal为仿真时间，假设省略则由系统默认；T为仿真时间向量，T=T0:dt:Tfinal,dt为连续系统离散化的采样周期，T0为仿真开始时间。impulse(sys)脉冲响应仿真6U,T=gensig(Type,Tau)信号发生器u,t=gensig(square,5,10)说明：Type为信号类型：sin、square、pulseTau为信号周期；U为信号值向量；T为与U对应的时间向量。lsim(sys,U,T)任意输入相应函数initialsys,x0用于计算零输入状态下，由初始状态x

14、0所引起的响应，只能用于状态空间模型。9T,y=ode45(f,tspan,yo,options)为一种单步显式，变步长RK-45算法,用于求解非刚性微分方程；ode23也一样；ode113届于多步法，求解非刚性。10例：已知某系统运动方程及初始条件为：y1=y2y3;y2=-y1y3;y3=-2y1y2求时间区问t=0,20微分方程的解。解:(1)建立描述系统微分方程的m-函数文件functiondy=rigit(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(3);dy3)=-2*y(1)*y(2);编写调用函数rigit()的M文件，并执行

15、T,y=ode45(rigit,0,20,0,0.5.-0.5);Plot(T,y(:,1),r,T,y(:,2),b*,T,y(:3),k-.);legend(y1,y2,y3)11已知二阶微分方程y-(1-yA2)y+y=0y(0)=0,y(0)=1求时间区问t=0,20微分方程的解。解：1将微分方程表示成一阶微分方程组：y1=y;y1=y2;y2=(1-y1A2)y2-y1(2)建立描述系统微分方程的m-函数文件functiondy=vdp(t,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=(1-y(1)A2)*y(2)-y(1)(3)编写MATLABt程序，并执行T

16、,y=ode45(vdp,0,20,0,1);Plot(T,y(:,1),r-,T,y(:,2),b:);legend(y1,y2)第七章系统频率响应及其仿真1频率响应是指系统对谐波输入的稳态响应；频率特性是指系统在正弦信号作用下，稳态输出与输入之比对频率的关系特性0G(jw)=Xo(jw)/Xi(jw)=A(jw)eAj4(w)2绘制G(s)=11(s+1)/s(sA2+15s+4)的幅频特性曲线和相频特性曲线解：num=11,11;den=11540;w=0.05:0.01:0.5*pi;p=0wcwgsubplot(2,1,1),plot(w,mag);gridon;title(幅频特性

17、ylable(|G|subplot(2,1,2),plot(w,theta);gridon;title(相频特性系数向量X:自变量、a为传递函数分子分母系数向量Xlable(omega(rad/s),ylable(deg)3多项式函数Y=polyval(P,X)P:4系统频率响应函数freqs(b,a,w)bh=freqs(b,a,w0h,w=freqs(b,a)freqs(b,a,w)5nyquist(sys)nyquist(sys,w)nyquist(sys1,指定正实角频率，返回响应值；自动确定200个频率点，返回响应值和对应角频率绘制对指定正实角频率向量的幅值和相频特性曲线。绘制nyq

18、uist图指定频率范围w,绘制nyquist图sys2,，sysn,w)在同一坐标系绘制多个nyquistRe,Im,w=nyquist(sys)bode(sys)bodemag(sys)7margin(sys)返回频率响应的实部和虚部及对应w,不绘图6绘制bode图绘制Bode图，标出幅值与相位裕度。仅绘制幅频bode图kg,y,wg,Wc=marginsys返回幅值裕度kg,相位裕度y,相位穿越频率Wc和幅值穿越频率Wg不绘制Bode图。第八章控制系统的综合校正1控制系统的性能指标大体上可分为两类，即频域指标和时域指标。2已知系统传递函数为G(s)=5(s+1)/(sA3+4sA2+6s+

19、5),求最大超调量Mp,调整时问ts和峰值时间tp.解：num=55;den=1465;sys=tf(num,den);finalvalue=polyval(num,0)/polyval(den,0);y,t=step(sys);%y(00)=limsY(s)=limG(s)yp,k=max(y);tp=t(k);G(s)=Y(s)/(1/s)=sY(s)Mp=100*(yp-finalvalue)/finalvalue;len=length(t);while(y(len)0.98*finalvalue)&(y(len)1.02*finalvalue)len=len-1;end;ts=t(le

20、n);step(sys)3已知单位负反馈系统的开环传递函数为Gk=10/s(s+1)(s+5)，求其单位斜坡输入时，系统稳态误差。解：Gk=zpk(,0-1-5,10);Xi=zpk(,00,1);sys=1/(1+Gk);Es=sys*Xi;%dcgain(sys)=limsys(s)ess=dcgain(tf(10,1)*Es);%计算稳态偏差ss=limsE(s)t=0:0.05:10;xi=t;E(s)=Xi(s)/(1+G(s)H(s)y=lsim(sys*Gk,xi,t);E(s)=H(s)E1(s)plot(t,xi,r-.,t,y,t,xi-y,k:)legend(shuru,shuchu,wucha,0);xlable(t(s),ylable(fuzhi、wucha)4频域指标包含：1开环频域指标：主要指幅值穿越频率