2022年自动化自动控制课程设计报告 .pdf

《2022年自动化自动控制课程设计报告 .pdf》由会员分享，可在线阅读，更多相关《2022年自动化自动控制课程设计报告 .pdf（30页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、动控制课程设计报告班级：自动化 08-1 班学号：08051116 名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 30 页 -姓名：刘加伟2011.7.17 任务一、双容水箱的建模、仿真模拟、控制系统设计一、控制系统设计任务1、通过测量实际装置的尺寸，采集DCS 系统的数据建立二阶水箱液位对象模型。（先建立机理模型，并在某工作点进行线性化，求传递函数）2、根据建立二阶水箱液位对象模型，在计算机自动控制实验箱上利用电阻、电容、放大器的元件模拟二阶水箱液位对象。3、通过 NI USB-6008 数据采集卡采集模拟对象的数据，测试被控对象的开环特性，验证模拟对象的正确性。4、采用纯比

2、例控制，分析闭环控制系统随比例系数变化控制性能指标（超调量，上升时间，调节时间，稳态误差等）的变化。LC3#4#名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 30 页 -5、采用PI 控制器，利用根轨迹法判断系统的稳定性，使用Matlab 中SISOTOOLS 设计控制系统性能指标，并将控制器应用于实际模拟仿真系统，观测实际系统能否达到设计的性能指标。6、采用 PID 控制，分析不同参数下，控制系统的调节效果。7、通过串联超前滞后环节校正系统，使用Matlab 中 SISOTOOLS 设计控制系统性能指标，并将校正环节应用于实际模拟仿真系统，观测实际系统能否达到设计的性能指标。

