LTI系统的建模,设计,仿真毕业设计.docVIP

  • 5
  • 0
  • 约5.72千字
  • 约 28页
  • 2017-06-08 发布于陕西
  • 举报
LTI系统的建模、设计和仿真 一、系统分析 1. 绘制三轴搬运系统开环根轨迹图、Bode 图和奈奎斯特图,并判断稳定性; 2. 如要求系统的阶跃响应的超调量小于16%,在此条件下当控制器为比例控制器,确定K的合适取值,并估算系统的调节时间(按2%准则); 3.当控制器为超前校正网络,设计合适的超前校正网络,使得系统的调节时间小于3s,阶跃响应的超调量小于16%;(要求用根轨迹法和频率响应法两种方法设计控制器) 4.如要求系统阶跃响应的超调量小于13%,斜坡响应的稳态误差小于0.125, 请设计一个能满足要求的滞后校正网络,并估算校正后系统的调节时间(按2%准则); 5.如将校正网络改为特定形式的PI控制器,请设计能满足要求的PI控制器; 6.为改善系统的瞬态性能,可以考虑引入前置滤波器Gp。请讨论当控制器为超前控制和滞后控制两种情况下,前置滤波器对控制系统性能的影响; 7.如希望尽可能的缩短系统的调节时间,请设计相对应的控制器,并画出校正后系统的阶跃响应曲线。 二、概要设计1. 在输出端按要求能准确复现给定信号; 2. 具有良好的相对稳定性; 3. 对扰动信号具有充分的抑制能力。 为使系统的控制满足预定的性能指标而有目的地增添的元件或装置,称为控制系统的校正装置。 在控制系统中,我们经常采用两种校正方案:串联校正和反馈校正,其结构如图4.1、图4.2所示。而校正方案的选择主要取决于系统的结构、系统各点功率的大小、可供选择的元件、设计者的经验及经济条件等因素。、超调量、上升时间和调整时间等,这时采用根轨迹法进行设计是比较方便的。 若性能指标以频域量的形式给出时,例如所期望的相角裕度、幅值裕度、谐振峰值、剪切频率、谐振频率、带宽()及反映稳态指标的开环增益、稳态误差或误差系数等,此时虽然时域响应和频域响应之间有间接的关系,但是采用Bode图法进行设计比较方便。采用Bode法设计时,如果给定的是时域指标、,则必须先将其换算成频域指标和,方可进行相应的运算和设计。 三、详细设计Bode 图和奈奎斯特图,并判断稳定性; 开环 Gs=tf([1],[1,6,5,0]) rlocus(Gs) bode(Gs) nyquist(Gs) 闭环 Gf=feedback(Gs,1) rlocus(Gf) bode(Gf) nyquist(Gf) 开环传递函数的根轨迹图 开环传递函数的BODE图 开环传递函数的Nyquist图 闭环传递函数的根轨迹图 闭环传递函数的BODE图 闭环传递函数的Nyquist图 自动控制原理中对系统稳定性的判定,由于该系统传递函数的闭环根轨迹的极点均分布在负半平面,故系统稳定。 。 2: 如要求系统的阶跃响应的超调量小于16%,在此条件下当控制器为比例控制器,确定K的合适取值,并估算系统的调节时间(按2%准则); 解决方案:此题利用MATLAB的SISOTOOL工具能很方便的得出Gc的参数。SISOTOOL是典型单入单出系统的分析工具,只需将开环传递函数Gs导入,就能方便的对它的校正装置和反馈装置进行调节并直观的看出它的各项性能。 源代码: sisotool(Gs) Sisotool界面 经过测试,当Kc=3时能满足题设的条件,超调量为11.7%,调节时间(Settling Time)为7.88 s。 3:当控制器为超前校正网络,设计合适的超前校正网络,使得系统的调节时间小于3s,阶跃响应的超调量小于16%;(要求用根轨迹法和频率响应法两种方法设计控制器) 解决方案: 根轨迹法: 源代码: sita=0.16; ts=3; delta=0.02; pesi=(1-1/(1+(1/pi)*log10(1/sita)))^(1/2); wn=4.4/(ts*pesi); (此处略去了在草稿纸上的画图分析,只是用MATLAB进行了部分参数的计算) 算出零极点后利用Sisotool观察是否满足条件:sisotool(Gs); Sisotool运行界面 根轨迹法—超前校正网络的系统阶跃响应图 分析:由图可知,超调量为2.59%,调节时间为1.84 s,均符合题设要求。 (2)频率响应法 源代码: K=8;Ts=3; Gs=tf([1],[1,6,5,0]); [mag,phase,w]=bode(Gs); [Gm,Pm]=margin(Gs); Gama=45; Mr=1/sin(Gama); gama=Gama*pi/180; alfa=(1-sin(gama))/(1+sin(gama)); %adb=20*log10(mag); %am=10*log10(alfa); %wc=spline(adb,w,am); wc=(2+1.5*(Mr-1)+2.5*(Mr-1)^2)*pi/Ts; T=1/(

文档评论(0)

1亿VIP精品文档

相关文档