沙子的开采问题。
摘要。本文主要研究的是估算某一地区的含沙量，并做出开采沙子的收益分析。
针对问题一，我们根据题中**中已有的数据，利用线性插值的方法对相邻两点作线性外推，估算出了缺少的数据。然后通过matlab画出地表、沙层顶部以及沙层低部的等高线。为了更精确的计算出表土层与沙层的体积，用matlab对数据进行线性插值，将50米50米的方格精确到1米1米的小方块。
再利用体积公式近似计算出表土层与沙层的体积，进而得到该矿含沙量超过400万吨，表土体积小于沙子体积19%，可以对这项工程进行投资。
针对问题二，构造单位面积的利润函数，利用matlab软件画出利润等值线。由于贴现的缘故，施工应按照利润由大到小的顺序进行。所以施工顺序是按照等值线的分布进行的。
在确定每个阶段的施工时间后，分别计算利润，进而求得总利润。
关键词： 沙子的开采线性插值 matlab 贴现率。
一。问题简述。
沙子作为最重要建筑材料之一，一直有着稳定的市场。沙子在形成之后上面覆盖着表土，由于表土层结构松散，沙子只能露天开采即先将表土移走再挖沙子。现有一开采地，四周有几块沼泽地。
蓝色隧道公司为生产建筑材料委托一支地质队调查该开采地，在开采地上画的网格点处用空心管垂直钻入地下．利用进入空心管的厚度的开采地的面积可得出表土和沙子的体积，其中每个网格点中最上面一项是地表高度，中间一项是沙层顶部的高度，最下面一项是沙层底部的高度．由于地表存在沼泽地，部分网格点的数据没有采集，即所给表中的数据空白处．
现有该公司提供的相关材料，需要解决以下工作：
1）通过计算表土和沙子体积决定是否该对这项工程进行投资，并画出地表、沙层顶部以及沙层低部的等高线。
2）如果可以投资，根据已给数据和条件设计开采方案，做出收益分析作并求出相应的利润。
二。问题分析。
问题一：开采地四周的沼泽地也存在沙子，所以要补充完整**中的数据，即分别把表土高度、沙层顶部高度、沙层底部高度的数据估算出来。根据网格中已有的相关数据，取前后或上下两个连着的数据，利用线性插值法进行数据处理，然后用matlab编程运算分别得出表土高度、沙层顶部高度、沙层底部高度的等高线图。
题中给出的地表高度、沙层顶部和沙层底部高度的原始数据，相临的点相距50米。为了更精确地计算表层体积，对沙子的储量作精确的估算，用matlab软件对数据进行插值，使得相临的点彼此相距1米，每个网格是11的小方块，经过插值后的数据可更精确地反映沙层厚度分布情形。我们将所有的沙土汇集在单位网格，网格点补充完整后就根据网格点中的数据分别算出每个小区域的表土厚度（表土高度-沙层顶部高度）和沙层厚度（沙层顶部高度-沙层底部高度），利用体积公式累加求出沙层、土层的总体积。
问题二：为解决开采沙子的收益分析，首先要构造利润函数。在不考虑成本中每月10万元的地区管理费的情况下，单位面积利润函数(单位：
元/米2 )只存在沙和土的密度与厚度这俩个变量，沙的销售单价，沙和土的开挖费和运输、处理费都为固定值。根据所得利润函数利用matlab软件求得函数的等值线。由于每个月交纳的10万元地区管理费，再考虑到每月最多开挖量10万吨的限制，只作出了该函数大于等于0的范围，我们的施工设计局限于该范围。
由于10万吨的限制将施工区域分成几个部分，按阶段施工，分别求出各阶段的利润，求和得出总利润。
三。模型假设。
1. 假设题目中所给数据都是准确的。
2. 开采地沙层与表土层都是连续的，近几年不会发生地质以及自然灾害造成地质构造的变化。
3. 月贴现只针对利润，并不包括其他费用。
4. 沙子的单价，开采费用以及其他投入**不会发生改变。
5. 每月利润当月结算，所得收入在月底统一贴现。
四。符号说明。
五。模型的建立与求解。
4.1 问题一：
要解决表土和沙子的体积就要填充网格中所空白数据，首先对网格中的所给出的数据进行分析，利用二维线性插值法较精确地估算出**中因沼泽地原因无法测得的数据。二维线性插值法方法如下：
首先在 x方向进行线性插值，得到：
然后在 y方向进行线性插值，得到：
这样就得到所要的结果：
即得到如下数据**：
再根据上述所给的**，利用matalab软件分别作出地表高度、沙层顶部高度、沙层底部高度的等高线图（matlab数据程序见附录）插值后的地表高度结果见图5-1：
插值后的沙层顶部高度结果见图5-2：
插值后的沙层底部高度结果见图5-3：
对沙子体积的计算采用最基本的体积计算公式：
体积=底面积×高度。
题中给出的地表高度、沙层顶部和沙层底部高度的原始数据，相临的点彼此相距50米。为了更精确地反映地质结构，对沙子的储量作精确的估算，用数学软件包matalab对沙层厚度进行插值，使得相临的点彼此相距1米，每个网格是1×1(米)的小方块，经过插值后的数据可更精确地反映沙层厚度分布情形。由合理假设，该区域没有复杂的地质构造。
可以认为经过插值后的在1×1(米)网格点处的函数值，是以此点为中心的单位方块内沙层厚度的平均值。计算出表土的体积是867480立方米，沙子的体积是4754900立方米。（**见附录）
再利用质量最基本的计算公式：
质量=密度×体积。
计算沙子的质量。由上述计算出的体积以及密度的值得：
表土的质量：
m1==1350千克/立方米×867480立方米。
1171098吨。
沙子的质量：
m2===1620千克/立方米×4754900立方米。
7702938吨。
因为沙子的质量为7702938吨》400万吨，沙子的体积是4754900立方米，475490019％=903431立方米。表土的体积是867480立方米<903431立方米，所以该项目值得投资。
5.2 问题二：
构造利润函数来解决开采沙子的收益分析，在不考虑成本中每月10万磅的地区管理费的情况下，单位面积利润函数(单位：磅 /平方米 )其中的x i和yi分别是沙和土的密度与厚度，m是沙的销售单价，ai和bi分别是沙和土的开挖费和运输、处理费(i =1，2)。
单位面积利润函数只存在沙和土的密度与厚度这两个变量，沙的销售单价，沙和土的开挖费和运输、处理费都为固定值。根据所得利润函数利用matlab软件求得函数的等值线。由于每个月交纳的10万元地区管理费，再考虑到每月最多开挖量10万吨的限制， 该函数的等值线图如图5-4所示。
图2只作出了该函数大于0的范围，我们的施工设计也局限于该范围。由于等值线将施工区域分成几个部分，计划按阶段施工，各阶段数据如表1所示。
在计算贴现后的利润时，要用到月贴现率，在进行工程预算时，对每个月的利润经过现金贴现得到纯现值，再相加得总利润。在计算工程开始后第n个月利润的纯现值，则用公式。
其中i是月贴现率，p是这个月的利润值。
根据各阶段的利润，求和得出总利润为10189579.723磅。
六。模型的评价与改进。
6.1模型的优点。
本题选择的模型简单易懂，计算简便，同时也结合了实际生产情况，对沙子的开挖过程具有一定的指导意义。
问题一通过matlab软件对沙层厚度进行线性插值，使得相邻的点彼此距离为1m的，插值后的数据可更精确的地反映沙层厚度分布情况。
问题二巧妙地结合的问题一求等高线的要求，求解得到了实际开挖中每阶段可挖采的范围。
6.2模型的缺点。
为计算方便，我们认为开挖方式是一次挖到沙层的底部，可认为是垂直开采，不考虑力学上的要求，故最终计算的利润会存在一定误差。
为结合实际生产情况，我们选择相邻阶段的施工区域相邻，而不是优先选择利润最大的位置，这会导致所求利润不是理论上的最大值。
参考文献。1】姜启源，谢金星，叶俊。 数学模型（第三版）. 高等教育出版社，2003.8
附录。附录一等高线绘制程序。
用matlab读入储存在excel中的数据。
附录1.1 绘制地表等高线。
dibiao=data(1:3:27,:)
x,y]=meshgrid([1:15],[1:9]);
contour(dibiao);
c,h]=contour(dibiao);
clabel(c,h);
title('地表等高线')
附录1.2 绘制顶部等高线。
dingbu=data(2:3:27,:)
x,y]=meshgrid([1:15],[1:9]);
contour(dingbu);
c,h]=contour(dingbu);
clabel(c,h);
title('顶部等高线')
附录1.3 绘制底部等高线。
dibu=data(3:3:27,:)
x,y]=meshgrid([1:15],[1:9]);
contour(dibu);
c,h]=contour(dibu);
clabel(c,h);
title('底部等高线')
附录二求表土、沙子的体积。
x0=[0:50:700];
y0=[0:50:400];
h1=data(1:3:27,:)data(2:3:27,:)
h2=data(2:3:27,:)data(3:3:27,:)
xi=0:1:700加密数据点。
yi=0:1:400;
zi1=interp2(x0,y0,h1,xi',yi,'linear');
zi2=interp2(x0,y0,h2,xi',yi,'linear');
tj1=sum(zi1表土体积。
tj2=sum(zi2沙子体积。
鞋子里的沙子作文
各位读友大家好！你有你的木棉，我有我的文章，为了你的木棉，应读我的文章！若为比翼双飞鸟，定是人间有情人！若读此篇优秀文，必成天上比翼鸟！有位哲人说过 让你疲倦的不是远方的高耸的山峰，而是你鞋子中的一粒沙子。仿佛这听起来有些不可思议，也许一粒微小的沙粒在我们眼中是那么的不起眼，但它着实能让一个人在人生...
大三数学建模我们的建模
2013年第六届 认证杯 数学中国。数学建模网络挑战赛。题目护岸框架 关键词河道整治护岸固堤四面六边透水框架建模 摘要 四面六边透水框架是一种新型的江河透水护岸工程技术。通过室内试验，研究了四面六边透水框架固堤护岸的可行性 计算分析了透水框架的布设方式 结构尺寸对流水减速落淤的影响。本文针对江河岸堤...
《沙子的命运》读后感
沙子的命运 读后感。今天，我读了 沙子的命运 这个故事，深有感触，命运读后感。在很久很久以前，有一个养蚌人，想培育出一颗最大最美的珍珠。于是，他就开始寻找沙粒。可沙粒们一颗一颗都说不愿意，因为蚌壳内阴冷潮湿又缺少空气。终于，一粒沙子答应了，因为，它一直想成为一颗珍珠。几年过去了，那颗沙粒终于成了一颗...