实验四微分方程的数值解.docVIP

实验四 微分方程数值解及应用 　 　 一　引例　　 在18世纪末，英国人马尔萨斯在他出版的一本专著中，对人口数量的增长趋势进行了模拟，提出了人口的指数增长模式，导致最后会出现人口数量超过地球的承受容量，即人口爆炸问题。虽然马尔萨斯人口模拟忽略了一些人口增长的重要因素，但是他为以后人口模型的改进提供了基础，下面将介绍著名的马尔萨斯人口模型及改进。 马尔萨斯认为：单位时间内，人口的出生数量和死亡数量与人口总数成正比例，即人口出生率和死亡率都是常数，因此人口的净增长率为常数。 设时刻的人口数量为，人口出生率为，死亡率为，则有： ， 从而 ， 其中： 称为净增长率（常数）。因此马尔萨斯人口模型如下： ，。 该微积分方程初值问题的解析解： 。 用此模型预测某国1850年人口数量的误差为1% ，1900年人口数量的误差为31%，2000年人口数量的误差达567%，2050年将达到80多个亿。 该模型对于长期预测不合理。 针对马尔萨斯人口模型的不足，1837年荷兰生物学家Verhulst提出了如下改进： 由于资源的限制，人口存在最大值（极限） 。因此，人口增长率不应该是常数，假设增长率是随着人口数量接近 而线性递增： ， 从而得到改进后的人口模型为 ，， 称为逻辑斯谛增长模型。 刻画世界千变万化的规律，微分方程是最有力的工具。 二、微分方程的数值解法 讨论一阶微分方程的初值问题 ， （1） 或一阶微分方程组的初值问题 其中是向量函数，为自变量。 微分方程的数值解法的基本原理： 引入自变量的取值点列，定义，称为步长，常用等间距的步长 （与无关，记为），假设精确解为。为了寻求的近似值，根据一定的原理，结合当前得到的近似解，近似地表示该点或前一点的导数值，由此推出计算的迭代公式。 欧拉方法 欧拉方法是一种简单的求解初值问题的数值逼近方法，其基本思路为： 对方程（1），在小区间上，用差商 代替导数， 1）用左端点替换函数中的，得到方程（1）的近似表达式 设，则的近似值为 （2） 其中为初始点，（2）称为显示欧拉公式，也称向前欧拉公式。 向前欧拉法简单易于计算，但精度却不高，收敛速度慢。 2）用右端点替换函数中的， （3） 其中为初始点，（3）称为隐示欧拉公式，也称向后欧拉公式。 这是一个非线性方程，无法直接计算。在计算精度、收敛速度方面与向前欧拉算法相同。 3）将（2）式与（3）式加以平均，得到 （4） 称为梯形公式。与前两个方法相比，该方法的计算精度高、收敛速度快。但迭代计算与(3) 式一样繁！ 4）改进的欧拉公式：先由公式（2）算出的预测值，然后带入梯形公式（4）的右端，作为校正，即 （5） 称为改进的欧拉公式。还可以写为 人们常用的是向前欧拉公式和改进的欧拉公式。且欧拉法可以推广到求解微分方程组的情形。 练习1：用上述方法求解微分方程初值问题 要求：取步长 ，分别用三种数值解法求解，并结合其精确解，对求解误差进行分析比较。 首先用解析法得到其精确解。 其次用数值解法 向前欧拉算法 迭代公式为 向后欧拉算法 迭代公式为 龙格-库塔方法 龙格-库塔方法，简称方法，是利用泰勒展式将在处展开，并取其前面若干项来近似而得到公式 如果，则 的近似值为 若 则称以上迭代公式为阶公式，的大小决定误差的高低，高阶高精度。要得到一个阶公式，关键在于如何选取 使之满足界的要求。 2-阶公式 中点公式 改进的欧拉公式 3-阶公式 4-阶公式 在MATLAB软件中含有数值求解的系统函数，其实现原理就是龙格-库塔方法。 三、实验作业 1．求下列微分方程的解析解。 1)一阶线性方程 ； 2)高阶线性齐次方程 3)高阶线性非齐次方程 . 2. 用两种方法求解方程 3.(单摆运动的近似解析解)当单摆初始角度 较小时， 也较小，从 ，单摆运动微分方程可近似写为 求此方程的解析解，并与数值解进行比较。 4.应用题 鱼雷击舰问题 一敌舰在某海域内沿正北方向航行时，我方战舰