系统时间响应及其仿真

《系统时间响应及其仿真》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第七章系统时间响应及其仿真本章主要内容如下：7.1仿真算法7.2系统仿真的MATLAB函数系统的时间响应是指系统在输入信号或初始状态作用下，系统输出随时间变化的情况。系统的时间响应反映了系统的特征和性能，如快速性、稳定性等。对系统时间响应的分析是我们设计、校正系统的基础。第七章系统时间响应及其仿真7.1仿真算法对系统的时间响应进行动态仿真，采用什么样的仿真算法是一个至关重要的问题。对连续时间系统进行数字动态仿真，主要是两种方法：基于数值积分的仿真方法;基于离散相似法的仿真方法。由于后者涉及到离散控制系统理论，因此本节重点介绍基于数字积分的连续系统仿真方法。

2、7.1仿真算法7.1.1数值求解的基本概念设微分方程为则求解方程中函数y(t)问题，就是已知初值的常微分方程求解问题。所谓数值求解就是要在时间区间[a,b]中取若干离散点求出微分方程在这些时刻的近似值常微分方程数值求解的基本方法是数值积分法。7.1仿真算法7.1.2数值积分的基本原理积分区间的划分将区间[a,b]分成N个小区间,时间间隔h（)也称为积分步长,在第k个间隔t=[tk,tk+1]内积分：则可用yk（k=0,1,…N）作为解y(t)的近似值，如图所示。abtky0ykyt数值积分图解tk+17.1仿真算法7.1.2数值积分的基本原理数值积分的展开

3、式为避免(2)式中的积分项，将y在tk，以h为增量展开成Taylor级数：式(3)是一个递推公式。积分值与实际微分方程解的误差取决于步长h和计算所用的阶数，它是数值积分的基础。7.1仿真算法7.1.2数值积分的基本原理有关概念单步法和多步法单步法指计算yk+1值只需利用tk时刻的信息，也称为自启动算法；多步法在计算yk+1值时，则需利用tk，tk-1，…时刻的信息。显示法和隐式法显示法在计算yk+1时所需数据均已算出；隐式法在计算yk+1时需用到tk+1时刻的数据，该算法必须借助予估公式。定步长和变步长定步长为积分步长在仿真运行过程中始终不变；变步长指在仿

4、真运行过程中自动修改步长。7.1仿真算法7.1.3数值积分的几个算法欧拉算法在(3)式中取前两项：可得欧拉算法：t0t1t2t3hy(t)y0y1yt欧拉近似解欧拉法【说明】欧拉法是用一条过各点的切线取代曲线来逼近精确解。该算法简单，计算量小，但精度较低。7.1仿真算法7.1.3数值积分的几个算法梯度法梯度法是欧拉法的改进。与欧拉法相比，梯度法是用两个点(tm，ym)、(tm+1,ym+1)的斜率的平均值来确定下一点的y值。由于上式计算时需要用到ym+1的值，而ym+1不能预先知道，故梯度法需要和欧拉法结合使用,即用欧拉法对ym+1进行予估，再由梯度法计算

5、ym+17.1仿真算法7.1.3数值积分的几个算法龙格-库塔法龙格-库塔法的基本思想欧拉算法的精度较低，主要是其微分方程解y的Taylor展开式所取的项数太少。显然为了提高计算精度，应当取泰勒公式（3）更高阶项。虽然增加高阶项可提高计算精度，但也同时带来了需要计算高阶导数的困难。龙格-库塔法的关键是利用低阶导数构成的曲线去拟合含有高阶导数的曲线，从而避免了计算高阶导数的问题。7.1仿真算法7.1.3数值积分的几个算法龙格-库塔法二阶龙格-库塔（RK）法取(3)式的前三项，则有设原微分方程(1)式解具有以下形式:式中，a1，a2，b1，b2为待定系数。将(8

6、)式中K2按二元函数展开成泰勒级数，并取前三项将K1，K2代入(8)式：比较(6-10)、(6-7)式：显然由(11)式并不能唯一确定a1,a2,b1,b2，因为只有三个方程。因此对于同一种算法可以有不同的表现形式。【说明】由于该算法只取到泰勒展开式的二阶导数项，所以称为二阶龙格-库塔法。但由(8)~(12)式可知，算法并没有用y的二阶导数。若设a1=a2，则即二阶RK法公式为7.1仿真算法7.1.3数值积分的几个算法龙格-库塔法龙格-库塔（RK）法的一般形式式中，i为待定权系数，ai,bij为待定系数，r为使用Ki的个数(即级数)，Ki为所取各点导数f

7、的值。Ki的个数与yk+1泰勒展开式所取的项数有关(即RK算法的阶数)，同时还与计算区间内所取导数值的点数有关。7.1仿真算法7.1.3数值积分的几个算法龙格-库塔法四阶RK公式四阶RK公式用到了y的泰勒展开式的四阶导数。在RK算法的一般公式（13）中，取r=4可得：由于(14)式在同级的RK算法中，计算精度较高，计算量较少，而在系统仿真的数值积分中应用十分广泛。称之为四阶四级RK公式。7.1仿真算法7.1.3数值积分的几个算法Gear算法“病态”常微分方程(刚性方程)的系数矩阵A的特征值具有如下特征：则称为“病态”方程。7.1仿真算法7.1.3数值积分的

8、几个算法Gear算法控制系统仿真中的“病态”问题病态系统中绝对值最