计算机仿真技术实验三(范文).docVIP

计算机仿真技术实验三(范文) 第 PAGE \* Arabic 1 页第 PAGE \* Arabic 1 页 计算机仿真技术实验 实验报告 实验三 利用欧拉法、梯形数法和二阶显示Adams法对RLC串联电路的仿真 1.1实验目的 熟悉MATLAB的工作环境； 掌握在MATLAB命令窗口调试运行程序； 掌握M文件编写规则及在MATLAB命令窗口运行程序； 4）掌握利用数值积分法构造离散模型的方法。 1.2实验内容 电路如图１所示电路进行仿真试验。元件参数：，，，。初始值：，。输出量电容电压。 利用欧拉法、梯形数法和二阶显示Adams法建立图１电路的离散数需模型； 建立计算机仿真模型； 选择一组离散时间间隔值，进行仿真试验； 分析仿真结果，从仿真模型实现的难易性、模型的稳定性、模型的精度及离散时间间隔等方面，对比分析上述方法构造的离散系统模型的优缺点。 1.3实验原理与方法 1.3.1系统的数学模型 连续系统模型： 微分方程 传递函数 1.3.2系统的仿真模型 离散系统模型： 1、欧拉法 (1)、前向欧拉法 = 1 \* GB3 ① = 2 \* GB3 ② (2)、后向欧拉法 = 1 \* GB3 ① = 2 \* GB3 ② 2、梯形法 = 1 \* GB3 ① = 2 \* GB3 ② 3、二阶显示Adams = 1 \* GB3 ① = 2 \* GB3 ② 需要用梯形法启动 1.4实验步骤 根据所得的系统仿真差分方程模型，编写程序，绘出图像得到仿真结果 1.4.1前向欧拉法 function []=RLC(R,L,C,Us,t,h) R=10; L=0.01; C=1.0e-6; Us=1; t=0.01; h = 5.0e-5; NNN = fix(t/h); AA = [-R/L -1/L;1/C 0]; BB = [1/L;0]; CC = [0 1]; Nr = 2; for i=1:1:Nr xx(1:Nr,1) = 0; end for k=1:NNN xx(1:Nr,k+1) = xx(1:Nr,k) + (AA* xx(1:Nr,k)+BB)*h; end 1.4.2后向欧拉法（未声明变量同上1.4.1） for i=1:1:Nr xxx(1:Nr,1) = 0; end EE = [1 0;0 1]; AA1 = inv(EE-AA*h); for k=1:NNN xxx(1:Nr,k+1) = AA1*(xxx(1:Nr,k) + BB*h); 1.4.3梯形法 for i=1:1:Nr xxxx(1:Nr,1) = 0; end EE = [1 0;0 1]; AA1 = inv(EE-AA*h/2); for k=1:NNN xxxx(1:Nr,k+1) = AA1*( xxxx(1:Nr,k) + BB*h + AA*xxxx(1:Nr,k)*h/2); 1.4.4二阶显式Adams法 AA1 = (EE-AA*h/2); for i=1:1:Nr xxxxx(1:Nr,1) = 0; end EE = [1 0;0 1]; for k=1:2 xxxxx(1:Nr,k+1) = AA1*(xxxxx(1:Nr,k) + BB*h + AA*xxxxx(1:Nr,k)*h/2); end for k=3:NNN Fk = 23*(AA*xxxxx(1:Nr,k)+ BB); Fk1 = -16*(AA*xxxxx(1:Nr,k-1)+ BB); Fk2 = 5*(AA*xxxxx(1:Nr,k-2)+ BB); xxxxx(1:Nr,k+1) = xxxxx(1:Nr,k)+(Fk+Fk1+Fk2)*h/12; end 1.5实验结果与分析 H = 5*10-5s 1.5.1前向欧拉法仿真图形 1.5.2后向欧拉法仿真图形 1.5.3 梯形法 1.5.4二阶显式Adams法 H=1.0*10-6s 1.5.5仿真结果分析 1、通过单个模型的分别仿真，可以得出结论：二