连续系统时域数值解法序设计.docVIP

毕业设计 题 目 连续系统时域数值解法程序设计 所在院(系) 物 电 学 院 专业班级 电子信息科学与技术102班 连续系统时域数值解法程序设计 [摘要] 采用Mathematica软件的NDSlove求连续系统的数值解，并编写求连续系统数值解的实用程序，从实例展示程序的用法和高效率。 [关键词] Mathematica程序; 数值解; 时域 The program design of the continuous systems time domain numerical solution Abstract: Use the NDSlove of Mathematica software to solve continuous system’s numerical solution.And compile the utility program of solving continuous system’s numerical solution.From the living example to diaplay program’s usage and high efficency Keywords: Mathematica program;numerical solution; time domain 目录 引言 1 1 描述LTI连续系统的方法 1 1.1 微分方程与初值的方法描述 1 1.2 系统函数 2 2 连续系统数值解的Mathematica程序设计思路 3 3 复杂电路系统数值解的求解程序设计思路 4 3.1 RLC三阶电路 4 3.2数学模型的建立 5 3.2.1 列出时域连续系统电路的微分方程 5 3.2.2 连续系统电路的时域映射到复频域（S域） 5 3.3 NDslove方法求解连续系统的优缺点 6 4 连续系统数值解的可实用Mathematica程序 6 4.1 程序设计思路 6 4.2 程序清单 6 5 动态电路数值解程序应用示例 7 5.1 时域动态电路 7 5.2 程序应用展示 7 6 结语 9 参考文献 9 引言 连续系统的各种解析解法虽然便于理论分析系统响应的变化趋势的系统特性，但实际系统总是多输入多输出的高阶系统，它们的解析形式的响应求解极为困难、即便较低阶系统的解析响应能够得到，其函数表示也比较复杂[1]。而连续系统的区间数值解法本质上用的是迭代解法 ，总是能够方便、快速地的得到，之后如果企图观察其响应随时间演化的趋势，可用数值解画出其波形来观察，甚至必要时做数据拟合寻找区间解的拟合函数也是可能的，而且数值解法还可以求一定区间上的非线性问题。 。 令，。根据时域微分定理，y（t）及其各阶导数的拉普拉斯变换为 （） （1.1-2） 如果f(t)是t=0时接入的，则在时f(t)及其各阶导数均为零，即。因而f(t)及其各阶导数的拉普拉斯变换为 （1.1-3） 取式（1.1-1）的拉普拉斯变换并将式（1.1-2）、（1.1-3）代入得 即 （1.1-4） 由上式可解得 （1.1-5） 式中，是方程（1.1-1）的特征多项式；，多项式和的系数仅与微分方程的系数、有关；，它也是s的多项式，其系数与和响应的各初始状态有关而与激励无关。 由式（1.1-5）可以看出，其第一项仅与初始状态有关而与输入无关，因而是零输入响应的象函数，记为；其第二项仅与激励有关而与初始状态无关，因而是零状态响应的象函数，记为[4]。于是式（1.1-5）可写为 （1.1-6） 式中，，。取上式逆变换，得系统的全响应 （1.1-7） 1.2系统函数 描述n阶LTI系统的微分方程一般可写为 （1.2-1） 设是时接入的，则其零状态响应的象函数为