计算机工与程应用复习题-2011-5.5 .ppt

非线性方程组 若利用牛顿迭代法求解，可转化成如下线性方程组求 对于截断的泰勒级数,从考察截断泰勒级数所引起的误差开始，当级数在包含(x-a)n项以后被截断时,表达f(x)的上式对f(x)的误差不大于 ?(x)=?(a)+?(a)(x–a)+?(a)/2!(x–a)2+ ?(a)/3!(x–a) 3+……+?(a)/n!(x–a)n 解:因为ex=e(0)+xe(0)+x2/2!e(0)+x3/3!e(0)+x4/4!e(0)+…… 所以ex=1+x+x2/2!+x3/3!+x4/4!+…… 若x=0.5,则 e0.5=1+0.5+0.52/2!+0.53/3!+0.54/4!+…… 若截至于e(0.5)3则: e0.5=1+0.5+0.52/2! =1.625 用泰勒法将ex在x=0附近展开,求e0.5及e(0.5)3 |d3(ex)/dx3|max(0.53/3!)=|ex|max(0.0208333) 这里的max表示在 区间上的最大值。 |ex|max=e0.5=1.6487213 所以，误差不大于 e(x)3=(1.6487213)*(0.0208333)=0.0343831 以二阶向后差分为例，说明以下递推公式成立: Δn ?j = Δ (Δn -1 ?j ) ▽n ?j = ▽(▽x –1 ? j ) ▽ 2? j= ▽ (▽? j )= ▽(? j - ? j –1) =▽ ? j - ▽ ? j –1 = ? j -? j –1- ? j –1+ ? j –2= ? j -2? j –1+ ? j –2 求解微分方程是一个复杂的过程，须事先附加一些条件，这就是对自变量的某些值规定出函数和（或）其导数的值 。 若在自变量为零的点上给出这些条件，这个问题就叫做“ ”。相应的条件叫做“ ”。 初值问题 初始条件 若在自变量不为零的值上做这种限定，这个问题便叫做“ ”，相应的限定条件叫做“ ”。 边值问题 边界条件 欧拉法求解初值问题的思路是： 根据泰勒级数，按初值条件展开： y(x0+h)=y(x0)+hy’(x0)+h2y(x0)/2+…… 若h相当小，h2及h的更高次幂就更小，略去h的高次项： y(x0+h)=y(x0)+hy’(x0) 据初始条件可从微分方程中得到y ’(x0)的值，于是在距离初始点x0只有一步h处得到了y(x)的近似值y(x0+h)。这个过程可以按下述关系持续下去： yn+1=yn+h?(xn,yn) n=1,2…… 改进欧拉法的递推公式为 yn+1=yn+1/2*h[?(xn?yn)+?(xn+1, y*n+1)] 在求解初值问题的常微分方程时，下列哪种方法最准确（） A、改进欧拉法 B、欧拉法 C、三阶的龙格—库塔法 D、四阶的龙格—库塔法 利用有限差分法求解已知边值问题的二阶微分方程的思路： 为简化讨论，以二阶方程： d2y/dx2=?(x,y,y’) 和边值条件 y(a)=A y(b)=B 使用有限差分法可以把求解边值问题转化为求解代数方程组问题，为了求解两点间边值问题： y’’=?(x?y?y) y (a)=A y(b)=B 可把边界a到b 的区间划分为“n”等分 Xi=x0+ih i=1,2,……,n 这里x0=a,xn=b, 而且 h=（b-a）/n y(xi)=1/2h(yi+1-yi-1) y’’(xi)=1/h2(yi+1-2yi+yi-1) 就能把微分方程写成“差分方程” 。如果对i=1,2,……，n写出了差分方程和两个边界条件。这个问题就化为n-1个未知数yi的n-1个方程的代数方程组。 在实际工程问题中，我们进行了大量的试验，从中得到了一个关于某个现象出现规律的描述方程，或者说