gusucode.com > ​DSISoft是由加拿大地质调查局发布的用于垂直地震剖面(VSP)数据处理的免费软件包 > dsisoftv3/demos2/demoproc4.m

    %Called by demo2 to perform f-k filtering

load june16reject.mat; %loads in the rejection filter for the down going waves
demodata=fkfilt(demodata,june16reject,300,10,1,5,1); %fkfilt the 1st record

state=4;
disp('f-k filtering has been applied to the data.');
disp('Please apply Bandpass Filtering next.');