基于VC++、MATLAB和EXCEL对RL暂态电路的数值分析.pdfVIP

第35卷(2007)第9期 计算机与数字工程 189 对RL暂态电路的数值分析‘ 刘琪1’ 巴力登’杨仁春2’ (新疆大学电气工程学院1’乌鲁木齐830046)(新疆大学化学化工学院2’乌鲁木齐830046) 单暂态电路的i～t进行数值分析。三程序计算结果并无差异，但Vc++在调用函数时需加人相应头文件，且其程序较 MAmB繁琐，ExcEL需做逐步单元格的自动填充实现，M删B则最为简单。 MATLAB 美键词vc++6．0 ExcEL龙格一库塔法数值分析RL暂态电路 中图分类号唧09 计算开关S闭合 1 引言 后电流i(t)在区 在工程计算和设计中，由于存在微分方程初值 间[O，O．01]上的 问题不能求得解析解，因而采用近似法求解，通常 电流值。 采用的近似方法为数值分析法，通过其可以解出一 则由电路理 些离散点上的近似值。运用经典四阶龙格一库塔 论可得i—t的常 法，在三种程序下实现了一阶RL简单暂态电路中 微分方程： i—t的数值分析。 J面2T”““一百‘f祟=譬sin“一譬t (2) 2 四阶龙格一库塔法简述 【i(O)=O 龙格一库塔方法是间接地使用泰勒级数法的 代入相关参数得数学表达式： 一种应用。其计算精度较高，数学描述如下： 』塞。622sill314z一207 (3) y(i+1)=y(i)+h}(kl+2‘k2+2’k3+k4)／6 【y(0)=0 (1) 3．1 vc++6．O算法 kl=f(xI，y，) 首先创建一个新文件，步骤如下：在New对话 k2=f(xI+h／2，yj+h+hl／2) console 框的工程选项卡中选择Win32 ApPlication k3=f(x．+h／2，y。“‘k2／2) 类型，输人工程名称及路径。选择AsimpleAppH— k4=f(x。+h，yI+h‘k3) cation后单击完成，系统自动生成一个工程项目。 其中：h表示计算过程中选取的步长；k。表示 毛点处的斜率；k：表示利用k。求得的(毛+h／2)点 输入如下程序代码： 处的斜率；k3表示利用k：求得的(x；+h／2)点处的 荆nclude“sld止．b” 斜率；k。表示利用k3求得的(x．+h)点处的斜率。 僦nclude“iostr朗m．h” 衔nclude“m