上海海事大学研究生数值分析复习大题.pdfVIP

1.求系数a,b，使求解常微分方程初值问题 y  f(x,y) 的数值解公式 yxa  s y  y  h(ay  by  )的局部误差为y 3 n1 n n n1 (xn1)  yn1  O(h ) 解：设步长 ，且y  y(x ) ，y  y(x ) 。 h n n n1 n1 因    2 ，故  2  3 y n1 y n hy n O(h ) y n1 y n (a b)hy n bh y n O(h )  1 2  3 3 1 又y(x )  y(x ) hy (x )  h y (x ) O(h ) ，比较得a  ，b   n1 n n n 2 2 2 2.对初值问题  写出以梯形公式所得近似解 的表达式。且当步长 y  y  0 y(0) 1 y n 时，讨论 的收敛性。讨论绝对稳定性对步长的限制 h 0 yn h 解：梯形公式为：yn1  y  [ f (x , y )  f (x , y )] n 2 n n n1 n1 n1 h 2 h  2 h  故：yn1  yn  [yn  yn1 ] ，整理得： yn1   yn    y0 2 2  h  2  h  n  2 h  所以：yn    。 2  h  n 2 h  x 又：设任给x  0 ，nh  x ，于是有lim yn 