实验六拟合与插值问题.pptx

《实验六拟合与插值问题.pptx》由会员分享，可在线阅读，更多相关《实验六拟合与插值问题.pptx（47页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、2023/3/161进一步巩固、加强进一步巩固、加强Matlab的应用能力、的应用能力、学会用学会用MATLAB软件进行数据拟合软件进行数据拟合了解在最小二乘意义下数据拟合的理了解在最小二乘意义下数据拟合的理论和方法论和方法.通过对实际问题的分析和研究，初步通过对实际问题的分析和研究，初步掌握建立数据拟合数学模型的方法掌握建立数据拟合数学模型的方法插值的基本原理插值的基本原理二、实验目的二、实验目的第1页/共47页2023/3/162三、实验内容与步三、实验内容与步三、实验内容与步三、实验内容与步骤骤骤骤1 1、建模实例：传染、建模实例：传染病模型病模型2 2、MATLABMATLAB求函求函

2、数的拟合与插值数的拟合与插值第2页/共47页2023/3/163 目目 的的1、传染病模型、传染病模型 描述传染病的传播过程描述传染病的传播过程 分析受感染人数的变化规律分析受感染人数的变化规律 预报传染病高潮到来的时刻预报传染病高潮到来的时刻 预防传染病蔓延的手段预防传染病蔓延的手段 按照传播过程的一般规律，按照传播过程的一般规律，用机理分析方法建立模型用机理分析方法建立模型第3页/共47页2023/3/164问题重述问题重述n问题:有一种传染病（如SARS、甲型H1N1）正在流行。现在希望建立适当的数学模型，利用已经掌握的一些数据资料对该传染病进行有效地研究，以期对其传播蔓延进行必要的控制

3、，减少人民生命财产的损失。考虑如下的几个问题，建立适当的数学模型，并进行一定的比较分析和评价展望。n1、不考虑环境的限制，设单位时间内感染人数的增长率是常数，建立模型求t时刻的感染人数。n2、假设环境条件下所允许的最大可感染人数为 。单位时间内感染人数的增长率是感染人数的线性函数，最大感染时的增长率为零。建立模型求t时刻的感染人数。第4页/共47页2023/3/1653、现有卫生防疫部门采集到的某地区一定时间内一定间隔区间的感染人数数据（见下表），利用该数据确定上述两个模型中的相关参数，并将它们的预测值与实际数据进行比较分析（计算仿真偏差）并对两个模型进行适当的评价。（注：该问题中，设最大可感

4、染人数为2000人）第5页/共47页2023/3/166实验问题实验问题第6页/共47页2023/3/167模型（一）模型（一）模型（二）模型（二）Matlab-Matlab-微分方程的求解微分方程的求解已解决已解决今天我们研究问今天我们研究问题题4 4的方法的方法拟合拟合插值插值第7页/共47页2023/3/168据人口统计年鉴，知我国从据人口统计年鉴，知我国从19491949年至年至19941994年人年人口数据资料如下：口数据资料如下：(人口数单位为：百万人口数单位为：百万)（1 1）在直角坐标系上作出人口数的图象。）在直角坐标系上作出人口数的图象。（2 2）建立人口数与年份的函数关系，

5、并估算）建立人口数与年份的函数关系，并估算19991999年年的人口数。的人口数。拟拟 合合 问问 题题 实实 例例 1年份年份19491954 1959 1964 1969人口数人口数 541.67602.66 672.09 704.99 806.71 年份年份 1974 1979 1984 1989 1994人口数人口数 908.59 975.42 1034.751106.761176.74 第8页/共47页2023/3/169如何确定如何确定线性模型线性模型第9页/共47页2023/3/16101 曲线拟合问题的提法曲线拟合问题的提法:已知一组（二维）数据，即平面上的已知一组（二维）数据

6、，即平面上的n个点个点),(iiyx，ixni,2,1L=互不相同，寻求一个函数（曲线）互不相同，寻求一个函数（曲线）)(xfy=，使使)(xf在某种准则下与所有数据点最为接近，即曲线拟合在某种准则下与所有数据点最为接近，即曲线拟合得最好，如图得最好，如图:xy0+一、曲线拟合一、曲线拟合确定确定f(x)使得使得 达到最小达到最小 最小二乘准则最小二乘准则 第10页/共47页2023/3/1611.用什么样的曲线拟合已知数据用什么样的曲线拟合已知数据?常用的曲线函数系类型：常用的曲线函数系类型：画图观察；画图观察；理论分析理论分析指数曲线：指数曲线：双曲线（一支双曲线（一支):):多项式：多项

7、式：直线：直线：第11页/共47页2023/3/1612线性最小二乘拟合线性最小二乘拟合 f(x)=a1r1(x)+amrm(x)中中函数函数 r1 1(x),),rm(x)的选取的选取 1.1.通过机理分析建立数学模型来确定通过机理分析建立数学模型来确定 f(x)；+f=a1+a2xf=a1+a2x+a3x2f=a1+a2x+a3x2f=a1+a2/xf=aebxf=ae-bx 2.2.将数据将数据 (xi,yi)i=1,n 作图，通过直观判断确定作图，通过直观判断确定 f(x)：第12页/共47页2023/3/1613 拟合函数组中系数的确定拟合函数组中系数的确定第13页/共47页2023

