数学建模实验9.docx

《数学建模实验》实验报告学号： 姓名： 实验九：插值与拟合2用给定的多项式，如,产生一组数据(xi,yi),i=1,2,..,m,再在yi上添加随机干扰(可用rand产生[0,1]区间均匀分布随机数，或用randn产生N(0,1)分布随机数)，然后用xi和添加了随机干扰的yi作3次多项式拟合，与原系数比较。如果作2或4次多项式拟合，结果如何？解答：Matlab程序如下 (1)先建立M文件fun.mfunction f=fun(x)f=x^3-6*x^2+5*x-3主程序NiheNum.m:for n=1:10 fun(n);end(2)在yi上添加随机干扰后产生的数据：建立M文件addDisturb.mfor n=1:10 y=fun(n)+randend用所得的数据进行多项式拟合：x=1:10;y=[-2.940 -8.9682 -14.7231 -14.9538 -2.9029 27.8235 81.6948 165.3171 285.9502 447.1988];A=polyfit(x,y,3)z=polyval(A,x);plot(x,y,k+,x,z,r)（3）作2次多项式拟合：x=1:10;y=[-2.940 -8.9682 -14.7231 -14.9538 -2.9029 27.8235 81.6948 165.3171 285.9502 447.1988];A=polyfit(x,y,2)z=polyval(A,x);plot(x,y,k+,x,z,r)（4）作4次多项式拟合：x=1:10;y=[-2.940 -8.9682 -14.7231 -14.9538 -2.9029 27.8235 81.6948 165.3171 285.9502 447.1988];A=polyfit(x,y,4)z=polyval(A,x);plot(x,y,k+,x,z,r)结果分析：以给定的多项式求出一组值 f=[-3 -9 -15 -15 -3 27 81 165 285 477]对每个xi求得的相应yi的值的列表如下yi=[-3 -9 -15 -15 -3 27 81 165 285 477]在yi上添加随即干扰后产生的数据：yi= [-2.940 -8.9682 -14.7231 -14.9538 -2.9029 27.8235 81.6948 165.3171 285.9502 447.1988];拟合得到运行结果：A= 0.9928 -5.8919 4.6375 -2.6296与原函数1 -6 5 -3很接近。3次拟合图像如图： 2次拟合得到2次多项式系数A= 10.4887 -70.9118 82.5494图像如图：4次拟合得到4次多项式系数A=-0.0012 1.0192 -6.0854 5.1663 -3.0421图像如图：由图像可以看出，3次拟合后的系数与元系数很接近，从2次和4次多项式拟合的图像上可以看出，词数越高，图像的拟合效果越好；反之，次数越低，图像的拟合越差。2.用最小二乘法求一形如y=aebx的经验公式拟合下表中的数据。xi15.320.527.436.649.165.687.87117.6解答：matlab程序如下：x=1:1:8;%a1=b,a2=log(a);lnc=lna+b*xc=[15.3 20.5 27.4 36.6 49.1 65.6 87.87 117.6];y=log(c);m=polyfit(x,y,1)b=m(1)a=exp(m(2))运行程序得到结果：a=11.4358 b=0.2913即3.(水箱水流量问题)许多供水单位由于没有测量流入或流出水箱流量的设备，因此只能测量水箱中的水位。试通过测得的某时刻水箱中水位的数据，估计在任意时刻t(包括水泵灌水期间)流出水箱的流量f(t).给出原始数据表(见下表)，其中长度单位为(1E=30.24cm)。水箱为圆柱体，其直径为57E.假设：影响水箱流量 的唯一因素是该区公众对水的普通需要;水泵的灌水速度为常数；从水箱中流出的水的最大流速小于水泵的灌水速度；每天的用水量分布都是相似的；水箱的流水速度可用光滑曲线来近似；当水箱的水容量达到514×103g时，开始泵水，达到 677.6×103g时，便停止泵水。时间/s水位/10-2E时间/s水位/10-2E031754463633503316311049953326066353054539363167106192994572543087139372947605743012179212892645542927212402850685352842252232795718542767285432752750212697