数值计算方法讲稿11.pptVIP

数值积分 §1. 插值型求积公式 (二) 抛物型求积公式 (三.1) Newton-Cotes求积公式 §2. 梯形,抛物线公式的误差估计 n=偶数时Newton-Cotes 求积公式的代数精确度 梯形公式的截断误差 定理证明 §3. 复化公式及其误差估计 误差公式: 区间越小, 误差更小——复化。 [a, b]2n 等分, 可得 T2n : 复化抛物型公式 ∵ 抛物线公式用到区间中点， ∴ 将区间看作等分偶数份 令 n=2m, m 是整数, 在每个 [?2k-2, ?2K] 上 用抛物线公式： 自动选步长计算 由误差要求可定出n ，但事先难估计：应边算边估计进而加密， 自动分半的 Simpson 公式 §4. Richardson 外推算法 §5. Romberg求积法 R(x)公式证明 复化公式误差 复化求积例 只要利用公式不断计算新分点之函数值， S2 ， S4 …， Sn ， S2n … 复化求积例（3） 如何仅通过构造 F1 的线性组合产生更高阶逼近函数 F2(h) 呢？ 通过适当线性组合就可以明显提高逼近阶! 同样可继续推广(Richardson 外推算法)： 两个步长分别为 h , 0.5h 的低级计算值的线性组合产生一个高级计算值，因此不断地分半计算是必须的： Romberg求积法的收敛性 Romberg 算法 §6. 高斯公式 牛顿—柯特斯型求积公式是封闭型的（区间[a,b]的两端点a, b均是求积节点）而且要求求积节点是等距的,受此限制，牛顿—柯特斯型求积公式的代数精确度只能是n(n为奇数)或n+1(n为偶数)。而如果对求积节点也适当的选取,即在求积公式中不仅Ak而且xk也加以选取,这就可以增加自由度，从而可提高求积公式的代数精确度。 高斯－勒让德公式 带权的高斯公式 * 第六章、数值微分与数值积分 数 值 微 分 x-h x x+h B C A T f(x) 差商型求导公式的余项 从截断误差的角度看，步长越小，计算结果越准确； 从舍入误差的角度来看，步长不宜太小。 2、插值型求导公式 两 点 公 式 三 点 公 式 同样，针对m也可扩展，如五点插值求积公式。 3、 样 条 求 导 (三.2)计算系数 Ai : 变换 x=a+th (三.3) Newton-Cotes系数Ci(n) 8 7 6 5 4 3 2 1 n Newton－Cotes公式的稳定性 衡量插值型求积公式的精度,可以用多项式的次数作为标准. 梯形求积公式的代数精确度 Newton-Cotes求积公式的代数精确度 抛物求积公式(Simpson)的截断误差(1) 证明思路: 1,将f(x)用插值多项式表示而且与抛物公式值相同 抛物求积公式误差证明(1) 抛物求积公式误差证明(2) 复化梯形公式的分半加密算法。 *