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

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

下地幔： α = 13.7 km/s, β = 7.2 km/s, ρ = 5.5 g/cm³

地核： α₂ = 8.0 km/s, β₂ = 0.0 km/s, ρ₂ = 9.9 g/cm³

一、实验原理

根据透射系数，反射系数，位移比和能流比公式计算反射系数，透射系数和位移比，和能流比。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