20115394王福临数学实验作业(拟合).docVIP

重 庆 大 学 学 生 实 验 报 告 实验课程名称 数学实验 开课实验室 DS1421 学 院 计算机学院 年级 大二 专业 计科三班 学 生 姓 名 王福临 学 号 开 课 时 间 至 学年第 学期 总 成 绩 教师签名 数 理 学 院 制 开课学院、实验室：数理学院DS1401 实验时间 ：2013 年 4 月 22 日 课程 名称 数学实验 实验项目 名 称 拟合 实验项目类型 验证 演示 综合 设计 其他 指导 教师 何光辉 成 绩 一、实验目的及意义 [1] 了解最小二乘拟合的原理，掌握用MATLAB作最小二乘拟合的方法。 [2] 通过实例学习如何用拟合方法解决实际问题，注意与插值方法的区,通过实例理解参数辨识的方法。 [3] 通过实验体验用函数拟合解决实际问题的全过程。 二、实验内容 1.人口问题预测 下表给出了近两个世纪某国人口的统计数据，试预测2010年该国的人口。 2.旧车价格预测 某年美国旧车价格的调查资料如下表其中xi表示轿车的使用年数，yi表示相应的平均价格。试分析用什么形式的曲线来拟合上述的数据，并预测使用4.5年后轿车的平均价格大致为多少？ 三、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 3．保存文件并运行； 4．观察运行结果(数值或图形)； 5．写出实验报告，并浅谈学习心得体会。 实验题目： 1.人口问题预测 下表给出了近两个世纪某国人口的统计数据，试预测2010年该国的人口。 马尔萨斯模型的基本假设是人口的增长率为常数，记为r。记时刻t的人口为，且初始时刻的人口为x0，于是得到如下微分方程 需要先求微分方程的解，再用数据拟合模型中的参数。 一、分析 有这个方程很容易解出 r0时，是表示人口箭杆指数规律随时间无限增长，称为指数增长模型。 将上式取对数，可得y=rt+a，y=lnx(t)，a=lnx0 二、用matlab 编码 t=1790:10:1980; x=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76.0 92 106.5 123.2 131.7 150.7 179.3 204.0 226.5]; p=polyfit(t,log(x),1); r=p(1) x0=exp(p(2)) x1=x0.*exp(r.*t); plot(t,x,r,t,x1,b) 三、结果和图像 exp(Polyval(p,2010)); ans = 622.4063;所以2010年人口为622.4063人. 2.旧车价格预测 某年美国旧车价格的调查资料如下表其中xi表示轿车的使用年数，yi表示相应的平均价格。试分析用什么形式的曲线来拟合上述的数据，并预测使用4.5年后轿车的平均价格大致为多少？ 一、分析 用matlab编码绘制出点图，预测图像大致形状。 分析图形呈现对数形状。 用函数来拟合以上数据。 二、用matlab编码 function f=curvefunl(k,xdata) f=k(1)*log(k(2)*xdata)+k(3); %???Dk(1)=a,k(2)=b,k(3)=c; xdata=1:1:10; ydata=[2615 1943 1494 1087 765 538 484 290 226 204]; k0=[1,1,0]; k=lsqcurvefit(curvefunl,k0,xdata,ydata) xi=linspace(0,10,100); yi=curvefunl(k,xi); plot(xdata,ydata,ro,xi,yi) xlable(x),ylable(y=f(x)), title(nonlinear curve fitting) Y=curvefun1(k,4.5); Y=971.6645 三、结果和图像 a=1.0e+003 * -1.1101 b=1.0e+003 *0.0005 c= 1.0e+003 *1.8307 Y=971.6645 所以第4.5年的平均价格大致为971.6645; 总结与体会 1.由本次实验的联系，将理论与实践充分结合，学习和理解了matlab中拟合问题的求解方法。 2.由应用问题的求解