使用C语言解常微分方程CODE.docx

文档来源为:从网络收集整理 文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持. PAGE PAGE # 解常微分方程 姓名：Vincent 年级：2010,学号：1033****,组号：5 （小组），4 （大组） 数值方法： 我们的实验目标是解常微分方程，其中包括几类问题。一阶常微分初值问题，高阶 常微分初值冋题，常微分方程组初值问题，二阶常微分方程边值问题,二阶线性常 微分方程边值问题。 对待上面的几类问题，我们分别使用不同的方法。 ?初值问题 使用龙格?库塔来处理 ?边值问题 用打靶法来处理 ?线性边值问题 有限差分法 初值问题 我们分别使用 ?二阶龙格-库塔方法 4阶龙格-库塔方法 来处理一阶常微分方程。 理论如下： 对于这样一个方程 当h很小时.我们用泰勒展开， 当我们选择正确的参数a国皿叩皿之后，就可以近似认为这就是泰勒展开。 对于二阶，我们有： 其中 经过前人的计算机经验，我们有. 选择则P=LQ=L于是又如下形式，这也使休恩方法的表达式。 所以我们称其为 龙格（库塔）休恩方法 对于4阶龙格方法，我们有类似的想法， 我们使用前人经验的出的系数，有如下公式 对于高阶微分方程及微分方程组 我们用 4阶龙格?库塔方法来解 对于一个如下的微分方程组 我们可以认为是-个一阶向量微分方程，所以可以用龙格?库塔方法的向虽形式解八 对于一个高阶的微分方程，形式如下