水道测量的设计与解法.docVIP

海床轮廓的构建及航道的确定论文 37组 杨艳林 刘汇川 周旭斌 整个问题，简言之，就是利用题目所给的14个点求出海床深度分布，再标记出深度小于5码的区域以禁止通行。针对这个问题，论文提出了三个基本假设： 所给数据是准确的； 讨论区与海地区面是光滑的，即曲面有一阶、二阶导数且都是连续的； 某位置的水深与其周围的水深相互作用，且这种作用随着距离的增大而减小。 根据假设，对于海底曲面光滑这一特性，可以用某个光滑曲面去逼近所给的14个节点。进而产生了两种想法：（1）、利用14个已给出的点，由假设三给所有网格计算出一个水深值，用网格及对应的水深值画出曲面；（2）、利用曲面光滑特性，直接由已给出的水深值直接求出所有格点的水深值。最后只需要对拟合的曲面求出深度小于5英尺的区域即达到了要求。针对以上两种不同的想法，论文提出了两种不同的方法： 反距离平方格式与双三次样条拟合 首先对确定X、Y坐标轴的范围，由题给的14个数据确定范围是X：（75,200）；Y：（-100,150）。在这一范围，对14个数据点的X、Y坐标分别进行分割。在这一区域内形成了一个不等距的16个节点的分割。 由假设三知道，给出的14个数据点对所有的格点都有影响，而且距离近的影响较大，距离远的影响较小。基于以上实事，论文采用了加权平均格式，用带距离的平方的逆权来表示对近距离的已知点影响更大，于是得出了最初的深度算法：。但是一个加权平均并不能完全反映曲面的走势。例如，再简单的一维情况下，当互相间隔4米的三个点A、B、C，若A点深8米，B点深4米，那么按照上面的算法就应该得到C点深度为4.8米，而这个4.8米显然与我们日常的认识：C点深度应该小于4米，是相悖的。从这个例子可以看出远点A对C点的影响大于B点而致使C点计算得到的深度大于4米。所以需要在加入一项，使得近点的影响更大，于是引入趋向项以增强近点对于所要求的点的影响。然而在二维平面上，G点并不在Q1、Q2 的连线上，于是由G点向Q1、Q2 的连线作垂线，设垂足为P，则Q1、Q2 对G的影响可以在P上体现出来。于是，再次引入关于P点的反距离加权影响，最终得表达式：。 根据以上表达式可推测出： 每个已知点对未知点的影响表现在两者间距离的平方分之一； 任意一已知点对于位置点都有影响，即，未知点的水深随着任一已知点的变化而变化； 任意两点对于未知点的影响取决于其中的三个距离； 若要求一个已知点A对某一未知点的影响，那必须把除这个点之外的所有已知点与A两两组合求出对未知点的总影响。 （不足之处：我们认为Z的左边与等式右边的式子一定是成正比的关系，但不一定是相等，即比例系数=1。应当存在比例系数K，最后应当利用已知点，代入公式求出K。） 以上所叙述的影响由以下几点原因产生： 两个已知点与未知点的距离； 未知点与两已知点连线的垂直距离； 两已知点之间的距离。 有以上的出的Z的坐标结合之前得到的X、Y坐标，经过双三次样条插值法便可拟合出曲面。 （2）极小曲率差分迭代法 在这一方法中，论文直接对区域作细密分割，只要求出每个节点上的水深值，就能直接描绘海底曲面。以下皆为求出海水深度的讨论： 将分割好的区域中不在网格上的点作微小平移，论文认为，这样做不会影响最后结果。将处理过的坐标代入微分方程： 。 这时有边界条件：1、分割矩形区域四条边上作用力为零；2、切线零混合矩；3、u(xn,yn)=wn。 离散化海底曲面总平方曲率得：。当涉及二维时，可以用简单曲率逼近表达式：。最后让uij出的导数为零，再代入调整后的数据，即可解出uij，最后将坐标画在坐标纸上便可得到海底大致曲面。 （3）IDW解法 将X、Y的区域作等距离分割，由假设三，我们可以设每个点对于未知点的影响表现为距离的p次方，则 ； 其中，，Di即为调整后的格点与未知点间的距离。 在上式中，指数p反映了未知点对于已知点的离散情况，p越大，说明两者离得越近，影响越大。在得出了Z的坐标之后，再通过双三次样条插值法便可得出所要求图形。 程序与图形： 反距离平方格式与双三次样条拟合法程序： #includestdio.h #includemath.h float dis(float x1,float y1,float x2,float y2) { float a=0; a=(float)sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1)); return(a); } float XP(float x1,float x2,float x3,float y1,float y2,float y3) { float xu, xd, x; xu=x3*(x1-x2)*(x1-x2)+(y1-y2)*(x2*(y1-y3