x=f2;
fs=100;?%?采樣率
N=length（x）;?%?數(shù)據(jù)長(zhǎng)度
n=0:N-1;
t=n/fs;?%?時(shí)間系列
y=cumtrapz（tx）;?%?積分變換
sz=[0;0];?%W.A?地震儀傳遞函數(shù)的零點(diǎn)
sp=[-5.4192+5.6848i;-5.4192-5.6848i];?%W.A?地震儀傳遞函數(shù)的極點(diǎn)
sk=0.999375;?%W.A?地震儀傳遞函數(shù)的增益
[zzzpzk]=bilinear（szspskfs）;