求助：如何模拟下地震记录，进行fk滤波？

我需要对地震数据sac文件进行fk滤波，在此之前先进行模拟地震记录，对其进行fk滤波，同时问一下fk滤波程序matlab咋写？