地震断层P轴和T轴计算

本次试验使用matlab

C-2. (a) 编写程序：用C-1的结果计算P轴和T轴的方向.

(b) 用该程序计算以下问题的P轴和T轴的方向

(1) φ = 330°, δ = 65°, λ = 70°

(2) φ = 280°, δ = 60°, λ = 270°

(3) φ = 280°, δ = 60°, λ = 90°

(4) φ = 40°, δ = 80°, λ = 20°

(5) φ = 40°, δ = 80°, λ = 200°

1. 实验原理
   

   根据法向矢量和滑动矢量可以计算P轴和T轴的方向，公式如下
   

   
   故可以利用上一问程序来计算滑动矢量和法向矢量，
   

   进而计算P轴和T轴的方向。
   

二、程序源代码

1.主函数

a1=input('inputa1');

a2=input('inputa2');

a3=input('inputa3');

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

t=n+d;

p=n-d;

b=cross(n,d);

t

p

b

2.子函数

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

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)];

三 、结果及分析

(1) φ = 330°, δ = 65°, λ = 70°

T轴方向[0.5508,-0.2699,1.2743]

P轴方向[0.3555,-1.2998,-0.4290]

b = [-0.8861,-0.3447,0.3100]

(2) φ = 280°, δ = 60°, λ = 270°

T轴方向[1.3453,-0.2372,-0.3660]

P轴方向[0.3605,-0.0636,1.3660]

b = [0.1736,0.9848 ,0]

(3) φ = 280°, δ = 60°, λ = 90°

T轴方向[0.3605,-0.0636,1.3660]

P轴方向[1.3453,-0.2372, -0.3660]

b =[-0.1736,-0.9848, 0]

(4) φ = 40°, δ = 80°, λ = 20°

T轴方向[0.1250,-1.3129,0.5105]

P轴方向[-1.3910,-0.1959,-0.1632]

b = [-0.1571,0.3448, 0.9254]

(5) φ = 40°, δ = 80°, λ = 200°

T轴方向[-1.3910,-0.1959,-0.1632]

P轴方向[0.1250,-1.3129, 0.5105]

b =[ 0.1571, -0.3448, -0.9254]