数学建模实验9课案.docxVIP

《数学建模实验》实验报告 学号： 姓名： 实验九：插值与拟合2 用给定的多项式，如,产生一组数据(xi,yi),i=1,2,..,m,再在yi上添加随机干扰(可用rand产生[0,1]区间均匀分布随机数，或用randn产生N(0,1)分布随机数)，然后用xi和添加了随机干扰的yi作3次多项式拟合，与原系数比较。如果作2或4次多项式拟合，结果如何？ 解答： Matlab程序如下 (1)先建立M文件fun.m function 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.m for n=1:10 y=fun(n)+rand end 用所得的数据进行多项式拟合： 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的经验公式拟合下表中的数据。 xi12345678yi15.320.527.436.649.165.687.87117.6 解答： matlab程序如下： x=1:1:8;%a1=b,a2=log(a);lnc=lna+b*x c=[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 即y=11.4358e0.2913 3.(水箱水流量问题)许多供水单位由于没有测量流入或流出水箱流量的设备，因此只能测量水箱中的水位。试通过测得的某时刻水箱中水位的数据，估计在任意时刻t(包括水泵灌水期间)流出水箱的流量f(t). 给出原始数据表(见下表)，其中长度单位为(1E=30.24cm)。水箱为圆柱体，其直径为57E. 假设： 影响水箱流量 的唯一因素是该区公众对水的普通需要; 水泵的灌水速度为常数； 从水箱中流出的水的最大流速小于水泵的灌水速度； 每天的用水量分布都是相似的； 水箱的流水速度可用光滑曲线来近似； 当水箱的水容量达到514×103g时，开始泵水，达到 677.6×103g时，便停止泵水。 时间/s水位/10-2E时间/s水位/10-2E03175446363350