总结:从连续到离散，从时域到ｓ域，到Ｚ域，到差分方程，到程序

f(t)通过傅里叶变换变成F(s)，或者通过时域微分方程推出，微分算子df/dt = s,积分算子1/s.然后s域到z域,s = (2./T) * (1 + z-1) / (1 - z-1).Tustin双线性法.,最后变成差分方程,f[k] = z-1，这个公式忘了，z变换里面的。