8、/3/1614用用MATLAB作线性最小二乘拟合作线性最小二乘拟合1.1.作多项式作多项式f(x)=a1xm+amx+am+1拟合拟合,可利用已有程序可利用已有程序:a=polyfit(x,y,m)2.2.多项式在多项式在x处的值处的值y可用以下命令计算：可用以下命令计算：y=polyval（a，x）输出拟合多项式系数输出拟合多项式系数a=a1,am,am+1(数组）数组）输入同长度输入同长度的数组的数组x，y拟合多项拟合多项式次数式次数第14页/共47页2023/3/1615人口预测线性模型人口预测线性模型对于开始提出的实验问题对于开始提出的实验问题,代如数据，计算得代如数据，计算得从而得到

9、人口数与年份的函数关系为从而得到人口数与年份的函数关系为把把x=1999代如，估算出代如，估算出1999年的人口数为年的人口数为 y（百万）亿（百万）亿1999年实际人口数量为年实际人口数量为.亿。亿。线性预测模型线性预测模型第15页/共47页2023/3/1616syms 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+-)人口模型的解人口

10、模型的解第16页/共47页2023/3/1617syms xx=0:1:14y=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 onplot(x,y1,b+-)传染病模型传染病模型3 3的解的解第17页/共47页2023/3/1618二、插二、插 值值 2.面面 对对 一一 个个 实实 际际 问问 题，应题，应 该该 用用 插插 值，还值，还 是是 拟拟 合。合。1.插值的基本原理；三种插值方法：拉格朗日插值的基本原理；三种插值方法

11、：拉格朗日插值，分段线性插值，三次样条插值。插值，分段线性插值，三次样条插值。第18页/共47页2023/3/1619插插 值值 问问 题题 实实 例例 1 1xy机翼下轮廓线第19页/共47页2023/3/1620插插 值值 问问 题题 的的 提提 法法已知已知 n+1个节点个节点其中其中互不相同，不妨设互不相同，不妨设求任一插值点求任一插值点处的插值处的插值节点可视为由节点可视为由产生产生,表达式复杂表达式复杂或无封闭形式或无封闭形式,，或未知或未知.。第20页/共47页2023/3/1621求求 解解 插插 值值 问问 题题 的的 基基 本本 思思 路路 构造一个构造一个(相对简单的相对

12、简单的)函数函数通过全部节点通过全部节点,即即再用再用计算插值，即计算插值，即第21页/共47页2023/3/1622插值的基本原理插值的基本原理选讲PPT（24-30）第22页/共47页2023/3/16231.1.拉格朗日(Lagrange)多项式插值1.1 1.1 插值多项式插值多项式有唯一解有唯一解第23页/共47页2023/3/16241.1.拉格朗日(Lagrange)多项式插值1.2 拉格朗插值多项式又（2）有唯一解，故（3）与（1）相同。第24页/共47页2023/3/16251.1.拉格朗日(Lagrange)多项式插值1.3 1.3 误差估计误差估计第25页/共47页202

13、3/3/16261.1.拉格朗日(Lagrange)多项式插值1.4 1.4 例例 将0,/2 n等分，用g(x)=cos(x)产生n+1个节点，作Ln(x)（取n=1,2),计算cos(/6),估计误差。解解:n=1,(x0,y0)=(0,1),(x1,y1)=(/2,0),L1(x)=y0l0+y1l1=1-2x/,n=2,(x0,y0)=(0,1),(x1,y1)=(/4,0.7071),(x2,y2)=(/2,0),L2(x)=y0l0+y1l1+y2l2=8(x-/4)(x-/2)/2-16x(x-/2)0.7071/2cos(cos(/6)=L/6)=L2 2(/6)=0.8508

14、 /6)=0.8508 精确值：精确值：coscos第26页/共47页2023/3/16271.1.拉格朗日(Lagrange)多项式插值 拉格朗日插值多项式的振荡To MATLAB(runge)Runge现象现象:第27页/共47页2023/3/16282.2.分段线性插值分段线性插值xjxj-1xj+1x0 xn计算量与n无关;n越大，误差越小.第28页/共47页2023/3/16293.3.三次样条插值三次样条插值第29页/共47页2023/3/16301.1.拉格朗日插值拉格朗日插值:自编程序自编程序,如名为如名为 的的M文件，文件，第一行为第一行为 function y=lagr1(

15、x0,y0,x)输入输入:节点节点x0,y0,插值点插值点x(均为数组，长度自定义均为数组，长度自定义)）；）；输出输出:插值插值y(与与x同长度数组同长度数组)）。）。应用时输入应用时输入x0,y0,x后后,运行运行 y=lagr1(x0,y0,x)2.2.分段线性插值分段线性插值:已有程序已有程序 y=interp1(x0,y0,x)3.3.三次样条插值三次样条插值:已有程序已有程序 y=interp1(x0,y0,x,spline)或或 y=spline(x0,y0,x)用MATLAB作插值计算第30页/共47页2023/3/1631对表格给出的函数，求出没有给出的函数值 例例2表中是待

