动控制课程设计报告班级:自动化08-1班学号:********姓名:***2018.7.17任务一、双容水箱的建模、仿真模拟、控制系统设计一、控制系统设计任务1、通过测量实际装置的尺寸,采集DCS系统的数据建立二阶水箱液位对象模型。
<先建立机理模型,并在某工作点进行线性化,求传递函数)2、根据建立二阶水箱液位对象模型,在计算机自动控制实验箱上利用电阻、电容、放大器的元件模拟二阶水箱液位对象。
3、通过NI USB-6008数据采集卡采集模拟对象的数据,测试被控对象的开环特性,验证模拟对象的正确性。
4、采用纯比例控制,分析闭环控制系统随比例系数变化控制性能指标<超调量,上升时间,调节时间,稳态误差等)的变化。
5、采用PI控制器,利用根轨迹法判断系统的稳定性,使用Matlab中SISOTOOLS设计控制系统性能指标,并将控制器应用于实际模拟仿真系统,观测实际系统能否达到设计的性能指标。
6、采用PID控制,分析不同参数下,控制系统的调节效果。
7、通过串联超前滞后环节校正系统,使用Matlab中SISOTOOLS设计控制系统性能指标,并将校正环节应用于实际模拟仿真系统,观测实际系统能否达到设计的性能指标。
(一)建立模型(二)实验模型及改变阶跃后曲线:1.取阶跃曲线按照以下模型建立系统辨识模型:一般取为0.4和0.8计算上行阶跃各参数:T1=171.26 T2=50.50 K=160.47 t1=141 t2=338建立传递函数为:G(s>=计算下行阶跃各参数:T1=84.20 T2=48.67 K=148.08 t1=89 t2=198建立传递函数为:G(s>=2.建立机理模型Q1=k1*u1;Q2=k2*u2*;Q=k3*u3*;k1=10 ;k2=1.9;k3=1.65;阀门开度u1=50; u2=52 ; u3=51 ;水箱面积A1=1050 ;A2=600理论传递函数G(s>=;取辨识传递函数G<s)=根据建立的二阶水箱液位对象模型,在计算机自动控制实验箱上利用电阻、电容、放大器的元件模拟二阶水箱液位对象。
根据传递函数可得:R2/R1=1.5 R2=300k(200k+100k>,R1=200kR2C=7.02 C=23.4ufR5/R4=1 R5=R4=200kR5C=4.06 C=20.3uF(四)通过NI USB-6008数据采集卡采集模拟对象的数据,测试被控对象的开环特性,验证模拟对象的正确性系统的开环特性曲线为<实际曲线理论曲线):(五)采用纯比例控制,分析闭环控制系统随比例系数变化时控制性能指标改变Kp得到实际与理论曲线:Kp=1 : KP=2:Kp=3:结论:由上图和系统指标分析可得,Kp越大,系统响应速度越快,上升时间越短,调节时间峰值时间也相应减少,且稳态误差减小,但超调量增大,系统振荡加剧,Kp过大时会对实际的系统造成破坏,所以,系统应选择合适的Kp。
采用PI控制器,利用根轨迹法判断系统的稳定性,使用Matlab中SISOTOOLS设计控制系统性能指标,并将控制器应用于实际模拟仿真系统,观测实际系统能否达到设计的性能指标Kp=2Ti=5:Kp=2 Ti=4:Kp=2 Ti=6 :由以上曲线可以看出,用根轨迹设计的设计的性能指标在实际系统中可以达到较好的效果,理论曲线与实际曲线较吻合。
(七)采用PID控制,分析不同参数下,控制系统的调节效果Kp=2 Ti=10 Td=1Kp=2 Ti=10 Td=2Kp=2 Ti=9 Td=2 Kp=3 Ti=9Td=2由上面四个图对比分析可知:(1>Ti,Td一定时,Kp增大,加快系统的响应,系统的超调量增大,调节时间变小,上升时间减小,减小余差;(2> Kp,Ti一定时,Td增大,系统的峰值时间减小,系统的超调量减小,振荡减小,调节时间减小。
(3>Kp,Td一定时,Ti增大,系统的超调量减小,减小振荡,使系统更加稳定,但余差消除的速度随之减慢。
以上各曲线参数列表如下:Kp Ti Td Ts Tpσ%Tr199******** 5.872999990141014.365399999012821.02 3.7240421049.7 4.8250311040 4.9260251032.6 5.12928.513 3.2 6.439211.510 5.4421011412 6.2 6.12102914 1.3 6.5为被控对象设计串联校正环节,使用Matlab中 SISOTOOLS设计控制系统性能指标,并将校正环节应用于实际模拟仿真系统,观测实际系统能否达到设计的性能指标<1)不加校正环节由图可知相角裕度为112deg,截至频率。
单位阶跃响应得上升时间,调节时间,峰值时间,超调量6.37 %,稳态值为0.597。
<2)由于原系统的相角裕度较大,为了使校正效果明显,给原系统加入积分环节1/s,于是其相角裕度为,截至频率,且系统已经发散,故要进行串联校正。
<3)由第一个拐点处读得。
假设调节时间,校正后的相角裕度,则:由可得由图,为了选择校正网络衰减因子,要保证已校正系统的截止频率为所选的,则等式成立,因,则可得到<4)校正前:相应校正后:于是,相角裕度:解得:所以:其中是装置的滞后部分,是校正装置的超前部分。
超前校正主要是利用超前网络的相角超前特性,减小系统的截止频率,而滞后校正则是利用滞后网络的高频幅值衰减特性,加强系统的抗干扰能力。
取合适的超调和调节时间,可取K=0.743,因而可得C(s>校正后的广义传递函数为:G(s>C(s>加入串联校正相角裕度为45.3deg,截止频率为0.0227rad/sec。
加入超前滞后环节后的阶跃响应,上升时间为48.7s,超调量为34%,峰值时间为131s,调节时间为404s,稳态值为1。
电路仿真曲线与理论曲线基本吻合,满足系统的要求。
(九)改变Kp及Ts对系统的影响由上图可以看出,Kp越大,系统响应速度越快,上升时间越短,调节时间峰值时间也相应减少,且稳态误差减小,但超调量增大,系统振荡加剧,Kp过大时会对实际的系统造成破坏。
由上图可以看出,Ts增加使系统调节时间变长,超调量增大,调节精度下降,Ts过大时使系统震动加剧,破坏系统的稳定性。
因此Ts应选取较小的值。
(十)为被控对象设计最小拍无差控制器,并进行实验分析传递函数 G(z>=零极点模型=用matlab作出结构图:仿真波形为:由上图可以看出,系统在Ts=0.2s时达到稳定,达到最小拍无差控制器的控制要求。
任务二、基于状态空间法单级倒立摆的控制系统设计一、单级倒立摆介绍倒立摆系统具有高阶次、不稳定、多变量、非线性和强耦合等特性,是控制理论的典型研究对象。
如机器人行走过程中的平衡控制、火箭发射中垂直度控制和卫星飞行中的姿态控制等均涉及到倒置问题对倒立摆系统的研究在理论上和方法论上均有着深远意义。
单级倒立摆系统的原理图,如图1所示。
假设已知摆的长度为2l,质量为m,用铰链安装在质量为M的小车上。
小车由一台直流电动机拖动,在水平方向对小车施加控制力u,相对参考系差生的位移s。
若不给小车实施控制力,则倒置摆会向左或向右倾倒,因此,它是个不稳定的系统。
控制的目的是通过控制力u的变化,使小车在水平方向上运动,达到设定的位置,并将倒置摆保持在垂直位置上。
⒈建立单级倒立摆的状态空间数学模型。
取状态变量。
测试系统的开环特性。
图是系统中小车和摆杆的受力分析图。
其中,N和P为小车与摆杆相互作用力的水平和垂直方向的分量。
已知单级倒立摆的各项数据如下所示:水平方向;垂直方向:根据力矩平衡方程:;因为很小,且<<1,所以=,=1。
设系统状态空间方程为:方程组消去P,N变量,对解代数方程,解得如下:整理得到系统将已知的代入上面所得的状态空间方程,得单级倒立摆的位移s及角度开环特性曲线:Mc=ctrb(a,b>Mc =0 0.4878 0 0.11660.4878 0 0.1166 00 -0.4878 0 -4.8970-0.4878 0 -4.8970 0>> rank(Mc>即能控性rank(Mc>得4,矩阵Mc的秩为4,满秩,所以系统能控。
能观性:能观性判别矩阵用MATLAB实现如下:ans =4>> No=obsv(a,c>No =1 0 0 00 0 1 00 1 0 00 0 0 10 0 -0.239 00 0 10.0390 00 0 0 -0.23900 0 0 10.039>> rank(No>ans =4rank(No>得4,即能观性矩阵的秩为4,满秩,所以系统能观。
稳定性:特征向量和特征根求解如下:[v,x]=eig(a>v =1.0000 -1 -0.0072 0.00720 0 -0.0227 -0.02270 0 0.3009 -0.30090 0 0.9534 0.9534x =0 0 0 00 0 0 00 0 3.1684 00 0 0 -3.1684由特征根X可知:状态矩阵A的特征值为0,0,3.1684,-3.1684。
平衡状态渐近稳定的充要条件是矩阵A的所有特征值均具有负实部。
此系统不满足条件,所以,系统不稳定。
另由开环特性曲线也可知系统不稳定。
通过状态反馈配置改变闭环系统极点。
闭环极点自行决定。
采用极点配置后,闭环系统的响应指标满足如下要求为:●摆杆角度和小车位移的稳定时间小于5秒●位移的上升时间小于2s●角度的超调量小于20度●位移的稳态误差小于2%。
因系统为4阶系统,配置极点可采用主导极点加两个副极点的方式。
调节时间<1)上升时间<2)由<1)得,在欠阻尼响应曲线中阻尼比越小,超调量越大,上升时间越短,一般,此时超调量适度,调节时间较短。
越大,响应时间越快,所以暂取,追求快速性;=2.5。
将其代入(2>进行验算,得Tr=0.8652,满足系统要求。
角度超调量和位移稳态误差暂时不好求取,可通过对极点配置后的系统进行观察来确定是否满足条件。
可得,主导极点: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]极点配置系统结构图位移角度曲线由图可知,前面增益选取K=-131,使位移稳定在1,且稳态误差小于2%,上升时间小于2秒,调节时间小于5秒,角度的最大超调量为0.3445rad,转换成角度为19.748小于20度,所以此极点配置满足要求。