推导一下这个一阶前向差分法的公式。
首先根据拉普拉斯变换的时域微分性质，有 d x ( t ) d t ↔ s X ( s ) \frac{dx(t)}{dt} \leftrightarrow sX(s) dtdx(t)↔sX(s)（网上有很多写的是 d x ( t ) d t ↔ s X ( s ) − f ( 0 ) \frac{dx(t)}{dt} \leftrightarrow sX(s)-f(0) dtdx(t)↔sX(s)−f(0)，不过计算机中的音频数据一般是因果信号，所以f(0) = 0）,我们利用差分来近似微分：
d x ( t ) d t = y [ n ] ≈ x [ n + 1 ] − x [ n ] T s \frac{dx(t)}{dt}=y[n] ≈ \frac{x[n+1]-x[n]}{T_s} dtdx(t)=y[n]≈Tsx[n+1]−x[n]
其中 T s T_s Ts是采样时间，即两个采样点之间时间间隔（对于连续系统，采样时间间隔为0）；故在T足够小时，将上式由时域变换到z域，并且按照 Z变换的时移性质 ： x [ n − n 0 ] ↔ z − n 0 X ( z ) x[n-n0] \leftrightarrow z^{-n_0}X(z) x[n−n0]↔z−n0X(z)，其中 n 0 n_0 n0=-1,即 x [ i + 1 ] ↔ z 1 X ( z ) x[i+1] \leftrightarrow z^{1}X(z) x[i+1]↔z1X(z),得：
Y [ z ] = z X [ z ] − X [ z ] T s Y[z] = \frac{zX[z]-X[z]}{T_s} Y[z]=TszX[z]−X[z]
即 Y [ z ] = ( z − 1 ) X [ z ] T s Y[z] = \frac{(z-1)X[z]}{T_s} Y[z]=Ts(z−1)X[z]
传递函数就是：
H ( z ) = ( z − 1 ) T s H(z) = \frac{(z-1)}{T_s} H(z)=Ts(z−1)
为了近似微分，用s替换H(z):
s = ( z − 1 ) T s s = \frac{(z-1)}{T_s} s=Ts(z−1)
即可得出公式。
这里其实会有让人疑惑的地方，为什么H(z)这么自然而然变为了s?,其实就是时域微分性质，
下面的推导会更简单：根据公式将s域离散化：
z
=
e
s
T
s
z = e^{sT_s}
z=esTs（这里类似于时域采样，只是是在复平面按采样时间间隔
T
s
T_s
Ts进行采样），当
s
T
s
sT_s
sTs很小时，根据泰勒展开式（或者等价无穷小）：
e
s
T
s
≈
1
+
s
T
s
e^{sT_s} ≈ 1 + sT_s
esTs≈1+sTs,得
z
≈
1
+
s
T
s
z ≈ 1 + sT_s
z≈1+sTs,整理上式，得
s = ( z − 1 ) T s s = \frac{(z-1)}{T_s} s=Ts(z−1)
1885
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
