MATLAB基础知识考试复习总结

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

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

3、的变量 exit、quit 关闭、退出 MATLAB6 c=3+5ic=3+5*i(a=3;b=5;c=a+b*i)1.3e-4 2.78e23 A.'转置 S.*B S./B B.S A.An7Comma nHistory (历史指令)窗记录着用户在命令指令窗中所输入过的所有指 令行，且所有这些被记录的指令行都能被复制，并送到指令窗中再运行。8Workspace Browser (工作空间浏览器)也叫内存浏览器，他保存了指令窗所使 用过的全部变量(除非有意删除),可通过该浏览器对内存变量进行操作。10点击MATLAB面工具条上的？图标，或选择下拉菜单项【Help】，都能提供 帮助；M

4、ATLAB还提供现场帮助，用鼠标点亮指令并点击右键，在弹出的菜 单中选择【Help On Selection】。第二章数据及其运算1简单数组生成方法：逐个元素输入法；冒号生成法x=a:i ntc:b;x=li nspace(a,b, n)=a:(b-a)/( n-1):b;logspace(w1,w2, n)2>> diag(3,3,3)产生对角形数组ans =300030003>> eye(3)产生单位数组ans =1000 1 00 0 1>> magic(3)产生魔方数组ans =81635492>> on es(3)ans =11 11

5、1>> rand(3)ans =0.81470.90580.63240.12700.0975>> zeros(3)ans =00 00 0产生全1数组产生0、1间随机数组0.27850.91340.54690.9575产生全0数组>> zeros(3,2)ans =000 00 03一位数组寻访>> x=2.0000 1.0472 1.7321 3.0000+5.0000i;>> x1=x(3)x1 =1.7321>> x2=x(1 2 4)x2 =2.0000>> x3=x(2:e nd)x3 =1.0472

6、>> x4=x(4:-1:1)x4 =3.0000 + 5.0000i 1.73211.04721.73213.0000 + 5.0000i3.0000 + 5.0000i1.04722.00004二维数组寻访：A (:, j)表示A矩阵的第j列全部元素；A 阵第i行全部元素；A( 1:3，2:4 )表示对A矩阵取第13 行, 兀素。例：>> A=rand(3,5)A =0.11900.34040.49840.58530.95970.2238>> A1=A(1,:)A1 =0.1190 0.3404 0.7513 0.6991 0.5472>>

7、A2=A(1:2,2:5)A2 =0.3404 0.7513 0.6991 0.5472i ,:）表示A矩 24列中所有0.75130.25510.50600.69910.89090.95930.54720.13860.14930.58530.25510.89090.1386>> A3=A(1,3,2,5)A3 =0.34040.54720.22380.14935非数(Not a Number)指的是0/0、x/x、0x之类的运算,NaN或 nan 表示。非数参与运算所得结果也是非数，具有传递性；非数没有大小，不能比较两 个非数的大小。6在MATLA中，空数组除了用表示外，某维或若

8、干维长度均为 0的数 组都是空数组7逻辑运算符：o&与1 例 t=(0:pi/50:2*pi)' k=0.4:0.1:1;X=cos(t)*k; plot(X)； Plot(t,X)|或 非第二草数据和数组的可视化%生成（101x1）的时间采样列向量%生成（1x7）的行向量生成（101x7）的矩阵绘制曲线，横坐标为每列元素对应的下标以t为横坐标，X为纵坐标，按丫的列数绘制曲线线型符号1-:-.含义丁实线虚线点划线双划线色彩符号b|grcmykw含义二蓝1蓝绿红青品红黄黑八、白%不画分格线使坐标开启在（x,y）处写字符注释 纵坐标轴名2 grid on box on title(

9、s) xlable(sgrid off box off text(x,y,s ) ylable(s)画出分格线使坐标圭寸闭书写图名)横坐标轴名多层叠绘axis ( 0,1,-1,17)坐标范围hold on ,hold offalpha(0)完全透明，alpha(0.5)半透明，alpha完全不透明 plotyy(X1,Y1,X2,Y2 )双纵坐标3 view(a z, el)View(vx,xy,xz)shadi ngflat,i nterpfaceted通过方位角、俯视角设置视点，例view （-82,58）通过直角坐标设置视点3 M文件包括M脚本文件和M函数文件；M脚本文件的执行过程与在

