Nghiên cứu bộ lọc thích nghi và ứng dụng trong khử nhiễu tín hiệu

1

ĐẠI HỌC THÁI NGUYÊN
TRƯỜNG ĐẠI HỌC CNTT VÀ TRUYỀN THÔNG

Lê Thị Uyên

NGHIÊN CỨU BỘ LỌC THÍCH NGHI VÀ ỨNG
DỤNG TRONG KHỬ NHIỄU TÍN HIỆU

LUẬN VĂN THẠC SĨ: KHOA HỌC MÁY TÍNH

Thái Nguyên, tháng 9 năm 2015

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

2

ĐẠI HỌC THÁI NGUYÊN
TRƯỜNG ĐẠI HỌC CNTT VÀ TRUYỀN THÔNG

Lê Thị Uyên

NGHIÊN CỨU BỘ LỌC THÍCH NGHI VÀ ỨNG
DỤNG TRONG KHỬ NHIỄU TÍN HIỆU
Chuyên ngành: Khoa học máy tính
Mã số:
60 48 01

LUẬN VĂN THẠC SĨ: KHOA HỌC MÁY TÍNH

NGƯỜI HƯỚNG DẪN: PGS.TS Ngô Quốc Tạo

Thái Nguyên, tháng 9 năm 2015

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

3

LỜI CẢM ƠN
Em xin bày tỏ lòng biết ơn sâu sắc đến thầy giáo: PGS.TS Ngô Quốc Tạo,
Viện Công nghệ thông tin – Viện Khoa học và Công nghệ Việt Nam. Trong suốt
thời gian học và làm luận văn tốt nghiệp, thầy đã dành rất nhiều thời gian quý báu
để tận tình chỉ bảo, hướng dẫn, định hướng cho em trong việc nghiên cứu, thực hiện
luận văn.
Em xin được bày tỏ lòng biết ơn chân thành đến Ban giám hiệu, Phòng sau
đại học, các thầy cô giáo Trường Đại học Công nghệ thông tin và Truyền thông –
Đại học Thái Nguyên đã giảng dạy trong quá trình học tập, thực hành, làm bài tập,
đọc và nhận xét luận văn của em, giúp em hiểu thấu đáo hơn lĩnh vực mà em nghiên
cứu và những hạn chế cần khắc phục trong việc học tập, nghiên cứu và thực hiện
bản luận văn này.
Tôi xin cảm ơn chân thành đến Ban giám hiệu cùng các thầy cô khoa Hệ
thống thông tin quản lý – Trường Đại học Hải Dương đã tạo điều kiện giúp đỡ tôi
trong suốt quá trình học tập.
Cuối cùng tôi xin bày tỏ lòng biết ơn chân thành tới gia đình, người thân, tới
bạn bè, đồng nghiệp...., đã tận tình giúp đỡ và tạo điều kiện cho tôi trong suốt quá
trình học tập và hoàn thành luận văn.

Mặc dù đã có nhiều cố gắng, song Luận văn không tránh khỏi những thiếu
sót. Rất mong nhận được những ý kiến đóng góp quý báu của thầy cô và các bạn.
Thái Nguyên tháng 09 năm 2015
Tác giả

LÊ THỊ UYÊN

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

4

DANH MỤC MỘT SỐ THUẬT NGỮ THƯỜNG DÙNG TRONG
BỘ LỌC THÍCH NGHI XỬ LÝ TÍN HIỆU

- Analog filter: Bộ lọc tương tự
- Digital filte: Bộ lọc số
- Digital Signal Processing: Xử lý tín hiệu số
- Dynamic Range: Dải động
- Linear Time Invariable: Tuyến tính - bất biến thời gian
- Finite Impulse Reponse (FIR): Đáp ứng xung chiều dài hữu hạn
- Infinite - duration Impulse Reponse (IIR): Đáp ứng xung chiều dài vô hạn
- Mean Square Error (MSE): Sai số bình phương nhỏ nhất
- Lattice Filter Structures: Cấu trúc mạch lọc ô mạng
- Adaptive Filter: Bộ lọc thích nghi
- Least Mean Square (LMS): Thuật toán xử lý tín hiệu thích nghi
- Normalized Least Mean Square (NLMS): Thuật toán xử lý tín hiệu thích
nghi chuẩn hóa
- Adaptive Filter: Cấu trúc lọc

- Performance Evaluation: Đánh giá hiệu suất
- Adaptive Algorithm: Thuật toán thích nghi

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

5

DANH MỤC HÌNH VẼ
Số hiệu
(a)
Tên hình
Trang
hình
Hình 1.1
Sơ đồ dòng tín hiệu
5
Hình 1.2
Sơ đồ khối của một hệ thống thực hiện được về mặt vật lý
6
(b)
Sơ đồ mô tả cấu trúc trực tiếp mạch lọc IIR bậc
Hình 1.3
9
N
Hình 1.4
Sơ đồ cấu trúc dạng trực tiếp II của mạch lọc IIR bậc N
9
Section I.2 Sơ đồ cấu trúc dạng trực tiếp I của mạch lọc

