空洞测量的几种解法.doc

空洞定位的几种解法 向为 王瑛 伍微 国防科技大学 四院一队 摘要：本文提出了图解排除法，及基于线性方程组的空洞定位方法（主要方法）。对第二种方案，先网格化得到超定的线性方程组，后用最小二乘法解之。在试图减少波源和接收器时，我们确定了若干基本原则，为机器判定作了准备。 问题描述 一块２４０×２４０的平板（如图），在AB边等距地放置７个波源Pi（i=1…7）,CD边对等地安放７个接收器Qj（j=1…7）,记录由Pi发出的弹性波达到Qj的时间t（秒）；在AD边等距地放置７个波源Ri（i=1…7）,BC边对等地安放７个接收器Sj（j=1…7）,记录由Ri发出的弹性波达到Sj的时间τ（秒）（如图1）。已知弹性波在介质和空气中的传播速度分别为２８８０（米/秒）和３２０（米/秒），且弹性波沿板边缘的传播速度于在介质中的传播速度相同。 确定该平板内空洞的位置。 仅根据由Pi发出的弹性波达到Qj的时间t（i=1…7），能确定空洞的位置吗？讨论在同样能够确定空洞位置的前提下，减少波源和接收器的方法。 图1 二．基本假设 1．假设所测时间数据真实无错误，但是有一定的误差。 2．假设平板可划分化为网格，空洞定位于每个网格单元内。 三．变量描述 WaveLine: 第条波线，即以PiQj为第=(i-1)*7+j条波线，以RiSj为第=(i-1)*7+j+49 条波线；（=1…98） Time: 测得的弹性波沿WaveLine的传播时间； Length: WaveLine的长度； Hole: 沿WaveLine的空洞的长度； MeshgridUnit: 第k个网格单元，具体表示如图2； (k=1…36) 图2 X: MeshgridUnit的空洞权值，X=0时，表示MeshgridUnit内无空洞，X=1时，表示MeshgridUnit内有空洞； WaveLineInUnit: WaveLine被MeshgridUnit所截得的长度，即WaveLine在 MeshgridUnit内的那部分长度。 四．思路及解法 数据初步处理： 由弹性波的物理传播方式可得： 　　　　 （=1…98） 运行Ｃ程序GetHoleLength.cpp(见附录)可求得Hole，矩阵Hole见附录。 由于真实的Hole≥0，分析数据得出误差限为9.0120米。 若算出空洞长度小于误差限的可认为没有空洞。 方案一：图解排除法 　　若Hole≤9.0120m, 则用直线连接与对应的PiQj或RiSj，得图3如下： 图3 由基本假设2.知：直线通过的网格单元内不存在空洞， 无直线通过的网格单元内（图3中阴影区）可能存在空洞。 由于R5S6通过一个阴影网格单元，故MeshgridUnit内有空洞； 由于P5Q6通过一个阴影网格单元，故MeshgridUnit内有空洞； 由于P4Q5通过一个阴影网格单元，故MeshgridUnit内有空洞； 由于P1Q4通过一个阴影网格单元，故MeshgridUnit内有空洞； 由于P3Q7通过一个阴影网格单元，故MeshgridUnit内有空洞； 由于P1Q7间空洞长度为160.8296米，且通过三个阴影网格单元，故MeshgridUnit, MeshgridUnit内有空洞； 由于R4S5间空洞长度为86.5502米，且通过二个阴影网格单元，故MeshgridUnit内有空洞。 根据已求出的Hole可知，洞径均在40米左右，可认为以上8个方格内几乎充满空洞。 注：本方案C源程序见附录。 方案二：基于线性方程组的空洞定位法 准备阶段： 从经典的线性方程组定位方法（如CT图像重建等）出发，本题的波线似乎应该有一定的线宽。但是，考虑到实际的波线线宽远小于山体等自然物，同时，弹性波的传播速度不会因为线宽的改变而改变。所以，弹性波的波线线宽可忽略不计。（参考书目[1]） 由于波线线宽忽略不计，实际中应该这样处理：左右网格单元邻边属于右侧网格单元，上下网格单元的邻边属于上侧网格单元。 运行Ｃ程序LineInHole.c可求得WaveLineInUnit。 前面已经求出了Hole。 算法实施： 记 （空洞权值矩阵） 就本题而言，m=98, n=36。 因为对WaveLine空洞长度守恒，有： ……………………………………………（＊） 这是一个超定方程组，理论上可用最小二乘法解决，即求 的 实际中，我们使用MATLAB软件的mldivide（\）和lsqnonneg两个函数分别求解，得： （注：A的秩