3、(一)建立模型(二)实验模型及改变阶跃后曲线：1 取阶跃曲线按照以下模型建立系统辨识模型：一般取为 0.4 和 0.8 计算上行阶跃各参数：T1=171.26 T2=50.50 K=160.47 t1=141 t2=338 建立传递函数为：()(0)yyKu12t/Tt/T*121221TTy(t)1eeTTTT*(t)y12212121212TT(tt)2.16T T(TT)1.74tt0.55名师资料总结-精品资料欢迎下载-名师精心整理-第 3 页，共 30 页 -G(s)=)150.50)(126.171(60.1ss计算下行阶跃各参数：T1=84.20 T2=48.67 K=148.0

4、8 t1=89 t2=198 建立传递函数为：G(s)=)167.48)(120.84(48.1ss2 建立机理模型Q1=k1*u1；Q2=k2*u2*1H；Q=k3*u3*2H；k1=10；k2=1.9；k3=1.65；阀门开度u1=50；u2=52；u3=51；水箱面积 A1=1050；A2=600 理论传递函数G(s)=)16.84)(12.61(41.1ss；取辨识传递函数G（s）=)106.4)(102.7(48.1ss(三)根据建立的二阶水箱液位对象模型，在计算机自动控制实验箱上利用电阻、电容、放大器的元件模拟二阶水箱液位对象。根据传递函数可得：R2/R1=1.5 R2=300k(

5、200k+100k),R1=200k R2C=7.02 C=23.4uf R5/R4=1 R5=R4=200k R5C=4.06 C=20.3uF(四)通过 NI USB-6008 数据采集卡采集模拟对象的数据，测试被控对象的开环特性，验证模拟对象的正确性实验A/DNI USB-6008OPC Server数据采集服务器OPC ClientMatlab1/*1/)1)(1(25451212210sCRRRsCRRRsTsTKUUi名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 30 页 -系统的开环特性曲线为（实际曲线理论曲线）：(五)采用纯比例控制，分析闭环控制系统随比例系数变

6、化时控制性能指标改变 Kp 得到实际与理论曲线：Kp=1：KP=2：Kp=3：名师资料总结-精品资料欢迎下载-名师精心整理-第 5 页，共 30 页 -结论：由上图和系统指标分析可得，Kp越大，系统响应速度越快，上升时间越短，调节时间峰值时间也相应减少，且稳态误差减小，但超调量增大，系统振荡加剧，Kp过大时会对实际的系统造成破坏，所以，系统应选择合适的Kp。(六)采用 PI 控制器，利用根轨迹法判断系统的稳定性，使用Matlab 中 SISOTOOLS设计控制系统性能指标，并将控制器应用于实际模拟仿真系统，观测实际系统能否达到设计的性能指标Kp=2 Ti=5：Kp=2 Ti=4：名师资料总结-

7、精品资料欢迎下载-名师精心整理-第 6 页，共 30 页 -Kp=2 Ti=6 :由以上曲线可以看出，用根轨迹设计的设计的性能指标在实际系统中可以达到较好的效果，理论曲线与实际曲线较吻合。(七)采用 PID 控制，分析不同参数下，控制系统的调节效果Kp=2 Ti=10 Td=1 Kp=2 Ti=10 Td=2 Kp=2 Ti=9 Td=2 Kp=3 Ti=9 Td=2 名师资料总结-精品资料欢迎下载-名师精心整理-第 7 页，共 30 页 -由上面四个图对比分析可知:(1)Ti，Td 一定时，Kp 增大，加快系统的响应，系统的超调量增大，调节时间变小，上升时间减小，减小余差；(2)Kp，Ti

8、一定时，Td增大，系统的峰值时间减小，系统的超调量减小，振荡减小，调节时间减小。(3)Kp，Td一定时，Ti 增大，系统的超调量减小，减小振荡，使系统更加稳定，但余差消除的速度随之减慢。以上各曲线参数列表如下：KpTiTdTsTp%Tr199999017145.872999990141014.365399999012821.023.7240421049.74.82503110404.9260251032.65.12928.5133.26.439211.5105.44210114126.26.121029141.36.5(八)为被控对象设计串联校正环节，使用 Matlab 中 SISOTOOLS

9、设计控制系统性能指标，并将校正环节应用于实际模拟仿真系统，观测实际系统能否达到设计的性能指标（1）不加校正环节名师资料总结-精品资料欢迎下载-名师精心整理-第 8 页，共 30 页 -由图可知相角裕度为112deg，截至频率sec/123.0/radc。单位阶跃响应得上升时间sTr1.7，调节时间sTS4.20，峰值时间sTp2.14，超调量6.37%，稳态值为 0.597。（2）由于原系统的相角裕度较大，为了使校正效果明显，给原系统加入积分环节 1/s，于是)06.41)(02.71(48.11*108.115.2848.1)(2sssssssG其相角裕度为0sec/1.31 rad，截至频

10、率sec/338.0/radc,且系统已经发散，故要进行串联校正。（3）由第一个拐点处读得sec/142.002.7/1radb。假设调节时间sTs500，校正后的相角裕度45,则：414.1sin451sin1Mr049.3)1M(5.2)1M(5.12K2rr0由C0sKt可得0192.0500/14.3*049.3c名师资料总结-精品资料欢迎下载-名师精心整理-第 9 页，共 30 页 -由图36dB20lgT)(Lcbc，为了选择校正网络衰减因子/1，要保证已校正系统的截止频率为所选的c，则等式020lgT)(L20lgcbc成 立，因7.0 2Tb，则 可 得 到095.63108.

11、1（4）校正前：)/1)(/1()/1)(/1()(babasssssC)96.8/1)(/0957.631()142.0/1)(/1(ssssaa相应校正后：)96.8/1)(/0957.631)(06.41()/1(48.1)()(ssssssCsGaa于是，相角裕度：)06.4arctan(90)/arctan(180cac)96.8/arctan()/0957.63arctan(cac45解得：sec/0191.0rada所以：)/1)(/1()/1)(/1()(babasssssC)96.8/1)(0191.0/095.631()142.0/1)(0191.0/1(ssss其中019

12、1.0/095.6310191.0/1ss是装置的滞后部分，96.8/1142.0/1ss是校正装置的超前部分。超前校正主要是利用超前网络的相角超前特性，减小系统的截止频率，而滞后校正则是利用滞后网络的高频幅值衰减特性，加强系统的抗干扰能力。取合适的超调和调节时间，可取K=0.743，因而名师资料总结-精品资料欢迎下载-名师精心整理-第 10 页，共 30 页 -可得 C(s)96.8/1)(0191.0/095.631()142.0/1)(0191.0/1(*743.0ssss校正后的广义传递函数为：G(s)C(s)106.4)(102.7(48.1*)96.8/1)(0191.0/095.

13、631()142.0/1)(0191.0/1(*743.0sssssss加入串联校正相角裕度为 45.3deg，截止频率为 0.0227rad/sec。加入超前滞后环节后的阶跃响应，上升时间为 48.7s，超调量为 34%，峰值时间为 131s，调节时间为 404s，稳态值为 1。电路仿真曲线与理论曲线基本吻合，满足系统的要求。(九)改变 Kp 及 Ts 对系统的影响名师资料总结-精品资料欢迎下载-名师精心整理-第 11 页，共 30 页 -由上图可以看出，Kp越大，系统响应速度越快，上升时间越短，调节时间峰值时间也相应减少，且稳态误差减小，但超调量增大，系统振荡加剧，Kp 过大时会对实际的系

14、统造成破坏。由上图可以看出，Ts 增加使系统调节时间变长，超调量增大，调节精度下降，Ts 过大时使系统震动加剧，破坏系统的稳定性。因此Ts 应选取较小的值。(十)为被控对象设计最小拍无差控制器，并进行实验分析传递函数G(z)=20.001012 z 0.000981.9246.220 925zz零极点模型()G z=988.0686 z0.9719z0.9519z1z0.9744用 matlab 作出结构图：名师资料总结-精品资料欢迎下载-名师精心整理-第 12 页，共 30 页 -仿真波形为：由上图可以看出，系统在Ts=0.2s 时达到稳定，达到最小拍无差控制器的控制要求。名师资料总结-精品

15、资料欢迎下载-名师精心整理-第 13 页，共 30 页 -任务二、基于状态空间法单级倒立摆的控制系统设计一、单级倒立摆介绍倒立摆系统具有高阶次、不稳定、多变量、非线性和强耦合等特性，是控制理论的典型研究对象。如机器人行走过程中的平衡控制、火箭发射中垂直度控制和卫星飞行中的姿态控制等均涉及到倒置问题对倒立摆系统的研究在理论上和方法论上均有着深远意义。单级倒立摆系统的原理图，如图 1 所示。假设已知摆的长度为2l，质量为 m，用铰链安装在质量为M 的小车上。小车由一台直流电动机拖动，在水平方向对小车施加控制力u，相对参考系差生的位移s。若不给小车实施控制力，则倒置摆会向左或向右倾倒，因此，它是个不

16、稳定的系统。控制的目的是通过控制力u的变化，使小车在水平方向上运动，达到设定的位置，并将倒置摆保持在垂直位置上。建立单级倒立摆的状态空间数学模型。取状态变量Tssx。测试系统的开环特性。图是系统中小车和摆杆的受力分析图。其中，N 和 P 为小车与摆杆相互作用力的水平和垂直方向的分量。已知单级倒立摆的各项数据如下所示：,5.0,1.0,2mlkgmkgMgmgkgmI/8.9,025.02名师资料总结-精品资料欢迎下载-名师精心整理-第 14 页，共 30 页 -水平方向NusM；)sin(22lsdtdmN垂直方向：cos22ldtdmmgP根据力矩平衡方程：INlPlcossin；因为很小，

17、且 rank(Mc)名师资料总结-精品资料欢迎下载-名师精心整理-第 16 页，共 30 页 -即能控性 rank(Mc)得 4，矩阵 Mc的秩为 4，满秩，所以系统能控。能观性：能观性判别矩阵用MATLAB 实现如下：ans=4 No=obsv(a,c)No=1 0 0 0 0 0 1 0 0 1 0 0 0 0 0 1 0 0-0.239 0 0 0 10.0390 0 0 0 0-0.2390 0 0 0 10.039 rank(No)ans=4 rank(No)得 4，即能观性矩阵的秩为4，满秩，所以系统能观。稳定性：特征向量和特征根求解如下：v,x=eig(a)v=1.0000-1-

18、0.0072 0.0072 0 0-0.0227-0.0227 0 0 0.3009-0.3009 0 0 0.9534 0.9534 x=0 0 0 0 0 0 0 0 0 0 3.1684 0 0 0 0-3.1684 由特征根 X可知：状态矩阵A的特征值为 0，0，3.1684,-3.1684。平衡状态渐近稳定的充要条件是矩阵A 的所有特征值均具有负实部。此系统不满足条件，所以，系统不稳定。另由开环特性曲线也可知系统不稳定。名师资料总结-精品资料欢迎下载-名师精心整理-第 17 页，共 30 页 -通过状态反馈配置改变闭环系统极点。闭环极点自行决定。采用极点配置后，闭环系统的响应指标满足

19、如下要求为：摆杆角度和小车位移的稳定时间小于5 秒位移的上升时间小于2s 角度的超调量小于20 度位移的稳态误差小于2%。因系统为 4 阶系统，配置极点可采用主导极点加两个副极点的方式。调节时间sTsn54.4（1）上升时间sTrn21arccos2（2）由（1）得88.0n，在欠阻尼响应曲线中阻尼比越小，超调量越大，上升时间越短，一般为宜8.04.0，此时超调量适度，调节时间较短。n越大，响应时间越快，所以暂取4.0，追求快速性；n=2.5。将其代入(2)进行验算，得 Tr=0.8652，满足系统要求。角度超调量和位移稳态误差暂时不好求取，可通过对极点配置后的系统进行观察来确定是否满足条件。

20、)3.21)(3.21(25.622222jsjsssssnn可得，主导极点：u1=-1+2.3j，u2=-1-2.3j 选取副极点为-10+0.01j和-10-0.01j，由此可得，极点矩阵 P=-10+0.01*j -10-0.01*j -1+2.3*j -1-2.3*j。反馈矩阵 K=place（A，B，P）=-131.5780 -68.1528 -452.0558 -113.2532 名师资料总结-精品资料欢迎下载-名师精心整理-第 18 页，共 30 页 -极点配置系统结构图位移角度曲线由图可知，前面增益选取K=-131，使位移稳定在 1，且稳态误差小于2%，上升时间小于 2 秒，调

21、节时间小于5 秒，角度的最大超调量为0.3445rad，转换成角度为 19.748 小于 20 度，所以此极点配置满足要求。假设系统的状态Tssx均无法测量，为实现上述控制方案建立系统的全维观测器，观测器极点自行决定。采用带有观察器极点配置后，闭环系统的响应指标满足如下要求为：摆杆角度和小车位移的稳定时间小于5 秒位移的上升时间小于2 秒角度的超调量小于20 度位移的稳态误差小于2%。名师资料总结-精品资料欢迎下载-名师精心整理-第 19 页，共 30 页 -由全维观测器方程)(/yyGBuAxxGyBuxGCA)(其中 G为输出误差反馈矩阵。改变极点的配置，P=-15+0.01*j -15-

22、0.01*j -1+2.3*j -1-2.3*j；G可由 Matlab 语句求出。G=place(A,C,P),得 G=3036.245509.129618.158507.07559.946317.153424.60382.16全维观测器模拟结构图全维观测器位移 s 和角度响应曲线：名师资料总结-精品资料欢迎下载-名师精心整理-第 20 页，共 30 页 -位移对比曲线：角度对比曲线：加入阶跃后各参量对比曲线如下：位移响应曲线：名师资料总结-精品资料欢迎下载-名师精心整理-第 21 页，共 30 页 -角度响应曲线：设置初始值后各参量对比曲线如下：位移初始值为 0.1：角度初始值为 0.05：

23、名师资料总结-精品资料欢迎下载-名师精心整理-第 22 页，共 30 页 -比较以上各图可以看出，此全维观测器的响应曲线和系统的响应曲线基本吻合，且调节时间小于 5 秒，上升时间小于 2 秒，角度超调量 0.3443rad，即 19.73度，稳态误差小于2%，均满足要求，所以，此全维观测器满足系统要求。假设系统的状态Tssx中，只用位移 s 可以测量，其他状态变量均无法测量，为实现极点配置，建立系统的降维观测器，观测器极点自行决定。采用带有观察器极点配置后，闭环系统的响应指标满足如下要求为摆杆角度和小车位移的稳定时间小于5 秒位移的上升时间小于2 秒角度的超调量小于20 度位移的稳态误差小于2

24、%系统为DuCxyBuAxx，D=0，即CxyBuAxx系统能观，且 rankC=1，则必存在线性变换xTx，使：222112111AAAAATTA，211BBBTB，ICTC0令0C=010000100001，所以1T=CC0=1000010000100001名师资料总结-精品资料欢迎下载-名师精心整理-第 23 页，共 30 页 -取100001000010000110CCT所以，00100039.100010000239.0001ATTA04878.004878.01BTB，1000CTCrankC=1，4-1=3，所以需设计 3 维降维观测器，由上可知0,001,000,0039.10

25、01000239.0022211211AAAA主导极点 P=-1+2.3j-1-2.3j -3;G=(place(A11,A21,P);9749.2884268.935G带入式：uBGByAGAGAGAwAGAw)()()(?)(21221221112111，wx?+yG得：ywxuyww9749.2884268.935?4878.004878.09631.506159.178671.2?0039.109749.288104268.930239.05反馈矩阵 K=place（A，B，P）=-131.5780 -68.1528 -452.0558 -113.2532 名师资料总结-精品资料欢迎下

26、载-名师精心整理-第 24 页，共 30 页 -加入降维观测器后各参量对比曲线：角速度对比曲线：速度对比曲线：角度对比曲线：名师资料总结-精品资料欢迎下载-名师精心整理-第 25 页，共 30 页 -设置速度初始值为0.5 设置角度初始值为0.01 设置角速度初始值为0.1 名师资料总结-精品资料欢迎下载-名师精心整理-第 26 页，共 30 页 -加入阶跃后各参量响应曲线：角速度对比曲线：速度对比曲线：角度对比曲线：名师资料总结-精品资料欢迎下载-名师精心整理-第 27 页，共 30 页 -由位移角度响应曲线可知，稳定时间小于 5 秒，上升时间小于 2 秒，角度超调量等于 0.3445rad

27、，转换成角度为19.74 度，稳态误差小于2%，系统满足要求。名师资料总结-精品资料欢迎下载-名师精心整理-第 28 页，共 30 页 -课程设计总结：为期半个月的实习今天结束了，这半个月的实习是我获益匪浅，我对自动控制原理和现代控制理论又有了一个全面性的理解，实习之前我所掌握的知识只是课本上的文字，这次实习之后我懂得了如何把这些知识应用于实习工作之中，既巩固了所学的知识，又掌握了将这些知识应用于实际生产中的方法，还增加了我团队合作的意识，这次实习是很有意义的一次实习。为了达到系统预定的性能指标，我们每个步骤都做了很多组试验，很多会因为仪器的故障一级参数的设置出现各种问题，但我们都一一解决了这

28、些问题，得到了较为理想的结论，使系统的性能达到最好。在这些过程中，我们掌握了各种器材以及软件的使用方法，掌握了自动控制原理的时域分析法，根轨迹法，频域分析法，以及各种校正装置的作用及用法，能够利用不同的分析法对给定系统进行性能分析，能根据不同的系统性能指标要求进行合理的系统设计，并调试满足系统的指标。能够熟练使用MATLAB语言,SISOTOOL以及 Simulink动态仿真工具来搭建系统结构图，并进行系统仿真与调试。学会使用软硬件结合的方法对系统进行模拟仿真以便确定系统的性能指标达到系统稳定性的要求。本次实习分为两部分：第一部分是二阶水箱液位控制系统的设计，第二部分是基于状态空间法单级倒立摆

29、的控制系统设计。在第一部分中，首先，通过实验，取得实际水箱系统的阶跃曲线，从而得到其数学模型，通过对理论系统的分析，建立电路模型，借用 MATLAB 这一强大的数学工具对其分析和仿真设计；接着通过 P,PI,PID 对其控制，并分析对比控制效果；并得出结论：(1)Ti，Td 一定时，Kp 增大，加快系统的响应，系统的超调量增大，调节时间变小，上升时间减小，减小余差；(2)Kp，Ti 一定时，Td 增大，系统的峰值时间减小，名师资料总结-精品资料欢迎下载-名师精心整理-第 29 页，共 30 页 -系统的超调量减小，振荡减小，调节时间减小。(3)Kp，Td 一定时，Ti 增大，系统的超调量减小，

30、减小振荡，使系统更加稳定，但余差消除的速度随之减慢。最后用串联超前滞后校正对其补偿，改善系统的稳态性能和动态性能，其中，超前校正主要是利用超前网络的相角超前特性，减小系统的截止频率，而滞后校正则是利用滞后网络的高频幅值衰减特性，加强系统的抗干扰能力。第二部分中，主要是应用现代控制理论的知识，首先查阅资料，建立倒立摆的模型，然后分别分析其能控性，能观性，稳定性；并分别通过极点配置，全维观测器，降维观测器对性能进行控制分析；了解了系统进行极点配置的方法，掌握了状态观测器的存在条件和实现条件，其存在的充分必要是系统的不能观子系统为渐进稳定；可实现的条件是系统完全能观。降维观测器与全维观测器相比，主要是维数的降低，可以减少控制器设计的成本，但是计算起来比较麻烦。通过本次实习对这些知识有了充分的了解，并且熟练的掌握了方法和要领，在以后的学习中一定会熟练的运用这些知识。名师资料总结-精品资料欢迎下载-名师精心整理-第 30 页，共 30 页 -