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）;