地震法向矢量和滑动向量计算

2016年2月12日 362 次阅读 0 条评论 0 人点赞

本次试验使用matlab,如下为我的实验报告

 C-1.(a) 编写子函数:给定三个断层角度,计算断层的法向矢量和滑动向量

(b) 用该子函数计算震源机制如

φ = 280°, δ = 60°, λ = 270°;

φ = 40°, δ = 80°, λ = 20°

的断层的法向矢量和滑动向量

(c) 检验断层的法向矢量和滑动向量正交

一 、实验原理

断层的法向矢量和滑动矢量可以用走向倾向和滑动角求出,公式为:


二、程序源代码

function [ n,d ] = dz4111(a1,a2,a3)

a1=input('inputa1');

a2=input('inputa2');

a3=input('inputa3');

n=[-sind(a2)*sind(a1);-sind(a2)*cosd(a1);cosd(a2)];%法向矢量

d=[cosd(a3)*cosd(a1)+sind(a3)*cosd(a2)*sind(a1);

-cosd(a3)*sind(a1)+sind(a3)*cosd(a2)*cosd(a1);

sind(a3)*sind(a2)];%滑动矢量

s=dot(n,d);

n

d

if dot(n,d)==0

fprintf('是正交');%判断正交

end

三、 结果及分析

1.φ = 280°, δ = 60°, λ = 270°

法向矢量n = [0.8529, -0.1504, 0.5000];

滑动矢量d = [0.4924, -0.0868, -0.8660];是正交

2.φ = 40°, δ = 80°, λ = 20°

法向矢量n =[-0.6330,-0.7544,0.1736]

滑动矢量d =[0.7580,-0.5585,0.3368];不是正交

菜鸟

文章评论(0)