高阶系统的时域分析
- 格式:doc
- 大小:345.50 KB
- 文档页数:18
课程设计任务书学生姓名: 专业班级: 自动化1002班 指导教师: 肖纯 工作单位: 自动化学院 题 目: 高阶系统的时域分析 初始条件:设单位系统的开环传递函数为))(105()()(2a s s s s b s K s G ++++=要求完成的主要任务: (包括课程设计工作量及其技术要求,以及说明书撰写等具体要求)(1) 当K=10,a=1,b=5时用劳斯判据判断系统的稳定性。
(2) 如稳定,则求取系统的单位阶跃响应、单位斜坡响应和单位加速度响应,用Matlab 绘制相应的曲线,并计算单位阶跃响应的动态性能指标和稳态性能指标,计算单位斜坡响应和单位加速度响应的稳态性能指标。
(3) 如不稳定,则计算系统稳定时K 、a 和b 的取值范围,在稳定范围内任取一值重复第2个要求。
(4) 绘制稳定时系统的根轨迹(在稳定范围内任取a 、b 值)。
分析K 变化对系统性能的影响。
时间安排:指导教师签名: 年 月 日系主任(或责任教师)签名: 年 月 日目录1 高阶系统的数学模型 (1)2 系统稳定性分析 (2)3 高阶系统的时域分析 (5)3.1 单位阶跃响应 (5)3.1.1 单位阶跃响应 (5)3.1.2 单位阶跃响应动态性能 (7)3.1.3 单位阶跃响应稳态性能 (8)3.2 单位斜坡响应 (9)3.2.1 单位斜坡响应 (9)3.2.2 单位斜坡响应稳态性能 (10)3.3 单位加速度响应 (11)3.3.1 单位加速度响应 (11)3.3.2 单位加速度响应稳态性能 (12)4 系统根轨迹 (13)5 设计心得体会 (14)参考文献 (14)高阶系统的时域分析1 高阶系统的数学模型一个高阶系统的闭环传递函数的一般形式为:10111011()(),()m m m mn n n nb s b s b s b C s s m n R s a s a s a a ----++++Φ==≤++++ 对分子、分母进行因式分解,得到零极点形式:11()()()()()mi i njj K s z C s s R s s p ==-Φ==-∏∏ (1)式(1)中,K=b 0/a 0;z i ,p j 分别为系统闭环零、极点。
本设计给定的单位反馈系统的开环传递函数为)a s )(10s 5s (s b)s ()s (2++++=K G P (2)则其闭环传递函数为(假设为负反馈):)3()10()510()5()()())(105()()(2342Kbs K a s a s a s b s K b s K a s s s s b s K s ++++++++=++++++=φ2 系统稳定性分析任何系统在扰动作用下都会偏离原平衡状态,产生初始偏差。
所谓稳定性, 是指系统在扰动消失之后,由初始偏差状态恢复到原平衡状态的性能。
线性系统 的稳定性仅取决于系统自身的固有特性,而与外界条件无关。
线性系统稳定的充 分必要条件是:闭环系统特征方程的所有根均具有负实部;或者说,闭环传递函 数的极点均位于 S 左半平面。
若求出闭环系统特征方程的所有根,就可判定系统的稳定性。
但对于高阶系统来说,求特征方程根很困难,并且不易对参数进行分析。
现使用一种不用求解特征根来判别系统稳定性的方法—劳斯稳定判据。
设系统的特征方程为10110()0,0n n n n D s a s a s a s a a --=++++=>,则可列出劳斯表如表1所示。
表1 劳斯表按照劳斯稳定判据,系统稳定的充分必要条件为:劳斯表中第一列各值均为正。
否则系统不稳定,且第一列各系数符号改变次数即为特征方程正实部根的数目。
当K=10,a=1,b=5时,代入式(3)得到系统闭环传递函数50201565010)(234+++++=s s s s s s φ则系统的闭环特征方程为:D(s)=s 4+6s 3+15s 2+20s+50=0. 按劳斯判据可列出如下劳斯表:由于劳斯表第一列数值符号有两次变化,故系统不稳定,且存在2个正实部根。
用 MATLAB 求出全部特征根如下:>> y=roots([1 6 15 20 50]) y =-3.1534 + 1.7836i-3.1534 - 1.7836i 0.1534 + 1.9458i 0.1534 - 1.9458i现继续用劳斯稳定判据求原给定系统稳定时K ,a ,b 的取值范围。
原给定系统的闭环特征方程为:D(s)=s 4+(5+a)s 3+(10+5a)s 2+(10a+K)s+Kb=0,按劳斯判据可列出如下劳斯表:根据劳斯稳定判据,令劳斯表中第一列各元素为正,即:⎪⎪⎪⎩⎪⎪⎪⎨⎧>>+-+++-++-++>++-++>+00)10()510)(5()5()10)](10()510)(5[(05)10()510)(5(052Kb K a a a a Kb K a K a a a a K a a a a (4) 即K 、a 和b 必须满足:⎪⎪⎩⎪⎪⎨⎧>>--+-++-++>-++>+002550)1015500()5250(50050255052232Kb K Kb K a Kb K a Kb K a K a a a (5) 系统才稳定。
3 高阶系统的时域分析取K=12,a=b=3时,此时系统由四阶变为三阶,系统开环传递函数为 )10s 5s (s 12)a s )(10s 5s (s b)s ()s (22++=++++=K G P (6)系统闭环传递函数为1210512)())(105()()(232+++=++++++=s s s b s K a s s s s b s K s φ (7) 经分析可知,此时 K 、a 、b 的值满足要求,系统稳定。
3.1 单位阶跃响应 3.1.1 单位阶跃响应单位阶跃响指的是系统在单位阶跃信号 r(t)=1(t)作用下的响应。
取其拉氏变换即 R(s)=1/s 。
此时,系统输出为:12105121)())(105()(1)()()(232+++=++++++==s s s s b s K a s s s s b s K s s s R s C φ对上式进行部分分式展开:425714.24286.035714.01)(2+++-+-=s s s s s s C 对部分分式进行拉普拉斯反变换,并设初始条件全部为零,得系统的单位阶跃响 应:)7321.1sin(2371.1)7321.1cos(4286.05714.01)]([)(31t e t e e s C L t c t t t -------==(8)对于一般的高阶系统来说,用这种方法来求取单位阶跃响应都比较麻烦,有时候甚至很难完成。
但利用 MATLAB 软件则可以很方便的得到响应,并绘制出响应曲线。
MATLAB 中tf2zp()函数能将传递函数模型转化为零极点模型,residue()函数可以直接求出传递函数部分分式展开,由这些结果可以直接写出系统的输出解析解。
另外,利用step()函数还能准确绘制系统单位阶跃响应曲线。
式(7)所表示系统可以用下面的MATLAB语句求解系统单位阶跃响应。
>> num=[12 36];den=[1 8 25 42 36] %描述闭环系统传递函数的分子、分母多项式sys=tf(num,den); %高阶系统建模[z,p,k]=tf2zp(num,den); %对传递函数进行因式分解zpk(z,p,k) %给出闭环传递函数的零极点形式[r,p,k]=residue(num,[den,0]) %对C(s)部分分式展开%在分母多项式后补零相当于乘以s step(sys) %绘制高阶系统的单位阶跃响应曲线grid %添加栅格title(‘单位阶跃响应’); %标注标题xlabel(‘t’); ylabel(‘c(t)’); %标注横、纵坐标轴绘制的单位阶跃响应曲线如图1所示。
图1 单位阶跃响应曲线由(8)式单位阶跃响应时域表达式可知系统闭环稳定时,单位阶跃响应的指数项和阻尼正弦余弦项均趋近于零,稳态输出为常数项1,这与用MATLAB绘制的响应曲线相符。
3.1.2 单位阶跃响应动态性能动态性能指标是指稳定的系统在单位阶跃函数作用下,动态过程随时间t的变化状况的指标,体现系统动态过程特征。
用解析法求解高阶系统的动态性能指标很困难,这里用MATLAB编程求解。
调用单位阶跃响应函数step(),获得系统的单位阶跃响应,当采用[y,t]=step(sys)的调用格式时,将返回值y及相应的时间t,通过对y和t进行计算,可以得到高阶系统各项动态性能指标。
利用MATLAB编程求取系统动态性能指标程序如下:sys=tf([12 36],[1 8 25 42 36])%系统建模%计算峰值时间tp和对应最大超调量MpC=dcgain(sys) %取系统终值[y,t]=step(sys); %求取单位阶跃响应,返回变量输出y和时间t[Y,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;endts=t(i)输出结果为: C =1 tp =2.2557 Mp = 0.1216 tr =1.6787 ts =3.2000即上升时间为 1.6787s ,峰值时间为 2.2557s ,最大超调量为12.16%,并且系统在3.2000s 后进入稳态。
3.1.3 单位阶跃响应稳态性能稳态性能是系统在典型输入作用下,当时间t 趋于无穷大时,系统输出量的最终复现输入量的程度。
稳态性能分析主要是指稳态误差的计算。
稳态误差是系统控制精度或抗干扰能力的一种度量。
现采用静态误差系数法计算单位阶跃响应稳态误差。
将K=12,a=3,b=3代入(2)式,得待分析系统的开环传递函数为:)105(12)(2++=s s s s G p 其静态位置误差系数为:∞=++==→→)105(12lim)(lim 2s s s s G K s p s p 所以单位阶跃输入作用下系统的稳态误差为:1()01ss pe K ∞==+3.2 单位斜坡响应3.2.1 单位斜坡响应单位斜坡输入21(),()r t t R s s==,此时 )42)(3(12)()()(22+++=Φ⋅=s s s s s s R s C 展开为部分分式:428571.06429.031905.08333.01)(22++++++-=s s s s s s s C对部分分式进行拉普拉斯反变换,并设初始条件全部为零,得系统的单位斜坡响应:)7321.1sin(1237.0)7321.1cos(6429.01905.08333.0)(3t e t e e t t c t t t ---++-= (9)利用 MATLAB 软件绘制该系统在单位斜坡响应曲线。