山野莽夫

  • 归档
    • 随笔
    • 建站资源
    • 分享
    • 代码
  • 地球物理学
    • 专业课
    • 概念解释
  • 计算机
  • 互联网
  • 教程
  • 规划
  • 实验室
    • 珍藏的软件
    • 贴吧云签到
    • A1账号自助申请
山野莽夫
小学生的挣扎的点点滴滴
  1. 首页
  2. 地球物理学
  3. 专业课
  4. 地震学
  5. 正文

实验一 基本地震学理论(斯奈尔定律)

2016年1月25日 3566点热度 0人点赞 0条评论

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('请输入入射角 a:');

vp1=input('上层p波速度 vp1:');

vp2=input('上层S波速度 vp2:');%如果作为子函数使用最好把这一块输入去掉,把输入调整到主函数中%

vs1=input('下层P波速度 vs1:');

vs2=input('下层S博速度 vs2:');%..........................................................%

if i==1 %P波

a1=a;

a2=asind(sind(a)*vs1/vp1);

b1=asind(sind(a)*vp2/vp1);

b2=asind(sind(a)*vs2/vp1);

fprintf('%f\t%f\t%f\t%f\n',a1,a2,b1,b2);

if vp2>vp1

r1=asind(vp1/vp2);

fprintf('p波临界角为:%f\n',r1);

else fprintf('p波没有临界角\n');

end

if vs2>vp1

r2=asind(vp1/vs2);

fprintf('sv波临界角为:%f\n',r2);

else fprintf('sv波没有临界角\n');

end

end

if i==2; %S波

a2=asind(sind(a)*vp1/vs1);

a1=asind(sind(a)*vs1/vs1);

b2=asind(sind(a)*vp2/vs1);

b1=asind(sind(a)*vs2/vs1);

fprintf('%f\t%f\t%f\t%f\n',a1,a2,b1,b2);

if vp2>vs1

r2=asind(vs1/vp2);

fprintf('p波临界角为:%f\n',r2);

else fprintf('p波没有临界角\n');

end

if vs2>vs1

r1=asind(vs1/vs2);

fprintf('sv波临界角为:%f\n',r1);

else fprintf('sv波没有临界角\n');

end

end

if i==3 %SH波

a2=asind(sind(a)*vs1/vs1);

b2=asind(sind(a)*vs2/vs1);

fprintf('%f\n%f\n',a2,b2);

end

end

三、 结果及分析

这个小问题并没有什么结果,特采用P波以30度角测试一下,选取模型,VP1=6.8km/s,VP2=8.0 km/s,VS1=3.9 km/s,VS2=4.8km/s。

30.000000 16.664362 36.031879 20.667316

p波临界角为:58.211669

sv波没有临界角

标签: 代码 地震学实验 地震学程序 基本理论 斯奈尔定律
最后更新:2016年1月25日

小菜菜

菜鸟

打赏 点赞
< 上一篇
下一篇 >

文章评论

razz evil exclaim smile redface biggrin eek confused idea lol mad twisted rolleyes wink cool arrow neutral cry mrgreen drooling persevering
取消回复

此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据。

联系方式

QQ群 | TG群 | 邮箱

最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
认识测井软件--CifLog Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享 实验三、化向地磁极 申请google voice 方法-纪念我撸到的第一个gv号 wordpress大学wpdx1.3主题分享 数据结构(C语言版)第2版参考答案
标签聚合
地震学程序 虚拟机 c语言 wordpress ppt 宝塔面板 模板 onedrive
最近评论
小菜菜 发布于 6 个月前(11月24日) 这玩意已经废了,成收割工具了,不能再用了。
eamon 发布于 7 个月前(11月07日) 我一年不用了才发现这个休眠管理费每月15,一共扣了我135元,然后我消费还消费不了,我宁愿消费掉也不...
magic 发布于 11 个月前(07月03日) 请问账号不注销会有什么影响吗?
magic 发布于 11 个月前(07月01日) 我想问一下 如果不注销账号就留着会怎么样
qwp6601 发布于 12 个月前(06月04日) 有没有方法改为bing

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2