《数控技术matlab在自控制原理中的应用》由会员分享，可在线阅读，更多相关《数控技术matlab在自控制原理中的应用（103页珍藏版）》请在金锄头文库上搜索。
1、第第8章章 Matlab在自动控在自动控制原理中的应用制原理中的应用n研究线形系统中的时不变系统u系统复杂，具有各种数学描述方法，有各个数学模型求出总的系统数学模型，研究系统特性的同时，注意研究系统中个部件的状态。u控制系统多为闭环反馈系统，一般研究其开环系统如何反映和影响系统的闭环特性u由于控制系统的应用广泛性，在Matlab的控制工具箱有很多，如鲁棒控制工具箱，神经网络工具箱等概述8.1 控制工具箱中的LTI对象n线性系统可以采用四种不同方法来进行描述，每种方法又有几个参数矩阵，基于面向对象的设计思想，建立专用的数据结构类型，将各种模型封装成统一的LTI对象。在一个名字下包含了改系统的所有
2、属性。n8.1.1 LTI 对象的类型和属性对象的类型和属性控制工具箱中一般有三个对象ss tf zpk共同属性p238 表8.1 Ts0，Ts1 Td 输入时延 特有属性p239 表8.28.1.2 LTI模型的建立n ndss(a,b,c,d)生成描述状态空间模型n nfilt(num,den)生成dsp形式的离散传递函数(z-1)n nss(a,b,c,d)生成状态空间模型n ntf(num,den,)生成传递函数模型n nzpk(z,p,k)生成零极增益模型8.1.3对象属性的获取和修改n n对象属性的获取uuget uuset n n模型的参数转换和提取uu dssdatauu ss
3、datauu tfdatauu zpkdatan n输出反馈n n反馈联接n n系统串联联接n n系统并联联接n n系统增广联接系统模型的连接输出反馈系统的结构图下图所示。可以使用函数cloop得到闭环系统的数学模型。格式1：numc,denc=cloop(num,den,sign)输入开环系统的传递函数，左变量返回系统的闭环参数，numc为分子多项式系数向量，denc为分母多项式系数向量；右变量中的sign=1为正反馈，sign=-1为负反馈。sign的缺省值为。+-uy输出反馈格式2：Ac,Bc,Cc,Dc=cloop(A,B,C,D,sign)MIMO系统时，状态空间模型的全输出反馈。格
4、式3： Ac,Ac,BcBc,Cc,Dc=,Cc,Dc=cloopcloop(A,B,C,D,outputs,inputs)(A,B,C,D,outputs,inputs)MIMO系统时，状态空间模型的选择输出反馈。向量outputs，inputs分别指明输出输入端口号数，正反馈时inputs取正值，负反馈时inputs取负值。例：单位反馈系统结构图如图所示，求闭环系统的数学模型。+-uynum=10;den=1,2,0;nc,dc=cloop(num,den,-1);printsys(nc,dc)num/den = 10 - s2 + 2 s + 10例：双输入、双输出系统结构图如下图所示，
5、求闭环系统的数学模型。a a=0 0;0 0;b=10,0;0,2;c=1 0;0 1;d=0 0;0 0;=0 0;0 0;b=10,0;0,2;c=1 0;0 1;d=0 0;0 0; ac,ac,bcbc,cc,dc=cloop(a,b,c,d,1,2,1,2,-1,-2,-1,-2),cc,dc=cloop(a,b,c,d,1,2,1,2,-1,-2,-1,-2)+-u1y1+-u2y2-ac = -20 0 0 -4bc = 10 0 0 2cc = 1 0 0 1dc = 0 0 0 0一般反馈系统的结构图如图所示。可以使用函数feedback得到闭环系统的数学模型。格式1： nu
6、mcnumc, ,dencdenc=feedback(num1,den1,num2,den2,si=feedback(num1,den1,num2,den2,sign)gn)系统1和系统2均为传递函数模式时，左变量为返回的闭环系统参数；右变量中，num1 和den1为系统1参数向量，num2和den2为系统2参数向量，符号变量sign=1为正反馈，+-uy反馈联接符号向量符号向量sign=-1sign=-1为负反馈。符号向量为负反馈。符号向量signsign的的缺省缺省值为值为-1-1。格式格式2 2： Ac,Ac,BcBc,Cc,Dc=feedback(A1,B1,C1,D1,A2,B2,C
7、,Cc,Dc=feedback(A1,B1,C1,D1,A2,B2,C2,D2,sign)2,D2,sign)系统系统1 1与系统与系统2 2均为状态空间模式时的使用格式。均为状态空间模式时的使用格式。格式格式3 3： Ac,Ac,BcBc,Cc,Dc=feedback(A1,B1,C1,D1,A2,B2,C,Cc,Dc=feedback(A1,B1,C1,D1,A2,B2,C2,D2,outputs1,inputs1)2,D2,outputs1,inputs1)MIMOMIMO系统时，系统系统时，系统1 1与系统与系统2 2为状态空间模型的为状态空间模型的反馈联接。向量反馈联接。向量outp
8、uts1outputs1和和inputs1inputs1分别指明当系分别指明当系统统2 2作为反馈模块时，联接系统作为反馈模块时，联接系统1 1的输出输入端口的输出输入端口号数，正反馈时号数，正反馈时inputs1inputs1取正值，负反馈取正值，负反馈时时inputs1inputs1取负值。取负值。例：反馈控制系统结构图如下图所示，求例：反馈控制系统结构图如下图所示，求闭环系统的数学模型。闭环系统的数学模型。n1=10;d1=1,2,0;n1=10;d1=1,2,0;n2n2=0.2,1;d2=0.01,1;=0.2,1;d2=0.01,1; n,d=feedback(n1,d1,n2,d
9、2,-1);n,d=feedback(n1,d1,n2,d2,-1);printsys(n,dprintsys(n,d) )num/den =num/den = 0.1 s + 10 0.1 s + 10 - - 0.01 s3 + 1.02 s2 + 4 s + 10 0.01 s3 + 1.02 s2 + 4 s + 10 +-uy系统串联联接的结构图如下图所示。在使用过程中可用函数series得到串联系统的数学模型。格式1：num,den=series(num1,den1,num2,den2)系统1和系统2均为多项式模型时，左变量为返回的闭环系统参数。右变量中num1和den1为系统1参
10、数向量；右变量中num2和den2为系统2参数向量。uy系统串联联接格式2： A,B,C,D=series(A1,B1,C1,D1,A2,B2,C2,D2)A,B,C,D=series(A1,B1,C1,D1,A2,B2,C2,D2)系统1与系统2均为状态空间模式时的使用格式。格式3：A,B,C,D=series(A1,B1,C1,D1,A2,B2,C2,D2,outputs1,inputs2)MIMO系统时，系统1与系统2之间的串联联接方式，向量outputs1指明系统1的输出号数，inputs2指明系统2和输出号数。例：已知系统串联校正结构图如下所示，求系统的传递函数。n1=0.5,1;d
11、1=0.1,1;n2=10;d2=1,2,0;no,do=series(n1,d1,n2,d2);nc,dc=cloop(no,do,-1);printsys(nc,dc)num/den = 5 s + 10 - 0.1 s3 + 1.2 s2 + 7 s + 10uy+-系统并联联接的结构图如图所示，可以使用函数parallel得到闭环系统的数学模型。具体使用格式如下：格式1： num,den=parallel(num1,den1,num2,den2)num,den=parallel(num1,den1,num2,den2)系统1和系统2均为传递函数时，左变量为返回的闭环系统参数；右变量中，
12、num1,den1为系统1参数向量，num2,den2为系统2参数变量。uy系统并联联接格式2：A,B,C,D=parallel(A1,B1,C1,D1,A2,B2,C2,D2)系统1与系统2均为状态空间模式时的使用格式。格式3： A,B,C,D=parallel(A1,B1,C1,D1,A2,B2,C2,D2,out1,out2)MIMO系统时，系统1与系统2之间的并联联接方式，向量out1指明系统1的输出号数，out2指明系统2的输入号数。例：系统结构图如图所示，求并联系统的数学模型。n1=2;d1=1 2;n2=5;d2=1,3;n,d=parallel(n1,d1,n2,d2)prin
13、tsys(n,d)num/den = 7 s + 16 - s2 + 5 s + 6uy系统增广联接的结构图如图所示。其中，系统1和系统2均为状态空间表达式，在使用中可利用函数append()得到增广联系统状态空间模型。使用格式如下： A,B,C,D=append(A1,B1,C1,D1,A2,B2,C2,D2)A,B,C,D=append(A1,B1,C1,D1,A2,B2,C2,D2)A1,B1,C1,D1为系统1的状态空间参数矩阵；A2,B2,C2,D2为系统2的状态空间参数矩阵；系统1系统2y1y2u2u1系统的增广联接A,B,C,D为增广系统的状态空间参数矩阵。例：系统结构图如图所示
14、，利用系统增广函数求闭环系统的数学模型。+-u1y1+-u2y2-ns1=10;ds1=1 0;a1,b1,c1,d1=tf2ss(ns1,ds1);ns2=2;ds2=1,0;a2,b2,c2,d2=tf2ss(ns2,ds2); as,as,bsbs, ,cscs, ,dsds=append(a1,b1,c1,d1,a2,b2,c2,d2);=append(a1,b1,c1,d1,a2,b2,c2,d2); a,b,c,d=cloop(as,bs,cs,ds,1 2 1 2,-1,-2,-1,-2)a,b,c,d=cloop(as,bs,cs,ds,1 2 1 2,-1,-2,-1,-2)
15、a = b = -20 0 1 0 0 -4 0 1c = d = 10 0 0 0 0 2 0 0 返回二、控制系统分析uu控制系统时域分析uu控制系统根轨迹分析uu控制系统频域分析uu控制系统稳定性分析返回n n控制系统时域响应分析阶跃响应函数功能：给定系统数学模型，求系统的单位阶跃响应。step(num,den) or step(A,B,C,D)step(num,den,t) or step(A,B,C,D,t)y,x=step(num,den) or y,x=step(A,B,C,D)y,x,t=step(num,den) or y,x=step(A,B,C,D,t)其中，num和de
16、n对应TF的系数向量。A,B,C,D对应SS模型的系数矩阵。格式1：给定num,den，求系统的阶跃响应并作图，时间向量t的范围自动设定。给定A, B, C,D求系统的阶跃响应并作图，时间向量的范围自动设定。格式2：时间向量t的范围由人工设定，等间隔。例：t=0:0.1:3。格式3：返回变量格式。返回输出变量y、状态变量x，不作图。格式4：返回变量格式。包括或不包括时间向量t，不作图。例：系统的传递函数为：求阶跃响应并作系统性能分析。绘制阶跃响应曲线num=4;den=1 1 4;step(num,den)计算峰值y,x,t=step(num,den);max(y)ans =1.4441求峰值
17、时间tp=spline(y,t,max(y)tp = 1.6062脉冲响应函数函数功能：给定系统数学模型，求系统的单位脉冲响应。impulse(num,den) or impulse(A,B,C,D)impulse(num,den,t) or impulse(A,B,C,D,t)y,x=impulse(num,den) or y,x=impulse(A,B,C,D)y,x,t=impulse(num,den) or y,x=impulse(A,B,C,D,t)其中，num和den对应TF的系数向量。A,B,C,D对应SS模型的系数矩阵。格式1：给定num,den，求系统的脉冲响应并作图，时间向
18、量t的范围自动设定。给定A, B, C,D求系统的脉冲响应并作图，时间向量t的范围自动设定。格式2：时间向量t的范围由人工设定，等间隔。例：t=0:0.1:3。格式3：返回变量格式。返回输出变量y、状态变量x至MATLAB窗口，不作图。格式4：返回变量格式。包括或不包括时间向量t，不作图。例：系统的传递函数为：例：系统的传递函数为：求阶跃响应并作系统性能分析。求阶跃响应并作系统性能分析。绘制脉冲响应num=4;den=1 1 4;impulse(num,den)计算误差面积y,x,t=impulse(num,den);trapz(t,y)ans =0.9983返回n n根轨迹作图与系统根轨迹分
19、析函数功能：给定系统的开环模型，绘制根轨迹图或计算绘图变量。用法：rlocus(num,den) or rlocus(num,den,k)r=rlocus(num,den,k) or r,k=rlocus(num,den)rlocus(a,b,c,d) or rlocus(a,b,c,d,k)r=rlocus(a,b,c,d,k) or r,k=rlocus(a,b,c,d)+-uy格式1：根据传递函数模型绘制系统的根轨迹图。格式rlocus(num,den)中，k为机器自适应产生的从0的增益向量；格式rlocus(num,den,k)中，k为人工给定的增益向量。格式2：返回变量格式，不作图。
20、k为返回增益值，r为返回的闭环根矩阵，矩阵行数为length(k)，其列数为(length(den)-1)。格式3：绘制SISO系统的根轨迹图；格式rlocus(a,b,c,d)中，k为机器自适应产生的从0的增益向量；格式rlocus(a,b,c,d,k)中，k为人工给定的增益向量。格式4：状态空间模型时返回变量格式。例：控制系统为绘制系统的根轨迹图，并确定系统临界稳定时系统参数。绘制系统的根轨迹图num=1;den=1,2,10,0;rlocus(num,den)求临界稳定时的系统参数r,k=rlocus(num,den);kg=spline(real(r(10:30,2),k(10:30)
21、,0) kg=20.0000例：非最小相位系统的开环传递函数为试绘制该系统的根轨迹，并确定使用闭环系统稳定的Kg范围。绘制根轨迹n=1 1;d1=1 0;d2=1 -1;d3=1 4 16;d=conv(d1,conv(d2,d3);rlocus(n,d)k1,p=rlocfind(n,d)Select a point in the graphics windowselected_point = -0.0046 + 1.5439ik1 =23.1109p1 =-1.5058 + 2.7146i -1.5058 - 2.7146i 0.0058 + 1.5486i 0.0058 - 1.5486
22、ik2,p2=rlocfind(n,d)Select a point in the graphics windowselected_point =-0.0046 + 2.5965ik2 =36.0909p2 = 0.0138 + 2.5855i 0.0138 - 2.5855i -1.5138 + 1.7627i -1.5138 - 1.7627i通过上述分析，可以得出系统稳定的Kg范围为23.1109Kg36.0909函数rlocfind()功能：用鼠标确定根轨迹上某一点的增益值和该点对应的n个闭环根。k,poles=rlocfind(num,den) or r,poles=rlocfind
23、(a,b,c,d)k,poles=rlocfind(num,den,p) or k,poles=rlocfind(a,b,c,d,p)格式1：在已经作出的根轨迹图上，用鼠标选择闭环极点的位置后，返回对应闭环极点的根轨迹增益值k和对应的n个闭环根，并在图上标注十字。使用该命令前，需要用rlocus()指令作出根轨迹图。格式2：给定右变量p为根轨迹上某点坐标值，返回对应该点的根轨迹增益值k和对应的n个闭环根。函数pzmap()功能：给定系统数学模型，作出零极点位置图。pzmap(num,den) or pzmap(a,b,c,d) or pzmap(p,z)p,z=pzmap(num,den) o
24、r p,z=pzmap(a,b,c,d)格式1：零极点绘图指令。零点标记为，极点标记为。格式2：返回零极点值，不作图。返回n n控制系统频域响应分析MATLAB提供了多个用于系统频域分析的函数。波特函数Bodebode(num,den) or bode(a,b,c,d,iu)bode(num,den,w) or bode(a,b,c,d,iu,w)m,p,w=bode(num,den) or m,p,w=bode(a,b,c,d,iu)函数功能：对数频率特性作图函数，即波特图绘图。格式1：给定开环系统的数学模型作波特图，频率向量w自动给出。当开环模型以状态空间模型给出时须指定第几个输入。格式2
25、：给定开环系统的数学模型作波特图，频率向量w由人工给出。w的单位为rad/s，可以由命令logspace得到对数等分的w值。格式3：返回变量式，不作图。其中，m为频率特性G(jw)的幅值变量，m=|G(jw)|；p为频率特性G(jw)的幅角向量，p=argG(jw)，单位为()；w为频率向量，单位为rad/s。例：控制系统的开环传递函数为作波特图，并确定谐振峰值的大小Mr与谐振频率wr。n=10;d=1 2 10;bode(n,d);m,p,w=bode(n,d);mr=max(m)mr = 1.6667wr=spline(m,w,mr)wr = 2.8284稳定裕度函数margin函数功能：
26、计算系统的稳定裕度，包括相位裕度Gm和幅值裕度Pm。margin(num,den) or margin(a,b,c,d)Gm,Pm,wg,wp=margin(num,den) or Gm,Pm,wg,wp=margin(a,b,c,d)Gm,Pm,wg,wp=margin(m,p,w)格式1：给定开环系统的数学模型作波特图，并在图上标注幅值裕度Gm和对应的频率wg，相位裕度Pm和对应的频率wp。幅值裕度Gm定义为：单位为dB(分贝)。相位裕度Pm定义为单位为度()。格式2：返回变量格式，不作图。返回幅值裕度Gm和对应的频率wg，相位裕度Pm和对应的频率wp。格式3：给定频率特性的参数向量、幅值
27、m、相位p和频率w，由插值法计算幅度裕度幅值裕度Gm和对应的频率wg，相位裕度Pm和对应的频率wp。例：控制系统结构图如图所示。求系统的稳定裕度，并分别采用格式2与格式3比较计算误差。n=2;d=1 3 20;margin(n,d)gm1,pm1,wg1,wp1=margin(n,d);m,p,w=bode(n,d);gm2,pm2,wg2,wp2=margin(m,p,w);gm1,gm2;pm1,pm2;wg1,wg2;wp1,wp2+-uy奈奎斯特函数nyquist()函数功能：奈奎斯特轨线作图命令，即极坐标图。nyquist(num,den) or nyquist(a,b,c,d)ny
28、quist(num,den,w) or nyquist(a,b,c,d,w)re,im,w=nyquist(num,den) or re,im,w=nyquist(a,b,c,d)格式1：给定开环系统的数学模型作波特图。频率w的范围自动给定。格式2：给定开环系统的数学模型作波特图。频率w的范围人工给定。格式3：返回极坐标图参数变量，不作图。其中re为复变函数G(jw)的实部变量，re=ReG(jw);im为复变函数G(jw)的虚部变量，im=ImG(jw);w为频率向量，单位为rad/s。例：二阶系统令n=1，分别作出=2,1,0.707,0.5时的奈奎斯特轨线。n=1;d1=1 4 1;d2
29、=1 2 1;d3=1 1.414 1;d4=1 1 1;nyquist(n,d1);hold onnyquist(n,d2);nyquist(n,d3);nyquist(n,d4);axis(equal);x=0:0.1:2*pi;plot(sin(x),cos(x),:)例：分别由例：分别由w w的自动变量与人工变量作下的自动变量与人工变量作下述系统的奈奎斯特轨线。述系统的奈奎斯特轨线。n n=1;=1;d d=1 1 0;=1 1 0;nyquist(n,dnyquist(n,d); );figure(2figure(2) )w=0.5:0.2:3w=0.5:0.2:3; ;nyquis
30、t(n,d,wnyquist(n,d,w); );axisaxis(equal);(equal);有时为了使曲线变化趋势看起来更明显，常有时为了使曲线变化趋势看起来更明显，常采用人工变量。采用人工变量。返回n n控制系统的稳定性分析借助于MATLAB语言工具，控制系统的稳定性分析采用直接方法，成为异常方便，异常简单的问题。因此，控制理论中经常使用的几种间接分析方法，如劳斯判据、赫尔维斯判据和奈奎斯特判据等，这里从时域和频域两个方面介绍MATLAB语言的应用。uu时域稳定性分析uu频域稳定性分析返回时域稳定性分析在数学模型的基础上，可采用直接求根法确定系统的稳定性。在前述课程的学习过程中，可以知
31、道系统的稳定性取决于系统特征根的正负，这样要求取系统的根，在一般的求解过程中比较繁琐，而MATLAB恰恰提供了这方面的便利。系统的传递函数模型一般表示为系统的特征方程为对于系统的传递函数模型，特征方程的根全部具有负实部，即对于特征方程其n个根为则系统是稳定的。而MATLAB提供了两个关于时域稳定性应用的函数。函数roots(d)函数功能：根据对应于特征方程式多项式向量求取系统对应的特征根。函数pzmap(num,den) or pzmap(a,b,c,d)函数功能：绘出零极点在s平面上的位置，以图形可视方式显示系统的稳定性。例：控制系统结构图如图所示，试分别确定例：控制系统结构图如图所示，试分
32、别确定k=2k=2，k=10k=10时系统的稳定性。时系统的稳定性。dzdz=0,-1,-2;=0,-1,-2;do=poly(do=poly(dzdz); );no1=2;no1=2;nc1,dc1=cloop(no1,do);nc1,dc1=cloop(no1,do);roots(dc1)roots(dc1)ansans = = -2.5214 -2.5214 -0.2393 + 0.8579i -0.2393 + 0.8579i -0.2393 - 0.8579i -0.2393 - 0.8579i+-uy因此k=2时，由于系统的闭环特征根全部具有负实部，系统是稳定的。no2=10;nc
33、2,dc2=cloop(no2,do);roots(dc2)ans = -3.3089 0.1545 + 1.7316i 0.1545 - 1.7316i因此，当k=10时，由于系统有一对共轭复数根的实部为正值，系统不稳定。利用函数pzmap(num,den)也可判断稳定性。dz=0,-1,-2;do=poly(dz);no1=2;nc1,dc1=cloop(no1,do);subplot(2,1,1)pzmap(nc1,dc1)no2=10;nc2,dc2=cloop(no2,do);subplot(2,1,2)pzmap(nc2,dc2)返回频域稳定性分析在频域稳定性分析中，主要使用波特图
34、上的在频域稳定性分析中，主要使用波特图上的nyquistnyquist稳定性判据，即稳定裕度命令稳定性判据，即稳定裕度命令margin()margin()。注注意，稳定裕度的定义是对最小相位系统而言的。意，稳定裕度的定义是对最小相位系统而言的。若系统为非最小相位系统，不能应用稳定裕度来若系统为非最小相位系统，不能应用稳定裕度来判别系统的稳定性。判别系统的稳定性。例：系统的开环传递函数为例：系统的开环传递函数为试分析系统的稳定性。试分析系统的稳定性。n=1;n=1;d1=1,0;d2=0.5,1;d3=0.1,1;d1=1,0;d2=0.5,1;d3=0.1,1;d=conv(d1,conv(d
35、2,d3);d=conv(d1,conv(d2,d3);margin(n,d)margin(n,d)如上图所示，当k=1时系统是稳定的。由插值函数spline确定系统稳定的临界增益。m,p,w=bode(n,d);wi=spline(p,w,-180)wi = 4.4720mi=spline(w,m,wi)mi = 0.08331/mians = 12.0036通过上述运算，得到系统临界增益值为k临=12，则有k0dB，相位裕度Pm0，系统是稳定的；当k12时，系统的幅值裕度Gm0dB，相位裕度Pm0，系统不稳定。例：非最小相位系统的开环传递函数为试在频域确定系统的稳定性。k=1;n=0.5,
36、1;no=k*n;do=1,-1,0;margin(no,do);从图中可以看出，系统的闭环是不稳定的，这时就要借助于nyquist判据来判稳。作奈氏曲线分析：作奈氏曲线分析：w=0.5:0.1:10;w=0.5:0.1:10;nyquistnyquist(no,do,w);(no,do,w);axis(equal)axis(equal)由插值法寻找临界增益k临：m,p,w=bode(no,do);wi=spline(p,w,-180)mi=spline(w,m,wi)1/miwi = 1.4142mi = 0.4999ans = 2.0004 因此，临界增益k临=2。由于p=1，根据奈氏判据
37、，该系统稳定条件为：角度增量=p=。所以k2时，在奈氏判图上，角度增量=，满足稳定条件，系统稳定。k=1和k=3的奈氏判据图如上页图所示。绘制程序如下：hold onno=no*3;nyquist(no,do);返回三、控制系统校正在控制分析的基础上对系统性能进行改善。为了保证一定的控制效果，按照传统方法，在原系统特性的基础上，将原特性加以修正称为系统的校正。例如改变原系统的波特图使之成为希望的形状即满足要求的性能指标就属于控制系统的校正内容。在本节主要介绍频域校正法，即基于波特图的校正方法。返回基于波特图的系统校正方法主要有两种。一种是相位超前校正，通过超前校正装置的相位超前特性使校正系统获
38、得希望的相位裕度；另一种是相位滞后校正，通过压缩频带宽度使校正系统获得希望的相位裕度。这两种校正方法都是在系统串联校正的基础上实现的。例：已知角位移随动系统的开环特性为要求：(1)r(t)=t时，ess0.1弧度； (2) c 4.41/s，c45。用频率法设计超前校正装置。为满足稳态性能，令K=10，作开环系统的波特图（见下页）。no1=10;do1=1,1,0;bode(no1,do1);hold ongm,pm,wg,wp=margin(no1,do1);gm,pm,wg,wpans = Inf 17.9642 NaN 3.0842相位裕度c=17.9645，wp=2.16642，K=3
39、0，满足设计要求。绘制系统校正前后的时域响应。程序程序1 1：将以上频域波特图江于一图。：将以上频域波特图江于一图。n=30;d=conv(1 0,conv(0.1 1,0.2 1);bode(n,d)gm,pm,wg,wp=margin(n,d);gm,pm,wg,wphold onnc=4,1;dc=50,1;bode(nc,dc);no2=conv(n,nc);do2=conv(d,dc);bode(no2,do2);gm,pm,wg,wp=margin(no2,do2);gm,pm,wg,wp返回四、应用举例例1：控制系统的传递函数分别为：作系统的脉冲响应曲线，并求解析解。num1=0
40、 0 1;num2=0 1 0;num3=1 0 0;den=1 2 10;impulse(num1,den);hold on;impulse(num2,den);impulse(num3,den);各传递函数对应的解析解分别为：syms sa=ilaplace(1/(s2+2*s+10);simple(a)-1/6*i*exp(-t+3*i*t)+1/6*i*exp(-t-3*i*t)vv1/3*exp(-t)*sin(3*t)b=ilaplace(s/(s2+2*s+10);simple(b)exp(-t)*cos(3*t)-1/3*exp(-t)*sin(3*t)c=ilaplace(s
41、2/(s2+2*s+10);simple(c)Dirac(t)-2*exp(-t)*cos(3*t)-8/3*exp(-t)*sin(3*t)例2：系统无零点与系统有零点表达式为试比较系统阶跃响应的判别。num1=10;num2=2 10;den=1 2 10;step(num1,den)hold on;step(num2,den)阶跃响应曲线如下页图所示！显然，系统有显然，系统有s=-5s=-5的零点时，其阶跃响应叠加了的零点时，其阶跃响应叠加了微分项，响应时间减小，但是超调量稍有增加。微分项，响应时间减小，但是超调量稍有增加。计算超调量增加幅度：计算超调量增加幅度：num1=10;num2
42、=2 10;den=1 2 10;step(num1,den)hold on;step(num2,den)y1,x1=step(num1,den);y2,x2=step(num2,den);max(y1),max(y2)ansans = = 1.3509 1.4348 1.3509 1.4348很明显，系统的超很明显，系统的超调量调量增加了增加了8%8%。例3：控制系统的开环传递函数为作波特图并确定系统的稳定裕度。num=15;den=1 2 10 0;margin(num,den)从下页图上可以看出系统的开环截止频率为wc=2.236，相位裕度为Pm=48.19，相位裕度Pm较大。但是由于二
43、阶环节谐振峰值的影响，幅值裕度Gm较小，且Gm=2.499dB，因此系统的动态性能较差。例4：给定系统作根轨迹图，并确定临界增益的值。num=1;den=1 2 2 0;subplot(1,2,1)rlocus(num,den);subplot(1,2,2)rlocus(num,den);k,r=rlocfind(num,den)Select a point in the graphics windowselected_point = 0.0038 + 1.4258ik = 4.0814r = -2.0135 0.0067 + 1.4237i 0.0067 - 1.4237i通过对图形和运算结果的分析，可以知道系统的临界增益约为k=4.0814。返回
