Linear Constant-Coefficient Difference Equations: M N X B K N y A
Linear Constant-Coefficient Difference Equations: M N X B K N y A
Difference Equations
N
k =0
a k y[n k ] =
bm x[n m]
for all n
m =0
b0
y[n]
TD
TD
b1
x[n-1]
a1
TD
TD
b2
x[n-2]
a2
y[n-2]
TD
TD
x[n-M]
y[n-1]
bM
aN
y[n-N]
y[n] =
bm x[n m]
m =0
Accumulator
Example
n
y[n] = x[k ]
k =
= x[n] +
n 1
k =
Example (continue)
y[n] =
M2
x[n k ]
M 2 +1
1
k =0
y[n] y[n 1] =
1
M 2 +1
y[n]
TD
x[n-1]
TD
x[n-2]
TD
x[n-M]
where b = 1/ (M2+1)
and TD denotes onesample delay
Second,
Solution of Difference
Equation
ak y[n k ] = 0
k =1
Solution of Difference
Equation
M
ak
bm
y[n] =
y[n k ] +
x[n m]
k =1 a0
m = 0 a0
M
ak
bk
y[n N ] =
y[n k ] +
x[n m]
k =1 a N
m =0 a N
LTI solution
If the input x[n] is zero for n less than some time n0, the
output y[n] is also zero for n less than n0.
The previous example does not satisfy this condition since
x[n] = 0 for n<0 but y[1] = c.
Frequency-Domain Representation of
Discrete-time Signals and Systems
y[n ] =
k =
h[k ]e jwk
k =
Eigenfunction of LTI
( ) h[k ]e
Let H e jw =
k =
[]
jwk
( )
we have y n = H e e
consequently, ejwn is the eigenfunction of the
system, and the associated eigenvalue is H(ejw).
We call H(ejw) the LTI systems frequency
response that consists of the real and
imaginary parts, H(ejw) = HR(ejw) + jHI(ejw), or in
terms of magnitude and phase,
jw
jwn
( ) = H (e )e
He
jw
jw
( )
jH e jw
Example of Frequency
Response
( )= e
He
jw
jwnd
( )
H e jw = 1,
( )
H e jw = wnd
Linear Combination
x[n] = k e jwk n
k
( )
A j jw0 n A j jw0 n
+ e e
x1 [n] + x2 [n]
x[n] = A cos(w0 n + ) = e e
2
2
( )
( )
( )
Difference to Continuous-time
System Response
( ) h[k ]e
H e jw =
k =
jwk
k =
M1 n M 2
h[n] = M 1 + M 2 + 1
0
otherwise
[ ]
m +1
= 1
k =n
m
, m>n
[ ]
H e jw
=
1
e jwM 1 e jw(M 2 +1)
=
M1 + M 2 + 1
1 e jw
1
e jw(M 1 + M 2 +1) / 2 e jw(M 1 + M 2 +1) / 2
M1 + M 2 + 1
1
=
M1 + M 2 + 1
1 e jw
e jw(M 1 + M 2 +1) / 2 e jw(M 1 + M 2 +1) / 2
e jw(M 2 M 1 +1) / 2
e jw(M 2 M 1 ) / 2
e jw / 2 e jw / 2
sin[ jw(M 1 + M 2 + 1) / 2] jw(M 2 M 1 ) / 2
1
e
=
M1 + M 2 + 1
sin[w / 2]
( ) exp
= He
jw
( )
jH e jw
Phase
response
2w
n
y[n] = h[k ]x[n k ] = h[k ]e jwk e jwn
k =
k
=0
n<0
n0
jwn
jwk
jwn
jwk
e
e
y[n] = h[k ]e
h[k ]e
k
k
n
=
0
=
+
1
jw jwn
jwk jwn
e
=H e e
h[k ]e
k = n +1
( )
[]
yt [n] =
( )
k = n +1
Transient response
yt [n] =
k = n +1
h[k ] Qn
k = n +1
Representation of Sequences by
Fourier Transforms
1
x[n] =
2
( )
X e jw e jwn dw
( ) = x[n]e
X e
jw
n =
jwn
Inverse Fourier
transform
Representation of Sequences by
Fourier Transforms (continue)
X (e jw ) = X (e jw )e jX (e )
jw
1
h[n] =
2
( )
H e jw e jwn dw
1
jwm jwn
x n
x[m]e
e dw
2 m =
1
=
x[m]
2
m =
1
Since
2
jw(n m )
jw(n m )
dw
sin (n m ) 1 m = n
dw =
=
= [n m]
(n m )
0 m n
x n =
x[m] [n m] = x[n]
m =
Absolutely summable
x[n] <
n =
Mean-square convergence:
lim
( )
( )
X e jw X M e jw
( ) x[n]e
dw = 0 , for X M e jw =
jwn
n=M
( ) 2 (w + 2r )
X e jw =
r =
( ) = 2 (w w
X e
+ 2r )
r =
1
x(n ) =
2
2 (w w0 )e jwn dw = e jw0 n
( ) = 1 ae
X e
jw
Linearity
x1[n] X1(ejw), x2[n] X2(ejw)
implies that
Time shifting
x[n] X(ejw)
implies that
( )
x[n nd ] e jwnd X e jw
Frequency shifting
x[n] X(ejw)
implies that
jw0 n
x[n] X e
j (w w0 )
Time reversal
x[n] X(ejw)
If the sequence is time reversed, then x[n] X(ejw)
Furthermore, if x[n] is real then x[n] X*(ejw), since X(ejw)
is conjugate symmetric.
Differentiation in frequency
x[n] X(ejw)
( )
implies that
dX e jw
nx[n] j
dw
Parsevals theorem
x[n] X(ejw)
implies that
E=
x[n]
n =
1
=
2
X (e )
jw
dw
( )
Ye
jw
1
=
2
X (e )W (e
j (w )
)d
a periodic
convolution
[n n0 ] e
jwn0
1 ( < n < )
2 [w + 2k ]
k =
a u[n] ( a < 1)
n
1
1 ae
jw
1
1 ae
jw
[w + 2k ]
k =
r n sin w p (n + 1)
sin w p
(1 ae )
jw 2
u[n] ( r < 1)
( )
1
1 2r cos w p e jw + r 2 e j 2 w
1
w < wc
sin wc n
jw
X e =
n
0 wc < w
e
sin[w / 2]
0 otherwise
e jw0 n
2 (w w0 + 2k )
k =
cos[w0 n + ]
j
j
(
)
w
+
2
k
+
e
(w + w0 + 2k )
0
k =
x2 [n] = x1 [n 5]
( )
1
1 ae
jw
( )
= X 1 e jw
( )
Thus, X 2 e jw = e j 5 w X 1 e jw =
( )
x[n] = a 5 x2 [n], so X e jw =
e j 5w
1 ae jw
a 5e j 5w
1 ae jw
Another Example
a u[n] ( a < 1)
n
1
1 ae
jw