实验4-1原理讲解.ppt

实验4 常微分方程的数值解 §1 常微分方程简介 如 一般地，凡表示未知函数，未知函数的导数与自变量之间的关系的方程叫做微分方程．未知函数是一元函数的，叫常微分方程;未知函数是多元函数的，叫做偏微分方程． 本章主要针对一阶常微分方程进行研究 一阶常微分方程 有唯一解的充分条件是 对 满足Lipschitz条件，即，存在 使 §2 数值解的欧拉方法 数值解的提法∶在求解函数表达式 处的函数值 较困难的情况下 转而求解 为了书写简便，设 的步长相等，为 在区间 上用 代替 针对 （一） 中 取 则 即 记 则 以此类推， 称为向前欧拉公式 （二） 中 取 则 即 则 称为向后欧拉公式 求解需要迭代 · · 向前欧拉公式 向后欧拉公式 以此类推， p0 P1 P2 （四） 将梯形公式的迭代过程简化为两步，得到 称为梯形公式的预估-校正格式 常写为 （三） 将向前欧拉公式和向后欧拉公式平均，得到 称为梯形公式 求解也需要迭代 预估： 校正： 例 利用向前欧拉公式和梯形公式的预估-校正法求方程 解：向前欧拉公式的matlab实现 h=.1; x=0:h:1; n=length(x); y(1)=1; for i=2:n y(i)=y(i-1)+h*(y(i-1)-2*x(i-1)/y(i-1)); end [x;y] xiangqian1.m x 0 0.1 0.2 0.3 0.4 0.5 y(x) 1.0000 1.1000 1.1918 1.2774 1.3582 1.4351 x 0.6 0.7 0.8 0.9 1.0 y(x) 1.5090 1.5803 1.6498 1.7178 1.7848 向前欧拉公式结果 x 0 0.1 0.2 0.3 0.4 0.5 y(x) 1.0000 1.0954 1.1832 1.2649 1.3416 1.4142 x 0.6 0.7 0.8 0.9 1.0 y(x) 1.4832 1.5492 1.6125 1.6733 1.7321 精确解 的对应值 h=.1; x=0:h:1; n=length(x); y=zeros(1,n); y(1)=1; for i=2:n k1=y(i-1)-2*x(i-1)/y(i-1); k2=(y(i-1)+h*k1)-2*x(i)/(y(i-1)+h*k1); y(i)=y(i-1)+h*(k1+k2)/2; end [x;y] 梯形的预估-校正公式的matlab实现 yugu1.m x 0 0.1 0.2 0.3 0.4 0.5 y(x) 1.0000 1.0959 1.1841 1.2662 1.3434 1.4164 x 0.6 0.7 0.8 0.9 1.0 y(x) 1.4860 1.5525 1.6165 1.6782 1.7379 改进欧拉公式结果 x 0 0.1 0.2 0.3 0.4 0.5 y(x) 1.0000 1.0954 1.1832 1.2649 1.3416 1.4142 x 0.6 0.7 0.8 0.9 1.0 y(x) 1.4832 1.5492 1.6125 1.6733 1.7321 精确解 的对应值 §2 数值解的龙格—库塔方法 龙格—库塔方法的一般形式 其中 经典