欧拉法改进欧拉法斐波那契法原理及流程图.docxVIP

1欧拉法求微分方程 方法说明 欧拉(Euler)法是解常微分方程初值问题 (4.1) 最简单的数值方法，其详细做法是，将区间[a，b]进行N平分： ，步长.并将式(4.1)写成等价的积分形 式 （4.2） 再对式(4.2)右端积分用矩形公式计算，则有 ,(4.3) 在式(4.3) 右端取 ，舍去余项。则得 , 作为 的近似值。 在式(4.3) 右端取 ，舍去余项，则得 精选文档 作为的近似值． 一般地，在式(4.3)右端取舍去余项，则得 (4.4) 作为的近似值．式(4.4)为欧拉法计算公式． 我们知道微分方程的解是平面上的一族积分曲线，这族曲线中 过点的积分曲线就是初值问题式(4.1)的解． 欧拉法的几何意义是，过点引斜率为的积分曲线的切线，此 切线与直线的交点为，再过点引以为斜率的切线 与直线的交点为，依此类推，从出发，作以 为斜率的切线，此切线与直线交点为．于是便获得过 点的一条折线，见图4.1．过的积分曲线则用此折线来代 替．因此，这种方法亦称折线法． 图4.1 —2 精选文档 例： 用欧拉法求微分方程 = -2x = = 区间为 [ ] y y 0,1 y 欧拉法流程图如下: 开始 读入x0,y0,b,h 计算N=fix((b-x0)/h) n=1 x0+h=x1 y0+h*f(x0,y0)=y1 输出x1,y1 n=N? 结束 欧拉法程序如下： clear; clc; x1=0; x2=1; h=0.1; x0=0; y0=1; N=(x2-x1)/h;%要计算的次数 x(1)=x0; y(1)=y0; forn=1:N x(n+1)=x(n)+h;  n=1+n x1=x0 y1=y0 —3 精选文档 y(n+1)=y(n)+h*(y(n)-2*x(n)/y(n)); end X=x Y=y 改良欧拉法求微分方程 方法说明 由于欧拉法采用矩形公式计算积分产生较大截断误差．改良欧拉法(又 称改良折线法)是采取梯形公式来计算式(4.3)右端积分，则有 （5.1） 在式(5.1)右端取，舍去余项，则得 将作为的近似值． 在式(5.1)右端再取，舍去余项，则得 将作为的近似值． 一般地，在式(5.1)右端取，舍去余项．则得 (5.2) —4 精选文档 将作为的近似值． 式(5.2)为改良欧拉法计算公式． 流程图如下： —5 精选文档 例：用改良欧拉法求微分方程y = y - 2x ,y(0) = 1,h = 0.1, [ ] 区间为 y 改良欧拉法程序如下： clear; clc; x1=0; x2=1; h=0.1; x0=0; y0=1; p(1)=0; N=(x2-x1)/h; x(1)=x0; y(1)=y0; forn=1:N x(n+1)=x(n)+h; y(n+1)=y(n)+h*(y(n)-2*x(n)/y(n)); p(n+1)=y(n)+h*(y(n+1)-2*x(n)/y(n+1)); y(n+1)=(y(n+1)+p(n+1))/2; end X=x Y=y 斐波那契法求极值 方法说明 斐波那契法原理近似于黄金切割法，只是搜寻区间的缩短率不再采用黄 金切割数0.618。如图7.1所示，只需在[a,b]内取两点x1,x2,并计算出 f(x1),f(x2),经过比较，可将区间[a,b]缩短为[a,x2]或[x1,b]。因为新的 区间内包含一个已经计算过函数值的点，所以再从其中取一个试点，又可将这个新区间再缩短一次，不断地重复这个过程，直至最终的区间长度缩短到知足预先给定的精准度为止。 —6 精选文档 图7.1 现在的问题是，怎样选用试点，在保证同样精准度的情况下使得计算f(x)函数值的次数最少？在计算函数值的次数一定的情况下，最初区间与最终区 间的长度之比可作为取点方式优劣的一个标准。计算n次函数值，怎样取点 使最终区间最小？或许最终区间长度为1，计算n次函数值，初始区间最多 为多长？为此，引入Fibonacci数列： 表7.1 —7 精选文档 系 所以当试点个数n确定之后，最初的两个试点分别选为: —8 精选文档 显然x1,x2对于区间[a,b]对称，即有x1-a=b-x2,如图7.2 所示 图7.2 —9