![控制系统计算机仿真](https://wfqqreader-1252317822.image.myqcloud.com/cover/532/678532/b_678532.jpg)
1.5 控制系统仿真的数学模型
从仿真的定义可知,建立数学模型是开展仿真的先决条件,因此在介绍各种仿真方法之前,必须先了解系统的数学模型。在自动控制中,常见的数学模型可以分为连续时间模型、离散时间模型和采样数据模型3大类。
1.5.1 连续时间系统的数学模型
连续时间模型通常采用下列3种数学形式描述连续系统。
1.输入、输出的微分方程
设系统的输入为u(t),输出为y(t),它们之间的关系可以表示为高阶微分方程
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0023_0001.jpg?sign=1739252861-AmS7AefYBk4JYCi1g9W5nRW85z6M9JOS-0-5008c03adcaecb8d4589226db2487aaf)
对于自动控制中最常见的线性定常系统,式(1.10)可以表示为
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0023_0002.jpg?sign=1739252861-1TlWxiKXoXZSipe32EYYfoa8RNv2gr7N-0-3ab21eba1c968b6fcad7f72699180e6f)
式中,ɑ0=1。
若引进微分算子,式(1.11)可以表示为
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0023_0004.jpg?sign=1739252861-J7ZsEFlsyPcVonjRgPdsi7IHTtDC85lV-0-5ee4da2a326b88f323d9f6d5b21d3f78)
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0023_0005.jpg?sign=1739252861-WJ2E3ZlG2cQ7O2QTWmZDuwVlIhufn2M7-0-bc92241ffdc4e76ff0b3d102e41eaf52)
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0023_0006.jpg?sign=1739252861-icYwY4oJLwTe5jYaW6JquatwNH4DFOrU-0-c0559b39f6046b00131efdb9e920f1a0)
则有
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0024_0001.jpg?sign=1739252861-m3197EhSLp6syT3Uzpn8Rwuu1w3l9nHx-0-45abd606b1e58ce6bb48dde3814fd612)
2.输入、输出的传递函数
对式(1.11)的两边取拉氏变换,若系统的输出y(t)和输入u(t)及其各阶导数的初始值都为零,则得到
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0024_0002.jpg?sign=1739252861-NOl0n1PslgZSoKv2jqjpnQamizvsMJO0-0-152236ed75c933f060c26ea9d899e23d)
式中,Y(s)和U(s)分别为y(t)和u(t)的拉氏变换。定义
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0024_0003.jpg?sign=1739252861-GfROMNPi6q3fEF4cYg9cJZuQrFz6xiUy-0-5fff602e84ffb3e0d0692048a93051f2)
称它为输入、输出的传递函数。由式(1.14),有
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0024_0004.jpg?sign=1739252861-1QMZt4I1QDn535hpKL4lgmylvwvM4ic0-0-887ad29e514fe9511c7f45bdf9cfd84e)
将式(1.13)与式(1.15)对比可知,在初值为零的情况下,p与s等价。
3.状态空间模型
式(1.11)和式(1.15)给出的数学模型仅仅描述了线性定常系统的外部特性,即仅仅确定了输入u(t)与输出y(t)之间的关系,故称为系统外部模型。
为了描述一个系统内部的特性,即确定组成系统的各个实体之间的相互作用而引起实体属性的变化,可以引进系统的内部变量——称为状态变量。实际系统中,真实的内部变量及数学上定义的内部变量可以是一致的,也可以是不一致的。例如,直流电动机的输入为电枢电压u(t),输出为转速ω(t),它的内部状态变量可以选为真实的物理量——电枢电流i(t),也可以选为其他变量。然而,无论状态变量如何选择,系统的数学描述总有以下形式
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0024_0005.jpg?sign=1739252861-PfGogng50g9eTMwWNhU3FdB2fTxVtGxF-0-af6332982eb551f34929cce23e654ed8)
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0024_0006.jpg?sign=1739252861-vZuhK06Mm6QpI7rpcE3GXzBXSktkaA20-0-bfdca68afdb758c5df8b4aa29be8f361)
其中,式(1.16)称为系统的状态方程,式(1.17)称为系统的输出方程。
式(1.16)和式(1.17)合称为系统的状态空间模型。由于它描述了系统的内部特性,所以又称为内部模型。
对于线性系统,上述3种数学模型的表示形式可以相互转换。
1.5.2 离散时间系统的数学模型
与1.5.1节中叙述的连续时间模型的表示形式相对应,离散时间模型通常也采用3种表示形式描述离散时间系统。
1.输入、输出的差分方程
设在离散时间点t0,t1,…,tk,…上控制系统的输入序列为{u(k)},输出序列为{y(k)},它们之间的关系可以表示为
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0001.jpg?sign=1739252861-iYa6WaeYRZKulY0RcodrOFmGpx1IpSFM-0-b3dabf261490d2f8b93a8d3fe285287a)
这是非线性差分方程的表达式,其中f是其变量的函数。如果可以从式(1.18)中解出y(n+k),则得到非线性递推形式的差分方程为
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0002.jpg?sign=1739252861-jHaiMQH1DCcdCmIta8toJmlzat1FBGPO-0-90fd1db0ef2d91a984ad077b9a02e583)
差分方程描述了离散控制系统在各个时刻输入、输出之间的相互关系。由式(1.18)可以看出,根据输入u(k),u(k+1),…,u(n+k)和以前的系统输出y(k),y(k+1),…,y(n+k-1),就可以递推出时刻tn+k处的输出值y(n+k)。
对于线性定常系统,式(1.18)可以表示为
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0003.jpg?sign=1739252861-UQx496eNsgpfwEWxWNL86EwX00A1Dg9p-0-0d825b12e686e651399b10bcb94d20d6)
若引进后移算子q-1,它定义为
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0004.jpg?sign=1739252861-NppMGlidRNEbyOG7lhL2A6GiDn4xB9uu-0-ae651491f6b80cbd05690292fa6c4a22)
则式(1.20)可以改写为
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0005.jpg?sign=1739252861-iMbGY5ZFDTeKgzWTucG3pfYaMYOJvW0R-0-bbcd8902d1d36993e354b0abe666058f)
式中,ɑ0=1。即有
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0006.jpg?sign=1739252861-mh6RNbIbcxgGE75yaTpGRkll4xs9argC-0-e50c23b44d1780d278232582440954aa)
式中
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0007.jpg?sign=1739252861-ZHJlIMMhcp5sMg7odwFl9ZB3qoOiPNdO-0-47ffca3653da1e74b4b6bffc5650e2fb)
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0008.jpg?sign=1739252861-3o6Vxgm8KO2JToPL2DCuhdPKRyb2cUgV-0-b8db80e1a0f826f92358fb86d6c83c33)
从而有
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0009.jpg?sign=1739252861-k6A2OuYFzaKJKbbxK0PzbGFKINZwqXdA-0-22b15a7612c0b5a55f5b772b54f5b487)
或者
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0010.jpg?sign=1739252861-9PSbVm3XYBHQ1Skk97T8zglCTlHtTjXl-0-61ad583c55e8676f27941871cdc03f54)
2.输入、输出的脉冲传递函数
对式(1.20)的两边取Z变换,并设系统的初始值及以前时刻的输入和输出都为零,即y(k)=u(k)=0,k≤0,则可得
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0011.jpg?sign=1739252861-r0DhchdjBAFBgpJWGB4sdk5HYx9ssak5-0-0d653f1354e5d05912829d98e4832412)
式中,Y(z)和U(z)分别为序列{y(k)}和序列{u(k)}的Z变换。定义
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0025_0012.jpg?sign=1739252861-F7N4s2m0FnBlgURLntpW0dn52DFpd3DW-0-bc405eb5cd8dd8399c365ef85d152ef4)
为线性差分方程式(1.20)描述的系统的脉冲传递函数,则有
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0026_0001.jpg?sign=1739252861-Hyx3v3uRpp0rQST0yoD7QyeA9FCmaygI-0-89ea8376dafe67b4596fed6f4b408ff7)
将式(1.23)与式(1.25)对比可知,在初始条件为零的情况下,q-1与z-1等价。
3.离散状态空间模型
式(1.20)和式(1.25)仅仅描述了线性定常离散系统的输入、输出之间的关系,故称为外部模型。如果引进系统的状态变量序列{x(k)},则可以构成离散状态空间模型
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0026_0002.jpg?sign=1739252861-hbC91sXDDQrbwteM6WRTjiQRr4FCDbSs-0-361329466bf01aa22b333790f8f007f8)
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0026_0003.jpg?sign=1739252861-TkcNZ9MBLhKvHhTnA5Ej8rOGDFiwxD4m-0-b95e703405eb7be71e4ba973d1c74cb5)
其中,式(1.26)称为系统的离散状态方程,式(1.27)称为系统的离散输出方程。离散状态空间模型属于内部模型。
对于线性定常系统,上述3种形式可以互相转换。
1.5.3 采样控制系统的数学模型
在采样控制系统中,控制器(由数字运算器构成)是对离散时间信号进行运算的部件,描述它的数学模型是离散时间模型;被控对象是连续过程,用连续时间模型描述。因此,这类系统需要用连续-离散时间混合模型,即采样数据模型来描述。计算机控制系统是典型的采样数据系统。如果忽略A/D和D/A转换器的转换时间及误差,则采样数据系统可由图1.8表示。
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0026_0004.jpg?sign=1739252861-GoL3oY5DHp5z7EOMDIIohCQ8UHpMjwNc-0-18ebe95845eafe76cb4395208e5170fb)
图1.8 采样数据系统结构
采样控制系统中,连续时间模型和离散时间模型都可以由任意一种形式给出。如果D(z)是控制器的脉冲传递函数,Gh(s)和G(s)分别是保持器和被控对象的传递函数,则采用传递函数形式描述的采样控制系统可由图1.9表示。
![](https://epubservercos.yuewen.com/DB7D15/3590428903440801/epubprivate/OEBPS/Images/figure_0026_0005.jpg?sign=1739252861-WE6bEzUAlnm4uEsUQPX0Dga8Fql6A8Js-0-86768213681c068fc8cc380a0094a0e4)
图1.9 采样控制系统的传递函数表示