高阶系统的时域分析课程设计

1、课程设计任务书题 目: 高阶系统的时域分析 初始条件：设单位系统的开环传递函数为要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）（1） 当K=10，a=1，b=4时用劳斯判据判断系统的稳定性。（2） 如稳定，则求取系统的单位阶跃响应、单位斜坡响应和单位加速度响应，用Matlab绘制相应的曲线，并计算单位阶跃响应的动态性能指标和稳态性能指标，计算单位斜坡响应和单位加速度响应的稳态性能指标。（3） 如不稳定，则计算系统稳定时K、a和b的取值范围，在稳定范围内任取一值重复第2个要求。（4） 绘制a=1，b=4时系统的根轨迹。时间安排：任 务时 间（天）指导老师下达任

2、务书，审题、查阅相关资料2分析、计算3编写程序2撰写报告2论文答辩1指导教师签名： 年 月 日系主任（或责任教师）签名： 年 月 日目 录1 高阶系统的数学模型12 系统稳定性分析13 高阶系统的时域分析33.1 单位阶跃响应43.1.1 求单位阶跃响应43.1.2 单位阶跃响应动态性能73.1.3 单位阶跃响应稳态性能93.2 单位斜坡响应103.2.1 求单位斜坡响应103.2.2 单位斜坡响应稳态性能113.3 单位加速度响应113.3.1 求单位加速度响应113.3.2 单位加速度响应稳态性能134 系统根轨迹135 设计心得体会15参考文献15高阶系统的时域分析1 高阶系统的数学模型

3、一个高阶系统的闭环传递函数的一般形式为：对分子、分母进行因式分解，得到零极点形式： (1)式(1)中，K=b0/a0；zi ,pj分别为系统闭环零、极点。本设计给定的单位反馈系统的开环传递函数为 (2)则其闭环传递函数为（假设为负反馈）： (3)2 系统稳定性分析线性系统稳定的充分必要条件为：闭环系统特征方程的所有根均具有负实部；或者说，闭环传递函数的极点均位于s左半平面。若求出闭环系统特征方程的所有根，就可判定系统的稳定性。但对于高阶系统来说，求特征方程根很困难，并且不易对参数进行分析。现使用一种不用求解特征根来判别系统稳定性的方法劳斯稳定判据。设系统的特征方程为，则可列出劳斯表如表1所示。

4、表1 劳斯表按照劳斯稳定判据，系统稳定的充分必要条件为：劳斯表中第一列各值均为正。否则系统不稳定，且第一列各系数符号改变次数即为特征方程正实部根的数目。当K=10，a=1，b=4时，代入式(3)得到系统闭环传递函数则系统的闭环特征方程为：D(s)=s4+5s3+12s2+18s+40=0. 按劳斯判据可列出如下劳斯表：由于劳斯表第一列数值符号有两次变化，故系统不稳定，且存在2个正实部根。现继续用劳斯稳定判据求原给定系统稳定时K，a，b的取值范围。原给定系统的闭环特征方程为：D(s)=s4+(4+a)s3+(8+4a)s2+(8a+K)s+Kb=0，按劳斯判据可列出如下劳斯表：根据劳斯稳定判据，

5、令劳斯表中第一列各元素为正，即：即K、a和b必须满足： (4)系统才稳定。3 高阶系统的时域分析取K=15，a=2，b=2时系统闭环传递函数 (5)分析，此时K、a、b的值满足不等式组(4)，系统稳定。3.1 单位阶跃响应 求单位阶跃响应单位阶跃输入r(t)=1(t)，R(s)=1/s.对于n（n3）阶系统先将系统闭环传递函数一般形式化成如(1)式所示零极点形式，则在单位阶跃输入作用下，系统输出可表示为（假设系统闭环极点均不相同）：将该式展开成部分分式的形式，响应可表示为式中，A0、Aj(j=1,2,q)、Bk和Ck(k=1,2,r)是由部分分式展开时获得的系数。对上式取拉普拉斯反变换得到系统

6、时域响应表达式：由上式可知，高阶系统的时域响应是由稳态值和一些惯性环节及振荡环节的瞬态响应分量组成。对于稳定系统，上式瞬态响应分量的指数衰减项和正弦衰减项均随响应时间t趋于无穷而趋于零，系统达到稳态值。各瞬态分量在过渡过程中所起作用的大小，将取决于它们的指数的值及相应系数项Aj、Bk、Ck的大小。在瞬态过程中，某衰减项的指数|pj|或的值越大，则该项衰减越快，反之亦然。而|pj|和就是系统的极点到虚轴的距离。因此，如果分布在s平面左半部分的极点离虚轴越远，则它对应的分量衰减越快。显然，对系统过渡过程影响最大的，是那些离虚轴最近的极点。各衰减项的系数不仅与相应的极点在s平面中的位置有关，而且还与

