资源描述：
《哈工大计算机仿真.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第3章连续系统的数字仿真通用算法3.1基于离散相似原理的数字仿真算法连续时间系统的状态方程形式的数学模型可以写成如下标准形式：(1)式(1)右端的为:对系统进行时域仿真分析，实际上就是要求解这类方程的“初值问题”，即求方程(2)(3)的解。求初值问题式(3)的解析解，意即求一个x(t),使之满足方程(3)。仿真是要利用计算机求数值解，求初值问题的数值解,意即将连续时间区间(0,t)离散为时间序列t=t1,t2,t3,…tn。然后求真实解在上述时刻的近似值，x1,x2,x3,…xn。事实上只有一些特殊类型的
2、状态方程才存在着解析解。大量从实际系统建立的状态方程的求解，主要依靠数值解。因此，在连续系统数值仿真中，首先应建立适于数值计算的数字仿真模型。离散相似原理将连续系统通过虚拟采样开关和信号保持器转换成离散系统，建立适用于数值计算的数值仿真模型。应用离散相似原理建立仿真模型时，根据对输入信号和信号保持器的不同选择，可以构造各种各样的算法，即相应的数字仿真模型。作为构造算法的例子，下面介绍欧拉法、梯形数值积分法、Adams方法及局部解析方法。3.1.1欧拉法(EulerMethod)(4)为了讨论方便，这里以一
3、阶连续系统为例。其状态方程为：在区间[tm,tm+1]上输出x(t)为：(5)为采样值Fi将代入式(5)就可以得到x(t)在区间[tm,tm+1]上的解析解。构成的简单函数。在区间[tm,tm+1]上，取信号保持器为零阶保持器，设其输出表达式为：(6)其中式(5)x(t)的解为:(7)当t=tm+1时，有即(9)(8)式(9)称作前向欧拉法的迭代公式，又称作显式欧拉法迭代公式。从式(9)可以看出，只要给出初始值x0就能一步一步地计算出变量x(t)在任意采样时间点的近似值。我们称这种迭代公式为自起步迭代公式
4、。当取零阶保持器的输出为：(10)式(11)称作后向欧拉法的迭代公式，又称作隐式欧拉法迭代公式。式(9)和式(11)两种迭代公式在计算xm+1的方法上有很大的差别,前者只需计算公式右端的函数f(x,u(t),t)在t=tm时的值便可得出xm+1，而后者则必须通过解代数方程才能求得xm+1。尽管后向欧拉法与前项欧拉法在计算xm+1的方法上有很大差别，而后向欧拉法在计算xm+1时，也只用到x(t)的一个过去值xm，因此只要给定初始值x0迭代就能进行，故此法也是自动起步的。欧拉法属于单步计算方法。则可得到xm的
5、迭代公式为：(11)图2前向欧拉法的几何意义图3后向欧拉法的几何意义几何意义：前向欧拉法的几何意义如图2所示。在区间(tm,tm+1)内，把真实解x(t)用一条直线来近似替代。该直线通过点xm，且斜率等于曲线x(t)在t=tm点的斜率。后向欧拉法的几何意义如图3所示。在区间(tm,tm+1)内，把真实解x(t)用一条直线来近似替代。该直线通过点xm，且斜率等于曲线x(t)在t=tm+1点的斜率。例:设系统方程为其数值解，取仿真步距。试用Euler法求解：原方程为则前向欧拉法的递推公式为后向欧拉法的递推公式
6、为显然，由后向欧拉法得到的递推公式为非线性代数方程。解此非线性方程可得t00.10.20.30.40.5…1.0yn10.91610.84470.78330.73000.6834…0.51653.1.2梯形法为了讨论方便，这里以一阶连续系统为例。其状态方程为：在区间[tm,tm+1]上输出x(t)为：取信号保持器为一阶保持器。一阶保持器是区间[tm,tm+1]上的一次多项式。其输出表达式为：其中a1和a0为常系数。(12)对一阶隐式保持器，曲线的插值点为(tm,Fm)和点(tm+1,Fm+1)因此，式(1
7、2)中的系数a0和a1可由下列方程确定：解方程(13)可得系数a0和a1为：(13)(14)将式(12)代入式(5)，则在区间[tm,tm+1]上的解为：其中h=tm+1-tm(15)当时，有(16)即:(17)梯形法的递推公式利用梯形法递推公式对系统进行仿真时，必须通过解代数方程才能求得xm+1。由于梯形法的推公式在计算xm+1时，只用到x(t)的过去值xm及，因此只要给定初始值x0递推就能进行。该方法为单步法，可以自起步计算。图4梯形法的几何意义几何意义：梯形法的几何意义如图4所示。在区间[tm,tm
8、+1]上，把真实解x(t)用折线来近似替代。段斜率为段的斜率为显然，梯形法较欧拉法具有更高的精度。，解：原方程为，例对例3.1给出的系统，试用梯形法求其数值解（取仿真步距h=0.1）则梯形法的递推公式为由梯形法得到的递推公式为非线性代数方程。解此非线性代数方程可得：t00.10.20.30.40.5…1.0yn10.90870.83280.76850.71360.6659…0.4994Euler法和梯形法的另一种解释——数值积
