火电厂仿真课件-第6章 离散相似法仿真.ppt

《火电厂仿真课件-第6章 离散相似法仿真.ppt》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第8章连续系统仿真的离散相似法“离散相似法”——将一个连续系统进行离散化处理,然后求得与它等价的离散模型(差分方程)的方法.2021/10/31获取离散相似模型的两个途径:（1）对传递函数作离散化处理得离散传递函数——称为“频域离散相似模型”；（2）基于状态方程离散化——称为“时域离散相似模型”；8.2状态方程的离散化8.1传递函数的离散化本章介绍:8.3典型环节的离散化模型8.4常用非线性环节的仿真8.5离散相似法仿真的特点2021/10/328.1传递函数的离散化对连续系统进行数字仿真可以先在系统加入虚拟的采样器

2、和保持器，如图8-1所示，然后利用Z变换的方法求出系统的脉冲传递函数，再从脉冲传递函数求出对应于系统G(s)的差分方程。保持器TTuy图8-1连续系统离散化结构图2021/10/33如此处理的理由是：（1）由采样定理可知：当采样频率和信号最大频率满足：（8-2）的条件时，可由采样后的信号唯一地确定原始信号。（2）把采样后的离散信号通过一个低通滤波器，即可实现信号的重构。常用的滤波器是零阶保持器和一阶保持器等。因此，研究系统G(s)，在一定的条件下可能等效地研究图8-1所示的系统。注意：图8-1所示系统的采样开关和保持

3、器实际上是不存在的，而是为了将(8-1)式离散化而虚构的。2021/10/34根据图8-1,有脉冲传递函数:其中Gh(s)是保持器的传递函数。若选择不同的保持器，则可得不同的G(z)，见表8-1。（8-1）保持器的传递函数Gh(s)脉冲传递函数G(z)零阶：一阶：三角形：表8-1不同保持器的G(z)2021/10/35表8-1中的G(z)还要求对或等进行Z变换，这样做有时不太方便，如果能对G(s)直接进行Z变换就方便多了。因此，对表8-1中的各式还可以进一步加以简化，得出二次加入虚拟采样器和保持器的脉冲传递函数。例

4、如:当采用零阶保持器时需要对进行Z变换，它相当于G(s)与1/s串联，如图8-2所示的模型。TT图8-2若在积分环节之前再加一个采样器和保持器，如图8-3所示:TT图8-3T2021/10/36那么就可以对G(s)及1/s分别求出脉冲传递函数，然后相乘，即:因此,二次加入采样器、零阶保持器后的脉冲传递函数为:（8-2）（8-3）从上述可看出，虚拟采样器和保持器可以不止一次地使用。在必要的地方加入，可为求脉冲传递函数带来方便。但须注意，由于保持器的频谱特性并非理想矩形，所以每加一次采样器和保持器都会带来误差。因此(8-

5、3)式较之表8-1中的式子误差要大些，所以要尽量减少虚拟采样器和保持器的使用。2021/10/37下面举例说明其使用方法。若:则根据表8-1，当加零阶保持器时，可得:所以得差分方程为:也可根据(8-3)式来求G(z)：其对应差分方程为:（8-4）（8-5）2021/10/38以上(8-4)式或(8-5)式是按一次和二次加入虚拟采样器和保持器时分别求得的差分方程。后者较前者误差要大些。有了上面的知识，不难将其推广应用到结构图表示的系统，求其等效的差分方程。在此不再多叙。2021/10/398.2状态方程的离散化假设连续

6、系统的状态方程为:若人为地在系统的输入端及输出端加上采样开关，同时为了使输入信号复原为原来的信号，在输入端还要加一个保持器，如图所示。保持器TTux假定为零阶保持器:——输入向量的所有分量在任意两个依次相连的采样瞬时为常值。例如:对第n个采样周期u(t)=u(nT)。T为采样间隔(周期)（8-6）图8-4采样控制系统结构图2021/10/310若对方程(8-6)式两边进行拉普拉斯变换，得：故对(8-7)式反变换可得：即：以(sI-A)-1左乘上式的两边可得：（8-7）考虑到：（8-8）（8-9）此为(8-6)式的连续

7、解，由此可推导出系统的离散解。注意：图8-4所示系统的采样开关和保持器同样是为了将(8-6)式离散化而虚构的。2021/10/311根据上式，n及n+1两个相连的采样瞬间，有：将(8-11)式减去(8-10)式后乘以eAT，得：（8-10）（8-11）（8-12）将(8-12)式右边积分进行变量代换，即令：则得：（8-13）（8-14）求解方程右边的卷积积分困难。2021/10/312但由图8-4可知：若系统采用零阶保持器时，则两个采样点之间输入量可看做常数，即u(nT+t)=u(nT)，这样(8-14)式可写为：

8、式中：（8-15）（8-16）（8-15）是在假定输入量U(t)在两个采样时刻之间保持不变（即如图8-5中矩形近似）的前提下导出的。而实际输入量在两个采样时刻之间是变化的，这样近似会引进误差。为减小误差，可假定两个采样时刻之间输入量U(t)为一斜坡函数，即如图8-5中梯形近似。u(t)t0kT(k+1)TT实际值梯形近似矩形近似图52021/1