计算机仿真技术实验报告-实验三[参照].pdfVIP

《仿真技术与应用》实验报告 计算机仿真技术实验报告 实验三 利用数值积分算法的仿真实验 - I - 《仿真技术与应用》实验报告 实验三 利用数值积分算法的仿真实验 一．实验目的 1) 熟悉 MATLAB 的工作环境； 2) 掌握 MATLAB 的 .M 文件编写规则，并在命令窗口调试和运行程序； 3) 掌握利用欧拉法、梯形法、二阶显式 Adams 法及四阶龙格库塔法构建系统 仿真模型的方法,并对仿真结果进行分析。 二．实验内容 系统电路如图 2.1 所示。电路元件参数：直流电压源E  1V ，电阻R  10 ，电感 L  0.01H ，电容C  1F 。电路元件初始值：电感电流iL (0)  0A ，电容电压 u (0)  0V 。系统输出量为电容电压u (t) 。 c c 连续系统输出响应u (t) 的解析解为： c u (t)  U (1eat (cost sin t a / )) （2-1 ） c s 2 R 1  R  其中，a  ，    。 2L L C 2L  R L iL (t)  u (t) DC E C c  图2.1 RLC 串联电路 三、要求 1）利用欧拉法、梯形法、二阶显式 Adams 法及显式四阶 Runge-Kutta 法构建系统仿真模 型，并求出离散系统的输出量响应曲线； - 2 - 《仿真技术与应用》实验报告 2 ）对比分析利用欧拉法、梯形法、二阶显式 Adams 法及显式四阶 Runge-Kutta 法构建系 统仿真模型的仿真精度与模型运行的稳定性问题； 3 ）分别编写欧拉法、梯形法、二阶显式 Adams 法及显式四阶 Runge-Kutta 法的.m 函数 文件，并存入磁盘中。.m 函数文件要求输入参数为系统状态方程的系数矩阵、仿真时间 及仿真步长。编写.m 命令文件，在该命令文件中调用已经编写完成的上述.m 函数文件， 完成仿真实验； 4 ）subplot 和 plot 函数将输出结果画在同一个窗口中，每个子图加上对应的标题。 四.实验原理 （1）连续系统解析解 连续系统输出响应u (t) 的解析解为： c u (t)  U (1eat (cost sin t x / )) c s 2 R 1  R  其中，x  ，    2L