Hình 1.5
10
IIR bậc N
Hình 1.6
Sơ đồ cấu trúc nối tiếp của mạch lọc IIR bậc 4
11
Hình 1.7
Sơ đồ cấu trúc dạng song song của mạch lọc IIR bậc 4
12
Hình 1.8
Sơ đồ cấu trúc dạng trực tiếp của mạch lọc FIR bậc M
12
Hình 1.9
Section I.3 Sơ đồ mạch lọc FIR bậc 6 ghép nối tiếp
13
Section I.4 Sơ đồ mạch lọc FIR pha tuyến tính với M =
Hình 1.10
13
6
Hình 1.11 Cấu trúc mạng mắc cáo toàn điểm cực
14
Hình 1.12 Cấu trúc thang - mắc cáo bậc N
15
Hình 2.1
Bộ lọc tuyến tính trong môi trường dừng
17
Hình 2.2
Sơ đồ khối của hệ thống xử lý tín hiệu thích nghi
18
Hình 2.3

Bộ lọc dàn hàng
19
Hình 2.4
Cấu trúc của bộ lọc FIR thích nghi dạng ngang.
22
Hình 2.5
Cấu trúc bộ lọc thích nghi dạng tổ hợp tuyến tính
23
Hình 2.6
Cấu trúc của bộ lọc IIR thích nghi
24
Hình 2.7
Cấu trúc cơ bản của bộ lọc thích nghi
25
Hình 2.8
Sơ đồ nhận dạng hệ thống
27
Hình 2.9
Nhận dạng hệ thống trong GSM
28
Hình 2.10
Khử nhiễu âm thanh
28
Hình 2.11 San bằng thích nghi
29
Hình 2.12 San bằng thích nghi trong thông tin di động
30
Hình 2.13 Hệ thống khử tiếng vang
31
Hình 2.14

Loại bỏ tiếng vang âm thanh
31
Hình 2.15 Khử tiếng vang trong mạng điện thoại
32
Hình 2.16 Cấu trúc của một bộ lọc thích nghi tổng quát
33
Hình 2.17 Sơ đồ biểu diễn thuật toán LMS
33
Hình 2.18 Sơ đồ cấu trúc bộ lọc ngang thích nghi bậc N
35
Hình 2.19 Sơ đồ khối mô phỏng thuật toán LMS
43
Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

6

Hình 3.1
Hình 3.2
Hình 3.3
Hình 3.4
Hình 3.5
Hình 3.6
Hình 3.7
Hình 3.8

Mô hình khử nhiễu thích nghi
Cấu trúc của bộ lọc FIR thích nghi dạng ngang.
Sơ đồ biểu diễn thuật toán LMS

Sơ đồ thuật toán
Giao diện của chương trình khi tạo tín hiệu đầu vào
Giao diện của chương trình khi tạo tín hiệu nhiễu
Giao diện của chương trình khi tạo các tín hiệu đầu vào
và nhiễu
Kết quả khi thực hiện quá trình lọc

Số hóa bởi Trung tâm Học liệu - ĐHTN

46
47
49
50
51
51
52
52

/>

7

MỞ ĐẦU

Ngày nay, với sự phát triển vượt bậc của khoa học kỹ thuật, các hệ
thống tương tự được thay thế dần bằng các hệ thống số. Các công nghệ mới
được ứng dụng rộng rãi cho xử lý tín hiệu. Bài toán loại bỏ nhiễu và tạp âm
luôn luôn là vấn đề lớn trong các hệ thống xử lý tín hiệu. Để loại bỏ can nhiễu
và tạp âm thường sử dụng các bộ lọc. Các bộ lọc kinh điển được thiết kế với
mục đích chọn lọc tần số (bộ lọc thông thấp, bộ lọc thông cao, bộ lọc thông

dải…) hay cực tiểu hóa bình phương trung bình của tín hiệu sai lệch. Tuy
nhiên những phương pháp này yêu cầu cần phải biết trước các đặc trưng
thống kê cơ bản của nhiễu như kỳ vọng, phương sai, hàm tương quan… giả
định nhiễu và tạp âm là những quá trình ngẫu nhiên dừng. Nhưng trong thực
tế, nhiễu và tạp âm là những quá trình ngẫu nhiên không dừng do đó các tham
số của nó thay đổi theo thời gian và do vậy việc thiết kế các bộ lọc theo
phương pháp kinh điển rất khó đạt được hiệu quả cao. Để phù hợp hơn với
điều kiện thực tế người ta đã đề xuất phương pháp xử lý tín hiệu thích nghi.
Mục đích của xử lý tín hiệu thích nghi là đạt được tín hiệu đầu ra tối ưu. Việc
nghiên cứu và xử lý tín hiệu trong môi trường không dừng dựa trên các thuật
toán xử lý thích nghi có một ý nghĩa thực tiễn rất lớn khi thiết kế các hệ thống
thông tin có độ chính xác cao.
Đề tài “Nghiên cứu bộ lọc thích nghi và ứng dụng trong khử nhiễu
tín hiệu” sẽ đi sâu vào nghiên cứu thuật toán xử lý tín hiệu thích nghi LMS
và các biến thể của nó. Từ đó thực hiện mạch xử lý tín hiệu thích nghi loại bỏ
can nhiễu.
Nội dung của luận văn bao gồm ba chương:
- Chương 1: Tổng quan về cấu trúc bộ lọc số
- Chương 2: Các bộ lọc thích nghi
- Chương 3: Mô phỏng ứng dụng khử nhiễu thích nghi

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

8

Trong quá trình làm luận văn, mặc dù đã có nhiều cố gắng, song không
thể tránh khỏi những sai sót, em rất mong nhận được sự góp ý tận tình của
Hội đồng bảo vệ để em có thể hoàn thiện luận văn hơn.

