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

实验六拟合与插值问题会计学第1页/共47页二、实验目的进一步巩固、加强Matlab的应用能力、学会用MATLAB软件进行数据拟合了解在最小二乘意义下数据拟合的理论和方法.通过对实际问题的分析和研究，初步掌握建立数据拟合数学模型的方法插值的基本原理第2页/共47页三、实验内容与步骤1、建模实例：传染病模型2、MATLAB求函数的拟合与插值 目 的第3页/共47页1、传染病模型 描述传染病的传播过程 分析受感染人数的变化规律 预报传染病高潮到来的时刻 预防传染病蔓延的手段 按照传播过程的一般规律，用机理分析方法建立模型第4页/共47页问题重述问题: 有一种传染病（如SARS、甲型H1N1）正在流行。现在希望建立适当的数学模型，利用已经掌握的一些数据资料对该传染病进行有效地研究，以期对其传播蔓延进行必要的控制，减少人民生命财产的损失。考虑如下的几个问题，建立适当的数学模型，并进行一定的比较分析和评价展望。1、不考虑环境的限制，设单位时间内感染人数的增长率是常数，建立模型求t时刻的感染人数。2、假设环境条件下所允许的最大可感染人数为 。单位时间内感染人数的增长率是感染人数的线性函数，最大感染时的增长率为零。建立模型求t时刻的感染人数。第5页/共47页3、现有卫生防疫部门采集到的某地区一定时间内一定间隔区间的感染人数数据（见下表），利用该数据确定上述两个模型中的相关参数，并将它们的预测值与实际数据进行比较分析（计算仿真偏差）并对两个模型进行适当的评价。（注：该问题中，设最大可感染人数为2000人）第6页/共47页实验问题第7页/共47页模型（一）模型（二）已解决Matlab---微分方程的求解拟合插值今天我们研究问题4的方法年份　1949　1954 1959 1964 1969　人口数　541.67　602.66 672.09 704.99 806.71 　年份1974 1979 1984 1989 1994人口数　908.59 975.42 1034.75　1106.76　1176.74　 第8页/共47页拟 合 问 题 实 例 1据人口统计年鉴，知我国从1949年至1994年人口数据资料如下： (人口数单位为：百万)（1）在直角坐标系上作出人口数的图象。（2）建立人口数与年份的函数关系，并估算1999年的人口数。第9页/共47页线性模型如何确定(x,y)n 已知一组（二维）数据，即平面上的个点，ii ==yf(x)i1,2,,n,xL互不相同，寻求一个函数（曲线），if(x)使在某种准则下与所有数据点最为接近，即曲线拟合: 得最好，如图第10页/共47页一、曲线拟合1 曲线拟合问题的提法: y+确定f(x)使得 ++++++达到最小 +0x最小二乘准则 第11页/共47页２. 用什么样的曲线拟合已知数据?１．画图观察；２．理论分析常用的曲线函数系类型：直线： 多项式： 双曲线（一支): 指数曲线： f=a1+a2xf=a1+a2x+a3x2f=a1+a2x+a3x2+++++++++++++++f=a1+a2/xf=aebxf=ae-bx+++++++++++++++第12页/共47页线性最小二乘拟合 f(x)=a1r1(x)+ …+amrm(x)中函数{r1(x),…,rm(x)}的选取 1. 通过机理分析建立数学模型来确定 f(x)； 2. 将数据 (xi,yi) i=1, …,n 作图，通过直观判断确定 f(x)：第13页/共47页３ 拟合函数组中系数的确定输出拟合多项式系数a=[a1, …,am , am+1] (数组）输入同长度的数组x，y拟合多项式次数第14页/共47页用MATLAB作线性最小二乘拟合1. 作多项式f(x)=a1xm+ …+amx+am+1拟合,可利用已有程序:a=polyfit(x,y,m)2.多项式在x处的值y可用以下命令计算： y=polyval（a，x）第15页/共47页人口预测线性模型对于开始提出的实验问题, 代如数据，计算得从而得到人口数与年份的函数关系为线性预测模型把x=1999代如，估算出1999年的人口数为 y（百万）＝亿1999年实际人口数量为１２.６亿。第16页/共47页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+-)人口模型的解第17页/共47页syms xx=0:1:14y=[39 53 72 96 129 171 232 314 386 502 629 760 920 1065