实验一 基本地震学理论(固固界面纵波和横波能流比反射系数透射系数数)

2016年1月25日 343 次阅读 0 条评论 0 人点赞

C-3. 编写程序:已知固-固界面两边介质速度和密度,求P波和S波垂直入射位移反射和透射系数、能流比 。用该程序计算核幔边界的这些量:

下地幔: α = 13.7 km/s, β = 7.2 km/s, ρ = 5.5 g/cm3

地核: α2 = 8.0 km/s, β2 = 0.0 km/s, ρ2 = 9.9 g/cm3

一、实验原理

根据透射系数,反射系数,位移比和能流比公式计算反射系数,透射系数和位移比,和能流比。P波公式为

同理S波公式和上述公式形式相同。

二、 程序源代码

v1=input('上层速度');

v2=input('下层速度');

p1=input('上层密度');

p2=input('下层密度');

j1=input('入射角');

p=sind(j1)/v1;%射线参数

j2=asin(p*v2);%利用斯奈尔定律计算透射角

T12=2*p1*v1/(p1*v1+p2*v2);%透射系数

R12=(p1*v1-p2*v2)/(p1*v1+p2*v2);%反射系数

Er=R12*R12;%反射能流比

Et=T12*T12*v2*p2*cosd(j2)/(v1*p1*cosd(j1));%透射能流比

fprintf('j2=%f,T12=%f,R12=%f,Eri=%f,Eti=%f',j2,T12,R12,Er,Et);

三、 结果及分析

P波: 反射系数-0.024911, 透射系数0.975089, 反射能流比0.000621,

透射能流比0.999379, 反射波位移比-0.024911, 透射波位移比0.975089

S波:反射系数1.000000, 透射系数2.000000, 反射能流比1.000000, 透射能流比0.000000, 反射波位移比1.00000, 透射波位移比2.00000

菜鸟

文章评论(0)