Em xin chân thành cảm ơn các thầy cô trong khoa, đặc biệt là thầy Ngô
Quốc Tạo đã tận tình hướng dẫn em hoàn thành luận văn này.
Thái Nguyên, ngày tháng năm 2015
Giáo viên hướng dẫn

Học viên thực hiện

PGS.TS Ngô Quốc Tạo

Lê Thị Uyên

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

9

Chương 1: TỔNG QUAN VỀ CẤU TRÚC BỘ LỌC SỐ

1.1. Khái niệm về bộ lọc số [2],[3]
1.1.1. Khái niệm về bộ lọc
Lọc số là quá trình rất quan trọng của xử lý tín hiệu số, vì chính những
khả năng phi thường của các bộ lọc số đã làm cho chúng trở nên rất phổ biến
như ngày nay. Các bộ lọc số gồm có hai công dụng chính: phân tích tín hiệu
và phục hồi tín hiệu. Phân tích tín hiệu được áp dụng khi tín hiệu mong muốn
bị giao thoa với các tín hiệu khác hay bị các loại nhiễu tác động vào nó. Còn
phục hồi tín hiệu là khi tín hiệu mà ta mong muốn hay cần để đánh giá, xét
nghiệm bị sai lệch đi bởi nhiều yếu tố của môi truờng tác động vào; làm cho
nó bị biến dạng gây ảnh hưởng đến kết quả đánh giá.
Có hai loại bộ lọc chính là bộ lọc tương tự và bộ lọc số. Hai loại bộ lọc

này hoàn toàn khác nhau về mặt thiết kế vật lý và cách thức hoạt động. Bộ lọc
tương tự (Analog filter) sử dụng các mạch điện tử tương tự được tạo ra từ các
thành phần như: điện trở, tụ điện và các bộ khuếch đại để tạo ra các hiệu ứng
lọc cần thiết trên các tín hiệu tương tự như điện áp, dòng điện,…Bộ lọc số
(Digital filter) sử dụng một bộ xử lý số để thực hiện các tính toán trên các tín
hiệu số. Bộ xử lý ở đây có thể là một máy tính đa năng chẳng hạn như một
máy tính hoặc một chip DSP (Digital Signal Processing) chuyên ngành.
Các bộ lọc tương tự có ưu điểm là giá thành rẻ, tác động nhanh, dải
động (Dynamic Range) về biên độ và tần số đều rộng và nói chung có thể
thực hiện được các công việc của một bộ lọc số. Tuy nhiên ngày nay bộ lọc số
đang ngày càng chiếm ưu thế.
1.1.2. Bộ lọc tần số
* Hệ thống tuyến tính bất biến - bộ lọc tần số
Các hệ thống LTI (Linear Time Invariable) hay các mạch lọc số có thể
được biểu thị dưới dạng các bộ nhân, bộ cộng và bộ trễ đơn vị (Hình 1.1) liên
y[n]=x[n]+s[n]

x[n]
Số hóa bởi
Trung tâm
+ Học liệu - ĐHTN
s[n]

/>

10

kết với nhau tạo thành sơ đồ dòng tín hiệu. Sơ đồ dòng tín hiệu thực hiện một
chức năng tính toán xác định, biểu thị bằng phương trình sai phân hoặc bằng
hàm truyền của một hệ thống hay của một mạch lọc số. Sơ đồ dòng tín hiệu

lại có nhiều dạng cấu trúc khác nhau, tuy nhiên ta luôn tìm được một cấu trúc
tối ưu hay còn gọi cấu trúc chính tắc. Đó là cấu trúc có các bộ nhân, bộ cộng
và bộ trễ đơn vị là ít nhất. Thiết lập cấu trúc là bước đầu tiên để thực thi phần
cứng và phần mềm cho mạch lọc số.

Hình 1.1: Sơ đồ dòng tín hiệu

Theo tính chất tự nhiên thì tín hiệu không tuần hoàn có năng lượng hữu
hạn có phổ bao gồm một dải liên tục. Hệ thống LTI thông qua hàm đáp ứng
tần số của nó sẽ làm suy giảm một số thành phần tần số nào đó của tín hiệu
vào đồng thời có thể khuếch đại các thành phần tần số khác. Hệ thống như
vậy có tác dụng như bộ lọc đối với tín hiệu đầu vào. Hiệu quả lọc đối với các
thành phần tần số khác nhau sẽ được thể hiện rõ thông qua đồ thị của H ( ω ) .
Mặt khác, góc pha của H ( ω ) sẽ xác định độ lệch pha của tín hiệu vào khi đi
qua hệ thống như một hàm của tần số.
Như vậy, tùy theo cách chọn các hệ số này chúng ta có thể thiết kế các bộ
lọc tần số cho phép truyền các tín hiệu với các thành phần tần số nằm trong

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

11

một dải nào đó trong khi sẽ làm suy giảm các tín hiệu có chứa các thành phần
tần số nằm trong một dải khác.
Trong trường hợp tổng quát, hệ thống tuyến tính bất biến theo thời gian
sẽ làm thay đổi phổ của tín hiệu đầu vào X ( ω ) tùy thuộc vào đáp ứng tần số
H ( ω ) để cho phổ đầu ra Y ( ω ) = H ( ω ) X ( ω ) . Với tác dụng này thì H ( ω ) có

