数学建模 拟合与插值.pptVIP

差分与差商概念 差分表 差商表 . xk f(xk) 一阶差商 二阶差商 三阶差商 …… n 阶差商 等距节点的差分性质 与拟合有关的MATLAB 函数 polyfit:多项式拟合 poly2sym:由多项式系数向量得多项式符号表达式 polyval:计算多项式函数在指定处的函数值 poly:计算过固定点的多项式 lsqcurvefit, lsqnonlin非线性最小二乘拟合 fit,fittype一般曲线拟合 Lagrange插值与Newton插值的异同点 Cubic Spline Interpolation Lagrange Interpolation 习题 1.某公交公司1路车过去20个季度内的客流量(单位：百万)如下表： * * 第5章作业 数学系《数学建模与数学实验》国家精品课程学习网站 拟合与插值 第5章 数值分析法建模 1. 拟合的基本原理；线性最小二乘拟合。 3. 面对一个实际问题，判断应该用插值，还是拟合。 2. 插值的基本原理； 三种插值方法：拉格朗日插值，Newton插值，三次样条插值。 曲 线 拟 合 问 题 已知一组（二维）数据，即平面上m个点（xi,yi) i=1,…m, 寻求一个函数（曲线）y=f(x), 使 f(x) 在某种准则下与所有数据点最为接近，即曲线拟合得最好。 + + + + + + + + + x y y=f(x) (xi,yi) ?i ?i 为点（xi,yi) 与曲线 y=f(x) 的距离 曲线拟合问题最常用的解法——线性最小二乘法的基本思路 先选定一组函数 r1(x), r2(x), …rn(x), nm, 令 f(x)=a1r1(x)+a2r2(x)+ …+anrn(x) （1） 其中 a1,a2, …an 为待定系数。 确定a1,a2, …an 的准则（最小二乘准则）： 使n个点（xi,yi) 与曲线 y=f(x) 的距离?i 的平方和最小 。 记 问题归结为，求 a1,a2, …an 使 J(a1,a2, …an) 最小。 线性最小二乘法的求解 记 当ATA可逆时，（4）有唯一解： 线性最小二乘拟合 f(x)=a1r1(x)+ …+anrn(x)中函数{r1(x), …rn(x)}的选取 1. 通过机理分析建立数学模型来确定 f(x)； + + + + + + + + + + + + + + + + + + + + y=a1+a2x (a20) 2. 作数据散点图，通过直观判断确定 f(x) + + + + + y=a1+a2x+a3x2（a30） + + + + + y=a1+a2x+a3x2(a30) y=a1+a2/x y=a1exp(a2x) （a20） y=a1exp(a2x) (a30) 1.曲线改直技巧 2.多项式阶数的确定 曲线拟合的几个问题 一阶前向差分 m阶前向差分 一阶差商 m阶差商 多项式拟合 差商与导数的关系 性质1(差分与函数值的关系) 各阶差分均可表示为函数值的线性组合 性质2(多项式的差分) 若f(x)∈Pn(n次多项式类), 则 性质3(差分与差商的关系): 性质4(差分与导数的关系） 1. 作多项式f(x)=anxn+ …+a1x+a0拟合,可利用Matlab命令: a=polyfit(x,y,n) 2. 对超定方程组 可得最小二乘意义下的解。 用 yy=polyval(a,xx)….计算多项式a在xx处的值 当精确函数 y = f(x) 非常复杂或未知时，在区间[a,b]上一系列节点 x0 … xm 处测得函数值 y0 = f(x1), …, ym = f(xm)，由此构造一个简单易算的 近似函数 g(x) ? f(x)，满足条件 g(xj) = f(xj) (j = 0, … m) (*) 这个问题称为“插值问题” 插值问题 g(x) 称为f(x) 的插值函数,一般取多项式函数。 x0 … xm称为插值节点, 条件(*)称为插值条件，区间[a,b]称为插值区间 x0 x0 x2 xm-1 xm x f(x) g(x) 精确函数 y = f(x) 非常复杂或未知 g(x) 称为f(x) 的插值函数 节点 x0 … xm称为插值节点 基本思想：在n次多项式空间Pn中找一组合适的基函数 ?0(x),?1(x),…, ?n(x),使 pn(x)=a0 ?0(x) +a1 ?1(x) +…+an ?n(x) 不同的基函数的