实验六、拟合和插值问题.ppt

数学模型;进一步巩固、加强Matlab的应用能力、 学会用MATLAB软件进行数据拟合 了解在最小二乘意义下数据拟合的理论和方法. 通过对实际问题的分析和研究，初步掌握建立数据拟合数学模型的方法 插值的基本原理;三、实验内容与步骤; 目 的;1.1问题重述;3、现有卫生防疫部门采集到的某地区一定时间内一定间隔区间的感染人数数据（见下表），利用该数据确定上述两个模型中的相关参数，并将它们的预测值与实际数据进行比较分析（计算仿真偏差）并对两个模型进行适当的评价。（注：该问题中，设最大可感染人数为2000人）;实验问题;模型（一） 模型（二）;据人口统计年鉴，知我国从1949年至1994年人口数据资料如下： (人口数单位为：百万);如何确定;1 曲线拟合问题的提法: ;２. 用什么样的曲线拟合已知数据?;线性最小二乘拟合 f(x)=a1r1(x)+ …+amrm(x)中函数{r1(x),…,rm(x)}的选取 ;３ 拟合函数组中系数的确定;用MATLAB作线性最小二乘拟合;人口预测线性模型;syms xx=1949:5:1994;y=[541.67 602.66 672.09 704.99 806.71 908.59 975.42 1034.75 1106.76 1176.74];plot(x,y,r*’)z=polyfit(x,y,1)y1=polyval(z,x)hold onplot(x,y1,b+-);syms x x=0:1:14 y=[39 53 72 96 129 171 232 314 386 502 629 760 920 1065 1232]; plot(x,y,r*) z=polyfit(x,y,2) y1=polyval(z,x) hold on plot(x,y1,b+-);二、插 值 ;插 值 问 题 实 例 1;插 值 问 题 的 提 法;?;插值的基本原理;1.拉格朗日(Lagrange)多项式插值;1.拉格朗日(Lagrange)多项式插值;1.拉格朗日(Lagrange)多项式插值;1.拉格朗日(Lagrange)多项式插值;1.拉格朗日(Lagrange)多项式插值;2.分段线性插值;3. 三次样条插值;1. 拉格朗日插值:自编程序,如名为 lagr1.m 的M文件， 第一行为 function y=lagr1(x0,y0,x) 输入:节点x0,y0, 插值点x (均为数组，长度自定义)）； 输出:插值y (与x同长度数组)）。 应用时输入x0,y0,x后,运行 y=lagr1(x0,y0,x);对表格给出的函数，求出没有给出的函数值 ;其中x0,y0是已知的节点坐标，是同维向量。y对应于x处的插值。y与x是同维向量。 解决上述问题,我们可分两步: 一 用原始数据绘图作为选用插值方法的参考. 二 确定插值方法进行插值计算;对y=1/(1+x2),-5≤x≤5，用n（=11）个节点（等分）作上述两种插值，用m（=21）个插值点（等分）作图，比较结果。 解：键入并运行如下命令 n=11;m=21;x=-5:10/(m-1):5;y=1./(1+x.^2); xo=-5:10/(n-1):5;yo=1./(1+xo.^2); y1=interp1(xo,yo,x); y2=spline(xo,yo,x); plot(x,y,r,x,y1,b,x,y2,k);例 3 一水库上游河段降暴雨.，根据预报测算上游流入水库的流量为Q(t) (102立方米/秒) ： t (时) 8 12 16 24 30 44 48 56 60 Q( t ) 36 54 78 92 101 35 25 16 13 利用这个预报值估计14：30 和 20：30 时上游流入水库的流量。 假设: 1 已知数据准确。 2 相邻两个时刻之间的流量没有突然的变化。 ;癸裁届克烈酬沼耳沽罢稚援晶棚贤课缉贵薄健佐招懈酵秆捡比闺曾昂褐竖实验六、拟合与插值问题实验六、拟合与插值问题;t=[8,12,16,24,30,44,48,56,60]; q=[36,54,78,92,101,35,25,16,13]; t1=8:0.5:60; q1=interp1(t,q,t1,linear); plot(t,q,b,t1,q1); hold on; q2=interp1(t,q,t1,spline); plot(t,q,b,t1,q1,t1,q2,r) q1 q2;化拖爆椒葬湿耸草词末庞埋戈渤募庆漾歪马猛萍缓吾糟趴昔划剧唆焊箔腾实验六、拟合与插值问题实验六、拟合与插值问题;拟合与插值的关系;最临近插