地震滤波1.m

x=f2;
fs=100; % 采样率
N=length（x）; % 数据长度
n=0:N-1;
t=n/fs; % 时间系列
y=cumtrapz（tx）; % 积分变换
sz=[0;0]; %W.A 地震仪传递函数的零点
sp=[-5.4192+5.6848i;-5.4192-5.6848i]; %W.A 地震仪传递函数的极点
sk=0.999375; %W.A 地震仪传递函数的增益
[zzzpzk]=bilinear（szspskfs）;