数值计算方法课件_xutao_update.ppt

四、龙贝格算法 得 四、龙贝格算法 2）龙贝格算法 四、龙贝格算法 或 第八章 常微分方程 一、引言 1、常微分方程最简单形式： 初始条件： 积分： 解得： 一、引言 2、常微分方程的一般形式 3、高阶微分方程 例如二阶微分方程： 令 则 二、欧拉法 二、欧拉法 二、欧拉法 三、改进欧拉法 1、观察积分项：若采用梯形求积公式 可得到改进欧拉法 三、改进欧拉法 三、改进欧拉法 三、改进欧拉法 2、截断误差分析 1）欧拉法截断误差 三、改进欧拉法 三、改进欧拉法 2）改进欧拉法截断误差 四、龙格-库塔法 一阶常微分方程 初始值： 四、龙格-库塔法 二、曲线拟合 整理得二次曲线的正规方程组 二、曲线拟合 3、m次多项式曲线拟合 设拟合曲线的m次多项式为： 拟合曲线的正规方程组： 第七章 数值积分 一、引言 给出一组数据 二、插值求积分 1、两点式（梯形公式） 二、插值求积分 2、三点式（辛卜生公式） 二、插值求积分 3、五点式（科特斯公式） 二、插值求积分 4、复化求积公式 1）复化梯形 在每个小区间上应用梯形公式 式中h为等距步长 二、插值求积分 2）复化辛卜生 在每两个小区间上应用辛卜生公式 式中h为等距节点 二、插值求积分 3）复化科特斯 在每四个小区间上应用科特斯公式 式中h为等距节点 三、求积公式的误差 的舍入对公式的影响； 求积公式的截断误差 。 1、舍入误差 三、求积公式的误差 三、求积公式的误差 三、求积公式的误差 2、截断误差 1）梯形公式 三、求积公式的误差 或： 三、求积公式的误差 2）辛卜生公式 三、求积公式的误差 三、求积公式的误差 三、求积公式的误差 3）科特斯公式 四、龙贝格算法 复化求积法 步长减小，精度提高 步长小，工作量增大 变步长求积法 为统一精度与计算量的矛盾，采用变步长，每次步长减半，直至迭代精度达到要求为止。 变步长梯形法 变步长辛卜生法 变步长科特斯法 尨贝格算法 四、龙贝格算法 1、变步长梯形法 四、龙贝格算法 四、龙贝格算法 2、变步长的辛卜生求积法 变步长梯形法算法简单，单精度低，收敛慢 四、龙贝格算法 四、龙贝格算法 3、龙贝格算法 1）变科特斯算法 五、埃特金插值 1、二次插值 五、埃特金插值 五、埃特金插值 2、三次插值 五、埃特金插值 五、埃特金插值 六、分段插值 六、分段插值 1、分段线性插值 2、分段抛物线插值 七、数值微分 1、插值型求导公式余项 2、一阶导数的二点公式 七、数值微分 3、一阶导数的三点公式 八、三次样条插值 高阶插值震荡不稳定； 分段插值不光滑； 飞机、船体放样要求在节点处要光滑，且不能产生震荡，这样求放样插值在节点处具有二阶连续导数； 样条：弹性木条，强迫它通过插值点（施加力），根据力学原理，样条所形成的曲线满足二阶连续导数； 三次插值多项式是具有二阶连续导数的函数。 八、三次样条插值 八、三次样条插值 八、三次样条插值 八、三次样条插值 八、三次样条插值 八、三次样条插值 八、三次样条插值 八、三次样条插值 八、三次样条插值 第六章 曲线拟合 一、最小二乘法 右图中 与测量值的误差，确定 使其与测量值误差的绝对值最小 或者 一、最小二乘法 对上式求偏导数 令 整理得 一、最小二乘法 解得 最后得到 二、曲线拟