10、指令窗中直 接输入指令的效果一样，但效率更高；M函数文件与脚本文件类似之处在于他们 都是一个扩展名为“ .m”的文件；M函数文件通常由以下几个部分组成：函数定 义行、H1行、函数帮助文本、函数体、注释。legend(s1,s2,)在图的右上角建立图例3 例：t=(0:0.02:2)*pi;x=si n( t);y=cos(t);z=cos(2*t);plot3(x,y,z,'b-',x,y, z,'bd')box onlegend('链','宝石')4 >>x=-4:4;y=x;X,Y=meshgrid(x,y);曲面

11、图网线图；box on %曲线图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) 5>> x,y=meshgrid(-3:0.1:3,-2:0.1:2);z=(x.A2+2*x).*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')hidde n off% 透视被叠压图形subplot

12、(1,2,2),mesh(x,y,z)>> title('xiao yin')hidden on%消隐被叠压图形>> axis(-3,3,-2,2,-0.5,1.0)第四章MATLAB编程1for x=array说明：循环体执行的次数由数组 array的列数决定(commands)例: for x=0.2:0.2:0.8endb=sqrt(1-xA2);end2while expressi on comma nds说明：当expression为真时，执行循环体直 到为假end3if expressi onif expressi oncomma ndsco

13、mma ndsendelseend4 ezplot(f )其中f=f(x,y)为用符号函数表示的隐函数 例：ezplot( ' xA2+x*y+yA2 -10')5x,y=ginput(n)该指令用鼠标从图形上获取n个点的坐标(x,y)6例：用冒泡法对数组 A=1 6 8 9 3 5 6进行升序排序解: >> r=length(A);for i=1:r-1for j=i+1:rif A(i)>A(j)%A(i)<A(j)temp=A(i);A(i)=A(j);A(j)=temp;endendendy=A时，进行降序排序第五章系统模型1计算机仿真过程：建模

14、；模型实现；仿真分析。2 num=5,3;den=1,6,11,6; %sys=tf( num,de n)sys=tf( nu m,de n)Tran sfer fun cti on:5 s + 3sA3 + 6 sA2 + 11 s + 6 sys=tf(5 3,1 6 11 6)Tran sfer fun cti on:5 s + 3 sA3 + 6 sA2 + 11 s + 6 z=-0.6;p=-3 -2 -1; sys=zpk(Z,P,K) k=5;sys=zpk(z,p,k)Zero/pole/ga in:5 (s+0.6)(s+3) (s+2) (s+1)>> sys

15、=zpk(-0.6,-3 -2 -1,5)Zero/pole/ga in:5 (s+0.6) (s+3) (s+2) (s+1)>> z,p,k=zpkdata(sys,'v')z =-0.6000p = -3-2-1k =5>> nu m,de n=tfdata(sys,'v')num =0053den =16 1163 pzmap(sys)注：该指令将在S平面上用符号“ O'表示零点，符号“ X”表示 极点。4 例：延时模型。G(s)=eA-0.5t(5s+3/sA3+6sA2+11s+6)num=5 3 ;sys=tf( n

16、u m,de n,' in putdelay',tao)den=1 6 11 6; sys=zpk(z,p,k,'i nputdelay',tao)sys=tf( nu m,de n,'i nputdelay',0.5)Tran sfer fun cti on:5 s + 3exp(-0.5*s) *sA3 + 6 sA2 + 11 s + 65 串联模型：sys=series(sys1,sys2)=sys1*sys2 并联模型：sys=parallel(sys1,sys2)=sys1+sys2 反馈连接：sys=feedback(sys1,sy

17、s2,sig n)第六章系统时间响应及其仿真1对连续时间系统进行数字动态仿真，主要有两种方法：基于数值积分的仿真方 法；基于离散相似法的仿真方法。2数值积分方法所得到的离散数值解只是精确解的近似，其误差来自于两个方面，即舍入误差和局部截断误差。舍入误差与计算步长h成反比，局部截断误差 与h成正比，显然选择一个合适的积分步长可使总误差达到最小。3数值积分常用方法：欧拉法、梯度法、龙格-库塔法、Gear法4 step(sys) step(sys,Tfi nal) step(sys,T)阶跃响应仿真说明：sys是系统模型；Tfinal为仿真时间，若省略则由系统默认；T为仿真时 间向量，T=TO:dt

