《实验二典型系统的时域响应分析实验仿真报告答案 (1)》由会员分享，可在线阅读，更多相关《实验二典型系统的时域响应分析实验仿真报告答案 (1)（4页珍藏版）》请在金锄头文库上搜索。
1、实验二典型系统的时域响应分析1. 实验目的1) 通过用MATLAB及SIMULINK对控制系统的时域分析有感性认识。2) 明确对于一阶系统，单位阶跃信号、单位斜坡信号以及单位脉冲信号的响应曲线图。3) 对于二阶系统阶跃信号的响应曲线图以及不同阻尼比、不同自然角频率取值范围的二阶系统曲线 比较图。4) 利用MATLAB软件来绘制高阶控制系统的零极点分布图，判断此系统是否有主导极点，能否用低阶 系统来近似，并将高阶系统与低阶系统的阶跃响应特性进行比较5) 编制简单的M文件程序。2. 实验仪器PC计算机一台，MATLAB软件1套3. 实验内容1) 一阶系统的响应(1) 一阶系统的单位阶跃响应在SIM
2、ULINK环境下搭建图1的模型，进行仿真，得出仿真曲线图。 理论分析：C (s) =1/s+1)由拉氏反变换得h(t)=1-e(-1/ (t=0) 由此得知，图形是一条单调上升的指数曲线，与理论分析相符。(2) 一阶系统的单位斜坡响应 在SIMULINK环境下搭建图2的模型，将示波器横轴终值修改为12进行仿真，得出仿真曲线图。 理论分析：C (s) =l/sT(4s+l)可求的一阶系统的单位斜坡响应为c(t)=(t-4)+4e%-t/4) e(t)=r(t)-c(t)=4-4e*t/4)当t=0时，e(t)=0,当趋于无穷时，误差趋于常值4.3) 一阶系统的单位脉冲响应在med it环境下，编
3、译一个.m文件，利用impulse ()函数可以得出仿真曲线图。此处注意分析在 SIMULINK 环境中可否得到该曲线图。理论分析：C (s) =5/+2)=(5/2)/+1)可求的g(t)i(-t/,是一个单调递减的函数。 两种环境下得到的曲线图不一致。2) 二阶系统的单位阶跃响应二阶系统的闭环传递函数标准形式为 其阶跃响应可以分以下情况解出当Q = 0时，系统阶跃响应为c(t) = l cos t)n当0 匚 1时，系统阶跃响应为c(t) = 1 -e-3r=sin( t + 9)Z 2 d其中9= tg-%1Q2/C ， = v 1C2dn当C = 1时，系统阶跃响应为c(t) = 1
4、(1 + w t)en当0 1时，系统阶跃响应为c(t) =1 =2 匚2 1入2其中九 1 二2 1 ,九2 二-匚 + ；匚 2 - 1自然角频率n二1二阶系统0 对系统响应的影响阻尼比:=0系统状态无阻尼状态欠阻尼状态临界阻尼状态过阻尼状态对系统响应 的影响系统的暂态响 应是恒定振幅 的周期函数系统的暂态响 应是振幅随时 间按指数规律 衰减的周期函 数，阻尼比越 大，振幅衰减 的越快系统的单位阶 跃响应随时间 的推移单调增 长，在时间趋 于无穷大时， 系统响应的最 大超调量为0暂态响应随时 间按指数规律 单调衰减。系 统无超调，但 过程缓慢。分析：当wn 一定时，越小，振荡越厉害，当增大
5、到1以后，曲线变为单调上升。阻尼比:=0.5选取不同自然角频率 =,，用MATLAB得到二阶系统阶跃响应曲线，并分析比较不同自然角频 n率对应的系统输出的情况。本题采用第三种，在SIMULINK环境下搭建图1的模型，进行仿真，二阶系统阶跃响应曲线。1.20.80.60.40.2268910%ime oifset: 0分析：当一定时，且处于欠阻尼状态时，wn越大，则系统达到稳定时，所需要的时间越短。（3）系统动态性能分析对于G=k必表示的二阶系统上升时间（S）峰值时间（S）最大超调量调整时间（s）曲线图12%公式计算12%解：wn=、；20 =25,二*5/4.可知系统处于欠阻尼状态，由课本上的
6、计算公式可得tr=, tp=,Mp=*100%,因为 0，所以 ts=.结论：通过比较得知，tp,Mp,ts，的理论值与图片中的值基本一致。3）高阶系统的单位阶跃响应已知高阶系统的闭环传递函数为用下式低阶系统近似原系统解: pl=-5,p2=+,p3=由于闭环极点与系统的原闭环极点传递函数之极点相同，零点则不同。对于高阶系统，极点均为负实数，而且无零点，则系统的暂态响应一定是非振荡的，响应主要取决于据虚 轴最近的极点，。若其他极点比最近极点的最大距离大5 倍以上，则可以忽略前者对系统暂态过程的影响。P1 距 p2 没有 5 倍以上，而 p3 和 p2 不能看成一对偶极子，由于 p4 和 p5 离原点很近，所以影响也不能忽 略。所以不能被低阶系统代替。（2）利用单位阶跃响应st ep（ ）、figure（）和hold on（）等函数和指令，在med it环境下，编译一个.m 文件，能够将原系统和降阶系统的单位阶跃响应绘制在一个图中，记录它们的响应曲线和暂态性能指标 （上升时间、峰值时间、超调量以及调整时间） ，进行比较分析。num=45;den=1,51,45;G=tf（num,den）;step（G）;figure（1）hold on num1=1; den1=1,1;G1=tf(num1,den1); step(G1);hold off
