三分量地震波分析及matlab代码

三分量地震波分析是指对地震信号的三个方向（东西向、南北向、垂直向）进行分析，以获取更为准确的地震信息。这种分析方法可以用于地震勘探、地震监测等领域。

以下是三分量地震波分析的matlab代码：

1. 读取数据

filename = 'data.txt';
data = load(filename);

2. 数据预处理

dt = 0.01; % 采样时间间隔
t = 0:dt:(size(data,1)-1)*dt; % 时间序列
data1 = detrend(data(:,1)); % 去趋势
data2 = detrend(data(:,2));
data3 = detrend(data(:,3));

3. 滤波

fc = 5; % 截止频率
fs = 1/dt; % 采样频率
[b,a] = butter(4,fc/(fs/2),'low'); % 低通滤波器
data1f = filtfilt(b,a,data1);
data2f = filtfilt(b,a,data2);
data3f = filtfilt(b,a,data3);

4. 计算速度

v = 4.0; % 声速
dx = 1.0; % 接收器间距
t0 = 0.0; % 起始时间
x = 0:dx:(size(data,2)-1)*dx; % 接收器位置
t1 = sqrt((x(1)-x).^2+v^2*t0^2)/v; % 到达时间
t2 = sqrt((x(2)-x).^2+v^2*t0^2)/v;
t3 = sqrt((x(3)-x).^2+v^2*t0^2)/v;

5. 叠加

dataf = data1f+data2f+data3f;

6. 可视化

figure;
subplot(2,1,1);
plot(t,data1f,'r',t,data2f,'g',t,data3f,'b');
xlabel('Time(s)');
ylabel('Amplitude');
legend('E-W','N-S','Vertical');
subplot(2,1,2);
plot(x,dataf);
xlabel('Distance(m)');
ylabel('Amplitude');