快速数字仿真算法优秀讲义.pptVIP

2.4 快速数字仿真算法 * 替换法 根匹配法 数值积分算法和离散相似算法的计算量比较大，因此仿真速度受到了一些限制。在一些要求尽可能快的获得仿真结果的场合，需要一些快速的数字仿真算法。 算法的精度和速度是一对矛盾，以前主要考虑的是精度，而在快速数字仿真算法中，考虑在满足计算稳定性及工程精度要求的情况下，尽可能提高仿真计算速度。 * * 一般快速数字仿真算法有如下两点基本要求： 每步计算量要小 算法具有良好的稳定性，允许采用较大的步长，同时又能保证必要的计算精度 为了提高速度，要尽量减少仿真中附加信息的计算。例如，如果只对输出感兴趣，就没必要将系统转换成状态空间描述，直接转换成脉冲传递函数有利于快速仿真。 2.4.2 替换法 相匹配的含义：指原系统的数学模型与仿真模型的动态和静态特性是一致的。其含义是：若被仿真系统的数学模型是稳定的，则其仿真模型也是稳定的，且动态、静态特性一致，即若G(s)是稳定的，则G(z)也应是稳定的。对于同一输入信号，由G(z)求得的输出y(kT)与由G(s)求得的y(t)应具有相一致的频率特性。 思路：如果能把一个高阶连续系统数学模型G(s)转换成与之相匹配、每步计算量较小、允许采用较大步长、具有合理精度的仿真模型G(z)，那么就可利用G(z)求得的差分方程进行数字仿真，以达到快速仿真的效果。 1 双线性变换法 连续系统的传递函数G(s)转化成脉冲传递函数G(z)一般使用部分分式展开或者级数求和法： * 但是对于高阶系统，这两种方法都比较困难。 替换为 线性映射 思路： 双线性变换法： 将s平面上的传递函数转换为z平面上的脉冲传递函数 双线性变换算法的稳定性 * S平面 z平面 因此如果原传递函数G(s)是稳定的，则通过双线性变换后得到的脉冲传递函数G(z)也必然是稳定的。 * 其次，双线性变换后的稳态增益不变。 设连续系统的传递函数为 　 　 其稳态增益为bn/an 。若将G(s)进行双线性变换，有 　　　　 显然，G(z)的稳态增益仍然是bn/an。 综上所述，双线性变换法符合相匹配原理。 双线性变换法适合于以分子分母多项式形式描述的G(s)。 * 例：已知连续系统的传递函数为 试采用双线性变换法求出相应的脉冲传递函数和差分方程，并分析脉冲传递函数的稳定性。 根据 得脉冲传递函数 离散系统特征方程的特征根的模为 因此脉冲传递函数G(z)是稳定的 * 双线性变换法步骤 设连续系统的传递函数为 　 　　 　 采用双线性变换法求取仿真模型的步骤如下： 将G(s)中的s替换得到G(z)，并整理成有理分式形式 　　 　　 　　 　　 根据得到的G(z)，由Y(z)=G(z)U(z)，两边取Z反变换，得到便于计算机递推计算的差分方程: 根匹配法的基本思想是要使离散化模型的瞬态特性和稳态特性与原连续系统模型的对应特性保持一致，也就是使离散化后所得脉冲传递函数的零点和极点与原连续系统传递函数的零点和极点相匹配。所以这种算法又被称为零极点匹配法。 离散系统的瞬态特性和稳态特性完全由脉冲传递函数G(z)的零、极点位置及增益大小决定。因此直接由G(s)的零点和极点确定与之对应的G(z)的零点和极点，再根据一定的原则确定G(z)的增益，则可很方便的由G(s)求出G(z)，使二者的瞬态特性和稳态特性保持一致。因此根匹配法适合于以零极点形式描述的G(s)。 * 2.4.4 根匹配法 * 使用根匹配法对连续系统离散化的步骤为： 1、将连续系统传递函数G(s)写成零极点形式 2、利用式 在z平面上一一对应的确定零极点的位置 3、在G(z)的分子上附加(n-m)个零点，使G(z)的分子多项式和分母多项式的阶次相同（在s平面的无穷远处还存在n-m个零点） * 4. 根据G(s)的特征确定Kz，一般有两种方法： 一种是选择适当的输入信号，应用终值定理，分别确定连续系统模型G(s)和离散化仿真模型G(z)的终值，然后根据终值（不等于零或无穷大）相等的原则确定G(z)的增益Kz。 另一种是使由G(s)和G(z)分别得到的频率特性在某一临界频率处相同，从而确定Kz。一般对于低通滤波器应使其低频增益相同，而对于带通滤波器则应使中频增益相同。 5. 根据得到的G(z)，由Y(z)=G(z)U(z)，两边取Z反变换，得到便于计算机递推计算的差分