7、零点的位置有关。极点的位置距原点越远，则相应分量的系数越小，该分量对系统过渡过程的影响就越小。如果某极点与零点很靠近，则相应分量的系数也很小，这时零极点对系统过渡过程的影响也将很小。因此，高阶系统的瞬态响应特性主要由系统传递函数中那些靠近虚轴而又远离零点的极点来决定。如果高阶系统有一个极点（或一对共轭复数极点）离虚轴最近，且其附近又无零点存在，而其他所有极点与虚轴的距离都在此极点与虚轴的距离的五倍以上，则可近似地认为系统的瞬态响应特性由这个（或这对）极点来确定，而其它极点的影响可以忽略不计，这个（或这对）极点就称为高阶系统的主导极点。下面以选取的系统进行分析，其闭环传递函数为上面的式(5)。在

8、单位阶跃输入作用下，系统输出为对上式进行部分分式展开：对部分分式进行拉普拉斯反变换，并设初始条件全部为零，得系统的单位阶跃响应： (6)对于高阶系统，用上述解析法求解系统单位阶跃响应比较复杂，若借助MATLAB软件将十分简单。MATLAB中tf2zp()函数能将传递函数模型转化为零极点模型，residue()函数可以直接求出传递函数部分分式展开，由这些结果可以直接写出系统的输出解析解。另外，利用step()函数还能准确绘制系统单位阶跃响应曲线。式(5)所表示系统可以用下面的MATLAB语句求解系统单位阶跃响应。num=15 30;den=1 6 16 31 30; %描述闭环系统传递函数的分子

