空间探测模型.pptVIP

空间探测模型 问题： 山体、隧洞、坝体等的某些内部结构可用弹性波测量确定。一个简化问题可描述为，一块均匀介质构成的矩形平板内有一些充满空气的空洞，在平板的两个邻边分别等距的设置若干波源，在它们的对边对等地安放同样多的接 收器，记录弹性波由每个波 源到达对边每个接收器的时 间，根据弹性波在介质中和 在空气中不同的传播速度， 来确定板内空间的位置。 现考虑如下具体问题： 一块240m×240m的平板(如图),在AB边等距的设置7个波源Pi (i=1,…,7),CD边对等地安放7个接收器Qj (j=1,…,7),记录由Pi发出的弹性波到达Qj的时间tij(s);在AD边等距的设置7个波源Ri (i=1,…,7),BC边对等地安放7个接收器Sj (j=1,…,7),记录由Ri 发出的 弹性波到达Sj的时间kij(s).已 知弹性波在介质和空气中的 传播速度分别为2880(m/s)和 320(m/s),且弹性波在铁板边 缘的传播速度与在介质中传 播速度相同。由此确定平板 内空洞的位置 分析： 模型建立： 求解： 求解使 最小的x的最优取值 模型的优点和缺点 附：MATLAB编程得到方程组 * * A B C D A B C D Pi Qj Ri Sj 0.1042 0.0688 0.0760 0.1954 0.3566 0.3397 0.4311 P7 0.1031 0.0939 0.2217 0.3064 0.2364 0.3177 0.3807 P6 0.1810 0.1768 0.0839 0.1917 0.2263 0.4529 0.3490 P5 0.2122 0.0740 0.1789 0.0738 0.4040 0.4453 0.3221 P4 0.1919 0.3563 0.4156 0.4153 0.0598 0.4131 0.3052 P3 0.3805 0.5242 0.4470 0.4318 0.4413 0.0592 0.0989 P2 0.5646 0.4923 0.4181 0.2032 0.1996 0.0895 0.0611 P1 Q7 Q6 Q5 Q4 Q3 Q2 Q1 tij 部分数据 A B C D Pi Qj Ri Sj A B C D Pi Qj Ri Sj 空洞位置 传播时间 正问题 反问题 空洞位置如何用数学方式表述？ 假设1 弹性波在空洞内传播速度一致，在其它介质中速度也一致 假设2 图中的36块小区域中，要么全是空气，要么全是介质 A B C D R3 2 1 4 3 5 7 6 9 8 11 10 15 12 33 27 21 S2 d1 d2 d3 d4 d5 d6 第i个小方块内全是介质 第i个小方块内全是空气 波线R3S2在介质中走的距离 在空气中走的距离 整理得 同理，利用计算机编程可得一个含36个未知数，84个方程的方程组 A B C D R3 2 1 4 3 5 7 6 9 8 11 10 15 12 33 27 21 S2 d1 d4 d2 d3 d5 d6 A B C D R3 S2 d1 d2 d3 d4 d5 d6 另： 设最终得到的方程组为 （A，b）极有可能列满秩，即无解 不符合实际数据的假设 二次规划 利用MATLAB求解可得 0.033 0.048 0.051 0.046 0.002 0.084 0.004 0.064 0.065 0.005 0.914 0.047 0.059 0.033 0.997 0.982 0.009 0.113 0.037 0.916 0.057 0.959 0.951 0.058 0.018 0.001 0.002 0.926 0.985 0.006 0.001 0.038 0.025 0.073 0.013 0.021 结论：表中红线区域内出现大面积空洞的可能性很大 进一步的问题 （2）增加细分程度，能否得到更精确的结果？ A B C D 可以得到含144个未知数，84个方程的方程组 （1）在保证探测结果的前提下，能否适当减少波源和接收器个数？ 0 0.002 0.377 0 0 0 0 0.007 0 0 0 0 0 0 0.014 0 0 0 0 0 0 0.227 0 0 0 0 0.059 0 0 0 1 0.626 0.998 0.688 0 0 0 0 0 0 0 0.036 0.818 0.998 0.798 0.940 0 0.109 0.057 0.042 0.682 1 0 0.071 0.753 1 0.861 0.818 0 0 0 0 0.954 0.675 0 0.00