vai trò như một hàm trọng số hoặc hàm dạng phổ đối với các thành phần tần
số khác nhau của tín hiệu đầu vào. Như vậy, khi ta chỉ quan tâm đến vai trò
này của hệ thống tuyến tính bất biến theo thời gian thì hệ thống có thể xem
như bộ lọc tần số và vì vậy thuật ngữ “Hệ thống bất biến theo thời gian” và
thuật ngữ “bộ lọc” là tương tự nhau và thông thường có thể thay đổi cho nhau.
Một hệ thống muốn thực hiện được về mặt vật lý thì nó phải nhân quả
và ổn định. Sơ đồ khối của hệ thống này được mô tả trên Hình 1.2
x(n)

Hệ thống tuyến tính bất biến
nhân quả và ổn định h(n)

y(n)

Hình 1.2: Sơ đồ khối của một hệ thống thực hiện được về mặt vật lý
Quan hệ giữa đầu vào và đầu ra của hệ thống này phải thỏa mãn điều
kiện sau đây:
y ( n) = h ( n) * x ( n) =

∞

∑ h( k) x( n −k)

k =0

L  h ( n )  = [ 0, ∞ ]
∞

∑ h ( n)

<∞

n =0

(1.1)
(1.2)
(1.3)

Các quan hệ này cho thấy chiều dài của đáp ứng xung h(n) là rất quan
trọng, các hệ số của h(n) là đặc trưng cho hệ thống. Vì thế chúng ta có thể

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

12

phân loại các hệ thống thành hai loại lớn tùy theo chiều dài của đáp ứng xung
h(n). Hai loại này như sau:
- Loại thứ nhất: Hệ thống được đặc trưng bởi đáp ứng xung có chiều
dài hữu hạn. Nó được gọi là hệ thống có đáp ứng xung chiều dài hữu hạn
(FIR), tức là h(n) chỉ khác không trong một khoảng có chiều dài hữu hạn N
(từ 0 đến N - 1).
- Loại thứ hai: Hệ thống được đặc trưng bởi đáp ứng xung có chiều dài
vô hạn. Nó được gọi là hệ thống có đáp ứng xung chiều dài vô hạn (IIR), tức
là h(n) khác không trong một khoảng vô hạn (từ 0 đến ∞ ).
Việc lựa chọn hệ thống FIR hay IIR trong khi thiết kế tùy thuộc vào
yêu cầu của bài toán và đặc tính kỹ thuật của đáp ứng tần số mong muốn.
Trên thực tế, bộ lọc FIR được sử dụng trong các bài toán lọc mà ở đó đòi hỏi
sự tuyến tính của đặc tuyến pha trong dải thông của bộ lọc. Nếu sự tuyến tính

này là không cần thiết thì có thể sử dụng bộ lọc FIR hay IIR. Tuy vậy, về
nguyên tắc chung, bộ lọc IIR có các thùy bên (búp phụ) trong dải chắn thấp
hơn so với bộ lọc FIR với cùng tham số. Do nguyên nhân này, nếu sự méo về
pha có thể chấp nhận được hoặc là không quan trọng thì khi đó bộ lọc IIR
thường được chọn do việc thiết kế bộ lọc dạng này thường đòi hỏi số lượng
tham số ít hơn, bộ nhớ cần sử dụng cũng ít hơn và độ phức tạp tính toán cũng
thấp hơn.
Việc thiết kế bộ lọc số FIR và IIR thường có độ mềm dẻo cao hơn do
khả năng thực hiện được bằng phần mềm. Tuy vậy, một điểm quan trọng cần
lưu ý khi thiết kế là phải chọn được bộ lọc phù hợp với ứng dụng và thỏa mãn
các yêu cầu của việc thiết kế.
Chúng ta đều biết rằng để có thể thực hiện bộ lọc số dễ dàng hơn, với
giá thành rẻ hơn thì số phần tử của bộ lọc phải ít nhất. Và để giảm thời gian
tính toán trong quá trình lọc thì một bộ lọc tốt nhất là bộ lọc có bậc nhỏ nhất.

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

13

Theo tiêu chí này nhiệm vụ của quá trình thiết kế bộ lọc là phải xác định các
hệ số { ak } và { bk } trong đặc tính đáp ứng tần số được đưa bởi (1.10) để có thể
đạt được các giá trị gần đúng nhất so với các tham số được yêu cầu. Bậc của
H ( ω ) được sử dụng để xấp xỉ hóa các tham số này một phần sẽ phụ thuộc vào

các tiêu chuẩn được sử dụng trong việc chọn lựa các hệ số { ak } và { bk } cũng
như số lượng (M, N) của các hệ số này.
Lọc tuyến tính và bất biến thời gian (LTI) được đặc trưng bởi đáp ứng
xung h(n) của nó. Đáp ứng đối với tín hiệu vào bất kỳ x(n) sẽ tổng chập của

h(n) với x(n). Tuy nhiên nhiều khi ta liên hệ trực tiếp tín hiệu ra và vào bằng
phương trình hiệu. Xét phương trình hiệu hay cấu trúc mạch lọc người ta chia
ra làm hai loại lớn đó là lọc phi đệ quy FIR và lọc đệ quy IIR. Trong chương
này chúng ta lần lượt xét các cấu trúc của hai loại mạch lọc đó.
1.2. Cấu trúc của bộ lọc IIR[4],[5]
Bộ lọc số IIR có đáp ứng xung dài vô hạn và được mô tả bằng phương
trình sai phân như sau:
y[ n ] =

M

N

