过程控制-第一到三章-作业

第一章作业1.1 常用的评价控制系统动态性能的单项性能指标有哪些？它与误差积分指标各有何特点？答：(1)衰减率、超调量、稳态误差ess、调节时间ts、振荡频率；(2)单项指标用若干特征参数评价系统优劣,积分指标用误差积分综合评价系统优劣。1.2 什么是对象的动态特性？为什么要研究对象的动态特性？答：(1)指被控对象的输入发生变化时,其输出（被调量）随时间变化的规律；(2)实现生产过程自动化时,对象的动态特性可以为控制工程师设计出合理的控制系统满足要求提高主要依据。1.3 通常描述对象动态特性的方法有哪些？答：微分方程或传递函数。1.4 过程控制中被控对象动态特性有哪些特点？答：无振荡、稳定或中性稳定、有惯性或迟延、非线性但在工作点附近可线性化。1.11 某水槽水位阶跃响应实验为：t /s0 10 20 40 60 80 100 150 200 300 400h /mm0 9.5 18 33 45 55 63 78 86 95 98其中阶跃扰动量µ=20%。(1)画出水位的阶跃响应曲线；(2)若该水位对象用一阶惯性环节近似,试确定其增益K 和时间常数T 。解：MATLAB编程如下：%作出标幺后的响应曲线t= 0 10 20 40 60 80 100 150 200 300 400 ;h= 0 9.5 18 33 45 55 63 78 86 95 98 ;x=0:0.01:400; y=interp1(t,h,x,'spline'); %三次样条函数据己知的t、h插出x的值 yy=y/y(end); %输出标幺plot(x,yy,'k'); xlabel('t/s');ylabel('h/mm');title('阶跃响应曲线','fontsize',10);grid;%找出最接近0.39和0.63的点less1=find(yy<=0.39); more1=find(yy>=0.39); front1=less1(1,end);behind1=more1(1,1);cha11=0.39-yy(1,front1);cha12=yy(1,behind1)-0.39;if cha11<=cha12 t1=x(1,front1)else t1=x(1,behind1)endless2=find(yy<=0.63); more2=find(yy>=0.63); front2=less2(1,end);behind2=more2(1,1);cha21=0.63-yy(1,front2);cha22=yy(1,behind2)-0.63;if cha21<=cha22 t2=x(1,front2)else t2=x(1,behind2)end %求增益K和时间常数T K=y(end)/20T=2*(t2-t1)(1)水位的阶跃响应曲线如图：(2) 计算结果如下：>> GK1_11t1 = 48.1700t2 = 96.5900K = 4.9T = 96.8400则该水位对象用一阶惯性环节近似后,得其增益K=4.9,时间常数T96.84。1.14 已知温度对象阶跃响应实验结果如下表：t /s0102030405060708090100150 /00.160.6511.51.521.751.881.941.971.992.02.0其中阶跃扰动量q=1t/h。试用二阶或n阶惯性环节写出它的传递函数。解：由表可知,t =30 s时 =11.5,显然存在粗大误差,故将此数据剔除。MATLAB编程如下：%作出标幺后的响应曲线t= 0 10 20 40 50 60 70 80 90 100 150 ;theta= 0 0.16 0.65 1.52 1.75 1.88 1.94 1.97 1.99 2.00 2.00 ;x=0:0.01:150; y=interp1(t,theta,x,'pchip'); %三次函数据己知的t,theta插出x的值 yy=y/y(end); %输出标幺plot(x,yy,'k'); xlabel('t/s'); ylabel('/');title('阶跃响应曲线','fontsize',10); grid on;%找出最接近0.4和0.8的点less1=find(yy<=0.4); more1=find(yy>=0.4); front1=less1(1,end); behind1=more1(1,1);cha11=0.4-yy(1,front1); cha12=yy(1,behind1)-0.4;if cha11<=cha12 t1=x(1,front1)else t1=x(1,behind1)endless2=find(yy<=0.8); more2=find(yy>=0.8); front2=less2(1,end); behind2=more2(1,1);cha21=0.8-yy(1,front2); cha22=yy(1,behind2)-0.8;if cha21<=cha22 t2=x(1,front2)else t2=x(1,behind2)end %求增益K和时间常数Ttt=t1/t2K=y(end)/1T=(t2+t1)/2.16/4温度的阶跃响应曲线如图：计算结果如下：>> GK1_14t1 = 23.1600t2 = 42.9200tt = 0.5396K = 2T = 7.6481因为tt=t1/ t20.5396,据表1-1取n=4。则该温度对象用四阶惯性环节近似后,得其增益K=2,时间常数T7.648,传递函数为第二章作业2.1 试确定题图中各系统调节器的正、反作用方式。答：(a)加热炉温度控制系统：燃料调节阀是气开式,则Kv>0；燃料增加时温度上升,则Gp(s)的K>0；前两者为正,则Gc(s)的Kc>0,即调节器为反作用方式。(b)锅炉汽包液位控制系统：给水调节阀是气关式,则Kv<0；给水增加时液位上升,则Gp(s)的K>0；前两者一负一正,则Gc(s)的Kc<0,即调节器为正作用方式。2.2 什么是调节器的动作规律？P、I、D控制规律各有何特点？答：(1)指调节器在动态中输入变化量e与输出变化量u之间的关系；(2) P控制规律：调节作用及时,偏差一旦产生,调节器立即使被控参数朝着减小偏差的方向变化,Kc增大则控制作用增强；有差调节,Kc增大则静差减小；但Kc过大会引起振荡,导致系统不稳定；I控制规律：调节作用不及时,TI减小则控制作用增强；无差调节,只要偏差存在,调节器就一直调整输出直至偏差为零；积分作用使系统稳定性变差；D控制规律：输出量与输入量的变化速度成正比,调节过程结束后偏差的变化速度为零,无论偏差多大调节器都不作用。2.3 某电动比例调节器的测量范围为100-200,其输出为0-10mA。当温度从140变化到160时,测得调节器的输出从3mA变化到7mA。求该调节器比例带。解： 2.4 某水槽液位控制系统如题图。已知：F=1000cm2,R=0.03s/cm2,调节阀为气关式,其静态增益|Kv|=28cm3/s·mA,液位变送器静态增益Km=1mA/cm。(1) 画出该系统的传递方框图；(2) 调节器为比例调节器,其比例带=40%,分别求出扰动Qd=56cm3/s以及定值扰动r=0.5mA时,被调量h的残差；(3) 若改为120%,其他条件不变,h的残差又是多少？比较(2)(3)的结果,总结对系统残差的影响；(4) 液位调节系统改用PI调节器后,h的残差又是多少？解：(1)其中,T=CR=FR=1000*0.03=30s,K=R=0.03s/cm2 Gp(s)=K/(Ts+1)=0.03/(30s+1)(2) 调节阀是气关式,则Kv<0；Gp(s)的K>0且Km>0,则Gc(s)的Kc<0(3)类似(2)的计算总结：比例带增大会导致残差增大。(4)h的残差为零,因为PI调节器的积分作用为无差调节。2.8被控对象传递函数为G(s)=K/s(Ts+1),如采用积分调节器,证明：积分速度So无论为何值,系统均不能稳定。证：令1+Gc(s)G(s)=1+So /sK/s(Ts+1)=0 Ts3+s2+KSo=0劳斯判据表：s3 T 0 s2 1 KSo s1 -TKSo s0 KSo根的分布：左半平面1个,右半平面2个所以可知,So无论为何值,系统均不能稳定。2.10 一个自动控制系统,在比例控制的基础上分别增加:适当的积分作用；适当的微分作用。试问：(1)这两种情况对系统的稳定性、最大动态偏差、残差分别有什么影响？(2)为了得到相同的系统稳定性,应如何调整调节器的比例带？并说明理由。答：(1)I：稳定性变差,最大动态偏差变大,消除残差；D：稳定性变好,最大动态偏差减小,不能消除残差。(2)I：因积分调节使系统稳定性变差,故可适当增大比例带,减弱比例调节作用；D：因为微分调节使系统稳定性变好,故可适当减小比例带,增强比例调节作用。2.11 比例微分控制系统的残差为什么比纯比例控制系统的小？答：微分调节总是试图抑制被调量振荡,可以提高系统的稳定性。在保持衰减率不变的情况下,适度引入微分作用后,可以允许减小比例带。而比例带减小,静差则减小。2.15 微分动作规律对克服被控对象的纯迟延和容积迟延的效果如何？答：纯迟延对象：在延迟时间段由于对象不发生变化,故微分作用在此阶段不起作用。此段时间过后对象的变化速率一定,故微分作用起作用,且强度保持不变；容积迟延对象：对象从一开始便一直变化,且变化速率不定,故此阶段内微分作用存在且随时间变化。最后对象不再变化时,微分作用消失,不再起作用。第三章作业3.1 为什么要对控制系统进行整定？整定的实质是什么？答：(1)不同的被控对象对调节器的特性要求不同,系统能否在最佳状态下工作,主要取决于控制器各参数的设置是否得当；(2)通过调整控制器的这些参数,使其特性与被控对象特性相匹配,以达到最佳的控制效果。3.2 正确选择系统整定的最佳性能指标有何意义？目前常用性能指标有哪些？答：(1)能够综合反映系统控制质量,而且便于分析和计算；(2)衰减率、超调量、调节时间ts、振荡频率。3.3 在简单控制系统中,调节器为比例动作。广义被控对象的传递函数如下,用衰减频率特性法求：=0.75(m=0.221)和=0.90(m=0.366)时,调节器的整定参数。(1) (2)3.5 某温度控制系统对象阶跃响应中,测得：K=10,T=2min,=0.1min,应用动态特性参数法设计PID调节器整定参数。解：=K/T=5/min =0.5 (1)一次P调节：=0.5即P=KC=2。 衰减振荡过程,但不满足=0.75的要求。 二次P调节：=0.48即P=Kc2.08。0.751,如下图：(2)一次PI调节：=1.1=0.528即P=KC1.894。 TI=3.3=0.33即I=KC/TI5.739。 衰减振荡过程,但不满足=0.75的要求。 二次PI调节：=0.485即P=KC2.06。 TI=2.06即I=KC/TI=1。 0.7514,=0.527,如下图：(3)一次PID调节：=0.7730.375即P=KC2.667。 TI=0.606TI1.248即I=KC/TI0.801。 TD=0.25TI=0.312即D=KCTD0.832。 振荡发散过程。 二次PID调节：=0.435即P=KC2.30。 TI=2.255即I=KC/TI1.020。 TD=0.024即D=KCTD0.055。 0.75,0.5,如下图：3.7已知被控对象阶跃响应曲线数据如下表,调节量阶跃变化 u=5。时间/min0510152025被调量0.6500.6510.6520.6680.7350.817时间/min303540455055被调量0.8810.9791.0751.1511.2131.239时间/min606570758085被调量1.2621.3111.3291.3381.3501.351(1)用一阶惯性环节加纯迟延近似对象,求出K、T、值；(2)应用动态特性参数法选择PID调节器参数。解：(1)MATLAB编程如下：%作出标幺后的响应曲线t=0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85;h=0.65 0.651 0.652 0.668 0.735 0.817 0.881 0.979 1.075 1.151 1.213 1.239 1.262 1.311 1.329 1.338 1.35 1.351;x=0:0.1:85; y=interp1(t,h,x,'spline'); %三次样条函数根据己知的t,h插出x的值 yyy=y-y(1); xlabel('时间/min'); ylabel('被调量');title('阶跃响应曲线','fontsize',10);plot(x,yyy,'k'); grid on; hold on;%作出x=38(即x0)处的切线x0=x(381); y0=yyy(381); %取切点x1=x(380); y1=yyy(380); %取临近点x2=x(382); y2=yyy(382);k=(y2-y1)/(x2-x1); b=y0-k*x0; %计算斜率和截距f=(x) k*x+b; %切线方程yk=f(x); plot(x,yk);阶跃响应曲线如下图： 计算结果：K=(1.351-0.65)/5=0.1402,T=54.5-17.5=37min,=17.5-10=7.5min。(2)=K/T0.0038/min =0.0285 (A)一次P调节：=0.0285即P=KC35.0877。 衰减振荡过程,但不满足=0.75的要求。 二次P调节：=0.0265即P=Kc37.736。0.7515,如下图：(B)一次PI调节：=1.10.029即P=KC34.483。 TI=3.3=24.75即I=KC/TI1.393。 衰减振荡过程,但不满足=0.75的要求。 二次PI调节：=0.027即P=KC37.037。 TI=46.3即I=KC/TI0.8。 0.7530,=0.494,如下图：(C)一次PID调节：=0.7730.021即P=KC47.619。 TI=0.606TI28.058即I=KC/TI1.697。 TD=0.25TI=7.0145即D=KCTD334.023。 振荡发散过程。 二次PID调节：=0.025即P=KC=40。 TI=40即I=KC/TI=1。 TD=0.6即D=KCTD=24。 0.7495,=0.519,如下图：3.11 对题图所示的控制系统中的调节器,试用稳定边界法整定参数。(a)当P=KC=0.328时,系统出现等幅振荡(如下图),得： cr3.049且Tcr=(42.73-2.973)/49.94(1)一次P调节：=2cr=6.098即P=KC0.164。 衰减振荡过程,但不满足=0.75的要求。 二次P调节：=29即P=KC0.0345。0.7486,如下图：(2)一次PI调节：=1.1=31.9即P=KC0.0313。 TI=0.85Tcr=8.449即I=KC/TI0.0037。 振荡发散过程。 二次PI调节：=83.333即P=KC0.012。 TI=24.490即I=KC/TI0.00049。 0.7510,1.036,如下图：(b)当P=KC=2.06724时,系统出现等幅振荡(如下图),得：cr0.488且Tcr=(74.16-62.08)/5=2.416(1)一次P调节：=2cr=0.976即P=KC1.025。 衰减振荡过程,但不满足=0.75的要求。 二次P调节：=1即P=KC=1。0.7519,如下图：(2)一次PI调节：=1.1=1.1即P=KC0.909。 TI=0.85Tcr=2.0536即I=KC/TI0.443。 单调稳定过程,无衰减振荡。 二次PI调节：=0.808即P=KC1.238。 TI=1.125即I=KC/TI1.100。 0.7493,=0.1755,如下图：(3)一次PID调节：=0.7590.613即P=KC1.631。 TI=0.588TI=0.6615即I=KC/TI2.466。 TD=0.25TI0.165即D=KCTD0.269。 振荡发散过程。 二次PID调节：=1.087即P=KC0.920。 TI=0.46即I=KC/TI=2。 TD=0.001即D=KCTD0.0009。 0.7502,0.5025,如下图：科教兴国13 / 13