这里主要有四个参数需要设置
wigb.m地震数据画图函数下载地址
a: seismic data 地震数据
% scale: multiple data by scale 扫描频率【1-3】范围内任选一个数字
% x: x-axis; 所有道的位置信息
% z: vertical axis (time or depth) 时间
% amx:这个参数可以不用管它
这里有一个时序信号;数据一行为所有道的采集数据,列为道数。
%加载地震数据的,也可以使用fopen()函数加载数据
load earthquak_data.mat %加载地震数据
rec=earthquak_data(1:1:1000,13:36);
[M,N]=size(rec);%地震数据是按照时间序列存放的,比如第一次采样时,所有道的数据是一行,列为道数
dt=0.002; %这里为采样时间间隔
t=0:dt:(M-1)*dt; %这里需要计算正确数据的时间
%下面这句是用来设置每道的道间距
for i=1:N
often_offset(i)=2*(i-1);
end
figure
xlabel('Traces');
ylabel('Time (s)');
wigb(rec,1,often_offset,t);