数学建模土地平整问题..doc

土地平整问题 摘要 随着社会的快速发展，对于山区的土地整理显得格外重要。为了更好地合理利用土地，选择合适的土地整理地点，达到节约成本的目的。我们会建立模型，根据几个重要的制约因素，使得土石方量最小，而且成本最低，在实际的操作过程中更接近于科学实施。 针对问题一：利用软件MATLAB中的绘图功能，对数据进行处理，得到这片土地的三维图形与等高线图。 针对问题二：采用二重积分的意义，即曲顶柱体的体积。利用分割、取值、求和、取极限的思想，将山地分割成无数个小曲顶柱体，然后对这些小曲顶柱体进行求和，并结合土石方量最小，挖填土时费用使用的最少限制条件，从而求得在什么海拔高度进行平整。对于平面位置的确定，利用枚举法，总费用最少，求得在什么位置，什么海拔高度平整，已达到最优的标准。 针对问题三：类比于问题二，采用相同的方法，在不同的海拔高度上平整，已满足二层的台阶状地块的要求，达到挖填土石方量最小的目标。 最后，对所建立的模型和求解方法的优缺点给出了客观的评价，并指出了改进的方法。 关键词：土石方量最小；二重积分；枚举法； 一、问题重述与分析 1.1问题的重述 十堰市是一个山区城市，向山要地是十堰市发展的一个必然的选择，但是如何在一片山地之中选择合适的方位与开挖深度，从而使总的土石方量最小，就是一个有意义的命题了！ 某工厂为了在一片长度为1500米，宽度为1000米的山地之中，开挖出一个800米×600米平坦连续的长方形地块作为工厂的厂房地基，前期已经在本块土地上测量出长、宽每隔30米的网格的对应网格点的海拔高度。 问题： ①用附件（1）中的数据画出工厂的这片土地的三维图形与等高线图； ②从什么地方，什么海拔高度平整一块800米×600米的连片土地能使总的土石方量最小？ ③如果允许平整出来的土地为二层的台阶状地块，要求各地块的长、宽不少于60米，又将从什么地方、什么海拔高度分别开挖，能使总的土石方量最小？ 1.2问题的分析 1.2.1问题一 我们可以利用MATLAB软件直接调用surf函数和contour函数分别画出这片土地的三维图形和等高线图。 1.2.2问题二 考虑要平整一块800米×600米的连片土地，且要使土石方量最小，我们需要确定这块800米×600米的连片土地的底面投影区域，以及海拔高度，使得其土石方量最小。让底面投影区域在长度为1500米，宽度为900米内枚举，此时计算出对应的土石方量最小的是，在这些所有的最小值中取得最小值的那块投影区域即为所求。我们采用此种方法求解。在平整土地的过程中，凸出的地方是要挖土的，凹下的地方是可以填土的，这样挖下来的土可以填充到凹下去的部分。通过计算体积的方式计算挖土和填土大小，求出最小值对应的海拔高度，即可确定整连片土地的具体高度。 二、基本假设 1、在平整土地的过程中，有些地方是要挖山的，但有些地方是要填土的，假设填土的每立方米所需的费用为挖山的每立方米土石方所需费用的1/3。 2、假设除了挖土和填土以外，在平整土地的过程中其他作业（如登高）产生的费用都与800米×600米的连片土地所处的位置方向和海拔高度均无关。这样我们只需考虑挖土和填土的土石方量及费用，据此来考虑连片土地所处的位置方向和海拔高度。 3、假设计算体积的过程中，分割的小曲顶柱体不能达到无穷小，而产生的误差忽略不计。实际计算中我们取一个很小的步长去划分，使其划分尽可能的小，产生的误差忽略不计。 三、符号约定 符号 说明 平整连片土地块的长 平整连片土地块的宽 平整连片土地块的底面所在闭区域 平整连片土地块底面划分的子区域 子区域的面积 平整土地所需的挖土量 平整土地所需的填土量 平整土地挖土和填土所需要的总费用 四、模型的建立 问题一求解 利用MATLAB7.0软件的三维绘图功能，画出工厂这片土地的三维图形如图4.1.1所示： 图4.1.1 山地的三维图形 这片土地的三维等高线图如图4.1.2 所示，二维等高线图4.1.3所示： 图4.1.2 山地的三维等高线图 图4.1.3 山地的二维等高线图 问题二模型建立 平整块的海拔高度的确定 通过土石方量费用最小的原则，确定平整土地海拔的开挖高度。平整土地快所对应的的底面区域记为，显然是一个800米×600米的矩形区域。设顶部海拔高度是一个非负连续函数，。假设 的位置已经确定，下面我们来讨论在什么样的海拔高度下，土石方量费用最小。 考虑底面区域对应的山体是一个曲顶柱体，把区域分成个小区域，，…，，为了简化运算，我们分成个等大的区域，每个区域都是边长为的正方形。再分别以这些小区域的边界为准线，作母线平行于轴的柱面，这些柱面将原曲顶柱体分成个细的曲顶柱体。每个区域上