数值分析习复题课.pptVIP

例12 对如下数据作形如 y = a eb x 的拟合曲线 解: 由于函数集合Φ={a eb x | a,b ∈R } 不成为线性空间，因此直接作拟合曲线是困难的。 在函数 y = a eb x 两端分别取对数得到 这时，需要将原函数表进行转换如下 令 z= ln y , A = ln a , B=b, 则 z=A+Bx ln y = ln a+bx xi 1 2 3 4 5 6 7 8 yi 15.3 20.5 27.4 36.6 49.1 65.6 87.8 117.6 第三十页，共八十页，2022年，8月28日 对 z=A+Bx 作线性拟合曲线，取 这时 xi 1 2 3 4 5 6 7 8 yi 15.3 20.5 27.4 36.6 49.1 65.6 87.8 117.6 xi 1 2 3 4 5 6 7 8 zi 2.72 3.02 3.31 3.60 3.89 4.18 4.48 4.77 第三十一页，共八十页，2022年，8月28日 得正则方程组 解得 于是有 拟合曲线为： 第三十二页，共八十页，2022年，8月28日 第三章数值积分 插值型积分公式 Newton-Cotes 型求积公式 复化求积公式 Romberg算法 Gauss积分 数值微分 第三十三页，共八十页，2022年，8月28日 Simpson公式n=2 梯形公式 n=1 第三十四页，共八十页，2022年，8月28日 (§1 ,2)需要掌握： 1. 等距节点(Newton-Cotes)的积分公式是如何构造的？ 第三十五页，共八十页，2022年，8月28日 2. N点等距节点的积分公式及其误差式怎么表示？ 3. 如何由上式给出梯形公式、Simpson公式及其误差？ 第三十六页，共八十页，2022年，8月28日 复化梯形公式及其误差 例 3-0 用四点复化梯形公式计算 0 1 第三十七页，共八十页，2022年，8月28日 例 3-1 用复化梯形公式计算积分 ，应将区间[0,1]多少等分，才可以使其截断误差不超过 解：复化梯形公式的误差为 而 从而 令 于是，只要将区间至少68等分,就可以达到需要的精度要求。 第三十八页，共八十页，2022年，8月28日 复化Simpson公式及其误差为： 第三十九页，共八十页，2022年，8月28日 要求： 了解各种积分公式的原理，构造方法, 会利用公式计算积分， （复化）梯形公式，(复化)Simpson公式及余项表达式，求解代数精度 会利用代数精度构造积分公式，并用构造的积分公式计算相应积分值 Romberg算法的实现原理，计算，外推加速技术； 数值微分公式的构造方法 第四十页，共八十页，2022年，8月28日 一、确定数值积分公式或数值微分公式，并推出余项 根据代数精度的概念 对Guass型求积公式，可借助Guass点与求积系数的关系确定参数 推导余项时，可设 对于数值微分公式，可构造适当的插值多项式或应用Taylor展开式推导 第四十一页，共八十页，2022年，8月28日 第四十二页，共八十页，2022年，8月28日 第四十三页，共八十页，2022年，8月28日 第四十四页，共八十页，2022年，8月28日 第四十五页，共八十页，2022年，8月28日 第四十六页，共八十页，2022年，8月28日 第四十七页，共八十页，2022年，8月28日 p72 方法二、p87例2 利用变量替换,将[a,b]转化为[-1,1]区间 第四十八页，共八十页，2022年，8月28日 第四十九页，共八十页，2022年，8月28日 余项表达式 第五十页，共八十页，2022年，8月28日 二、计算定积分和函数的导数的近似值 对于给定的被积函数与求导函数，应用指定的数值积分公式或数值微分公式计算， t9,t12,t13,t18,t19，t25，t26等 明确积分公式与微分公式 三、确定复化求积公式和数值微分公式的步长或节点数，使计算结果满足所给精度要求 根据复化求积公式和数值微分公式的余项或截断误差表达式，对满足精度要求解一个相应的不等式，即可确定所需的步长或节点数 第五十一页，共八十页，2022年，8月28日 n=213 第五十二页，共八十页，2022年，8月28日 n=4 N169.2 n=170 第五十三页，共八十页，2022年，8月28日 n=4 第五十四页，共八十页，2022年，8月28日 常微分方程数值解 复习 常微分方程初值问题 Euler法（显式Euler公式，隐式Euler公式，梯形公式，改进Euler公式，变形Euler公式）基本公式 Runge-Kutta