16、加工零件下轮廓线的一组数据，现需要得到x坐标每改变时所对应的y的坐标.x03571112131415y01.21.72.12.01.81.21.01.6下面是关于插值的两条命令(专门用来解决这类问题)：y=interp1(x0,y0,x)分段线性插值y=spline(x0,y0,x)三次样条插值 第31页/共47页2023/3/1632其中x0,y0是已知的节点坐标，是同维向量。y对应于x处的插值。y与x是同维向量。解决上述问题,我们可分两步:一 用原始数据绘图作为选用插值方法的参考.二 确定插值方法进行插值计算x0=0,3,5,7,9,11,12,13,14,15;y0=0,1.2,1.7,

17、2.0,2.1,2.0,1.8,1.2,1.0,1.6plot(x0,y0)%完成第一步工作x=0:0.1:15;y=interp1(x0,y0,x);%用分段线性插值完成第二步工作plot(x,y)y=spline(x0,y0,x);plot(x,y)%用三次样条插值完成第二步工作第32页/共47页2023/3/1633对y=1/(1+x2),-5x5，用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+

18、xo.2);y1=interp1(xo,yo,x);y2=spline(xo,yo,x);plot(x,y,r,x,y1,b,x,y2,k)练习：第33页/共47页2023/3/1634例 3 一水库上游河段降暴雨.，根据预报测算上游流入水库的流量为Q(t)(102立方米/秒)：t(时)8 12 16 24 30 44 48 56 60Q(t)36 54 78 92 101 35 25 16 13 利用这个预报值估计14：30 和 20：30 时上游流入水库的流量。假设:1 已知数据准确。2 相邻两个时刻之间的流量没有突然的变化。第34页/共47页2023/3/1635第35页/共47页202

19、3/3/1636t=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)q1q2第36页/共47页2023/3/1637第37页/共47页2023/3/1638拟合与插值的关系拟合与插值的关系 函数插值与曲线拟合都是要根据一组数据构造一个函数作函数插值与曲线拟合都是要根据一组数据构造一个函数作为近似，由于近似

20、的要求不同，二者在数学方法上是完全不同为近似，由于近似的要求不同，二者在数学方法上是完全不同的的 实例：实例：下面数据是某次实验所得，希望得到X和 f之间的关系？问题：问题：给定一批数据点，需确定满足特定要求的曲线或曲面解决方案：解决方案：若不要求曲线（面）通过所有数据点，而是要求它反映对象整体的变化趋势，这就是数据拟合数据拟合，又称曲线拟合或曲面拟合若要求所求曲线（面）通过所给所有数据点，就是插值问题插值问题；第38页/共47页2023/3/1639最临近插值、线性插值、样条插值与曲线拟合结果：最临近插值、线性插值、样条插值与曲线拟合结果：第39页/共47页2023/3/1640即要求即要求

21、 出二次多项式出二次多项式:中中 的的使得使得:1、对下面一组数据作二次多项式拟合、对下面一组数据作二次多项式拟合上机操作（上机操作（1、2拟合；拟合；3、4用插值法）用插值法）第40页/共47页2023/3/16412、最佳营销策略问题：某公司有一批以每桶2元购进的彩漆为了获得较高的利润希望以较高的价格卖出但价格越高,售出量就越少,二者之间的关系由表一给出.于是打算用作广告的办法来促销.而广告费与销售量的关系可由销售增长因子来描述.例如,投入3万元的广告费,销售因子为1.85,意味着做广告后的销售量将是未做广告销售量的倍.根据经验,广告费与销售因子的关系如表2,现请你作出决策：投入多少广告费

22、和售价为多少时所获得的利润最大?第41页/共47页2023/3/1642售价 2.002.503.003.504.004.505.005.506.00预期销售量（千桶）413834322928252220表1广告费（千元 0102030405070销售增长因子 1.001.401.701.851.952.001.80表2第42页/共47页2023/3/1643提示用描点法画出预期销售量售价；销售增长因子广告费的关系图.建立上述两个函数关系.设x为售价,y为广告费,P为所得利润,建立P关于x和y的函数关系.用多元函数的极值理论计算并回答此问题.第43页/共47页2023/3/16443、在某处测得海洋不同深度处水温如下：求深度为500、1000、1500米处的水温。深度44671495014221636水温7.044.283.402.542.13第44页/共47页2023/3/1645xy4、用MATLAB作插值计算机翼下轮廓线已知下轮廓线上数据如下，求已知下轮廓线上数据如下，求x每改变时的每改变时的y值。值。第45页/共47页2023/3/1646第46页/共47页2023/3/1647感谢您的观看！第47页/共47页