∑b x[ n − m] − ∑ a y[ n − k ]

m =o

m

k =1

k

(1.4)

Đối với bộ lọc LTI nhân quả thì M ≤ N, ở đây N là bậc của mạch lọc số.
1.2.1. Cấu trúc dạng trực tiếp của bộ lọc IIR
Lấy biến đổi z phương trình sai phân (1.5) và sau đó lập tỉ số Y(z)/X(z)
ta sẽ thu được hàm truyền của mạch lọc này:
M

H(z) =

∑b

m =0
N

m

z −m

1 + ∑a k z −k

(1.5)

k =1

Nếu đặt H(z) = H1(z).H2(z)

Số hóa bởi Trung tâm Học liệu - ĐHTN

(1.6)

/>

14

M

Trong đó: H1(z) =

∑b

m =o

m

z −m

(1.7)

1
N

H2(z) = 1 + a z − k
∑ k

(1.8)

k =1

Thì chúng ta sẽ thu được cấu trúc trực tiếp của bộ lọc IIR. Đó là cấu
trúc gồm đủ các bộ cộng, bộ nhân và bộ trể đơn vị như trong phương trình sai
phân.

b0

x[n]

y[n]

z −1

b1

− a1

z −1

z −1

b2

− a2

z −1

z −1

bM −1

− a N −1

z −1

z −1

bM

− aN

z −1

Hình 1.3. Sơ đồ mô tả cấu trúc trực tiếp mạch lọc IIR bậc N
Bây giờ nếu thực thi hàm truyền H2(z) trước, sau đó đến H1(z) có nghĩa
là thực hiện hàm truyền H(z) dưới dạng sau:
H(z) = H2(z).H1(z)

(1.9)

Thì ta sẽ thu được cấu trúc dạng trực tiếp 2.
b0

x[n]

−a1
−a2

−a N −1
−a N

z

−1

z

−1

y[n]

b1
b2

z −1
bN −1
z −1b

N

Hình 1.4: Sơ đồ cấu trúc dạng trực tiếp II của mạch lọc IIR bậc N
Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

15

Sau khi thu được cấu trúc dạng trực tiếp 2, ta thực hiện phép chuyển vị
sao cho H1(z) thực hiện trước rồi mới đến H2(z), có nghĩa là ta thực thi hàm
truyền H(z) dưới dạng:
H(z) = H1(z).H2(z)

(1.10)

Như vậy sẽ thu được cấu trúc dạng trực tiếp 1.
b0

x[n]

y[n]

b1 z

−1

−1
b2 z

z −1
bN −1
bN

z −1

− a1
− a2

− a N −1
− aN

Hình 1.5 Sơ đồ cấu trúc dạng trực tiếp I của mạch lọc IIR bậc N
1.2.2. Cấu trúc dạng nối tiếp của bộ lọc IIR
Trong nhiều ứng dụng thực tế, hàm truyền của mạch lọc bậc N thường
được khai triển thành tích các hàm truyền bậc hai. Do vậy khi thực thi dạng
khai triển này, sẽ thu được cấu trúc dưới dạng nối tiếp các hệ thống bậc hai.
Thật vậy, hàm truyền H(z) từ (1.9) có thể khai triển dưới dạng sau:
M

H(z) =

∑bm z −m

m =0
N

1 + ∑a k z −k

1 + b1i z − 1 + b2i z − 1
∏i=1 1 + a z −1 + a z −1 = b0
1i
2i
= b0
L

L

∏ H ( z)
i

i =1

(1.11)

k =1

 N + 1

Trong đó L = 
. Vậy nếu N lẻ thì ngoài các hệ thống bậc hai còn

 2 
có một hệ thống bậc nhất ghép nối tiếp, và

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

16

1 + b1i z −1 + b2i z −1
Hi(z) =
1 + a1i z −1 + a 2i z −1

(1.12a)

là hàm truyền của mạch lọc bậc hai thứ i.

x[n]

− a11
− a21

z

−1

b11

z−1 b

21

− a12

− a22

z

y[n]

−1

b12

z −1 b

22

Hình 1.6: Sơ đồ cấu trúc nối tiếp của mạch lọc IIR bậc 4
1.2.3. Cấu trúc dạng song song của bộ lọc IIR
Nếu phân tích các hàm truyền H(z) của mạch lọc IIR bậc N (1.2) thành
tổng các hàm truyền bậc hai như sau:
M

H(z) =

∑b

m =0
N

m

z −m

1 + ∑ ak z

−k

b0i + b1i z −1
= C0 + ∑
= C0 +
−1
+ a 2i z −1
i =1 1 + a1i z
L

L

∑ H ( z)
i =1

i

(1.12 b)

k =1

 N + 1

Trong đó L = 
, thì sẽ thu được cấu trúc gồm các hệ thống bậc
 2 
hai ghép song song với nhau như trên hình 1.8. Nếu N lẻ thì ngoài các hệ
thống bậc hai còn có các hệ thống bậc nhất ghép song song. Mạch lọc bậc hai
thứ i có hàm truyền dạng:
b0i + b1i z −1
Hi(z)=
1 + a1i z −1 + a 2i z −1

Số hóa bởi Trung tâm Học liệu - ĐHTN

(1.13)

/>

17

C0
b01
z −1

− a11

x[n]

y[n]

b11

z −1

− a 21

b02
z

− a12

−1

b12

z −1

− a 22

