地震矩张量计算

C-3. (a)用 C-1的结果编写程序计算地震矩张量

(b) 用该程序计算下列震源机制的地震矩张量

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

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

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

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

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

1. 实验原理

矩张量可以根据滑动矢量和正交矢量求得，公式为

二、程序源代码

1.主函数

a1=input(‘inputa1’);

a2=input(‘inputa2’);

a3=input(‘inputa3’);

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

M=M0*[2*n(1)*d(1),n(1)*d(2)+n(2)*d(1),n(1)*d(3)+n(3)*d(1);

n(2)*d(1)+n(1)*d(2),2*n(2)*d(2),n(2)*d(3)+n(3)*d(2);

n(3)*d(1)+n(1)*d(3),n(3)*d(2)+n(2)*d(3),2*n(3)*d(3)];

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°

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

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

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

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

结果分析：当走向和倾向相同，滑动角相差180度时，两者的矩张量相加等于零。