沙子的开采-数学建模.docVIP

沙子的开采问题 摘要 本文主要研究的是估算某一地区的含沙量，并做出开采沙子的收益分析。 针对问题一，我们根据题中表格中已有的数据，利用线性插值的方法对相邻两点作线性外推，估算出了缺少的数据。然后通过matlab画出地表、沙层顶部以及沙层低部的等高线。为了更精确的计算出表土层与沙层的体积，用matlab对数据进行线性插值，将50米50米的方格精确到1米1米的小方块。再利用体积公式近似计算出表土层与沙层的体积，进而得到该矿含沙量超过400万吨，表土体积小于沙子体积19%，可以对这项工程进行投资。 针对问题二，构造单位面积的利润函数，利用matlab软件画出利润等值线。由于贴现的缘故，施工应按照利润由大到小的顺序进行。所以施工顺序是按照等值线的分布进行的。在确定每个阶段的施工时间后，分别计算利润，进而求得总利润。 关键词： 沙子的开采 线性插值 matlab 贴现率 一.问题简述 沙子作为最重要建筑材料之一，一直有着稳定的市场。沙子在形成之后上面覆盖着表土，由于表土层结构松散，沙子只能露天开采即先将表土移走再挖沙子。现有一开采地，四周有几块沼泽地。蓝公司委托一支地质队的网格点处用空心管垂直钻入地．利用进入空心管，其中每个网格点中最上面一项是地表高度，中间一项是沙层顶部的高度，最下面一项是沙层底部的高度．由于地表，部分网格点的数据没有采集，即表中的空白处．现相关材料，：（1）（2）如果可以投资，并求出相应的利润 二.问题分析 问题一：要补充完整表格分别把表土高度、沙层顶部高度、沙层底部高度的数据出来根据网格中已有的相关数据，取数据利用插值法进行数据处理，然后用Matlab编程运算分别得出表土高度、沙层顶部高度、沙层底部高度的图。1的小方块，经过插值后的数据可更精确地反映沙层厚度分布情形。我们将所有的沙土汇集在单位网格网格点补充完整后就根据网格点中的数据分别算出每个小区域的表土厚度（表土高度-沙层顶部高度）和沙层厚度（沙层顶部高度-沙层底部高度）累加求出沙层、土层。10万元的地区管理费的情况下，单位面积利润函数(单位: 元/米2 )只存在沙和土的密度与厚度这俩个变量，沙的销售单价，沙和土的开挖费和运输、处理费都为固定值。根据所得利润函数利用matlab软件求得函数的等值线。由于每个月交纳的10万元地区管理费,再考虑到每月最多开挖量10万吨的限制,只作出了该函数大于等于0的范围,我们的施工设计局限于该范围。由于10万吨的限制将施工区域分成几个部分，按阶段施工，分别求出各阶段的利润，求和得出总利润。 三.模型假设 假设题目中所给数据都是准确的。 开采地沙层与表土层都是连续的，近几年不会发生地质以及自然灾害造成地质构造的变化。 月贴现只针对利润，并不包括其他费用。 沙子的单价，开采费用以及其他投入价格不会发生改变。 5. 每月利润当月结算，所得收入在月底统一贴现。 四.符号说明 符号 符号说明 单位 m1 m2 p 表土的质量 表土的密度 表土的体积 沙子的质量 沙子的密度 沙子的体积 单位面积利润函数 开采单位面积沙土所获利润 吨 吨/立方米 立方米 吨 吨/立方米 立方米 磅/平方米 磅 五.模型的建立与求解 问题一： 要解决表土和沙子的体积就要填充网格中所空白数据对网格中的所给出的数据进行分析，首先在 方向进行线性插值，得到 然后在 方向进行线性插值 这样就得到所要的结果 即得到如下数据表格： A B C D E F G H I J K L M N O 0 1 2 3 4 5 6 7 8 22.3 22.04 21.9 22.0 22.4 23.2 23.4 23.2 22.8 22.5 23.0 23.2 23.2 22.6 21.9 19.1 19.0 19.3 19.5 20.0 20.6 20.4 19.9 19.35 18.4 17.8 18.0 19.0 19.0 18.8 16.4 14.2 12.1 8.7 5.8 4.1 2.1 1.1 0.5 0.4 0.4 0.6 0.2 0.2 3.6 22.4 22.4 22.5 22.6 23.0 23.1 23.2 23.4 23.5 24.0 24.0