Hình 1.7: Sơ đồ cấu trúc dạng song song của mạch lọc IIR bậc 4
1.3. Cấu trúc các bộ lọc FIR [2],[3],[4]&[5]
Mạch lọc FIR bậc M có hàm truyền H(z) dạng sau:
M

H(z)= ∑ bm z

−m

(1.14)

m =0

Quan hệ giữa tín hiệu vào x[n] và tín hiệu lối ra y[n] được biểu thị bằng
phương trình sai phân bậc M:
M

y[n]= ∑ bm x[n − m]

(1.15)

m =0

Ở đây hệ số bm (m=0,1,...,M) cũng chính là đáp ứng xung đơn vị của
mạch lọc.
1.3.1. Cấu trúc dạng trực tiếp
Phương pháp thường dùng nhất để thực hiện bộ lọc FIR là phương pháp
dạng trực tiếp (direct form), phương pháp này sử dụng đường trì hoãn rẽ
nhánh (tapped delay line) được biểu thị trên hình 1.9.

z −1

x[n]

b0

b1

z −1
b2

z −1

bM − 2

z −1

z −1

bM − 1

bM

y[n]

Hình 1.8. Sơ đồ cấu trúc dạng trực tiếp của mạch lọc FIR bậc M
Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

18

Cấu trúc này yêu cầu M + 1 phép nhân, M phép cộng và M trì hoãn.
Tuy nhiên, nếu có các đối xứng trong đáp ứng xung đơn vị, ta có thể giảm bớt
số lượng phép nhân.
1.3.2. Cấu trúc nối tiếp
Nếu phân tích hàm truyền H(z) từ (1.14) thành tích các hàm truyền bậc
hai dưới dạng:
L

M

∑ bm z −m

H(z)= m =0
Trong đó L= [

=b0

∏ (1 + b

1i

z −1 + b2i z − 2 )

(1.16)

i =1

M +1
] thì mạch lọc FIR được thực thi dưới dạng các
2

mạch lọc bậc hai ghép nối tiếp. Trường hợp M lẻ thì ngoài các mạch lọc bậc
hai còn có một mạch lọc bậc nhất ghép nối tiếp.
x[n]

z −1
z −1

z −1

b11

b12

z −1

b21

b22

z −1

b13

z −1

b23

y[n]

Hình 1.9: Sơ đồ mạch lọc FIR bậc 6 ghép nối tiếp
1.3.3. Cấu trúc mạch lọc FIR pha tuyến tính
Nếu đáp ứng xung của bộ lọc FIR thoả mãn điều kiện đối xứng:
h[n]= ± h[ M- n ], n=0,1,2,...,M-1

(1.17)

thì trong sơ đồ dòng tín hiệu có thể rút bớt được một nửa bộ nhân hệ số.

z −1

x[n]

z −1
z −1

z −1
y[n]

b0

b1

z −1

b2

z −1
b3

Hình 1.10: Sơ đồ mạch lọc FIR pha tuyến tính với M = 6
Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

19

1.4. Cấu trúc mắt cáo
Cấu trúc mạng mắt cáo hay còn gọi là cấu trúc mạch lọc ô mạng
(lattice filter structures) rất có ích trong xử lý tiếng nói và trong thực thi các
mạch lọc thích nghi dùng để tiên đoán tuyến tính. Trong xử lý tiếng nói, cấu

trúc này thường được sử dụng nhiều hơn các cấu trúc FIR và IIR, bởi vì trong
phân tích và tổng hợp tiếng nói chỉ cần một lượng nhỏ các hệ số cũng có thể
cho phép một số lượng lớn các formants được mô hình hoá theo thời gian
thực. Có hai cấu trúc mắt cáo chính là mạng toàn điểm không và mạng toàn
điểm cực. Mạng vừa có điểm không vừa có điểm cực thì được gọi là mạng
bậc thang mắt cáo (ladder).
Cấu trúc dạng mắt cáo có thể được dùng để thực hiện các bộ lọc số FIR
và cả IIR. Trong mục này chỉ mô tả đối với bộ lọc IIR.
Mạch lọc IIR toàn điểm cực có hàm truyền H(z) dạng:
1

H(z) = 1 +

N

∑a
k =1

k

(1.18)

z −k

Sẽ có sơ đồ dòng tín hiệu cho trên hình 1.11
x[n]

f M [n]
− kM
kM

g M [n]

f M −1 [n]

f 1 [ n]

− k M −1

− k1
k1

k M −1
z −1

f 0 [n]

g M −1[n]

y[n]

z −1 g [n]
1

z −1 g [ n]
0

Hình 1.11. Cấu trúc mạng mắt cáo toàn điểm cực
Biến đổi z giữa tín hiệu fm[n] và fm-1[n] liên hệ với nhau bằng hệ thức
sau:
Fm ( z ) − k m z − m Fm ( z −1 )

Fm-1(z) =
1 − k m2

Số hóa bởi Trung tâm Học liệu - ĐHTN

, m = M, M-1, ..., 1

(1.19)

/>

20

Phương trình (1.19) cho phép tính đa thức bậc thấp hơn Fm-1(z) từ Fm(z).
Do đó, phương pháp này còn có tên là phương pháp hạ cấp, bắt đầu từ m và
lui dần tới m=1.
Các hệ số phản xạ km liên hệ với các hệ số ak của H(z) từ (1.18) bằng hệ
thức truy hồi:

am-1,i=

a mi − k m a m ,m −i

(1.20)

1 − k m2

