LTI系统的建模 - 设计 - 仿真
更新时间:2024-04-16 01:57:01 阅读量: 综合文库 文档下载
- Lti系统推荐度:
- 相关推荐
11001421 邱欢
LTI系统的建模、设计和仿真
一、系统分析
考虑如图所示的三轴搬运系统,要求机械臂能在三维空间中准确定位。
R(s)Gp(s)?G(s)Gcc(s)1s(s?1)(s?5)C(s)
1. 绘制三轴搬运系统开环根轨迹图、Bode 图和奈奎斯特图,并判断稳定性; 2. 如要求系统的阶跃响应的超调量小于16%,在此条件下当控制器Gc为比例控制器,确定K的合适取值,并估算系统的调节时间(按2%准则);
3.当控制器Gc为超前校正网络,设计合适的超前校正网络,使得系统的调节时间小于3s,阶跃响应的超调量小于16%;(要求用根轨迹法和频率响应法两种方法设计控制器)
4.如要求系统阶跃响应的超调量小于13%,斜坡响应的稳态误差小于0.125, 请设计一个能满足要求的滞后校正网络,并估算校正后系统的调节时间(按2%准则); 5.如将校正网络改为特定形式的PI控制器,请设计能满足要求的PI控制器; 6.为改善系统的瞬态性能,可以考虑引入前置滤波器Gp。请讨论当控制器为超前控制和滞后控制两种情况下,前置滤波器对控制系统性能的影响;
7.如希望尽可能的缩短系统的调节时间,请设计相对应的控制器,并画出校正后系统的阶跃响应曲线。
二、概要设计
控制算法设计思想概述:
单位负反馈控制系统具有良好的性能是指:
1. 在输出端按要求能准确复现给定信号; 2. 具有良好的相对稳定性;
3. 对扰动信号具有充分的抑制能力。
为使系统的控制满足预定的性能指标而有目的地增添的元件或装置,称为控制系统的校正装置。
在控制系统中,我们经常采用两种校正方案:串联校正和反馈校正,其结构如图4.1、图4.2所示。而校正方案的选择主要取决于系统的结构、系统各点功率的大小、可供选择的元件、设计者的经验及经济条件等因素。
在控制工程实践中,解决系统的校正问题,采用的设计方法一般依性能指标而定。若单位反馈控制系统的性能指标以时域量的形式给出时,例如期望的闭环主导极点的阻尼比ζ和无阻尼自振频率?n、超调量?00、上升时间tr和调整时间ts等,这时采用根轨迹法进行设计是比较方便的。 若性能指标以频域量的形式给出时,例如所期望的相角裕度?、幅值裕度Kg、谐振峰值Mr、剪切频率?c、谐振频率?r、带宽(0~?b)及反映稳态指标的开环增益K、稳态误差ess(?)或误差系数Ci(i?0,1,2)等,此时虽然时域响应和频域响应之间有间接的关系,但是采用Bode图法进行设计比较方便。采用Bode法设计时,如果给定的是时域指标?00、ts,则必须先将其换算成频域指标?和?c,方可进行相应的运算和设计。
三、详细设计与源代码,仿真结果与性能分析
1: 绘制三轴搬运系统开环根轨迹图、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%,在此条件下当控制器Gc为比例控制器,确定K的合适取值,并估算系统的调节时间(按2%准则);
解决方案:此题利用MATLAB的SISOTOOL工具能很方便的得出Gc的参数。SISOTOOL是典型单入单出系统的分析工具,只需将开环传递函数Gs导入,就能方便的对它的校正装置和反馈装置进行调节并直观的看出它的各项性能。 源代码: sisotool(Gs)
Sisotool界面
经过测试,当Kc=3时能满足题设的条件,超调量为11.7%,调节时间(Settling Time)为7.88 s。
3:当控制器Gc为超前校正网络,设计合适的超前校正网络,使得系统的调节时间小于3s,阶跃响应的超调量小于16%;(要求用根轨迹法和频率响应法两种方法设计控制器) 解决方案: (1) 根轨迹法:
源代码: sita=0.16; ts=3; delta=0.02;
pesi=(1-1/(1+(1/pi)*log10(1/sita)))^(1/2); wn=4.4/(ts*pesi);
算出零极点后利用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)); -b=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/(wc*sqrt(alfa)); alfat=alfa*T;
Gc=tf([T 1],[alfat 1]); figure(1) margin(K*Gc*Gs); figure(2)
step(feedback(K*Gc*Gs,1)); 运行结果:
结论:超调量:8.6% 调节时间:2.63 s 均满足题设要求。
4:如要求系统阶跃响应的超调量小于13%,斜坡响应的稳态误差小于0.125, 请设计一个能满足要求的滞后校正网络,并估算校正后系统的调节时间(按2%准则);
采用频率响应法进行设计 源代码: K=8;
Gs=tf([1],[1,6,5,0]); P0=57;
fic=-180+P0+5;
[mu,pu,w]=bode(Gs); wc2=spline(pu,w,fic);
d1=conv(conv([1 0],[1 1]),[1 5]); na=polyval(K,j*wc2); da=polyval(d1,j*wc2); G=na/da; g1=abs(G); L=20*log10(g1); alfa=10^(L/20); T=1/(0.1*wc2); alfat=alfa*T;
Gc=tf([T 1],[alfat 1]); figure(1) margin(K*Gc*Gs); figure(2)
step(feedback(K*Gc*Gs,1)); 运行结果
结论:超调量:12.9 调节时间:30.4 s。。
5:如将校正网络改为特定形式的PI控制器,请设计能满足要求的PI控制器 解决方案:
分别观察Kp,Ki,Kd对系统的影响: (1)观察Kp对系统阶跃响应的影响:
源代码: for Kp=3:1:5
Gp=feedback(Kp*Gs,1) step(Gp),hold on end
axis([0,10,0,1.5]); gtext('Kp=3'); gtext('Kp=4'); gtext('Kp=5'); 运行结果:
结论:Kp对超调量和上升时间有影响,Kp越大,超调量越大,上升时间越快。 (2)观察Ki对系统阶跃响应的影响: 源代码:
Gs=tf([1],[1,6,5,0]); Kp=3;
for Ki=0.3:0.1:0.5 Go=tf(Kp*[Ki 1],[Ki 0]); Gpi=feedback(Go*Gs,1); step(Gpi),hold on end
axis([0,10,-5,5]); gtext('Ki=0.3'); gtext('Ki=0.4'); gtext('Ki=0.5'); 运行结果:
结论:Ki对系统的稳态误差有影响,Ki越小,稳态误差越大
(3)观察Kd对系统阶跃响应的影响: 源代码:
>> Gs=tf([1],[1,6,5,0]); Kp=3;Ki=10 for Kd=1:1:3
Go=tf(Kp*[Ki*Kd Ki 1],[Ki 0]); Gpid=feedback(Go*Gs,1); step(Gpid),hold on end
gtext('Kd=1'); gtext('Kd=2'); gtext('Kd=3'); 运行结果:
结论:Kd对系统的上升时间和超调量有影响,Kd越大,上升时间越快,超调量越小。
(4)设计PI控制器: 源代码: s=tf('s');
Gs=tf([1],[1,6,5,0])
[Kc,Pm,Wc]=margin(Gs); %计算频域响应参数,增益裕量kc和剪切频率wc Tc=2*pi/Wc;
PIKp=0.4*Kc; %频域响应整定法计算并显示PI控制器 PITi=0.8*Tc;
PIGc=PIKp*(1+1/(PITi*s))
step(feedback(PIGc*Gs,1)),hold on [PIKp,PITi] gtext('PI'); 运行结果:
使用试探法来满足要求: Kp=2;Ki=12
Go=tf(Kp*[Ki 1],[Ki 0]); Gpi=feedback(Go*Gs,1); step(Gpi);
当Kp=2;Ki=12时,系统阶越响应比较满意
6:为改善系统的瞬态性能,可以考虑引入前置滤波器。请讨论当控制器为超前控制和滞后控制两种情况下,前置滤波器对控制系统性能的影响 (1) 当控制器为超前控制情况下,前置滤波器对控制系统性能的影响:
Gc=tf([1 0.9],[1 18.4]); G=tf([1],[1 6 5 0]); G1=series(308*Gc,G); Gp=tf([0.9],[1 0.9]); Gq=feedback(G1,1); Gh=series(Gq,Gp); step(Gq) hold on step(Gh)
引入滤波器前系统的超调量13.8%,调节时间为2.14s,引入滤波器后系统的超调量为0%,调节时间为4.74s,改善了系统的动态性能。
(2)当控制器为滞后控制情况下,前置滤波器对控制系统性能的影响: Gc=tf([1 0.01],[1 0.0002]); G=tf([1],[1 6 5 0]); G1=series(3*Gc,G);
Gp=tf([0.01],[1 0.01]);Gq=feedback(G1,1); Gh=series(Gq,Gp); step(Gq) hold on step(Gh)
引入滤波器前系统的超调量13.6%,调节时间为8.48s,引入滤波器后系统的超调量为0%,调节时间为386s,改善了系统的动态性能。
7:如希望尽可能的缩短系统的调节时间,请设计相对应的控制器,并画出校正后系统的阶跃响应曲线。 PID调节: s=tf('s'); Kp=2;Ki=8;Kd=1;
Go=Kp*(1+1/(Ki*s)+Kd*s/((Kd/10)*s+1)); Gpid=feedback(Go*Gs,1); step(Gpid); (1)调节P至临界震荡 >> Kp=30;Ki=0;Kd=0; Go=Kp
Gpid=feedback(Go*Gs,1); step(Gpid);
(2)加上I,D项根据经验公式联调 s=tf('s'); Kp=12;Ki=18;Kd=1;
Go=Kp*(1+1/(Ki*s)+Kd*s/((Kd/10)*s+1)); Gpid=feedback(Go*Gs,1); step(Gpid); 运行结果:
分析:PID所能达到的最佳调节时间仍比较长 采用超前校正:
修改系数,调整零点到0.87,极点调整到10 则得到增益K为10*10=100 代码:
numc=[1 0.87]; denc=[1 10];
gc=tf(numc*100,denc);%gc是校正网络
num=[1];
den=[1 6 5 0];
g=tf(num,den);%原系统 gs=g*gc %串联后系统 sisotool(gs) 整个系统为:
Transfer function: 100 s + 87
---------------------------- s^4 + 16 s^3 + 65 s^2 + 50 s
此时系统的调节时间是1.81秒,超调量为4.24%。
可以修正一下增益,用sisotool调节,则最佳增益是K=90 此时,调节时间为1.25秒,超调量为1.85%。
提高系统的时域指标,取极点在40 则得到增益K为40*10=400 numc=[1 0.87]; denc=[1 40];
gc=tf(numc*400,denc);%gc是校正网络 num=[1];
den=[1 6 5 0];
g=tf(num,den);%原系统 gs=g*gc %串联后系统 sisotool(gs)
整个系统为:
Transfer function: 400 s + 348
------------------------------ s^4 + 46 s^3 + 245 s^2 + 200 s
则此时没有超调量,系统的调节时间为1.29秒,时域指标较好。 修改增益,K为477 新的系统为:
Transfer function: 477 s + 415
------------------------------ s^4 + 46 s^3 + 245 s^2 + 200 s
没有超调量,调节时间为0.981秒,调节时间为1秒左右。
四、小结:
(1) 对自动控制原理和控制系统的设计方法有了更加深刻的认识; (2) 掌握了MATLAB在控制系统设计中的使用方法;
(3) 对各类控制器的设计方法(根轨迹法、频域响应法)有了全面的认识及理
解并能利用其解决实际问题;
(4) 对超前校正,滞后校正以及PID控制有了系统全面的认识;
五、参考文献:
(1) 《Matlab编程》---------------------- Stephen J. Chapman (科学出版社) (2) 《自动控制原理》--------------------胡寿松
(科学出版社)
没有超调量,调节时间为0.981秒,调节时间为1秒左右。
四、小结:
(1) 对自动控制原理和控制系统的设计方法有了更加深刻的认识; (2) 掌握了MATLAB在控制系统设计中的使用方法;
(3) 对各类控制器的设计方法(根轨迹法、频域响应法)有了全面的认识及理
解并能利用其解决实际问题;
(4) 对超前校正,滞后校正以及PID控制有了系统全面的认识;
五、参考文献:
(1) 《Matlab编程》---------------------- Stephen J. Chapman (科学出版社) (2) 《自动控制原理》--------------------胡寿松
(科学出版社)
正在阅读:
LTI系统的建模 - 设计 - 仿真04-16
苏州珠江钢琴|钢琴1-10级考级曲目大全,这下知道该练哪首了04-06
沪科版学第一学期七年级生物教学计划精修订04-24
君合律师事务所习题 - 图文04-02
太极拳期末考试试题样本04-17
黄华华省长闭门会议主持词(定稿)01-12
挖笋记作文700字06-19
护理部二季度护理质量考核情况总结06-06
- 多层物业服务方案
- (审判实务)习惯法与少数民族地区民间纠纷解决问题(孙 潋)
- 人教版新课标六年级下册语文全册教案
- 词语打卡
- photoshop实习报告
- 钢结构设计原理综合测试2
- 2014年期末练习题
- 高中数学中的逆向思维解题方法探讨
- 名师原创 全国通用2014-2015学年高二寒假作业 政治(一)Word版
- 北航《建筑结构检测鉴定与加固》在线作业三
- XX县卫生监督所工程建设项目可行性研究报告
- 小学四年级观察作文经典评语
- 浅谈110KV变电站电气一次设计-程泉焱(1)
- 安全员考试题库
- 国家电网公司变电运维管理规定(试行)
- 义务教育课程标准稿征求意见提纲
- 教学秘书面试技巧
- 钢结构工程施工组织设计
- 水利工程概论论文
- 09届九年级数学第四次模拟试卷
- 建模
- 仿真
- 设计
- 系统
- LTI
- 食品安全问题引发的企业诚信与道德的思考
- 郑淑媛申报省级教学团队 - 图文
- 第六期全国医师定期考核管理干部业务培训合格结业名单
- 3300多回路组合开关使用说明手册
- 检验检测机构资质认定管理办法条文释义
- 习作教学(较好)
- 浙江省温州中学2017届高三3月高考模拟数学试题
- 2009年高考数学压轴题系列训练含答案及解析详解
- 广东省揭阳一中2014届高三上学期开学摸底考试数学理试题
- 浅析小学英语有效课堂活动的设计与组织-最新教育资料
- 高效磨套铣工具的研制
- PACE 毛细管电泳操作手册-Operational Manual for MDQ with 32 K
- 应用Matlab对图像信号进行频谱分析及滤波
- 杜甫
- 2012年全国各地中考数学解析汇编 第二十一章 勾股定理(按章节考
- 易地(私庄)扶贫开发项目搬迁农户申请表1
- 学校薪酬管理制度范本-word范文(6页)
- 基于PLC的自动售货机毕业论文 - 图文
- 船舶压载水处理技术进展现状
- 有机课堂练习