（matlab代码）地震数据带通滤波

DATA_for_bandpass=data1; % 选择要分析的数据
f1=0.001; f2=1; f3=15; f4=20;
for j=1:size(DATA_for_bandpass,2)
data_bandpass(:,j)=bandpass_filter(DATA_for_bandpass(:,j),dt,f1,f2,f3,f4);
end
data3=data_bandpass;
figure
plot(data3,‘k’)
% axis([0,size(data3,1),-1,1]);