với km= amm ; m = M, M-1,..., 1 và i = 0, 1, 2,.., m-1 ;
km ≠ 1

•

Cấu trúc thang - mắt cáo(lattice- ladder)
Cấu trúc thang mắt cáo hay còn gọi là cấu trúc Gray- Markel được thực

thi đối với mạch lọc IIR tổng quát.
M

H(z) =

∑b

m=0
N

m

z −m

BM ( z )

1 + ∑ a k z −k

= A ( z)
M

(1.21)

k =1

Cấu trúc mày được thực thi theo hai bước. Bước đầu tiên là thực thi cấu
trúc mắt cáo toàn điểm cực như trong hình 1.12 với các hệ số phản xạ km, 1
≤ m ≤ N . Bước thứ hai là cộng các lối ra g m[n] với các hằng số Cm như trong

hình 1.12 để tạo ra cấu trúc thang- mắt cáo.
f M [n]

x[n]

f1 [ n]

− k M −1

− kM
kM
CM

f M −1[n]

− k1
k1

k M −1
g M [n]

z −1
C M −1

f 0 [ n]

g M −1[n]

z −1
C1

−1
g 0 [ n]
g1 [ n ] z
C0

y[n]

Hình 1.12. Cấu trúc thang - mắc cáo bậc N

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

21

Chương này nêu tổng quan về lọc số, các thông số của hệ thống ở miền
thời gian, ở miền tần số.
Phần này cũng đã nêu lên được một cách tổng quát hai cấu trúc của bộ
lọc số là: FIR và IIR.
Từ những cấu trúc của các bộ lọc số cơ bản, ta ứng dụng các thuật toán
thích nghi cải biến bộ lọc để cho ra một bộ lọc thích nghi hoàn chỉnh với
những tính năng rất thực tế và hiệu quả; làm nền tảng để tìm hiểu các bộ lọc
thích nghi. Chúng ta cùng qua chương 2 để tìm hiểu kỹ hơn về vấn đề này.

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

22

Chương 2: CÁC BỘ LỌC THÍCH NGHI

2.1. Giới thiệu lọc thích nghi
Trong các bộ lọc số quy ước (FIR và IIR), mọi thông số của quá trình
lọc dùng để xác định các đặc trưng của hệ thống coi như đã biết. Các thông số
này có thể biến đổi theo thời gian, trong một số bài toán thực tiễn cho thấy
một số thông số có độ bất ổn định cao và bản chất của sự biến thiên thì không
tiên đoán được. Để giải quyết vấn đề đó, người ta nghiên cứu thiết kế bộ lọc
sao cho có thể tự thích nghi với hoàn cảnh hiện hành, có nghĩa là nó có thể tự
điều chỉnh các hệ số trong bộ lọc để bù lại các thay đổi trong tín hiệu vào, tín
hiệu ra, hoặc trong thông số của hệ thống. Đó chính là bộ lọc thích nghi.
Các bộ lọc thích nghi (Adaptive Filter) được sử dụng tốt nhất trong các
loại bộ lọc, ở các tín hiệu có điều kiện hay các thông số hệ thống thay đổi rất
chậm và bộ lọc đã được điều chỉnh để bù cho sự thay đổi này. Thuật toán
LMS là một thuật toán dò tìm được sử dụng để cung cấp một kế hoạch quản
lý tốt việc điều chỉnh các hệ số bộ lọc, ngoài ra còn có một số thuật toán khác
cũng có khả năng thích nghi như: RLS, NLMS,… Mỗi thuật toán có các ưu,
khuyết điểm khác nhau; chúng ta sẽ tìm hiểu rõ điều đó qua các phần sau đây.
2.1.1. Một số khái niệm cơ bản
Xử lý tín hiệu thực chất là một quá trình lấy ra tín hiệu mong muốn từ
một tập tín hiệu có lẫn nhiễu tại đầu vào máy thu. Tín hiệu khi được truyền đi
trong môi trường bị biến dạng bởi các tác động của can nhiễu và tạp âm. Do
vậy tại thiết bị thu ta phải thiết kế như thế nào để càng giảm được tác động
của nhiễu càng nhiều càng tốt. Với mục đích nâng cao độ tin cậy cho thiết bị
thu thì các hệ thống thông tin cần phải tích hợp các khối xử lý để giảm ảnh

hưởng của nhiễu và tạp âm. Những khối này luôn tồn tại trong các hệ thống
thông tin tương tự cũng như các hệ thống thông tin số, chúng có thể qui về
các bộ lọc và các bộ san bằng. Một trong những ứng dụng quan trọng của các

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

23

bộ lọc là loại bỏ nhiễu và tạp âm. Các bộ lọc kinh điển được thiết kế với mục
đích chọn lọc tần số (bộ lọc thông thấp, bộ lọc thông cao, bộ lọc thông dải…)
hay cực tiểu hóa bình phương trung bình của tín hiệu sai lệch. Tuy nhiên
những phương pháp này yêu cầu cần phải biết trước các đặc trưng thống kê
cơ bản của nhiễu như kỳ vọng, phương sai, hàm tương quan…và giả định
nhiễu và tạp âm là những quá trình ngẫu nhiên dừng. Hình 2.1 mô tả cấu trúc
của một bộ lọc tuyến tính hoạt động trong môi trường dừng.
Đầu vào bộ lọc

Đầu vào bộ lọc
Bộ lọc

Sai số

Tham số bộ lọc

+

