实验三：数据插值与拟合实验报告.docVIP

广东金融学院实验报告 课程名称：MATLAB程序设计教程 实验编号 及实验名称 实验三：数据插值与拟合实验 系 别 应用数学系 姓 名 刘瑶琰 学 号 091614203 班 级 信息与计算科学二班 实验地点 新电-403 实验日期 2011-04-18 实验时数 8 指导教师 胡蓉 同组其他成员 无 成 绩 实验目的及要求 了解插值、最小二乘拟合的基本原理 掌握用MATLAB计算一维插值和两种二维插值的方法 掌握用MATLAB作最小二乘多项式拟合和曲线拟合的方法 实验环境及相关情况（包含使用软件、实验设备、主要仪器及材料等） WindowsXP MATLAB软件 实验内容及步骤（包含简要的实验步骤流程） 实验内容： （1）针对实际问题，试建立数学模型。用MATLAB计算一维插值和两种二维插值的方法求解 （2）用MATLAB中的函数作一元函数的多项式拟合与曲线拟合，作出误差图 （3）用MATLAB中的函数作二元函数的最小二乘拟合，作出误差图 （4）针对预测和确定参数的实际问题，建立数学模型，并求解 实验步骤： 开启软件平台——MATLAB，开启MATLAB编辑窗口 根据各种数值解法步骤编写M文件 保存文件并运行 观察运行结果(数值或图形) 根据观察到的结果写出实验报告，并浅谈学习心得体会 实验结果（包括程序或图表、结论陈述、数据记录及分析等，可附页） 见附页 实验总结（包括心得体会、问题回答及实验改进意见，可附页） 通过本次实验 （1）了解了数据插值与拟合的方法，相比之下，通过拟合方法画出的函数图像较光滑，而通过插值方法画出的图像仅仅是有限个离散的数据点连接起来的曲线图。 （2）运用数据的插值与拟合方法，可以建立实际问题的数学模型，对预测和确定实际问题的参数有重要作用。 教师评语 评语 评语等级 优 良 中 及格 不合格 1.实验态度认真，实验目的明确 2.实验方案、程序设计合理 3.实验过程（实验步骤详细，记录完整，数据合理） 4.实验结论正确，分析透彻 5.实验报告独立完成，无抄袭现象，并按时提交，格式规范，文字叙述流畅，逻辑性强 综合评定： 附页 1、（1） x=linspace(0,8.534,11); y=[0.914,5.06,7.772,8.717,9.083,9.144,9.083,8.992,8.687,7.376,2.073]; x1=-0:8.534/99:8.534; y1=interp1(x,y,x1,spline); plot(x,y,’*’,x1,y1,’-’) S=trapz(x1,y1) S = 65.2729 (2) x=0:400:2800; y=0:400:2400; z=[1430 1450 1470 1320 1280 1200 1080 940 1450 1480 1500 1550 1510 1430 1300 1200 1460 1500 1550 1600 1550 1600 1600 1600 1370 1500 1200 1100 1550 1600 1550 1380 1270 1500 1200 1100 1350 1450 1200 1150 1230 1390 1500 1500 1400 900 1100 1060 1180 1320 1450 1420 1400 1300 700 900]; Z=flipud(z); [X,Y]=meshgrid(x,y); mesh(X,Y,Z) [xi,yi]=meshgrid(0:50:2800,0:50:2400); zi=interp2(x,y,Z,xi,yi,cubic);mesh(xi,yi,zi) [xi,yi]=meshgrid(0:50:2800,0:50:2400); zi=interp2(x,y,Z,xi,yi,cubic); meshc (xi,yi,zi) 2、 x=dsolve(Dx=r*x,x(0)=x0,t) x =x0*exp(r*t) 将x=x0*exp(rt)两边取对数可得lnx=lnx0*exp(rt),(lnx=lnx0+rt,另y=lnx,a= lnx0,所以可得y= rt+a。利用1790-1900年的数据进行试验 t=1790:10:19; 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]; p=polyfit(t,log(x),1); r=p(1) x0=exp(p(2)) plot(t,x,+,t, x0*exp(r