DSP Formula
DSP Formula
04, 2006-04-25
sin(W t) 1, |Ω| < W Zeros: B(z) = 0; Poles: A(z) = 0; where H(z) = B(z)/A(z)
πt ↔ X(jΩ) =
0, |Ω| > W Frequency, magnitude/amplitude, phase response, z ← ejω
jω
δ(t) ↔ 1 H(ejω ) = |H(ejω )| ej∠H(e )
δ(t − tk ) ↔ ejΩtk jω
H(e ) = H(z)|z=ejω
1
e−at µ(t) ↔ a+jΩ , where Real{a} > 0 H[k] = H(ejω )|ω=2πk/N
d
Group delay τ (ω) = − dω ∠H(ejω )
Discrete-time Fourier-transform (DTFT):
1 jω Important transform pairs and properties:
) ejωn dω
R
x[n] = 2π 2π∞X(e
X(e ) = n=−∞ x[n] e−jωn , X(ejω ) periodic 2π
jω
P a δ[n − k] ↔ a e−jkω ↔ a z −k
an µ[n] ↔ 1/[1 − a e−jω ] ↔ 1/[1 − a z −1 ]
x[n − k] ↔ e−jkω X(ejω )
a x[n − k] ↔ a e−jkω X(ejω ) ↔ a z −k X(z)
ejωk n x[n] ↔ X(ej(ω−ωk ) )
y[n] = h[n] ⊛ x[n] ↔ Y (z) = H(z) · X(z)
x1 [n] ⊛ x2 [n] ↔RX1 (ejω ) · X2 (ejω )
1 jθ j(ω−θ) rectangular ↔ sinc, sinc ↔ rectangular
x1 [n] · x2 ↔ 2π 2π X1 (e )X2 (e ) dθ
d jω LTI filter design (synthesis)
nx[n] ↔ j dωP X(e )
ejω0 n ↔ 2π P l δ(ω − ω0 − 2πl)
Bilinear transform H(z) = H(s)|s and prewarping
cos(ω0 n) ↔ π Pl [δ(ω − ω0 − 2πl) + δ(ω + ω0 − 2πl) s = k · (1 − z −1 )/(1 + z −1 ), k = 1 or k = 2/T = 2fT
sin(ω0 n) ↔ jπ P l [δ(ω + ω0 − 2πl) − δ(ω − ω0 − 2πl) Ωprewarp,c = k · tan(ωc /2), k = 1 or k = 2/T = 2fT
x[n] = ( 1 ↔ 2π l δ(ω − 2πl) Spectral transformations, ω̂c desired cut-off
1, |n| ≤ N1 LP-LP z −1 = (ẑ −1 − α)/(1 − αẑ −1 ), where
x[n] = ↔ sin(ω(N 1 +0.5))
sin(ω/2) α = sin(0.5(ωc − ω̂c ))/ sin(0.5(ωc + ω̂c ))
0, |n| > N1
sin(W n)
LP-HP z −1 = −(ẑ −1 + α)/(1 + αẑ −1 ), where
W Wn
πn = πsinc(
( π ) ↔ ... α = − cos(0.5(ωc + ω̂c ))/ cos(0.5(ωc − ω̂c ))
1, 0 ≤ |ω| ≤ W Windowed(Fourier series method
. . . X(ejω ) = 1, |ω| < ωc
0, W < |ω| ≤ π H(ejω ) = ↔ ...
δ[n] ↔ 1 0, |ω| ≥ ωc
δ[n − k] ↔ e−jkω . . . h[n] = sin(ω c n)
= ωπc sinc( ωπc n )
πn
an µ[n] ↔ 1−ae1−jω , |a| < 1 hF IR [n] = hideal [n]
1
R π · w[n]
Discrete Fourier-transform (DFT): HF IR (ejω ) = 2π jθ
−π Hideal (e )W (e
j(ω−θ)
) dθ
Connection to DTFT: X[k] = X(ejω )|ω=2πk/N Fixed window functions, order N = 2M , −M ≤ n ≤ M :
WN = e−j2π/N Rectangular w[n] = 1
x[n] = N1
PN −1 −kn Hamming w[n] = 0.54 + 0.46 cos((2πn)/(2M ))
k=0 X[k]WN , 0≤k ≤N −1
PN −1 kn Hann w[n] = 0.5 · 1 + cos((2πn)/(2M ))
X[k] = n=0 x[n]WN , 0 ≤ k ≤ N − 1
Blackman w[n] = 0.42 + 0.5 cos((2πn)/(2M )) +
Laplace transform: 0.08 cos((4πn)/(2M ))
Convergence with a certain ROC (region of convergence). Bartlett w[n] = 1 − (|n|/M )
Connection to continuous-time Fourier-transform: s = jΩ
1
R σ+j∞ Multirate systems
x(t) = 2πj X(s)est ds
R ∞ σ−j∞ −st Upsampling with factor L, ↑ L
X(s) = −∞ x(t) e dt (
x[n/L], n = 0, ±L, ±2L, . . .
z-transform: xu [n] =
xu [n] = 0, otherwise
Convergence with a certain ROC (region of convergence).
ConnectionH to discrete-time Fourier-transform: z = ejω Xu (z) = X(z L ), Xu (ejω ) = X(ejωL )
1
x[n] = 2πj X(z)z n−1 dz Downsampling with factor M , ↓ M
P∞
X(z) = n=−∞ x[n]z −n xd [n] = x[nM ]
PM−1 −k
x[n − k] ↔ z −k X(z) Xd (z) = (1/M ) k=0 X(z 1/M WM ),
jω
PM−1 j(ω−2πk)/M
x1 [n] ⊛ x2 [n] ↔ X1 (z) · X2 (z) Xd (e ) = (1/M ) k=0 X(e )