梯形数值分析法.docVIP

第九章 解常微分方程初值问题 梯形法 在⑵右端积分用梯形求积公式近似，并用儿代替y(耳)，儿+|代替y(x〃+J 儿+产 儿 + ￡(/(兀“，儿)+/(丟+1，%+1)) (4) 梯形法的迭代公式为： =儿+力*/(兀,儿) /； r 元丁 =儿 + -[/Un，儿)+ f(xn+l，必] ⑸ (()」,2,…). 梯形方法也是隐式方法，耍迭代求解?迭代收敛的条件是处〈1? 2 梯形法的谋差估计梯形方法的局部截断误差是。(F)?整体截断误差是0 (F). 上机实验⑶ 上机题目：贬值求一阶常微分方程的初值问题 实验目的：掌握齐种Euler方法和梯形法。 进过计算结果來分析四种方法的优缺点，掌握规律。 分析结果。 实验要求：用不同的方法来解同一个例子。 上机前充分准备，复习有关内容，写出计算步骤，查对程 序。②梯形法在Matlab环境中运算，并分析岀最好的方法， 再给出它的流程图。③实验结朿后写出完整的实验报告。 算法说明：①经过所给出的方程组和初值初步改变方程组。 由以上四种方法的计算公式來逐步计算y的每个值。 最后为了方便比较列为表最适合。 上机例题：例1.后退Euler方法解初值问题,/7==().1 y = y?2x/y, 0(x(1 y(o)= 1 Mat lab 程序: %定义步长%计算x(j) %定义步长 %计算x(j) for j=l :n+l x(j)=a+h*(j-l); end y(l )=y0; for i=l:n %计算 y(i+l)y(i+1 %计算 y(i+l) y(i+l)=y(i)+h*(subs(subs(f；x(i)),y(i))+subs(subs(t；x(I+l)),y(i+l)))/2; %算法公式来计算 t二y； 运行结果： t=Tixing(,y-2*x/y,0,1,0,1,10) t =1.0000 1.0957 1.1836 1.2654 1.3423 1.4151 1.4843 1.5504 1.6139 1.6751 1.7341 TG 今[/(a?)] S?=S+/(x) xV=x+〃 xb? I输出耳 I输出耳 复习思考题 为什么耍研究微分方程数值解？本章主要研究的是怎样一类初值问题？ 何谓方法的阶与局部截断误差？ 何谓单步法、多步法？何谓显示公式、隐式公式？四阶龙格一库塔法属于哪一 类？ 线性多步法的优缺点是什么？ 预测一校正公式有何特点与优点？ 试述用带误差修正的四阶阿当姆斯预测一校正公式（或哈明公式）解初值问 题的计算步骤。 一阶微分方程组与高阶方程初值问题如何求解？导出下述三阶方程初值问题 的龙格一库塔公式。 什么叫高斯公式和高斯点？无伙= 1,2,…/）是高斯点的充耍条件是什么？怎 样求高斯点和高斯公式？高斯公式有何优点和缺点？ 什么叫待定系数法？怎样利用它确定数值积分公式？有什么优点和缺点？ 插值型求导公式是怎样形成的？其误差怎样估计？为什么使用插值型求导 公式要注意误差分析？样条插值求导有什么优点？