›
›

Bài giảng Xử lý số tín hiệu DSP - Chương 7: Thiết kế bộ lọc số

pdf 28 trang phuongnguyen 3150
Download
Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Xử lý số tín hiệu DSP - Chương 7: Thiết kế bộ lọc số", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pdfbai_giang_xu_ly_so_tin_hieu_dsp_chuong_7_thiet_ke_bo_loc_so.pdf

Nội dung text: Bài giảng Xử lý số tín hiệu DSP - Chương 7: Thiết kế bộ lọc số

  1. Xử lý số tín hiệu Chương 7: Thiết kế bộ lọc số
  2. 1. Các bước thực hiện bộ lọc số.  Bộ lọc : hệ thống thực hiện hiệu chỉnh tín hiệu ở một số thành phần tần số nào đó.  Bộ lọc chọn lọc tần số: cho qua một số thành phần tần số của tín hiệu và loại bỏ tất cả các thành phần còn lại.  3 bước thiết kế bộ lọc số: Yêu cầu của bộ lọc Phụ thuộc vào ứng dụng Thiết kế bộ lọc Thực hiện trên phần cứng Phụ thuộc vào phần cứng
  3. 1. Các bước thực hiện bộ lọc số (tt) Mục đích của thiết kế bộ lọc số: Xác định hàm truyền H(z)  Đối với bộ lọc IIR: M 1 b z i i 0 i H(z) M 2 1 a z i i 1 i →Xác định các vector tham số tử số b=[b ,b , ,b ] 0 1 M1 và mẫu số a=[1,a , ,a ] 1 M2  Đối với bộ lọc FIR: M H(z) 1 b z i i 0 i → Xác định vector b=[b ,b , ,b ], đây cũng chính là 0 1 M1 đáp ứng xung của bộ lọc
  4. 2. Các yêu cầu của bộ lọc.
  5. 2. Các yêu cầu của bộ lọc. ωp: Cạnh dải thông. ω s: Cạnh dải chắn. δp1, δp2: độ gợn dải thông. δs: độ gợn (suy hao) dải chắn.
  6. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ.  Thiết kế bộ lọc FIR dùng phương pháp cửa sổ dựa trên nguyên tắc xấp xỉ đáp ứng xung của bộ lọc mong muốn.  Giả sử ta có bộ lọc có đáp ứng tần số mong muốn: H () h (n)e jn d n d  Đáp ứng xung mong muốn: 1 jn hd (n) Hd ()e d 2  VD: tìm đáp ứng xung của bộ lọc thông thấp lý tưởng: H() 1,  c
  7. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt) Đáp ứng xung của một số bộ lọc thông dụng: sin( n)  Thông thấp: c hd (n) n sin(cn)  Thông cao: hd (n)  (n) n sin(bn) sin(an)  Thông dải: hd (n) n sin(bn) sin(an)  Chắn dải: hd (n)  (n) n Nhận xét: Các bộ lọc lý tưởng có đáp ứng tần số thay đổi đột ngột giữa dải thông và dải chắn → đáp ứng xung dài vô hạn và không nhân quả.
  8. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt)  Để xấp xỉ đáp ứng xung của các bộ lọc lý tưởng này, cách đơn giản nhất là xén bớt từ đáp ứng xung lý tưởng này.  Giả sử cần thiết kế bộ lọc bậc M(M chẵn), đáp ứng xung của hệ thống xấp xỉ: M M h(n) hd (n), - n 2 2  Hay ta có thể viết cách khác: h(n) hd (n)w(n)  Với M M 1 n w(n) 2 2 0 n khác
  9. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt)  Do h(n) sau khi nhân với hàm cửa sổ là không nhân quả, ta cần làm trễ đi M/2 mẫu để có đáp ứng xung nhân quả. h(n-M/2)=hd(n-M/2)w(n-M/2)
  10. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt) VD: Dùng phương pháp cửa số chữ nhật xác định đáp ứng xung nhân quả của bộ lọc thông thấp FIR bậc 10 có tần số cắt là ωc= /4 Giải: sin( n) h(n) c , 5 n 5 n 2 2 1 2 1 2 1 2 2 h(n) ,0, , , , , , , ,0, 10 6 2 2 4 2 2 6 10 Đáp ứng xung nhân quả: 2 2 1 2 1 2 1 2 2 h(n 5) ,0, , , , , , , ,0, 10 6 2 2 4 2 2 6 10
  11. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt) Tìm lại đáp ứng tần số từ đáp ứng xung này: 5 2 2 1 1 H() e j5 h(n)e jn e j5 cos(5) cos(3) cos() n 5 5 3 4 1 2 2 1 1 H() cos(5) cos(3) cos() 5 3 4 0.8 H() 5 :pha tuyến tính. 0.6 )|  Nhận xét: |H( _ Bộ lọc kết quả có độ dốc giữa dải thông và 0.4 dải chắn lài. _ Đáp ứng biên độ của bộ lọc kết quả có 0.2 những gợn trong cả dải thông và dải chắn. 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 /
  12. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt)  Hiện tượng Gibbs 1 H() Hd ()W( )d 2 M / 2 sin(M 1) / 2 W() e jn n M / 2 sin( / 2)
  13. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt)  Thay đổi chiều dài cửa sổ:  Tăng chiều dài cửa sổ: độ rộng búp chính của hàm cửa sổ giảm → khoảng chuyển tiếp nhỏ. Tuy nhiên, tần số các gợn biên tăng lên.  Giảm chiều dài cửa sổ: khoảng chuyển tiếp lớn.
  14. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt)  Để giảm bớt hiện tượng Gibbs cần thay đổi hình dạng cửa sổ: tránh thay đổi đột ngột tại cạnh hàm cửa sổ.  Cửa sổ Bartlett (cửa sổ tam giác) 2n / M 0 n M / 2 w(n) , M even 2 2n / M M / 2 n M  Cửa sổ Hann: w(n) 0.5 0.5cos(2 n/ M), 0 n M  Cửa sổ Hamming: w(n) 0.54 0.46cos(2 n/ M), 0 n M  Cửa sổ Blackman: w(n) 0.42 0.5cos(2 n/ M) 0.08cos(4 n/ M), 0 n M
  15. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt)
  16. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt) Ví dụ: Bộ lọc thông thấp dùng cửa sổ chữ nhật và cửa sổ Hamming với cùng chiều dài N=81.
  17. 3. Thiết kế bộ lọc FIR dùng phương pháp cửa sổ. (tt)  Thiết kế bộ lọc FIR bằng cửa sổ Kaiser: 2 n M / 2 I 1 0 M / 2 w(n) , 0 n M I0 ( ) Với  훼: hệ số hình dạng.  I0(x): hàm Bessel sửa đổi loại 1 bậc 0. Cửa sổ Bessel có 2 tham số giúp bù trừ độ rộng búp chính (độ rộng khoảng chuyển tiếp) và độ cao búp chính (tăng/giảm hiện tượng Gibbs) qua hệ số hình dạng và chiều dài bộ lọc ->linh động hơn khi thiết kế.
  18. 4. Thiết kế bộ lọc IIR từ các bộ lọc thời gian liên tục  Đối với bộ lọc IIR thì các phương pháp thiết kế thường bắt đầu bằng thiết kế bộ lọc tương tự, sau đó dùng các phép biến đổi để chuyển thành bộ lọc số.  Các phương pháp thiết kế bộ lọc IIR tương tự đã rất phát triển.  Các bộ lọc tương tự thường được xác định bằng các công thức cụ thể -> đơn giản hoá việc thiết kế.  Các phương pháp xấp xỉ toán học khi áp dụng cho bộ lọc số IIR thường rất phức tạp.
  19. 4. Thiết kế bộ lọc IIR từ các bộ lọc thời gian liên tục (tt)  Biến đổi bộ lọc tương tự thành bộ lọc số: H(s) → H(z)  Các yêu cầu của phép biến đổi: Miền s Miền z s  j z re j Trục tần số   Tính ổn định  0 r 1
  20. 4. Thiết kế bộ lọc IIR từ các bộ lọc thời gian liên tục (tt)  Phép biến đổi song tuyến tính: 1 z 1 s 1 1 z  Kiểm tra các yêu cầu của phép biến đổi: 1 s 1  j z 1 s 1  j  Chứng minh j   e j  cho s j  : 1 j z | z | 1: z e j 1 j 1 e j  Mối liên hệ giữa Ω và ω: s j j tan 1 e j 2  Như vậy:  tan hay  2arctan  2
  21. 4. Thiết kế bộ lọc IIR từ các bộ lọc thời gian liên tục (tt)  Chứng minh  0  r 1 1 s 1  j z 1 s 1  j Dễ thấy: Nếu σ 0: tử số lớn hơn mẫu số →|z|>1
  22. 4. Thiết kế bộ lọc IIR từ các bộ lọc thời gian liên tục (tt) Ví dụ: Cho các tiêu chuẩn thiết kế của một bộ lọc số như sau: 0.89125 H() 1, 0  0.2 H() 0.17783, 0.3  a) Thiết kế bộ lọc tương tự Butterworth thoả các điều kiện trên biết: 2 1 H () 2N  1 c b) Sử dụng phép biến đổi song tuyến tính tìm hàm truyền của bộ lọc IIR tương ứng.
  23. 4. Thiết kế bộ lọc IIR từ các bộ lọc thời gian liên tục (tt)  Nhận xét:  tan 2 ánh xạ trục tần số vô hạn vào vòng tròn đơn vị hữu hạn dẫn đến các tần số được ánh xạ không tuyến tính -> không áp dụng được cho bộ lọc có đáp ứng biên độ hay pha tuyến tính.
  24. 4. Thiết kế bộ lọc IIR từ các bộ lọc thời gian liên tục (tt) Bộ lọc 1 cực: 1 b0 b1z H(z) 1 1 a1z Hàm truyền chuẩn hoá trong miền s: H a (s) s Yêu cầu: thiết kế bộ lọc thông thấp với đáp ứng tần số 2 tại ωc là Gc
  25. 4. Thiết kế bộ lọc IIR từ các bộ lọc thời gian liên tục (tt) Ta có: 2 2 Ha () Ha () 2 2 j   tan( / 2) Tại c c : 2 2 G H ( ) G2 c  a c 2 2 c 2 c c 1 G 2 c Nếu Gc =1/2 thì 훼=Ωc Suy ra hàm truyền trong miền z: 1 z 1 H(z) H (s) 1 z 1 b a s 1 1 1 1 z 1 z s s 1 az 1 z 1 1 Với a , b 1 1
  26. 4. Thiết kế bộ lọc IIR từ các bộ lọc thời gian liên tục (tt) Ví dụ: Thiết kế bộ lọc thông thấp 1 cực với tần số cắt 3dB là 0.2 sử dụng phép biến đổi song tuyến tính áp dụng cho bộ lọc tương tự sau: H(s) s Giải:  Ωc =tan(0.1 )=0.325 2  Do |Ha(Ωc)| =1/2, 훼= Ωc 0.325  Hàm truyền trong miền s: H(s) s 0.325  Hàm truyền trong miền z: 1 j 1 z 1 e H(z) 0.245 1 H() 0.245 j 1 0.509z 1 0.509e  Tại ω=0.2 , |H(ω)|0.707
  27. 5. Bộ lọc có pha tuyến tính  Người ta chứng minh được rẳng bộ lọc có pha tuyến tính nếu đáp ứng xung của nó thoả điều kiện đối xứng hoặc phản đối xứng: h(n) h(M n)  H(z) z M H(z 1)  Đối với bộ lọc FIR: đáp ứng xung của các bộ lọc chọn lọc tần số và các hàm cửa sổ sẵn có dạng đối xứng -> pha tuyến tính có thể đạt được dễ dàng.  Đối với bộ lọc IIR: Từ phương trình của H(z), các cực xuất hiện thành cặp nghịch đảo. Nếu bộ lọc có 1 cực nằm trong vòng tròn đơn vị thì cũng có 1 cực nằm ngoài → không ổn định
  28. 6. Ưu/khuyết điểm của bộ lọc IIR so với bộ lọc FIR FIR IIR Ưu điểm: Ưu điểm:  Pha tuyến tính  Bậc thấp → chi phí tính toán  Ổn định (không có các cực) thấp.  Dễ thiết kế. Nhược điểm: Nhược điểm:  Chi phí tính toán lớn do để  Dễ trở nên bất ổn nếu quá có cùng một đáp ứng biên độ trình lượng tử tham số làm giống nhau thì bộ lọc FIR cho các cực bị đẩy ra khỏi luôn có bậc lớn hơn vòng tròn đơn vị  Không thể đạt pha tuyến tính trên toàn khoảng Nyquist