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

本次试验使用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]；不是正交