计算机仿验3数值积分算法.doc

实验三 利用数值积分算法的仿真实验 学号： 姓名： 学院： 一． 实验目的 熟悉MATLAB的工作环境； 掌握在MATLAB命令窗口调试运行程序； 掌握M文件编写规则及在MATLAB命令窗口运行程序； 掌握利用欧拉法、梯形法、二阶显式Adams法及四阶龙格库塔法构建系统仿真模型的方法,并对仿真结果进行分析。 二．实验内容 电路如图１所示电路进行仿真试验。电路元件参数： SKIPIF 1 0 ， SKIPIF 1 0 ， SKIPIF 1 0 ， SKIPIF 1 0 。电路元件初始值： SKIPIF 1 0 ， SKIPIF 1 0 。系统输出量为电容电压 SKIPIF 1 0 。 SKIPIF 1 0 三． 实验步骤 1. 求连续系统传递函数 根据所示电路图，我们利用电路原理建立系统的传递函数模型，根据系统的传递函数是在零初始条件下输出量的拉普拉斯变换与输入量的拉普拉斯变换之比，可得该系统的传递函数： SKIPIF 1 0 2. 离散系统仿真模型 在连续系统的数字仿真算法中,较常用的有欧拉法、梯形法、二阶显式Adams法及显式四阶Runge-Kutta法等。欧拉法、梯形法和二阶显式Adams法是利用离散相似原理构造的仿真算法，而显式四阶Runge-Kutta法是利用Taylor级数匹配原理构造的仿真算法。对于线性系统，其状态方程表达式为: SKIPIF 1 0 SKIPIF 1 0 式（3-1）中， SKIPIF 1 0 是系统的n维状态向量, SKIPIF 1 0 是系统的m维输入向量， SKIPIF 1 0 是系统的r维输出向量。A为 SKIPIF 1 0 阶参数矩阵，又称动态矩阵，B为 SKIPIF 1 0 阶输入矩阵，C为 SKIPIF 1 0 阶输出矩阵，D为 SKIPIF 1 0 阶交联矩阵。根据图所示电路，系统状态方程模型: SKIPIF 1 0 式中，状态变量 SKIPIF 1 0 ，输出变量 SKIPIF 1 0 ，系数矩阵为： SKIPIF 1 0 ， SKIPIF 1 0 ， SKIPIF 1 0 。 欧拉法 利用前向欧拉法构建线性系统的仿真模型为： SKIPIF 1 0 式中， SKIPIF 1 0 为积分步长， SKIPIF 1 0 为单位矩阵。利用后向欧拉法构建线性系统的仿真模型为： SKIPIF 1 0 对于前向欧拉法，系数矩阵为： SKIPIF 1 0 ， SKIPIF 1 0 ， SKIPIF 1 0 ，D=0。对于后向欧拉法，系数矩阵为： SKIPIF 1 0 ， SKIPIF 1 0 ， SKIPIF 1 0 。 梯形法 利用梯形法构建线性系统的仿真模型为： SKIPIF 1 0 对图所示的系统，利用梯形法构造的系统差分方程具有形式： SKIPIF 1 0 其系数矩阵为： SKIPIF 1 0 ，， SKIPIF 1 0 ， SKIPIF 1 0 ，D = 0。 （3）二阶显式Adams法 利用二阶显式Adams法构建线性系统的仿真模型为： SKIPIF 1 0 式中： SKIPIF 1 0 二阶显式Adams法为多步计算方法，利用多步计算方法对系统进行仿真时，需要与之具有相同计算精度的单步计算方法辅助计算。二阶显式Adams法的计算精度为二阶，可以采用梯形法或改进的Euler法等辅助计算。利用改进的Euler法构建线性系统的仿真模型为： SKIPIF 1 0 其中， SKIPIF 1 0 。由式计算出 SKIPIF 1 0 和 SKIPIF 1 0 后，便可以转入由二阶显式Adams法构造的离散系统模型计算，即系统差分方程。其计算方程为： SK