实用标准文档文案大全实验三：最小拍控制系统一、实验目的1．建立计算机最小拍控制系统的一般概念；2．掌握有纹波最小拍控制器的设计方法3．观察无纹波最小拍控制器的设计方法；4．了解最小拍控制器的优缺点；5．掌握最小拍控制系统的改进方法。二、实验仪器1.PC计算机一台三、实验内容()DZ1TseS5(1)SSREC图11.针对图一所示的计算机控制系统,考虑输入为单位速度信号时，进行计算控制算法D(Z)设计，编程实现最小拍有纹波系统；2.讨论纹波的生成原因，编程实现最小拍有纹波系统；纹波产生原因：控制信号u(t)的变化引起输出的波动无波纹最小拍系统的设计出了满足有波纹最小拍系统的三个条件外，还必须满足Φ(z）包括G(Z)的全部零点。3.讨论最小拍系统的特点，采取惯性因子法对最小拍控制器加以改进，并研究惯性因子对系统性能的影响。最小拍系统的特点：系统对应于典型的输入信号具有最快的响应速度，被控量在最短的时间达到设定值。局限性：1)仅适应与一种类型的输入信号2)系统参数变化引起系统极点位置产生偏移，系统动态性能明显变坏。实用标准文档文案大全3)采样周期的限制。惯性因子法：基本思想：牺牲有限拍的性质为代价，换取系统对不同输入类型性能皆能获得比较满意的控制效果。方法：引进惯性因子，改进系统的闭环脉冲传递函数，使系统对多种信号的输入信号都有满意的响应。相当于增加单位圆内实轴上的极点。设被控对象传递函数)1(5)(sszG，采样周期sTs1.0，在单位速度输入作用下，采用惯性因子法设计有限拍控制系统，编程仿真实现，并讨论惯性因子对系统性能的影响。四、实验步骤1.第一步：求广义脉冲传递函数>>Gs=tf([5],[1,1,0])；>>Gz=c2d(Gs,0.1,'zoh')//求解广义对象的脉冲传递函数Transferfunction:0.02419z+0.02339----------------------z^2-1.905z+0.9048Samplingtime:0.1第二步：确定闭环脉冲传递函数Φ(Z)由广义对象的闭环脉冲传递函数得延迟因子为z1，则Φ(Z)应包含该延迟因子>>[num,den]=c2dm([5],[1,1,0],0.1,'zoh')num=00.02420.0234den=1.0000-1.90480.9048实用标准文档文案大全>>tf2zpk(num)ans=-0.9672由此可得：G(z)没有单位圆上或圆外的零点>>tf2zpk(den)ans=1.00000.9048由此可得：G(z)包含一个单位圆上极点根据物理可实现条件和稳定性条件,Φ(Z)=z1F1(z)根据零稳态误差条件和稳定性条件,Φe(Z)=1-Φ(Z)=(1-z1)^2*F2(z)根据Φe(Z)=1-Φ(Z)F1(z)=2-z1F2(z)=1所以Φe(Z)=(1-z1)^2>>Qez=tf([1-21],[1,0,0],0.1)>>Qz=1-QezTransferfunction:2z-1-------z^2Samplingtime:0.1>>Dz=1/Gz*Qz/(1-Qz)Transferfunction:2z^5-4.81z^4+3.715z^3-0.9048z^2----------------------------------------------------0.02419z^5-0.02498z^4-0.0226z^3+0.02339z^2Samplingtime:0.1实用标准文档文案大全>>Qz=minreal(Dz*Gz/(1+Dz*Gz));>>Q=2;>>t=0:0.1:Q;u=t;plot(0:0.1:Q,u,'r*');holdon>>yt=lsim(Qz,u,t,0);plot(0:0.1:Q,yt);结果如图一2.根据物理可实现条件和稳定性条件以及无纹波条件Φ(Z)=z1*(1+0.9672z1)F1(z)根据零稳态误差条件和稳定性条件Φe(Z)=1-Φ(Z)=(1-z1)^2*F2(z)根据Φe(Z)=1-Φ(Z)F1(z)=1.2666-0.7583z1F2(z)=1+0.7334z1所以Φe(Z)=1-1.2666z1-0.4668z2-0.7334z3>>Qez=tf([1-1.2666-0.46680.7334],[1,0,0,0],0.1);>>Qz=1-QezTransferfunction:1.267z^2+0.4668z-0.7334-----------------------------z^3Samplingtime:0.1设计最小拍无纹波系统控制器>>Dz=1/Gz*Qz/(1-Qz)Transferfunction:1.267z^7-1.946z^6-0.4765z^5+1.819z^4-0.6636z^3---------------------------------------------------------------------0.02419z^7-0.007241z^6-0.04092z^5+0.006818z^4+0.01716z^3Samplingtime:0.1>>Qz=minreal(Dz*Gz/(1+Dz*Gz));实用标准文档文案大全>>Q=2;>>t=0:0.1:Q;u=t;plot(0:0.1:Q,u,'r*');holdon>>yt=lsim(Qz,u,t,0);plot(0:0.1:Q,yt);结果如图二3.改变α的值，观测实验结果α=0.8Qez=tf([1-1.2666-0.46680.7334],[1,0.8,0,0],0.1);Qz=1-Qez;Dz=1/Gz*Qz/(1-Qz);Qz=minreal(Dz*Gz/(1+Dz*Gz));Q=2;t=0:0.1:Q;u=t;plot(0:0.1:Q,u,'r*');holdonyt=lsim(Qz,u,t,0)...
