微分方程数值方法习题二.docxVIP

微分方程数值方法 常微分方程初值问题习题一 对初值问题 p/f= ax + b. [(0) = 0, 分别写LB Euler法和改进的Euler法的近似解％的表达式，并求 它们与真解 u(x) = ^-ax2 +bx 的差 u(xm)-um . J 取步长/? = 0.1 ,分别用Euler法和改进的Euler法求下列初值问 题的解，并与真解和比较. (1)0 x 1, (1) w(0) = 1, 真解 w(x) = Jl + 2x ; (2)心厂产I*, 讥 1) = 2, 真解 u(x) = x(8-31nx)3 ; 1 x 1 x 1.5, ⑶“心丟一乔 讥 1) = 1, 3 ] 真解 u(x) = (4x2 -3x2)\ 用Euler法计算[efldt在x = 0.1, 0.2的近似值. 取步长A = 0.2,用四阶Runge-Kutta法解 0= U + X. 0 X 1, M(0) = 1, 并与真解u(x) = 2ex - x-1相比较? 设 f(x.u)关于 u 满足 Lipschitz 条件，证明 N 级 Runge-Kutta 法中的增量函数(pZ,h)关于”也满足Lipschitz条件. 对初值问题 (u= u + x-l, t ?(0) = 1, 写岀四阶Taylor级数法和四阶Runge-Kutta法的计算公式，它们 是否和同. 7?证明改进的Euler法的绝对稳定区间是(-2, 0)? 8?证明：当h(=帥)满足 _2 _3 _4 + A+ — + — + — 1 6 24 时，四阶Runge-Kutta法绝对稳定. 用Tayor展开确定下面多步法中的系数，使其阶尽可能高，并求 局部截断误差的主项. (1)= aiUm +。2%1 + 力 00 九+1 ； (2 ) um+] = um + h(0\ fm + %九一】)； (3)%i =勺％_1 + W + 0| ftn + 02九+1)? 对初值问题 wH= f(x,u),u(xQ) = uQ.u\xQ) = u 确定求解公式 %+= 2知一％1 +力2(0()九_1 +01 九 +02九—1) 中的系数0°，02与局部截断误差主项. 求公式 h 2 2 %i =% +#/(?”，％ ) + 3/(? +~h,uin +-hf(xm,um))] 的阶和局部截断误差. 取步长A = 0.1,用四阶Adams方法的的预测■校正公式(72)解 初值问题 [u= u + x, 0 x 1, [(0) = 1, 并与真解u(x) = 2ex -x-l相比较. 讨论下列公式的相容性、稳定性、和绝对稳定性. %i = 一4陷 + 5%i + 加4九 + 2九); 1 3 心+】=-(9um - %_2)+ g 力(九+i + 2九一九-)? 求a ,使线性多步法 h , %+2 一 (1 + Q)%+i + a% = —[(3)九+1 一 (1 + Q)九] 是相容的和稳定的. 证明三阶Adams内插公式 h Um+l = Um + — (5九+1 + 8九- fin_i ) 的绝对稳定区间是(-6,0)? 16?证明中点公式(90)是A■稳定的. 求下列方程组的刚性比.若用四阶Adams内插公式求解吋，最 大步长应小于多少？ u= 一2000“ + 999.75V +1000.25, (1 ) v= u - V, h(0) = 0,v(0) = 2; (2)u= -0.1m -49.9v, v= -50v, 0=150—200° u(0) = 2,v(0) = l,^(0) = 2. (2) 把下列高阶方程化为一届方程组，并写出它们的Euler公式和四 阶 Runge-Kutta 公式. un= 2xu un= 2xul+4x,u(0) = 0,u(0) = l; w,!= —20+(—-— l)w, X w(l) = . 椭圆型方程的差分法习题三 用五点差分格式解下列椭圆型方程边值问题. Aw = 2, 0 x, y 1, (1 厂 w(0, y) = 0, u(l, y) = 1 + y, q 2 u(x90) = X ,w(xj) = X + X, 取 h = r =—(真解 u = x2 4- xy ); ? Aw = 0, 0 x, y 3, (2) w(0,y) = -y2,w(3,y) = 9-y2, m(x,0) = x2, w(x,3) = x2 -9, 取/? = r = 1 (真解“ =x2 - y2 ); 3 3 (3)-Aw = 9 sin(— 7ix) sin(—码J, ?|r2= 0, 2 ? (3) 取 h = r =-; 3 * 一 Aw = (x2 + y2)sin(xy), 0