数学建模实验 五要点解析.doc

西北农林科技大学实验报告 学院名称：理学院 专业年级：2013级信计1班 姓 名： 学 号： 课 程：数学模型与数学建模 报告日期：2015年12月8日 实验五 插值模型与误差估计 实验目的 插值模型是数据挖掘的另一类模型，插值的目的是根据能够获得的观测数据推测缺损的信息，此时观测数据被视为精确地基准数据，寻找一个至少满足连续条件的函数，使得，i=1,2,...,n。在本节我们强调的是插值模型的应用，而不是插值方法的构造。 一、一维插值 一元函数插值公式为 其中是满足的函數，依据插值的方式，如最近邻插值、线性插值、分段三次Hermite插值等，分别取阶梯函数、线性函数、三次多项式函数等，相应的数学表达可以查阅教材。下面先通过简单的MATLAB一维插值指令了解。 ?x=0:4:20;?y=[37?51?45?74?83?88];?xx=0:1:20;?y1=interp1(x,y,xx,nearest);?y2=interp1(x,y,xx);?y3=interp1(x,y,xx,cubic);?y4=interp1(x,y,xx,spline);?plot(y1,DisplayName,y1,YDataSource,y1);figure(gcf)?plot(y2,DisplayName,y2,YDataSource,y2);figure(gcf)?plot(y3,DisplayName,y3,YDataSource,y3);figure(gcf)?plot(y4,DisplayName,y4,YDataSource,y4);figure(gcf) 例一 水库库容量与高程 1 实验题目 设一水库将河道分为上、下游两个河段，降雨的开始时刻为8时，这时水位的高程为168m，水库容量为21.9×10^8m3.预测上游流量Q（t）（m3/s）d取值如表2.2.1所示 表2.2.1 上游流量Q（t）的预测 t/h 8 12 16 24 30 44 48 56 Q/m3s-1 3600 5400 7800 9200 10100 3500 2500 1600 已知水库中水的库容量V(108m3)与水位高程H(m)的数值关系为表2.2.2 表2.2.2 水库库容量与水位高程关系 V/108m3 23.93 24.06 24.12 24.33 24.47 24.6 24.75 H/m 168.75 168.8 168.85 168.9 168.95 169 169.05 如果从当日8时起,水库一直保持每秒1000立方米的泄流量请按所给数据1000立方米的泄流量, 请按所给数据,预报当日20时水库的库容量与水位高程. 2 实验内容 解 为了给出每小时的预报，需要补充每小时整点时刻上游流量的数据，以及相应库容量的水位高程数据 假设:? (a）已知数据准确 (b)相邻两个时刻之间的流量变化是线性的. (c)相邻两个水位高程之间的高程对水的库容量的变化也是线性的. 首先，利用MATLAB线性插值指令，确定每小时的上游流量q(t) ?T=[8,12,16,24,30,44,48,56];?Q=[3600,5400,7800,9200,10100,3500,2500,1600];?t=8:56;?q=interp1(T,Q,t,linear)?plot(q,DisplayName,q,YDataSource,q);figure(gcf) 然后确定每个时刻t的水库容量v(t).因为，水库容量=原库存量+流入量-泄流量（108m3/s）,即 这里我们遇到数值积分，被积函数q（s）没有解析表达式，只是由一个数列表示，qi表示在i时刻的流量。利用MATLAB得到水库容量v在每一时刻t∈[8,56]的值 ?v=21.9+36*10^(-6)*cumtrapz(t,q-10^3*ones(size(t)));?plot(v,DisplayName,v,YDataSource,v);figure(gcf) ?vmax=max(v); 最后确定每个时刻t水库的水位高程h(t).因为最大的水库容量已经远远超出已知数据范围。需要利用外插方法补充数据，确定水库高程对水库容量的依赖关系h=H（v），最后利用函数复合得到水位高程h(t)=H（v(t)） ?V=[21.9?23.93?24.06?24.12?24.33?24.47?24.6?24.75];?H=[168?168.75?168.8?168.85?168.9?168.95?169?169.05];?h=interp1(V,H,v,linear,extrap