发布时间 : 星期三 文章计算机控制实验指导书(LABVIEW)2016.1.22.更新完毕开始阅读9a4299f5185f312b3169a45177232f60ddcce7d3
计算机控制技术实验指导书
实验三 D(s)离散化方法的研究
一．实验目的
1．学习并掌握数字控制器的混合仿真实验研究方法。
2．熟悉常用的从连续化途径（先按连续系统设计再离散化）设计数字控制器的方法。 3．学习并掌握将模拟控制器D(s)离散化为数字控制器D(z)的方法。
4．通过混合仿真实验，对D(s)的各种离散化方法作比较研究，并对D(s)离散化前后闭环系统性能作比较研究，以加深对计算机控制系统特性的理解。
二．实验内容
1．按连续系统设计串联校正控制器D(s)，并利用实验设备测取该连续系统的动态特性。 2．利用实验设备，设计并构成用于混合仿真实验的计算机闭环控制系统。
3．采用保持器等价、匹配Z变换、一阶差分和双线性变换等方法离散化D(s)，从而得到控制算法，并加以实现。
4．研究采样控制周期变化时，不同离散化方法对闭环控制系统性能的影响。 5．对上述连续系统和计算机控制系统的动态性能作比较研究。
三．实验步骤
1．在已知二阶被控对象传递函数的条件下，用连续系统综合方法设计串联校正控制器D(s)；再利用实验设备，设计并连接用于模拟该闭环控制系统的电路，在上位机测取并记录其阶跃响应。
2．设计并连接模拟二阶被控对象的电路，并利用C8051F060构成的数据采集系统完成计算机控制系统的模拟量输入、输出通道的设计和连接。利用上位机的虚拟仪器功能对此模拟二阶被控对象的电路进行测试，根据测试结果调整电路参数，使它满足实验要求。
3．选定采样控制周期，采用保持器等价法离散化D(s)，再从D(z)推导控制算法，然后在上位机上完成该算法编程、调试和运行。以阶跃信号作为系统输入，观测系统输出的屏幕显示，并作记录。
4．以相同的采样控制周期，采用匹配Z变换法离散化D(s)，再从D(z)推导控制算法，然后在上位机上完成该算法编程、调试和运行。以阶跃信号作为系统输入，观测系统输出的屏幕显示，并作记录。
5．以相同的采样控制周期，采用一阶差分法离散化D(s)，再从D(z)推导控制算法，然后在上位机上完成该算法编程、调试和运行。以阶跃信号作为系统输入，观测系统输出的屏幕显示，并作记录。
6．以相同的采样控制周期，采用双线性变换法离散化D(s)，再从D(z)推导控制算法，然后在上位机上完成该算法编程、调试和运行。以阶跃信号作为系统输入，观测系统输出的屏幕显示，并作记录。
7．将采样控制周期减小，重复3－6，观测采样控制周期减小对系统阶跃响应的影响。 8．将采样控制周期不断增加，重复3－6，观测采样控制周期增大对系统阶跃响应的影响。如有系统不稳定情况发生，记下此时的采样控制周期和所采用的离散化方法。
9．对实验结果进行比较、分析和研究，写成实验报告。
9
计算机控制技术实验指导书
四．附录 1．按连续系统串联校正及其动态性能 实验系统被控对象的传递函数为 550 （3－1） ?0.2s(0.5s?1)s(s?2)具有串联校正控制器G(s)的线性连续系统的结构方块图如图3.1所示 G(s)?r(t)e(t)D(s)u(t)G(s)y(t)图3.1 按以下要求设计期望系统的开环对数幅频特性： （1）超调量Mp?20% （2）调节时间（过渡过程时间）ts?1s （3）校正后系统开环增益（静态速度误差系数）Kv?25% 从期望系统的开环对数幅频特性中，减去上述二阶被控对象的对数幅频特性，可以得到串联校正控制器的对数幅频特性，由此得到其传递函数 D(s)?0.5s?1 0.05s?1 （3－2） 已知G(s)和D(s)后，参阅图3.1所示系统结构，设计系统被控对象的模拟电路如图3.2所示。
R2200kR31u500k100k+-+++CR4R1-+r(t)200k200k-++1u-+200kc(t)图3.2 10 计算机控制技术实验指导书
2．计算机闭环控制系统的混合仿真
1u500k1u-++Om100k-++200kIn图3.3实验系统被控对象的传递函数见式（3－1），它可以用图3.3所示电路来模拟。 图中电路可利用U9 、U11等单元组成。其中Om端连U3单元O1，In端连U3单元I1，以便利用该单元的数据处理功能与上位机的虚拟仪器功能观测和记录系统动态过程，用于对比分析。 R(t)e(t)D(z)Z.0.Hy(t)G(s)图3.4 计算机控制系统的方框图如图3.4所示，除了虚线框内部分用电路模拟外，其余部分由上位机和数据处理系统完成。 图3.4中，Z.0.H为零阶保持器，其传递函数为1?e?Tss（以后实验中均同此义）。 3．D(s)的离散化方法 （1）保持器等价 1?e?Ts10?(e?20T?9)z?1D(z)?Z[D(s)]?
s1?e?20Tz?1u(k)?e?20Tu(k?1)?10e(k)?(e?20T?9)e(k?1)
（2）匹配Z变换
对应D(s)的 s1??20、s2??2，有z1?e?20T、z2?e?2T
K(1?e?2Tz?1)故D(z)? ?20T?11?ez1?e?20T从limD(s)?limD(z) 确定 K 得K? s?0z?11?e?2T?20T?2T所以 u(k)?eu(k?1)?Ke(k)?Kee(k?1)
（3）一阶后向差分
0.5s?11?z?1以s?代入D(s)?有
0.05s?1T 11
计算机控制技术实验指导书 0.5?T?0.5z?110.5?T?0.5z?1D(z)??? ?10.050.05?T?0.05z0.05?T1?z?10.05?T0.050.5?T0.5u(k?1)?e(k)?e(k?1) 所以 u(k)?0.05?T0.05?T0.05?T（4）双线性变换
0.5s?121?z?1D(s)?以s??代入有 ?10.05s?1T1?z1?T?(1?T)z?111?T?(1?T)z?1D(z)??? ?10.1?T0.1?T?(0.1?T)z0.1?T1?z?10.1?T0.1?T1?T1?Tu(k?1)?e(k)?e(k?1) 所以 u(k)?0.1?T0.1?T0.1?T 12