Đầu ra bộ lọc
Tham số bộ lọc

Tín hiệu
mong muốn
Đầu ra bộ lọc

+ tuyến tính trong môi trường
Hình 2.1: Bộ lọc
Sai số
Tín hiệu
dừng
Nhưng trong thực tế, nhiễu và tạp âm là mong muốn những quá trình
ngẫu nhiên không dừng do đó các tham số của nó thay đổi theo thời gian và
do vậy việc thiết kế các bộ lọc theo phương pháp kinh điển rất khó đạt được
hiệu quả cao. Để phù hợp hơn với điều kiện thực tế người ta đã đề xuất
phương pháp xử lý tín hiệu thích nghi. Mục đích của xử lý tín hiệu thích nghi
là đạt được tín hiệu đầu ra tối ưu theo nghĩa này hay nghĩa khác. Do không
biết trước được các tham số đặc trưng cho nhiễu hay tín hiệu có lẫn nhiễu tại
đầu vào máy thu nên các thuật toán xử lý tín hiệu thích nghi sẽ xử lý theo
từng mẫu dữ liệu thu được và sử dụng các mẫu đó để tìm các mẫu dữ liệu kế
tiếp theo phương pháp đệ quy. Mọi thuật toán xử lý tín hiệu thích nghi đều
xuất phát từ một tập điều kiện ban đầu. Điều kiện ban đầu chính là những gì

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

24

biết được về môi trường truyền dẫn. Trong môi trường dừng, ta sẽ tìm được
một giá trị tối ưu sau khi thực hiện một số chu kỳ xử lý thành công. Nhưng
ngược lại trong môi trường không dừng, không tồn tại một giải pháp tối ưu

duy nhất cho quá trình xử lý tín hiệu thích nghi. Để đảm bảo đạt được tín hiệu
thu tốt nhất thì các bộ lọc thích nghi vẫn phải thực hiện quá trình điều chỉnh
trọng số bộ lọc dù không biết trước được các tính chất thống kê của tín hiệu
vào. Nhưng thay vì phải đưa ra tất cả mọi thông tin về một quá trình nào đó
thì ta chỉ phải đưa ra một chuỗi mẫu tín hiệu tại thời điểm đó và sử dụng
phương pháp đệ quy để tìm các mẫu tín hiệu trong các thời điểm kế tiếp. Có
rất nhiều biện pháp để có thể tìm ra được tín hiệu mong muốn nhưng phương
pháp hiệu chỉnh theo sai số bình phương trung bình là phổ biến hơn cả. Sơ đồ
khối của hệ thống xử lý tín hiệu thích nghi được mô tả trên hình 2.2.
Đầu vào
Đầu vào bộ lọc

Bộ lọc

Tham số bộ lọc

Thích nghi

Đầu ra
-

+

-

+

Sai số
Hình 2.2: Sơ đồ khối của hệ

Tín hiệu
mong muốn
Đầu ra bộ lọc

thống xử lý tín hiệu thích
Tín
nghi
hiệu mong
Từ trên sơ đồ ta thấy rằng hoạt động của muốn
thuật toán xử lý
tín hiệu thích nghi gồm hai quá trình chính:


Sai số

Quá trình lọc: quá trình này thực hiện lấy tín hiệu ra từ tín hiệu đầu

vào.

Số hóa bởi Trung tâm Học liệu - ĐHTN

/>

25



Quá trình xử lý thích nghi: mục đích của quá trình này là điều

chỉnh hàm truyền đạt của hệ thống theo sự thay đổi của môi trường. Quá trình

thích nghi điều chỉnh theo một tín hiệu sai lệch. Thông thường thì giá trị bình
phương trung bình của tín hiệu sai lệch được sử dụng trong quá trình xử lý
thích nghi.
Hai quá trình trên xử lý luân phiên nhau. Do vậy việc lựa chọn một cấu
trúc cho quá trình lọc có ảnh hưởng rất lớn đến toàn bộ quá trình xử lý của
toàn bộ thuật toán. Có ba dạng cấu trúc bộ lọc hay được sử dụng trong thuật
toán xử lý tín hiệu thích nghi là:


Bộ lọc dàn hàng: bộ lọc này chỉ bao gồm ba phần tử cơ bản như

được chỉ ra trên Hình 2.3, đó là: bộ trễ, bộ nhân và bộ cộng. Số phần tử trễ
trong bộ lọc sẽ tương ứng với số đáp ứng xung hay bậc của bộ lọc. Tín hiệu
đầu ra của bộ lọc dàn hàng:
y ( n) =

M −1

∑ w x(n − k )
k =0

(2.1)

k

với: y (n) là tín hiệu đầu ra của bộ lọc.
wk là hệ số của bộ lọc hay còn gọi là trọng số lọc k = 1,2,...M − 1
x(n) là đầu vào của bộ lọc.
x(n)

w0

z-1

x(n-1)

z-1

x(n-2)

x(nM+2)

z-1

x(n-M+1)

w1

w2

wM-2

wM-1

∑

∑

∑

∑

y(n)

Hình 2.3. Bộ lọc dàn hàng


Bộ dự báo mắt cáo: Bộ dự báo này gồm các môđun riêng biệt,

mỗi môđun này chỉ xuất hiện trong một mắt lưới. Quan hệ giữa tín hiệu đầu
vào và tín hiệu đầu ra của bộ dự báo mắt cáo.
Số hóa bởi Trung tâm Học liệu - ĐHTN

/>