Xử lý số tín hiệu
Chương 4: Lọc FIR và tích chập
1. Các phương pháp xử lý khối
Khối vào gồm L mẫu:
x = [x
0
x
1
x
2
x
3
… x
L-1
]
Đáp ứng xung có chiều dài M+1: (bộ lọc FIR bậc M)
h = [h
0
h
1
h
2
h
3
… h
M
]
H
x
0
x
)()()(
)(
H
x
0
x
1
x
2
x
L-1
y
0
y
1
y
2
y
3
y
4
…
1. Các phương pháp xử lý khối
b. Dạng trực tiếp (Direct form)
Bộ lọc nhân quả FIR, bậc M: h = [h
0
h
1
h
2
+ L
h
- 1
mnxmhny
m
)(
1. Các phương pháp xử lý khối
0 ≤ m ≤ M (1)
0
≤ n – m ≤ L – 1 n – L + 1 ≤ m ≤ n (2)
(1) & (2)
=> max(0, n – L + 1)
≤ m ≤ min(n,M)
Công th
ức tích chập trực tiếp:
v
ới n = 0, 1, …, L + M – 1
),min(
)1,0max(
)(
Mn
Lnm
mnxmhny
h
0
x
2
h
0
x
3
h
0
x
4
h
1
h
1
x
0
h
1
x
1
h
1
x
2
h
1
x
3
x
0
h
3
x
1
h
3
x
2
h
3
x
3
h
3
x
4
1. Các phương pháp xử lý khối
Ví dụ: tính tích chập của
h = [1, 2, -1, 1] và x = [1, 1, 2, 1, 2, 2, 1, 1]
h x 1 1 2 1 2 2 1 1
1 1 1 2 1 2 2 1 1
2 2 2 4 2 4 4 2 2
-1 -1 -1 -2 -1 -2 -2 -1 -1
1 1 1 2 1 2 2 1 1
h x 1 1 2 1 2 2 1 1
1 1 2 1 2 2 1 1
2 2 2 4 2 4 4 2 2
-1 -1 -1 -2 -1 -2 -2 -1 -1
1 2 2 1 1
h x 1 1 2 1 2 2 1 1
1
2 2
-1 -1 -1
1 2 1 1
h x 1 1 2 1 2 2 1 1
1
2
-1 -1
1 1 1
h x 1 1 2 1 2 2 1 1
1
2
-1
1 1
y = [ 1 3 3 5 3 7 4 3 3 0 1]
1. Các phương pháp xử lý khối
d) Dạng tuyến tính bất biến theo thời gian (LTI)
x = [x
0
x
1
x
2
x
3
x
4
]
1. Các phương pháp xử lý khối
h
0`
h
1
h
2
h
3
h
4
x
0
.h
0
x
0
.h
1
x
0
.h
2
x
0
.h
3
x
0
.h
2
.h
2
x
2
.h
3
x
2
.h
4
x
3
.h
0
x
3
.h
1
x
3
.h
2
x
3
.h
3
x
3
.h
3
0 0 0 0
x
0
x
0
h
0
x
0
h
1
x
0
h
2
x
0
h
3
x
1
x
1
h
0
x
1
h
1
h
0
x
3
h
1
x
3
h
2
x
3
h
3
x
4
x
4
h
0
x
4
h
1
x
4
h
2
x
4
n
1 3 3 5 3 5 -1 2
1. Các phương pháp xử lý khối
e. Dạng ma trận
+ x là vector chiều dài L
y là vector chi
ều dài L + M
+ D
ạng ma trận: y = Hx với H: ma trận (M+L) x L, xác định
t
ừ đáp ứng xung h(n)
+ D
ễ dàng thấy
00
000
h
h
h
h
h
hh
hhh
hhhh
hhh
hh
h
H
1. Các phương pháp xử lý khối
+ Cũng có thể viết:
y = X.h v
ới X là ma trận xác định từ x như sau:
000
x
xx
xxx
xxxx
xxxx
xxx
xx
x
X
1. Các phương pháp xử lý khối
f. Dạng lật và trượt
yn = h
0
x
n
+ h
1
x
n-1
+ … + h
M
x
n-M
0 0 0 x
0
x
1
x
2
h
2
h
1
h
0
h
3
h
2
h
1
h
0
h
3
h
2
h
1
h
0
y
0
y
1
y
n
y
L-1+M
x(M) + h
1
x(M-1) + … + h
M-1
x(1) + h
M
x(0)
y(M+1) = h
0
x(M+1) + h
1
x(M) + … + h
M-1
x(2) + h
M
x(1)
…
y(L-1) = h
0
x(L-1) + h
1
x(L-2) + … + h
M
x(L-1-M)
=> kho
ảng thời gian [M; L-1]: trạng thái tĩnh (xác lập)
y(L) = h
1
x(L-1) + h
2
=
y
2
=
n = 0 n = L n = 2L n = 3L
1. Các phương pháp xử lý khối
Ví dụ :
Tính tích chập x = [1, 1, 2, 1, 2, 2, 1, 1] với h = [1, 2, -1, 1]
b
ằng pp cộng dồn khối, chọn L = 3
Giải :
Chia ngõ vào thành các khối nhỏ:
Ch
ập từng khối nhỏ với h, ở đây sử dụng bảng tích chập.
x =[ 1, 1, 2, 1, 2, 2, 1, 1, 0 ]
1. Các phương pháp xử lý khối
Block 0 Block 1 Block 2
h x 1 1 2 1 2 2 1 1 0
1 1 1 2 1 2 2 1 1 0
2 2 2 4 2 4 4 2 2 0
-1 -1 -1 -2 -1 -2 -2 -1 -1 0
1 1 1 2 1 2 2 1 1 1
n 0 1 2 3 4 5 6 7 8 9 10
y0 1 3 3 4 -1 2
y1 1 4 5 3 0 2
y2 1 3 1 0 1
y 1 3 3 5 3 7 4 3 3 0 1
2. Phương pháp xử lý mẫu
Các khối cơ bản của hệ thống DSP
Khối cộng:
(n) = x(n-1)
z
-1
z
-1
x(n) x(n-1)
2. Phương pháp xử lý mẫu
- Đối với bộ làm trễ bậc cao hơn: dùng nhiều biến trạng
thái n
ội (thanh ghi hơn).
T
ổng quát, khi trễ D đơn vị, nội dung các thanh ghi là ω
i
(n),
i = 1, 2, …, D. Ký hi
ệu ngõ vào ω
0
(n)
Phương trình I/O của bộ trễ D đơn vị:
y(n) =
ω
D
(n)
ω
0
(n) = x(n)
ω
i
(n+1) = ω
i-1
-1
z
-1
x(n – 2)
x(n)
x(n – 1)
y(n)
2. Phương pháp xử lý mẫu
Đặt các trạng thái nội:
Giải thuật xử lý mẫu: với mỗi mẫu vào x(n):
ω
0
= x
y = h
0
ω
0
+ h
1
ω
1
+ h
2
ω
2
+ h
3
ω
3
ω
2
(n)
y(n)
ω
0
(n)
ω
3
(n)
2. Phương pháp xử lý mẫu
2. Phương pháp xử lý mẫu
Ví dụ: Xác định thuật toán sử lý mẫu trực tiếp, với
h = [1, 2, -1, 1]
x = [1, 1, 2, 1, 2, 2, 1, 1]
S
ử dụng thuật toán để tính đáp ứng ngõ ra.
Giải:
Phương trình I/O của bộ lọc:
y(n) = x(n) + 2x(n – 1) – x(n – 2) + x(n – 3)
V
ới trạng thái nội ω
i
(n) = x(n – 1), i = 1, 2, 3 và đặt ω
0
(n) = x(n).