18、:Tfinal,dt为连续系统离散化的采样周期，T0为仿真开始时间。5 impulse(sys)脉冲响应仿真6U,T=gensig(Type,Tau) 信号发生器u,t=gensig( square' ,5,10)说明： Type为信号类型：sin '、square '、pulse ' Tau为信号周期；U为信号值向 量；T为与U对应的时间向量。7 lsim(sys,U,T)任意输入相应函数8 initial (sys,x0 )用于计算零输入状态下，由初始状态x0所引起的响应，只能用于状态空间模型。9T,y =ode45( f' ,tspan,yo,op

19、tions)为一种单步显式，变步长 RK-45算法,用于求解非刚性微分方程；ode23也一样；ode113属于多步法，求解非刚性。10例：已知某系统运动方程及初始条件为：y1'=y2y3; y2'=-y1y3; y3'=-2y1y2y1(0)=0; y2(0)=0.5; y3(0)=-0.5求时间区间t=0,20微分方程的解。解:(1)建立描述系统微分方程的 m-函数文件rigit.mfunction dy=rigit(t,y) dy=zeros(3,1);dy(1)=y(2)*y(3); dy(2)=-y(1)*y(3);dy(3)=-2*y(1)*y(2);编写调用

20、函数rigit() 的M文件，并执行T,y=ode45('rigit',0,20,0,0.5.-0.5);Plot(T,y(:,1),'r',T,y(:,2),'b*',T,y(:3),'k-.');lege nd('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

21、-y1A2)y2-y1(2)建立描述系统微分方程的m-函数文件vdp.mfun cti on dy=vdp(t,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=(1-y(1)A2)*y(2)-y(1)(3)编写MATLA主程序，并执行T,y=ode45('vdp',0,20,0,1);Plot(T,y(:,1),'r-',T,y(:,2),'b:');lege nd('y1','y2')第七章系统频率响应及其仿真1频率响应是指系统对谐波输入的稳态响应；频率特性是指系统在正弦信号作用下，稳态输出与输

22、入之比对频率的关系特性。G(jw)=Xo(jw)/Xi(jw)=A(jw)eAj© (w)2绘制G (s)=11(s+1)/s(sA2+15s+4)的幅频特性曲线和相频特性曲线解：nu m=11,11;de n=1 15 4 0;w=0.05:0.01:0.5*pi;Gw=polyval( nu m,j*w)./polyval(de n,j*w); mag=abs(Gw); %(p=0) wc<wg 系统稳定theta=a ngle(Gw); wc>wg不稳定subplot(2,1,1),plot(w,mag); grid on;title('幅频特性)ylabl

23、e('|G|'subplot(2,1,2),plot(w,theta); grid on ;title('相频特性)Xlable('omega(rad/s)'),ylable('deg')3多项式函数 Y=polyval(P,X) P:系数向量 X :自变量4系统频率响应函数freqs(b,a,w) ba为传递函数分子分母系数向量h=freqs(b,a,w0 h,w=freqs(b,a) freqs(b,a,w) 5nyq uist(sys) nyq uist(sys,w) nyq uist(sys1,指定正实角频率，返回响应值；自动确定

24、200个频率点，返回响应值和对应角频率 绘制对指定正实角频率向量的幅值和相频特性曲线。绘制nyquist图指定频率范围w,绘制nyquist图sys2,，sysn,w)在同一坐标系绘制多个nyquistRe,lm,w=nyquist(sys)返回频率响应的实部和虚部及对应w,不绘图bode(sys) bodemag(sys)绘制bode图仅绘制幅频bode图7 margin(sys)绘制Bode图，标出幅值与相位裕度。kg,y,wg,Wc=margin (sys)返回幅值裕度kg，相位裕度y,相位穿越频率 Wc 和幅值穿越频率Wg不绘制Bode图。第八章控制系统的综合校正1控制系统的性能指标大

25、体上可分为两类，即频域指标和时域指标。2已知系统传递函数为 G(s)=5(s+1)/(sA3+4sA2+6s+5),求最大超调量Mp,调整时间ts和峰值时间tp.解: >> num=5 5;de n=1 4 6 5;sys=tf( nu m,de n);fin alvalue=polyval( nu m,0)/polyval(de n,0);y,t=step(sys); %y(x )=lim sY(s)=limG(s)yp,k=max(y);tp=t(k); G(s)=Y(s)/(1/s)=sY(s)Mp=100*(yp-fi nalvalue)/fi nalvalue;len=l

26、e ngth(t);while(y(le n) >0.98*fi nalvalue)&(y(le n)<1.02*fi nalvalue) len=len-1;e nd;ts=t(le n);step(sys)3已知单位负反馈系统的开环传递函数为 入时，系统稳态误差。Gk=10/s(s+1)(s+5),求其单位斜坡输解：G k=zpk(,O -1-5,10);Xi=zpk(,O 0,1)；sys=1/(1+Gk);Es=sys*Xi;%dcga in( sys)=lim sys(s)ess=dcgain(tf(1 0,1)*Es);% 计算稳态偏差 £ ss=lim sE(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'