没有合适的资源?快使用搜索试试~ 我知道了~
地震学实验报告二(matlab)
需积分: 44 12 下载量 79 浏览量
2020-12-17
10:55:00
上传
评论 1
收藏 336KB DOCX 举报
温馨提示
试读
13页
编写程序计算入射P波和S波在固-固界面上垂直入射的位移反射和透射系数,能量比值。利用程序估算核-幔边界(尽管是固-液界面)处的这些物理量,假定下地幔p波速度为13.7km/s,S波速度为7.2km/s,密度为5.5g/cm³,地核的P波速度为8.0km/s,S波速度为0.0km/s,密度为 9.9g/cm³。
资源详情
资源评论
资源推荐
实验二
C-1 编写程序计算入射 P 波和 S 波在固-固界面上垂直入射的位移反射和透
射系数,能量比值。利用程序估算核-幔边界(尽管是固-液界面)处的这些
物理量,假定下地幔 p 波速度为 13.7km/s,S 波速度为 7.2km/s,密度
为 5.5g/cm³,地核的 P 波速度为 8.0km/s,S 波速度为 0.0km/s,密度
为 9.9g/cm³。
一. 实验原理:
根据公式计算。
R=
ρ
1
v
1
− ρ
2
v
2
ρ
1
v
1
+ ρ
2
v
2
T =
2 ρ
1
v
1
ρ
1
v
1
+ ρ
2
v
2
E
R
E
=R
2
二. 程序源代码:
function [R,T,QR,QT]=yu(p1,p2,v1,v2)
R=(p1*v1-p2*v2)/(p1*v1+p2*v2);
if v2==0;
T = 0;
else
T=2*p1*v1/(p1*v1+p2*v2);
end
QR=R*R;
QT=(T*T*p2*v2)/p1*v1;
end
clc;
close all;
%参数
vp1=13.7;
vp2=8.0;
p1=5.5;
vs1=7.2;
vs2=0;
p2=9.9;
%计算
[Rp,Tp,QRp,QTp]=yu(p1,p2,vp1,vp2);
[Rs,Ts,QRs,QTs]=yu(p1,p2,vs1,vs2)
三. 结果与分析:
P 波 S 波
C-2
一、实验原理
运用斯奈尔定律求出特定地层模型下,P 波或 S 波不同入射角下的反射波和透射波
的角度,根据波射线的角度绘制一个入射点所有波射线的图像。根据不同入射角划
分 x 轴为多个点,分别绘制不同入射角的波产生的所有射线。
二、程序源代码
i=input('请输入第几幅图 i');%图号
vp1=input('input vp1:');%上层 p 波速度
vp2=input('input vp2:');%下层 p 波速度
vs1=input('input vs1:');%上层 s 波速度
vs2=input('input vs2:');%下层 s 波速度
if i==1;%第一幅图 p 波从上方入射
x=[0 0];
y=[2 -1];
plot(x,y,'k')
hold on
for a=10:10:90
[ a1,a2,b1,b2 ]=snel( i,a,vp1,vp2,vs1,vs2 );
剩余12页未读,继续阅读
weixin_53668688
- 粉丝: 0
- 资源: 4
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功
评论0