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