第3章 节 系统仿真算法《电力电子电路仿真——MATLAB和PSpice应用》课件.pptVIP

第3章系统仿真算法 第3章系统仿真算法 3.1 数值积分法 在数学仿真中，从一个实际系统抽象出数学模型只是第一步，这一步将实际系统变成了数学模型，可以称之为系统建模或系统辨识，这是第一次模型化过程。这次模型化所得到的只是数学方程式，必须使用一定的仿真工具才能求解。将已获得的数学模型变成能在一定仿真工具中运算求解的仿真模型，这是第二次模型化过程。 数值解法 3.1 数值积分法（针对连续系统的微分方程形式。利用数值积分方法对常微分方程（组）建立离散化形式的数学模型——差分方程，并求其数值解 ） 3.1.1 数值积分法基本原理 把被仿真系统表示成一阶微分方程组或状态方程的形式 （3-1） 则上式在 在形式上的连续解为 通过上式的演化，连续系统的数值解就转化为相邻两个时间点上的数值积分问题。  3.1 数值积分法 所谓数值解法，就是寻求初值问题式（3-1）的真解在一系列离散点 上的近似数值解 相邻两个时间离散点的间隔hn = tn+1 - tn，称为计算步距或步长，通常取hn = h为定值。 数值积分法的主要问题归结为对函数F(t，Y)的数值积分问题 1. Euler法 对于初值问题 假定y(t)为其解析解，将展成泰勒级数： 则 写成差分方程的形式：  3.1 数值积分法 只要给定初始条件y0及步长h，就可根据f(t0，y0)算出y1，再由y1，算出y2，如此递推算出y3，y4，y5…。 【例3-1】 试用Euler法求其数值解。  3.1 数值积分法 将这种思想引申如果在每个积分步中多取几个点（如取r个点），分别求出其斜率 k1，k2，…，kr，然后取不同的权值，则公式为： 得出一系列龙格 — 库塔（Runge-Kutta）法 。 3. 几个基本概念 （1）算法自启动 （2）单步法与多步法 （3）显式与隐式 （4）截断误差 （5）舍入误差 （6）初始误差  3.1 数值积分法 3.1.2 龙格－库塔（Runge-Kutta）积分法 是用几个点上的的一阶导函数值的线性组合来近似代替在某一点的各阶导数，然后用Taylor级数展开式确定线性组合中各加权系数。这就是Runge-Kutta（简称RK）法的基本思想。 各阶RK公式及说明如下： 一阶RK公式为： 它就是Euler公式，也就是说，Euler公式是RK公式的特例 二阶RK法是每步取两个斜率加权。 第一斜率 第二斜率 加权后递推公式为  3.1 数值积分法 经过与泰勒级数对比可确定四个参数αβω1ω2,可得二阶RK法的计算公式： 依次取1/2,1/2,0,1时为 依次取1,1,1/2,1/2时为  3.1 数值积分法 三阶RK法公式： 四阶RK法公式：  3.1 数值积分法 四阶RK法程序框图  3.1 数值积分法 【例3-2】已知系统方程 取步长 试用RK4公式计算 时y的值。 解：将原系统方程转化成下列方程组： （1）计算所有变量的第一个RK系数 （ ， ）  3.1 数值积分法 （2）计算所有变量的第二个RK系数 （ ， ） （3）计算所有变量的第三个RK系数  3.1 数值积分法 （4）计算所有变量的第四个RK系数 （5）计算时的近似值 用同样的方法，在y11，y21的基础上可以求得y12和y22的值。即为原系统中在t=0.2时的近似值。  3.1 数值积分法 3.1.3 线性多步法 计算n+1时刻的值yn+1时，不只使用第n时刻的值yn和fn,而是充分利用一系列历史时刻的数据来求解yn+1,则可望既加快仿真速度又能获得较高的仿真精度。这就是多步法的具体出发点。线性多步法中以Adams方法