数学建模课件——数值计算方法要点.pptVIP

一、线性插值（n=1，一次插值） 二、抛物线插值（n=2，二次插值） 二维插值的定义 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? x y O 第一种（网格节点）： 例2 用下面一组数据拟合 中的参数a，b，k 该问题即解的最优化问题： 1）编写M-文件 curvefun1.m function f=curvefun1(x,tdata) f=x(1)+x(2)*exp(-0.02*x(3)*tdata) %其中 x(1)=a; x(2)=b；x(3)=k; 2）输入命令 tdata=100:100:1000 cdata=1e-03*[4.54,4.99,5.35,5.65,5.90,6.10,6.26,6.39, 6.50,6.59]; x0=[0.2,0.05,0.05]; x=lsqcurvefit (curvefun1,x0,tdata,cdata) f= curvefun1(x,tdata) F(x，tdata)= ，x=(a，b，k) 解法1. 用命令lsqcurvefit 3）运算结果： f =0.0043 0.0051 0.0056 0.0059 0.0061 0.0062 0.0062 0.0063 0.0063 0.0063 x =0.0063 -0.0034 0.2542 4）结论：a=0.0063, b=-0.0034, k=0.2542 1）编写M-文件 curvefun2.m function f=curvefun2(x) tdata=100:100:1000; cdata=1e-03*[4.54,4.99,5.35,5.65,5.90, 6.10,6.26,6.39,6.50,6.59]; f=x(1)+x(2)*exp(-0.02*x(3)*tdata)- cdata 2）输入命令: x0=[0.2,0.05,0.05]; x=lsqnonlin(curvefun2,x0) f= curvefun2(x) 函数curvefun2的自变量是x，cdata和tdata是已知参数，故应将cdata tdata的值写在curvefun2.m中 解法 2 用命令lsqnonlin x=（a，b，k） 3）运算结果为 f =1.0e-003 *（0.2322 -0.1243 -0.2495 -0.2413 -0.1668 -0.0724 0.0241 0.1159 0.2030 0.2792） x =0.0063 -0.0034 0.2542 可以看出，两个命令的计算结果是相同的。 4）结论：即拟合得a=0.0063 b=-0.0034 k=0.2542 设有一个年产240吨的食品加工厂, 需要统计其生产费用, 但由于该厂各项资料不全, 无法统计。在这种情况下, 统计部门收集了设备、生产能力和该厂大致相同的五个食品加工厂的产量与生产费用资料如下表： 如何根据已知五个食品加工厂的资料较准确地估计该厂的生产费用呢？ 引例 插值问题 （2）有的函数甚至没有表达式，只是一种表格函数， （1）如果函数表达式本身比较复杂，且需要多次 实际问题中经常要涉及到函数值的计算问题： 重复计算时，计算量会很大； 方便且表达简单的函数来近似代替，这就是插值问题。 实际对于这两种情况，我们都需要寻找一个计算 而我们需要的函数值可能不在该表格中。 问题：插值与拟合的区别？ 一 维 插 值 一、插值的定义 二、插值的方法 三、用Matlab解插值问题 拉格朗日插值 牛顿插值 反插值 三次样条插值 分段插值法 二维插值 一、二维插值定义 二、网格节点插值法 三、用Matlab解插值问题 最邻近插值 分片线性插值 双线性插值 网格节点数据的插值 散点数据的插值 一维插值的定义 已知 n+1个节点 其中 互不相同，不妨设 求任一插值点 处的插值 ? ? ? ? ? ? 构造一个(相对简单的)函数 通过全部节点, 即 再用 计算插值，即 ? ? ? ? ? ? 代数插值的唯一性 是惟一的。 定理 个不同节点，满足插值条件 的n次插值多项式 当选择代数多项式作为插值函数时，称为代数多 项式插值。 代数插值 求解 L1(x)=a1 x+a0 已知 使得 L1(xi) = yi