大江截流数学建模.docxVIP

大江截流1997年11月8日电视正在播放十分壮观的长江三峡工程大江截流的实况。截流从8：55开始，当时龙口水面宽40m，水深60m。到11：50时，播音员报告宽为34.4m。到13：00时，播音员又报告水面宽31m。这时电视机旁的小明说，现在可以估算下午几点合龙。从8：55到11：50，进展的速度为每小时宽度减少1.9m。从11：50到13：00，每小时宽度减少2.9m。小明认为回填速度是越来越快的，近似地每小时速度加快1m。从下午1：00起，大约要5个多小时，即到下午6点多才能合龙。但到了下午3点28分，电视里传来了振奋人心的消息：大江截流成功！小明后来想明白了，他估算的方法不好。现在请你根据上面的数据设计一种合理的估算方法（建立一种合理的数学模型）进行估算，使你的计算结果更切合实际。matlab程序： t=[0,175,245]; %已知的时间向量 l=[40,34.4,31]; %与时间向量相对应的江面宽度%根据公式1/2*a*t.^2=1200-3/4*l.^2，a为拟合参数，将1/2*t.^2看做x，1200-3/4*l.^2看做y x=1/2*t.^2; %设x=1/2*t^2，所得的横坐标点 y=1200-3/4*l.^2; % 设y=1200-3/4*l.^2，所得的纵坐标值 A=x; d=A\y’%解超定线性方程组Ab=y’d = 0.0169% 求得a值，即回填加速度， a=d; solve(1200-1/2*0.0169*t^2=0,t)%当江面宽度为0时，求得对应的时间ans = -376.84457581279663049188985764706 376.84457581279663049188985764706% 取正值，即从大江截流开始至完成的时间为376.8分钟[b,bint,r,rint,stats]=regress(y,A)%对数据进行拟合，并对拟合效果残差分析b = 0.0169%回填加速度a的值为0.0169bint = 0.0114 0.0224%a的置信区间为[0.0114,0.0224]r = 0 53.9270 -27.5138%相应的残差值rint = -260.4841 260.4841 53.9270 53.9270 -27.5138 -27.5138%残差的95%置信区间stats = 1.0e+003 * 0.0010 NaNNaN 1.8326%R^2决定系数为1.0e+003，残差平方和为1.8326 rcoplot(r,rint)%画出残差图附：拟合y=ax图像：残差分析图分析：由图可知，三个点中有两个点都为异常点，这可能是因为原始数据太少而导致的拟合精度不高，若将大坝合拢的时间点加入，则程序如下：t=[0,175,245,393]; l=[40,34.4,31,0]; x=1/2*t.^2; y=1200-3/4*l.^2; plot(x,y,*) A=x; [d,bint,r,rint,stats]=regress(y,A)d =0.0158% 此时所得的回填加速度值为0.0158bint = 0.0142 0.0174r = 0 71.2414 6.4224 -16.6222rint = -165.2558 165.2558 44.6706 97.8123 -147.3156 160.1605 -71.1022 37.8579stats = 1.0e+003 *0.0010 NaNNaN 1.7976%此时所得的残差平方和为1.7976，小于1.8326 solve(1200-1/2*0.0158*t^2=0,t) ans = -389.74188147697854242909705613752 389.74188147697854242909705613752 rcoplot(r,rint)此时拟合所得的残差图如下，显然，拟合精度有所提高。