沙层体积问题.docVIP

求沙层体积的数学模型 摘要：为了求出沙层体积，我们建立了一个回归模型。根据表中所给数据拟合出了沙尘的大致的图像，并将缺的数据用插值法拟和求得，近似求出沙尘的体积，在求体积问题上，我们用回归求出面积和高度的关系，把沙层看由多个小立方体组成，求出小立方体体积之和，即沙层的体积。 关键词：标高 含沙量 数学模型 1．问题的提出： 某地区有优质细沙埋在地下，某公司欲在该地开采钻探资料如下表：在每个格点上有3个数字，分别代表覆盖表层的标高，沙层顶部的标高，沙层底部的标高，阴影部分表示沼泽地地带，没有钻探资料，试估计整个矩形区域内的含沙量。 A B C D E F G H 0 22.4 * * * * 22.5 23.0 23.2 20.0 * * * * 18.4 17.8 17.8 5.8 * * * * 0.5 0.4 0.4 1 23.0 23.1 23.2 23.4 23.5 24.0 24.0 24.0 19.9 20.0 20.0 19.8 19.9 20.0 19.8 19.6 6.6 3.2 1.6 1.0 1.1 1.0 0.8 0.9 2 23.1 23.2 23.4 23.4 23.5 24.2 24.1 24.1 17.8 19.7 19.4 20.0 20.1 20.3 20.3 20.5 2,2 1.4 0.6 0.5 0.3 -0.2 -0.1 0.2 2．问题的分析 由于阴影的沙层没有数据记录，我们利用插值法，求出阴影部分的大致数据（见附录1），然后用MATLAB画出沙层示意图（如下），同时用所给的数据把高插值出来，接着我们把它切割成很多个小立方体，求出小立方体体积和，即沙层总体积。 在沙层体积的求解问题，为了使求出的小立方体更精确，取小立方体高(h)的值为四边高的平均值，长和宽均取步长为1，其体积=1*1*h 3、符号说明 v：沙层的总体积 4、模型的假设 1、沙层的分布是连续的； 5、模型的求解 用MATLAB软件编程（见附录2）循环累加小立方体体之和，即沙层的体积： 6.3063e+005 6附录1 插值求解阴影部分的大致数据： x=[0 250 300 350 ]; y=[22.4 23.5 23.8 23.2] y1=interp1(x,y,[50 100 150 200]) 解得 y1 = 22.6200 22.8400 23.0600 23.2800 插值求解： x=[0 250 300 350 ]; y=[20.0 18.4 17.0 18.0] y2=interp1(x,y,[50 100 150 200]) 解得 y2 = 19.6800 19.3600 19.0400 18.7200 插值求解： x=[0 250 300 350 ]; y=[5.8 0.5 0.4 0.4] y3=interp1(x,y,[50 100 150 200]) 解得：y3 = 4.7400 3.6800 2.6200 1.5600 7．附录2 x=[0 0 0 50 50 50 100 100 100 150 150 150 200 200 200 250 250 250 300 300 300 350 350 350]; y=[0 50 100 0 50 100 0 50 100 0 50 100 0 50 100 0 50 100 0 50 100 0 50 100]; z= [17.6000 13.9000 14.2000 18.3000 16.8000 15.3000 18.8000 18.4000 16.7000 19.5000 18.8000 17.8000 19.8000 18.8000 18.2000 20.5000 19.0000 17.9000 20.4000 18.0000 17.4000 20.5000 18.7000 17.6000]; cx=0:0.5:350; cy=0:0.5:100; cz=griddata(x,y,z,cx,cy,cubic); v=0; i=1; j=1; while i200; while j700; v1=0.5*0.5*(cz(i,j)+cz(i,j+1)+cz(i+1,j)+cz(i+1,j+1))/4; v=v+v1; j=j+1; end i=i+1; j=1; end y 运行结果： v = 6.3063e+005