EEE40003 Digital Signal and Image Processing: LAB 3: Discrete LTI Systems

EEE40003 Digital Signal and Image Processing
SWINBURNE UNIVERSITY OF TECHNOLOGY
EEE40003
Signal and
Processing
Digital
Image
EEE40003 Digital Signal and Image

Processing
Yew Wee Wong 10012429
LAB 3: Discrete LTI Systems
Page 0
Table of Contents
Aim........................................................................................................................ 1
Result and Analysis................................................................................................ 2
Part B: A Simple FIR System............................................................................... 2
Part C: A Simple IIR System................................................................................9
Conclusion........................................................................................................... 11
Aim
Familiarize with MATLAB

Examinnation of a simple FIR system
Examination of a simple IIR system
Page 1
Result and Analysis

Part B: A Simple FIR System
1)
By comparing both equations, the theoretical unit pulse response of the system is determined
as below:
1
1
1
1
1
h(0) , h(1) , h( 2) , h(3) , h(4)
5
5
5
5
5
h( r ) 0, r 0
Since the system shows characteristic such that
and
system is causal.
2)
M-file for maverage (x)
Page 2
1
h(r ) , r 0,1,2,3,4
5
, the

3)
Code
Result
Comment: The result of the unit pulse response complements the theoretical response in Part
B (1) where
1
h(r ) , r 0,1,2,3,4
5
and 0 for everywhere else.
Page 3

4)
Code
Result
Comment: From the result obtained, we can see that the first 5 elements of the output signal
are distorted because the number of input sequence is less than the window- averaging length
(N = 5). When the input sequences are 5, the reconstruction of the output signal is similar to
the input similar.
Page 4

5)
Code
Result
Comment: As expected, the first 5 elements of the output signal are distorted due to the
limitation of window averaging. The result above also shows that the noise of the original
signal is greatly reduced by using window averaging without altering the shape of the
original signal.
Page 5

6) and 7)
Code
Result
Comment: The result by using convolution to calculate the unit pulse response of the system
is the simular to the result obtained in step 3) above. However, there are a slight difference in
both techniques where the result by using convolution extend up to the 36th elements while
the result obtained in step 3) above only reaches up to the 32th elements.
8)
Code
Result
Comment: By using convolution function, the noise of the original BV waveform is greatly
reduced which is similar to the result obtained in step 5). However, the result obtained by
using convolution has a slight distortion at the last few elements.
Page 6

9)
a0 y (n) a1 y (n 1) a2 y (n 2) ... bo x(n) b1 x(n 1) b2 x (n 2) ... [5]
y ( n)
1
1
1
1
1
x (n) x(n 1) x(n 2) x( n 3) x (n 4) [3]
5
5
5
5
5
By comparing [3] with [5], we get:

a0 1, an 0
n0
for
1
b0 , b1 , b2 , b3 , b4 , bn 0
n otherwise
5
for
Code
Result
Comment: The result obtained is as expected as it is the same as the result obtained in step
3). The amplitude of the unit impulse response is 0.2 at 0 t 4.
Page 7

10)
a0 y (n) a1 y (n 1) a2 y (n 2) ... bo x(n) b1 x (n 1) ... b5 x(n 5) ... [5]
1
y (n) y ( n 1) [ x (n) x(n 5)] [6]
5
By comparing equation [6] with equation [5], we get:
a0 1, a1 1, an 0
n otherwise
for
1
b0 0.18, b5 , bn 0
n otherwise
5
for
Code
Result
Comment: Since the result of the unit impulse response for equation (6) is same as the result
obtained for equation (3), thus, the system given for equation (3) and (6) are the same. It is
surprising as both equations have different equation coefficients.
Question 1:
From Equation (3),
1
y (n) [( x(n) x(n 1) x( n 2) x (n 3) x(n 4)] [1]
5
1
y (n 1) [( x( n 1) x(n 2) x(n 3) x (n 4) x (n 5)] [2]
5
By subtracting [2] from [1], Eq [6] can be obtained:
1
y (n) y (n 1) [ x (n) x(n 5)] [6]
5
Since equation [3] can be modified to acquire equation [6], both system shows the same
characteristic, thus, both equations produce the same output.
Page 8
Page 9

Part C: A Simple IIR System
1)
y (n) y (n 1) (1 ) x( n)
0.82
By substituting
,
y (n) 0.82 y (n 1) (1 0.82) x(n)
Comparing the formula with:
a0 y (n) a1 y (n 1) a2 y (n 2) ... bo x( n) b1 x (n 1) b2 x(n 2) ...
Therefore,
a0 1, a1 0.82, an 0
for
b0 0.18, bn 0
for
n otherwise
n otherwise
2)
Code
Result
Code
Result
3)
Page 10
4)
Code
Original BV Wave
Filtered BV Wave
Question 2:
y (n) (1 ) x(n) y (n 1)
Page 11
(1 ) x(n)
y ( n 1)
y (n 1)
Comment: The diagram above proves that the equation in step (1) follows from the block
diagram given.
Question 3: This discrete system shows similarity with a RC circuit. The characteristic of the
system shown in step (2) demonstrates the discharging phase of the capacitors voltage while
step (3) shows the charging phase of the capacitors voltage.
Conclusion
In a nutshell, student successfully indulges in more familiarization with Matlab in Part A. In
Part B and C of the lab, student successfully examines the characteristic of a simple FIR
system and IIR system and provides relevant comment of their thoughts about their
examination.
Page 12

EEE40003 Digital Signal and Image Processing: LAB 3: Discrete LTI Systems

Uploaded by

EEE40003 Digital Signal and Image Processing: LAB 3: Discrete LTI Systems

Uploaded by

EEE40003 Digital Signal and Image Processing

SWINBURNE UNIVERSITY OF TECHNOLOGY