5.3.3 Phương pháp lấy mẫu tần số 5.3.3a Cơ sở của phương pháp lấy mẫy tần số

Phương pháp lấy mẫu tần số sử dụng phép biến đổi Fourier rời rạc (DFT). Cơ sở của phương pháp lấy mẫu tần số là xấp xỉ đặc tính biên độ tần số HN(ej) của bộ lọc số cần tổng hợp theo đặc tính biên độ tần số H(ej) của bộ lọc số lý tưởng cùng loại.

Việc xấp xỉ được thực hiện bằng cách lấy mẫu tần số qua DFT, tức là làm cho các mẫu của HN(ej) và H(ej) bằng nhau tại các tần số rời rạc k = k1 = (k.2/N) :

Hay : [5.3-9]

Bằng cách như vậy, tại các điểm tần số rời rạc k = k1 , sai số xấp xỉ giữa HN(ej) và H(ej) bằng 0, còn tại các tần số ở giữa khoảng k1 và (k + 1)1 thì sai số xấp xỉ là hữu hạn. Sai số xấp xỉ sẽ giảm nhỏ nếu giảm tần số lấy mẫu cơ bản 1 = (2/N), điều đó tương ứng với tăng độ dài N của đặc tính xung h(n)N của bộ lọc số được tổng hợp.

Trong miền k của DFT, biểu thức [5.3-9] có dạng :

Hay : Trong đó, H(k)N được lấy mẫu tần số từ đặc tính biên độ tần số H(ej) của bộ lọc số lý

tưởng cùng loại, tức là :

Hình 5.21 mô tả cách lấy mẫu đặc tính biên độ tần sốH(ej)của bộ lọc dải thông lý tưởng có các tần số cắt :

c1 = 3/10 = 0,94c2 = 8/10 = 2,51

Việc lấy mẫu tần số được thực hiện trong một chu kỳ 0 < 2 , ứng với 0 k 9 và N = 10.

Đặc tính biên độ tần số được lấy mẫu H(k)10trên hình 5.21 có

H(ej)

H(k)10

Hình 5.21 : Lấy mẫu H(ej) H(k)10.

dạng đối xứng khi k [1 , 9] . Trong dải thông của bộ lọc lý tưởng, H(k)10có sáu mẫu giá trị 1 , năm mẫu ở ngoài dải thông giá trị 0 :

Sau khi xác định được , theo [4.4-2] ở chương bốn có :[5.3-10]

Từ các mẫu DFT H(k)N có thể tìm được đặc tính tần số HN(ej) của bộ lọc số cần tổng hợp theo công thức nội suy [tl] :

[5.3-11]

Từ [5.3-11] , đối với bộ lọc số số FIR pha tuyến tính loại 1 và loại 2 có :

Trong đó :

Do đó có đặc tính tần số của bộ lọc FIR pha tuyến tính loại 1 và loại 2 :

216

3 82 97 101 4 60 5

k

1

1

[5.3-12]

Tương tự, đặc tính tần số của bộ lọc FIR pha tuyến tính loại 3 và loại 4 :

[5.3-13]

Từ [5.3-12] và [5.3-13] có biểu thức xác định đặc tính biên độ tần số của cả bốn loại bộ lọc số FIR pha tuyến tính cần tổng hợp :

[5.3-14]

Đặc tính pha () của bộ lọc FIR pha tuyến tính loại 1 và loại 2 :

; với [5.3-15]Đặc tính pha () của bộ lọc FIR pha tuyến tính loại 3 và loại 4 :

; với và [5.3-16]

Khi đặc tính tần số HN(ej) của bộ lọc cần tổng hợp thỏa mãn các chỉ tiêu kỹ thuật đã cho, thì bước tiếp theo có thể xác định đặc tính xung h(n)N của các bộ lọc số cần tổng hợp theo IDFT :

[5.3-19]Các bộ lọc số FIR pha tuyến tính loại 1 và loại 2 có đặc tính xung h(n)N đối xứng khi 0 n

(N - 1) , và (k) có dạng [4.4-12] và [4.4-16] :

[5.3-20]

Hơn nữa, bộ lọc số FIR pha tuyến tính loại 1 có N lẻ, theo mục 4.4.2 chương bốn thì A(k)N

