用PREM模型追踪震源地表的地震射线，表示出上地幔间断面和地核对射线的影响。

海洋地壳模型反射波时距曲线绘制以及两种公式汇出的时距曲线的对比。

C-4. 编写程序：给定介质两边速度，用C-2的结果生成如图2.6-11 所示的射线；用该程序计算核幔边界可能入射波的射线路径（问题C-3） 下地幔： α = 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波不同入射角下的反射波和透射波的角度，根据波射线的角度绘制一个入射点所有波射线的图像。根据不同入射角划分x轴为多个点，分别绘…

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('上层速度…

C-2. 编写子函数：对于固-固介质界面，已知P波和S波速度以及入射角，求出P波和S波反射及透射角。（注意：该子函数可以计算并列出入射波的临界角） 一、 实验原理 二、程序源代码 function [ a1,a2,b1,b2] = snel( i,a,vp1,vp2,vs1,vs2 )%本程序一些地方有中文，在某些版本的matlab中可能无法正确显示，但不影响使用。 i=input('请输入波形，p波请输入1，S波请输入2，SH波请输入3:');%........................% a=input(…

C-1. 编写cos (ω t − kx)函数，并画出以下结果： (a) 时间 t 从 0 到 10, 在 x = 1, 当 ω = 1, k = 1. (b) 时间 t 从 0 到 10, 在 x = 0, 当 ω = 4, k = 1. (c) 位置 x从 0 到 10, 在 t = 0, 当 ω = 1,k = 2. (d) 位置x从 0 到 10, 在 t = 0, 当 ω = 1,k = 4. 一、 实验原理 运用matlab编写程序，根据函数解析式 cos (ω t − kx) 分别改变函数角速度等参数…