插值和数值积分汇.pptVIP

大学数学实验;计算机会“算”吗？靠得住吗？;实验3的基本内容;什么是插值(Interpolation)？从查函数表说起;插值的基本原理;?;1.拉格朗日(Lagrange)多项式插值;1.1 拉格朗日插值多项式;*;*;*;;1.3 拉格朗日插值多项式的振荡;2.分段线性插值;MATLAB xch11，xch12，xch13， xch14．;比分段线性插值更光滑;;3. 三次样条插值;三次样条插值确定4n个系数需增加 2个条件;思考;例;三种插值方法小结;1. 拉格朗日插值:自编程序,如名为 lagr.m 的M文件， 第一行为 function y=lagr(x0,y0,x) 输入:节点x0,y0, 插值点x (均为数组，长度自定义)）； 输出:插值y (与x同长度数组)）。 应用时输入x0,y0,x后,运行 y=lagr(x0,y0,x);用MATLAB作插值计算;插值的应用;例：从1点到12点的11小时内，每隔1小时测量一次温度，测得的温度的数值依次为：5，8，9，15，25，29，31，30，22，25，27，24．试估计每隔1/10小时的温度值。;?;二维插值的定义; ;第二种（散乱节点）：;已知n个节点; 注意：最邻近插值一般不连续．具有连续性的最简单的插值是分片线性插值．; 将四个插值点（矩形的四个顶点）处的函数值依次简记为： ;插值函数为：; 双线性插值是一片一片的空间二次曲面构成． 双线性插值函数的形式如下：; 要求x0,y0单调；x，y可取为矩阵，或x取行向量，y取为列向量，x,y的值分别不能超出x0,y0的范围．;例：测得平板表面3×5网格点处的温度分别为： 82 81 80 82 84 79 63 61 65 81 84 84 82 85 86 试作出平板表面的温度分布曲面z=f(x,y)的图形．;再输入以下命令: xi=1:0.2:5; yi=1:0.2:3; zi=interp2(x,y,temps,xi,yi,cubic); mesh(xi,yi,zi) 画出插值后的温度分布曲面图. ; 通过此例对最近邻点插值、双线性插值方法和双三次插值方法的插值效果进行比较．; 插值函数griddata格式为:; 例 在某海域测得一些点(x,y)处的水深z由??表给出，船的吃水深度为5英尺，在矩形区域（75，200）×（-50，150）里的哪些地方船要避免进入．;;*;*;数 值 积 分 的 基 本 思 路;为什么要作数值积分;1.从矩形公式到梯形公式;2.辛普森(Simpson)公式 （抛物线公式）; 对k求和(共m段)，得（复合）辛普森公式：;梯形公式在每小段上是用线性插值函数T(x)代替 f(x); ;同理可得：;梯形公式和辛普森公式的收敛性;*;积分步长的自动选取;高斯(Gauss)求积公式;梯形公式的代数精度（考察T1）;高斯公式的思路;将f(x)代入计算得;;思路：将积分区间分小，在小区间上用n不太 大的 。而在节点加密一倍时能够利用原节点的函数值，可以把区间的端点作为固定节点。;用MATLAB 作数值积分;用MATLAB 作数值积分;矩形域上计算二重积分的命令：;用MATLAB 作数值积分;数值积分的应用;s1=439km, s2= 2384km, r=6371km;用梯形公式和辛普森公式计算;*;*;*;*;*;;布置实验