地震震源监测的数学建模分析.docVIP

《》使用数学非线性函数建模，模型中的参数非线性进行最佳逼近拟合 2008年5月12日发生在四川汶川的一场大地震给我们每一位中国人带来了巨大的伤痛，痛定思痛，我们应该为减少震后灾害做些事情。 当地震发生时，震中位置的快速确定对第一时间展开抗震救灾起到非常重要的作用，而震中位置可以通过多个地震观测站点接收到地震波的时间推算得到。这里假定地面是一个平面，在这个平面上建立坐标系见图1。图1中给出了10个地震观测站点（A—J）的坐标位置。 图1 地震观测站点示意图 ? 2009年4月1日某时在某一地点发生了一次地震，图1中10个地震观测站点均接收到了地震波，观测数据见表1。 地震灾区的地形、地貌对抗震救灾的进展会有很大影响，根据卫星遥感和飞机航拍得到的照片可以构建灾区地形、地貌图。所构建的灾区局部地形、地貌等高线图见图2。 注：（1）图2中等高线标注数字的单位为千米； （2）图1和图2不属于同一坐标系。 图2 灾区局部地形、地貌等高线图 ? 问题一：假定地震波在各种介质和各个方向的传播速度均相等，并且在传播过程中保持不变。请你根据表1中的数据确定这次地震的震中位置、震源深度以及地震发生的时间（不考虑时区因素，建议时间以分为单位）。 问题二：请你根据图2中的局部地形、地貌等高线图，建立数学模型，绘制出相应的三维地形、地貌曲面图。 问题三：查阅有关资料，了解地震波在各种介质和各个方向的传播速度问题，给出合理假设，根据表1中的数据确定这次地震的震中位置、震源深度以及地震发生的时间 表1 地震观测站坐标及接收地震波时间 地震观测站 横坐标x（千米） 纵坐标y（千米） 接收地震波时间 A 500 3300 4月1日9时21分9秒 B 300 200 4月1日9时19分29秒 C 800 1600 4月1日9时14分51秒 D 1400 2200 4月1日9时13分17秒 E 1700 700 4月1日9时11分46秒 F 2300 2800 4月1日9时14分47秒 G 2500 1900 4月1日9时10分14秒 H 2900 900 4月1日9时11分46秒 I 3200 3100 4月1日9时17分57秒 J 3400 100 4月1日9时16分49秒 概述： —① —② —③ —④ —⑤ —⑥ —⑦ — ⑧ —⑨ —⑩ 利用MATLAB编写程序如下： [h,yfun,key,c]=fsolve(earthqu,[1600;2000;100;1000;20])%非线性方程的求解主函数function q=earthqu(h)%非线性方程函数 x1=500;x2=300;x3=800;x4=1400;x5=1700;x6=2300;x7=2500;x8=2900;x9=3200;x10=3900;%十个点的x坐标 y1=3300;y2=200;y3=1600;y4=2200;y5=700;y6=2800;y7=1900;y8=900;y9=3100;y10=100;%十个点的y坐标 t1=1269;t2=1169;t3=891;t4=797;t5=706;t6=887;t7=614;t8=706;t9=1077;t10=1009;%以9：00为基准点的时间的变化量 q=zeros(10,1); %H(1)=x H(2)=y p(3)=z H(4)=v H(5)=t %(x-xi).^2+(y-yi).^2+z.^2-(ti-t).^2.*v.*v q=[(h(1)-x1).^2+(h(2)-y1).^2+h(3).^2-(t1-h(5)).^2.*h(4).*h(4);%A点非线性方程 (h(1)-x2).^2+(h(2)-y2).^2+h(3).^2-(t2-h(5)).^2.*h(4).*h(4);%B点非线性方程 (h(1)-x3).^2+(h(2)-y3).^2+h(3).^2-(t3-h(5)).^2.*h(4).*h(4);%C点非线性方程 (h(1)-x4).^2+(h(2)-y4).^2+h(3).^2-(t4-h(5)).^2.*h(4).*h(4);%D点非线性方程 (h(1)-x5).^2+(h(2)-y5).^2+h(3).^2-(t5-h(5)).^2.*h(4).*h(4);%E点非线性方程 (h(1)-x6).^2+(h(2)-y6).^2+h(3).^2-(t6-h(5)).^2.*h(4).*h(4);%F点非线性方程 (h(1)-x7).^2+(h(2)-y7).^2+h(3).^2-(t7-h(5)).^2.*h(4).*h(4);%G点非线性方程 (h(1)-x8).^2+(h(2)-y8).^2+h(