数学实验之2.微分方程20021.pptVIP

蓝色曲线 ——y（1）； （原方程解） 红色曲线 ——y（2）； (导函数曲线) 1.3855 1.3487 1.3679 1 1.3241 1.2874 1.3066 0.9 1.2665 1.2305 1.2493 0.8 1.2132 1.1783 1.1966 0.7 1.1645 1.1314 1.1488 0.6 1.1209 1.0905 1.1065 0.5 1.0830 1.0561 1.0703 0.4 1.0513 1.029 1.0408 0.3 1.0264 1.01 1.0187 0.2 1.0091 1 1.0048 0.1 1 1 1 0 向后欧拉 向前欧拉 精确解 x （1）步长h=0.1的数值解比较表 计算结果 （2）步长h=0.01的数值解比较表 1.3697 1.3660 1.3679 1 1.3084 1.3047 1.3066 0.9 1.2511 1.2475 1.2493 0.8 1.1983 1.1948 1.1966 0.7 1.1504 1.1472 1.1488 0.6 1.1080 1.1050 1.1065 0.5 1.0717 1.0690 1.0703 0.4 1.0419 1.0397 1.0408 0.3 1.0195 1.0179 1.0187 0.2 1.0053 1.0044 1.0048 0.1 1 1 1 0 向后欧拉 向前欧拉 精确解 x 结论：显然迭代步长h 的选取对精度有影响。 图形显示 有什么方法可以使精度提高？ 向后欧拉法 返 回 对方程y’=f(x,y), 两边由xi到xi+1积分，并利用梯形公式，有： 2、使用数值积分 即梯形法： 梯形公式 改进欧拉公式 yn + hf(xn,yn) 返 回 以例1为例，用改进欧拉公式编程计算，再与精确解的比较。 yn+1 = yn + (h/2)*[(-yn + xn+ 1) + (- yn +1+xn +1+1)] = yn + (h/2)* [(-yn + xn+ 1) -(yn+ h*(-yn + xn+ 1) )+xn +1+1] = yn + (h/2)* [ (1-h)*xn + xn+1 + 2-h + (h-2)*yn] died1.m 使用改进欧拉公式的例 1.3685 1.3855 1.3487 1.3679 1 1.3072 1.3241 1.2874 1.3066 0.9 1.2500 1.2665 1.2305 1.2493 0.8 1.1972 1.2132 1.1783 1.1966 0.7 1.1494 1.1645 1.1314 1.1488 0.6 1.1071 1.1209 1.0905 1.1065 0.5 1.0708 1.0830 1.0561 1.0703 0.4 1.0412 1.0513 1.029 1.0408 0.3 1.019 1.0264 1.01 1.0187 0.2 1.005 1.0091 1 1.0048 0.1 1 1 1 1 0 改进欧拉 向后欧拉 向前欧拉 精确解 x 步长 h= 0.1 的数值解比较表 使用改进欧拉公式的例 3、使用泰勒公式 以此方法为基础，有龙格-库塔法、线性多步法等方法。 4、数值公式的精度 当一个数值公式的截断误差可表示为O（hk+1）时（k为正整数，h为步长），称它是一个k阶公式。 k越大，则数值公式的精度越高。 欧拉法是一阶公式，改进的欧拉法是二阶公式。 龙格-库塔法有二阶公式和四阶公式。 线性多步法有四阶阿达姆斯外插公式和内插公式。 返 回 [t，x]=solver（’f’,ts,x0） ode23 ode45 由待解方程写成的m-函数文件 ts=[t0，tf]，t0、tf为自变量的初值和终值 函数的初值 ode23：组合的2/3阶龙格-库塔算法 ode45：运用组合的4/5阶龙格-库塔算法 自变量值 函数值 Matlab软件计算数值解 1）首先建立M-文件 （weif.m） function f = weif(x,y) f=-y+x+1; 2）求解：[x，y]=ode23(‘weif’, [0, 1], 1) 3) 作图形: plot(x, y, ‘r’); 4) 与精确解进行比较 hold on ezplot(‘x+exp(-x)’,[0, 1]) 例1 y’= - y+x+1，y(0) = 1 标准形式： y’= f(x ,