đối xứng trong khoảng 1 k (N - 1) và đặc tính xung h(n)N của bộ lọc số cần tổng hợp được xác định theo [4.4-13] :

[5.3-21] Còn bộ lọc

số FIR pha tuyến tính loại 2 có N chẵn, theo mục 4.4.3 chương bốn thì A(k)N phản đối xứng trong khoảng 1 k (N - 1) và đặc tính xung h(n)N của bộ lọc số cần tổng hợp được xác định theo [4.4-17] :

[5.3-22]

Các bộ lọc số FIR pha tuyến tính loại 3 và loại 4 có đặc tính xung h(n)N phản đối xứng khi 0 n (N - 1), (k) có dạng [4.4-20] và [4.4-24] :

[5.3-23]

Hơn nữa, bộ lọc số FIR pha tuyến tính loại 3 có N lẻ, theo mục 4.4.4 chương bốn thì A(k)N

phản đối xứng trong khoảng 1 k (N - 1) và đặc tính xung h(n)N của bộ lọc cần tổng hợp được xác định theo [4.4-21] :

[5.3-24] Còn bộ lọc

số FIR pha tuyến tính loại 4 có N chẵn, theo mục 4.4.5 chương bốn thì A(k)N đối xứng trong khoảng 1 k (N - 1) và đặc tính xung h(n)N của bộ lọc số cần tổng hợp được xác định theo [4.4-25] :

[5.3-25]

5.3.3b Các bước tổng hợp bộ lọc số theo phương pháp lấy mẫu tần sốCác bước tổng hợp đặc tính xung h(n)N của bộ lọc số FIR pha tuyến tính theo phương pháp

lấy mẫu tần số như sau :

217

Bước 1 : Chọn số điểm lấy mẫu N [chính là độ dài của đặc tính xung h(n)N ]. Thực hiện lấy mẫu đặc tính biên độ tần số H(ej) của bộ lọc lý tưởng cùng loại trong một chu kỳ 0 < 2 để nhận được đặc tính biên độ tần số rời rạc H(k)N của bộ lọc số FIR pha tuyến tính cần tổng hợp.

Bước 2 : Xác định đặc tính biên độ tần số HN(ej) của bộ lọc số FIR pha tuyến tính cần tổng hợp bằng biểu thức nội suy [5.3-14] .

Để tìm HN(ej) theo [5.3-14] , trước hết cần xác định A(k)N :- Bộ lọc loại 1 có A(k)N đối xứng trong khoảng 1 k (N - 1). - Bộ lọc loại 2 có A(k)N phản đối xứng trong khoảng 1 k (N - 1). - Bộ lọc loại 3 : A(k)N phản đối xứng trong khoảng 1 k (N - 1), A(0)N = 0.- Bộ lọc loại 4 có A(k)N đối xứng trong khoảng 1 k (N - 1), và A(0)N = 0 .

Bước 3 : Kiểm tra đặc tính biên độ tần số HN(ej) có đạt các chỉ tiêu kỹ thuật đã cho 1 , 2 , c

, hay không ? Nếu đạt tất cả các chỉ tiêu kỹ thuật đã cho thì giảm số điểm lấy mẫu N và thực hiện lại các

bước trên cho đến khi chọn được Nmin đảm bảo đạt tất cả các chỉ tiêu kỹ thuật đã cho. Nếu không đạt thì tăng số điểm lấy mẫu N và thực hiện lại các bước trên cho đến khi chọn

được Nmin để HN(ej) của bộ lọc cần tổng hợp đạt được tất cả các chỉ tiêu kỹ thuật đã cho.Bước 4 : Xác định đặc tính xung h(n)N của bộ lọc số FIR pha tuyến tính cần tổng hợp :

[5.3-26]- Đối với bộ lọc loại 1 , để tìm h(n)N có thể tính IDFT theo [5.3-21]. - Đối với bộ lọc loại 2 , để tìm h(n)N có thể tính IDFT theo [5.3-22].- Đối với bộ lọc loại 3 , để tìm h(n)N có thể tính IDFT theo [5.3-24].- Đối với bộ lọc loại 4 , để tìm h(n)N có thể tính IDFT theo [5.3-25].

Nếu N lớn, thì có thể sử dụng các thuật toán FFT để tính IDFT [5.3-26].

218