插值与数值积分.pptVIP

将f(x)代入计算得 用n个节点，Gn的代数精度可达2n-1, 但是需解 复杂的非线性方程组，实用价值不大。 常 用 的 高 斯 公 式 将(a,b)分小，把小区间变换为(-1，1), 再用G2 代数精度为3 节点加密时，原计算信息无法利用 思路：将积分区间分小，在小区间上用n不太 大的 。而在节点加密一倍时能够利用原节点的函数值，可以把区间的端点作为固定节点。 改进的高斯公式 Gauss-Lobatto求积公式 其中a, b为小区间的端点， 为2n-2个参数， 代数精度可达到2n-3 注意：实际计算中一般采用自适应方法确定步长 用MATLAB 作数值积分 矩形 公式 Sum(x) 输入数组x(即fk),输出x的和(数) cumsum(x) 输入数组x,输出x的依次累加和(数组) 梯形 公式 trapz(x) 输入数组x,输出按梯形公式x的积分(单位步长) trapz(x,y) 输入同长度数组 x,y,输出按梯形公式 y对x的积分(步长不一定相等) 用MATLAB 作数值积分 辛普森公式 quad(@fun,a,b,tol,trace) [I,fn]=quad(…) 用自适应辛普森公式计算 tol为绝对误差，缺省时为10-6 Gauss-Lobatto公式 quadl(@fun,a,b,tol,trace) [I,fn]=quadl(…) 用自适应Gauss-Lobatto公式计算 tol为绝对误差，缺省时为10-6 注意：fun.m中应以自变量为矩阵的形式输入(点运算) 矩形域上计算二重积分的命令： dblquad(@fun,xmin,xmax,ymin,ymax,tol) 广义积分、二重和三重积分 长方体上计算三重积分的命令： triplequad(@fun,xmin,xmax,ymin,ymax, zmin,zmax,tol) 注：fun是被积函数，本身可以有自己的参数 广义积分： 通过分析和控制误差，转换成普通积分 quadv(@fun,a,b,tol,trace) 向量值积分： 用MATLAB 作数值积分 例. 计算 1）矩形公式和梯形公式 将(0, ? /4)100等分 2）辛普森公式和Gauss-Lobatto公式 精确、方便 无法计算用数值给出的函数的积分 Jifen1a.m Jifen1b.m 精确值为 数值积分的应用 实例 人造卫星轨道长度 轨道长度 y x o 近地点s1=439km,远地点s2= 2384km s1 s2 地球半径r=6371km r 需要作数值积分 s1=439km, s2= 2384km, r=6371km y x o s1 s2 r s1 s2 y x o r a c b 数值积分实例 人造卫星轨道长度 用梯形公式和辛普森公式计算 只将区间5等分，梯形公式就给出很好的结果 轨道长度 L=4.8707?104千米 数值积分实例 人造卫星轨道长度 Jifen2.m 布置实验 目的 1、掌握用MATLAB计算拉格朗日、分段线性、三次样条三 种插值的方法，改变节点的数目，对三种插值结果进行初步分析。 2、 掌握用MATLAB及梯形公式、辛普森公式计算数值积分。 3、 通过实例学习用插值和数值积分解决实际问题。 内容 实验练习 10；12； 大学数学实验 Mathematical Experiments 实验3 插值与数值积分 计算机会“算”吗？靠得住吗？ 例：把4开n次方，再平方n次，结果是4？存在误差？ 英国著名数值分析学家 Higham (1998): Can you count on computers? 精确计算： 解析结果 (Analytical) 近似计算： 数值结果 (Numerical) n=55左右：结果变成1 计算功效= 计算工具* 计算方法(算法) 浮点运算：舍入误差 实验3的基本内容 3.数值积分的梯形公式、辛普森公式和高斯公式。 1.插值的基本原理； 三种插值方法：拉格朗日插 值，分段线性 插值，三次样条插值。 2.插值的 MATLAB 实现及插值的应用。 4.数值积分的 MATLAB 实现及数值积分的应用。 什么是插值(Interpolation)？从查函数表说起 查 函 数 表 标准正态分布函数表 求 ? (1.114) ?(1.114)=0.8665? (0.8686?0.8665)?0.4=0.8673 插值 插值在图像处理/数控加工/外观设计等领域有重要应用 插值的基本原理 插值问题的提法 已知 n+1个节点 其中 互不相同，不妨设 求任一插值点 处的插值 ? ? ? ? ? 节点可视为由 产生, g表达式复杂,甚至无表达式 ? ? ? ? ? ? 求解插值问题的基本思路 构造一个(相对简单的)函数 通过全部节