9、、分母多项式sys=tf(num,den); %高阶系统建模z,p,k=tf2zp(num,den);%对传递函数进行因式分解zpk(z,p,k) %给出闭环传递函数的零极点形式r,p,k=residue(num,den,0) %对C(s)部分分式展开 %在分母多项式后补零相当于乘以sstep(sys) %绘制高阶系统的单位阶跃响应曲线grid %添加栅格title(单位阶跃响应); %标注标题xlabel(t); ylabel(c(t); %标注横、纵坐标轴程序运行后得到系统零极点形式、部分分式展开式，这里不列出。绘制的单位阶跃响应曲线如图1所示。 图1 单位阶跃响应曲线由(6)式单位阶跃响

10、应时域表达式可知系统闭环稳定时，单位阶跃响应的指数项和阻尼正弦余弦项均趋近于零，稳态输出为常数项1，这与用MATLAB绘制的响应曲线相符。现将(6)式中三个瞬态分量曲线用MATLAB软件画出，如图2所示。其中曲线1为瞬态分量，曲线2为分量，曲线3为分量.由比较曲线可以看到，各分量的衰减速率和初始值都与相应的极点到虚轴的距离密切相关。与e-3t项相比，e-0.5t项具有慢得多的衰减速率。因此，对于除了t趋近于零以外的所有时间，e-3t项在合成的时域响应中的贡献可以忽略不计。所以可以说，e-0.5t项在响应中起着主导作用，相应地，s = -0.5±2.1794i是该系统的主导极点。图2

11、单位阶跃响应瞬态分量比较 单位阶跃响应动态性能动态性能指标是指稳定的系统在单位阶跃函数作用下，动态过程随时间t的变化状况的指标，体现系统动态过程特征。用解析法求解高阶系统的动态性能指标很困难，这里用MATLAB编程求解。调用单位阶跃响应函数step()，获得系统的单位阶跃响应，当采用y,t=step(sys)的调用格式时，将返回值y及相应的时间t，通过对y和t进行计算，可以得到高阶系统各项动态性能指标。利用MATLAB编程求取系统动态性能指标流程图如图3所示。图3 求取系统动态性能程序流程图利用MATLAB编程求取系统动态性能指标程序如下：sys=tf(15 30,1 6 16 31 30);

12、 %系统建模 %计算峰值时间tp和对应最大超调量MpC=dcgain(sys) %取系统终值y,t=step(sys); %求取单位阶跃响应，返回变量输出y和时间tY,k=max(y); %求输出响应的最大值Y（即峰值）和位置ktp=t(k) %取峰值时间Mp=(Y-C)/C %计算最大超调量 %计算上升时间trn=1;while y(n)<C %循环求取第一次到达终值时的时间 n=n+1;endtr=t(n) %计算调节时间(误差带取2%)i=length(t); %求取仿真时间t序列的长度while(y(i)>0.98*C)&(y(i)<1.02*C) i=i-1

13、;endts=t(i)程序运行后，输出结果为：C=1tp=1.8265Mp=0.3685tr=1.2645ts=7.5868即上升时间为1.2645s，峰值时间为1.8265s，最大超调量为36.85%，并且系统在7.5868s后进入稳态。 单位阶跃响应稳态性能稳态性能是系统在典型输入作用下，当时间t趋于无穷大时，系统输出量的最终复现输入量的程度。稳态性能分析主要是指稳态误差的计算。稳态误差是系统控制精度或抗干扰能力的一种度量。现采用静态误差系数法计算单位阶跃响应稳态误差。将K=15，a=2，b=2代入（2）式，得待分析系统的开环传递函数为，其静态位置误差系数为：所以单位阶跃输入作用下系统的稳

14、态误差为：3.2 单位斜坡响应 求单位斜坡响应单位斜坡输入，此时展开为部分分式：对部分分式进行拉普拉斯反变换，并设初始条件全部为零，得系统的单位斜坡响应： (7)用MATLAB绘制系统单位斜坡响应曲线使用lsim()函数，lsim()可以绘制线性定常系统在任意输入信号作用下的时间响应曲线，程序代码如下：sys=tf(15, conv(1 3,1 1 5); %系统建模t=0:0.01:10; %响应时间u=t; %单位斜坡输入lsim(sys,u,t) %单位斜坡响应gridxlabel(t); ylabel(c(t) %标注横、纵坐标轴title(单位斜坡响应); %标注标题程序运行后得到系

15、统单位斜坡响应曲线如图4所示。图4 单位斜坡响应曲线由(7)式单位斜坡响应时域表达式分析可知，本系统的单位斜坡响应的稳态分量为(t-0.5333)，系统稳态输出速度恰好与单位斜坡输入速度相同，即系统能跟踪斜坡输入，在位置上存在稳态跟踪误差，这与图4所示曲线相符合。 单位斜坡响应稳态性能待分析系统的开环传递函数为，其静态速度误差系数为所以系统在单位斜坡输入作用下的稳态误差为3.3 单位加速度响应 求单位加速度响应单位加速度输入，此时展开为部分分式形式：对部分分式进行拉普拉斯反变换，并设初始条件全部为零，得系统的单位加速度响应：(8)下面用MATLAB绘制系统单位加速度响应曲线，仍然使用lsim(

16、)函数。在MATLAB工作空间中输入如下程序代码：num=15; den=conv(1 3,1 1 5); sys=tf(num,den); %系统建模t=0:0.01:10; %响应时间序列u=0.5*t.2; %单位加速度输入lsim(sys,u,t) %绘制单位加速度响应曲线gridxlabel('t'); ylabel('c(t)');title('单位加速度响应');程序运行后，得到系统单位加速度响应曲线如图5所示。由(8)式单位加速度响应时域表达式分析可知，系统单位加速度响应的稳态输出为()，稳定时系统不能跟踪加速度输入，随响应时间t

17、的增大，稳态位置误差将越来越大，从图5所示单位加速度响应曲线也可以看出。图5 单位加速度响应曲线 单位加速度响应稳态性能待分析系统的开环传递函数为，其静态速度误差系数为：所以在单位加速度输入作用下的系统稳态误差为：4 系统根轨迹绘制a=1，b=4时系统的根轨迹图，将a、b的值代入式(2)，得到系统开环传递函数为 (9)MATLAB中提供了rlocus()函数，可以直接用于绘制开环系统的闭环根轨迹。先在MATLAB中输入系统开环传递函数，然后调用rlocus()函数就可以绘制出精确的根轨迹曲线，具体MATLAB程序代码如下：num=1 4; den=conv(1 0,conv(1 4 8,1 1

18、);sys=tf(num,den);rlocus(sys) %绘制根轨迹图title(根轨迹图)xlabel(实轴); ylabel(虚轴);程序运行后，得到系统根轨迹图如图6所示。图6 系统根轨迹图单击根轨迹上的点，则可以显示出该点处的增益值和其他相关信息。例如，若单击根轨迹和虚轴相交的点，则可以得出该点处增益的临界值为6.64，如图6所示。可以看出，当0<K<6.64时，系统闭环稳定。5 设计心得体会对高阶系统进行时域分析，运用经典解析方法，采用拉普拉斯反变换求解瞬态响应时域表达式比较复杂，要计算出各项动态性能指标也很困难。但对于许多高阶系统，利用主导极点法可以简化系统的分析和性能指标的估算。而利用MATLAB软件可以方便地对高阶系统时域响应进行准确分析。通过本次课程设计，加深了对所学自动控制原理课程知识的理解，特别是系统稳定性分析，系统各项动态性能指标，稳态误差以及系统根轨迹等相关知识的理解。设计时借助MATLAB软件进行控制系统分析，进一步熟悉了MATL