Đang tải... (xem toàn văn)
Với xu hướng số hóa các hệ thống thông tin hiện nay,việc xử lý tín hiệu số ngày càng trở nên quan trọng với khả năng xử lý thông tin một các ưu việt. Để có thể tiếp cận được lĩnh vực này, chúng ta cần có những kiến thức cơ bản về tín hiệu số và các phương pháp xử lý. Một trong những kiến thức quan trọng đó là thiết kế bộ lọc số hệ thống có thể làm thay đổi tín hiệu để phù hợp với mục đích của con người. Trong xử lý số tín hiệu, tồn tại nhiều bộ lọc số khác nhau như: bộ lọc thôngthấp, bộ lọc thông dải, bộ lọc vi phân,… Để thiết kế các bộ lọc thích hợp, trước hết phải xác định yêu cầu thực tế dựa trên các chỉ tiêu kỹ thuật cho trước, trên cơ sở đó định hình cấu trúc bộ lọc và phương pháp thiết kế tối ưu. Cấu trúc bộ lọc có thể là:cấu trúc FIR (bộ lọc số có đáp ứng xung chiều dài xác định) hoặc cấu trúc IIR (bộ lọc số có đáp ứng xung chiều dài không xác định). Ph ương pháp thiết kế có thể là: phương pháp cửa sổ, phương pháp lấy mẫu tần số, hoặc phương pháp xấp xỉ tối ưu,… Được sự phân công của thầy giáo , trên cơ sở những kiến thức đã học, tôi đã tìm hiểu bộ lọc FIR theo phương pháp lấy mẫu tần số. Tôi xin chân thành cám ơn thầy giáo , bạn bè cùng lớp đã tận tình hướng dẫn để tôi có thể hoàn thành đề tài này. Chắc chắn đề tài sẽ không tránh khỏi những thiếu sót rất mong được sự góp ý của quý thầy cô và các bạn.
ĐẠI HỌC QUY NHƠN oOo BÀI TẬP LỚN MÔN HỌC: XỬ LÝ SỐ TÍN HIỆU ĐỀ TÀI: THIẾT KẾ BỘ LỌC FIR THÔNG CAO BẰNG PHƯƠNG PHÁP LẤY MẪU TẦN SỐ Người hướng dẫn : Sinh viên thực hiện : Lớp : Quy Nhơn , tháng 5/2011 Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 2 LỜI NÓI ĐẦU Với xu hướng số hóa các hệ thống thông tin hiện nay,việc xử lý tín hiệu số ngày càng trở nên quan trọng với khả năng xử lý thông tin một các ưu việt. Để có thể tiếp cận được lĩnh vực này, chúng ta cần có những kiến thức cơ bản về tín hiệu số và các phương pháp xử lý. Một trong những kiến thức quan trọng đó là thiết kế bộ lọc số - hệ thống có thể làm thay đổi tín hiệu để phù hợp với mục đích của con người. Trong xử lý số tín hiệu, tồn tại nhiều bộ lọc số khác nhau như: bộ lọc thông thấp, bộ lọc thông dải, bộ lọc vi phân,… Để thiết kế các bộ lọc thích hợp, trước hết phải xác định yêu cầu thực tế dựa trên các chỉ tiêu kỹ thuật cho trước, trên cơ sở đó định hình cấu trúc bộ lọc và phương pháp thiết kế tối ưu. Cấu trúc bộ lọc có thể là: cấu trúc FIR (bộ lọc số có đáp ứng xung chiều dài xác định) hoặc cấu trúc IIR (bộ lọc số có đáp ứng xung chiều dài không xác định). Phương pháp thiết kế có thể là: phương pháp cửa sổ, phương pháp lấy mẫu tần số, hoặc phương pháp xấp xỉ tối ưu,… Được sự phân công của thầy giáo , trên cơ sở những kiến thức đã học, tôi đã tìm hiểu bộ lọc FIR theo phương pháp lấy mẫu tần số. Tôi xin chân thành cám ơn thầy giáo , bạn bè cùng lớp đã tận tình hướng dẫn để tôi có thể hoàn thành đề tài này. Chắc chắn đề tài sẽ không tránh khỏi những thiếu sót rất mong được sự góp ý của quý thầy cô và các bạn. Xin chân thành cảm ơn! Quy nhơn, tháng 5 năm 2011 Người thực hiện Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 3 MỤC LỤC LỜI NÓI ĐẦU 2 MỤC LỤC 3 Phần 1. CƠ SỞ LÝ THUYẾT 4 1.1. Dẫn nhập 4 1.2. Cấu trúc của bộ lọc FIR 6 a. Cấu trúc dạng trực tiếp 6 b. Cấu trúc dạng ghép tầng: 7 c. Cấu trúc dạng pha tuyến tính: 7 1.3. Các đặc tính của bộ lọc FIR pha tuyến tính 8 a. Đáp ứng xung h(n): 9 b. Đáp ứng tần số H(e j ): 11 1.4 Phương pháp thiết kế lấy mẫu tần số : 14 a. Phương pháp thiết kế đơn giản 15 b. Phương pháp thiết kế tối ưu 15 Phần 2. THIẾT KẾ LỌC FIR THÔNG CAO 16 2.1. Bài toán thiết kế 16 2.2. Phương pháp thiết kế 16 2.3. Thuật toán và chương trình Matlab 17 a. Lưu đồ thuật toán: 17 b) Chương trình 19 c) Kết quả 21 TÀI LIỆU THAM KHẢO 22 Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 4 Phần 1. CƠ SỞ LÝ THUYẾT 1.1. Dẫn nhập Việc thiết kế một bộ lọc số tiến hành theo 3 bước: Đưa ra các chỉ tiêu: Trước khi thiết kế một bộ lọc chúng ta cần xác định các chỉ tiêu. Các chỉ tiêu được xác định bởi các ứng dụng cụ thể khác nhau. Tìm các xấp xỉ: Một khi chỉ tiêu đã được xác định, ta sử dụng các khái niệm và công cụ toán học khác nhau để tiến tới biểu diễn và tính gần đúng cho bộ lọc với tập các chỉ tiêu đã cho. Và đây là chủ đề chính của việc thiết kế lọc số. Thực hiện bộ lọc: Kết quả của các bước trên cho mô tả dưới dạng một phương trình sai phân, hoặc một hàm hệ thống H(z), hoặc một đáp ứng xung h(n). Từ các mô tả này chúng ta có thể thi hành bộ lọc bằng phần cứng hoặc phần mềm mô phỏng trên máy tính. Trong những phần tiếp sau đây, để mô tả các chỉ tiêu của bộ lọc chúng ta sẽ xem xét việc thiết kế một bộ lọc thông thấp như một ví dụ cơ sở. Có 2 nhóm chỉ tiêu: Các chỉ tiêu tuyệt đối (Absolute Specifications) và các chỉ tiêu tương đối (Relative Specifications - DB). Trên hình 1.1 là mô tả các chỉ tiêu của bộ lọc FIR thông thấp (Low Pass Filter): 1 1+ 1 1- 1 Độ gợn dải thông Độ gợn dải chắn Dải chuyển tiếp 2 R p A s w p w s 0 0 |H(e jw )| Decibels Hình 1.1 Các chỉ tiêu của bộ lọc FIR: các chỉ tiêu tuyệt đối và tương đối Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 5 Trong đó: Band [0, w p ] được gọi là dải thông, và 1 là dung sai (gợn sóng) được chấp nhận trong đáp ứng dải thông lý tưởng. Band [w s , ] được gọi là dải chắn, và 2 là dung sai ở dải chắn. Band [w p , w s ] được gọi là dải chuyển tiếp, và không có ràng buộc nào về đáp ứng biên độ trong dải này Các chỉ tiêu tương đối gồm có: R p : Độ gợn sóng trong dải thông tính theo dB. A s : Suy hao trong dải chắn tính theo dB. Quan hệ giữa các chỉ tiêu này như sau: 0 1 1 1 1 10 log20 p R (0) (1.1) 0 1 1 2 10 log20 s A (>>1) (1.2) Các chỉ tiêu trên được đưa ra đối với bộ lọc FIR thông thấp, và tất nhiên đối với các bộ lọc khác như thông cao HPF (High Pass Filter), thông dải BPF (Band Pass Filter) đều có thể được định nghĩa tương tự. Tuy nhiên, các tham số thiết kế quan trọng nhất là các dung sai dải tần và các tần số cạnh-dải (tolerance or ripples and band-edge frequencies). Bởi vậy, trong phần 1 về cơ sở lý thuyết này chúng ta chỉ tập trung vào bộ lọc FIR thông thấp. Việc thiết kế cụ thể cho bộ lọc FIR thông dải bằng kỹ thuật cửa sổ sẽ được phát triển trên cơ sở lọc thông thấp và sẽ được mô tả chi tiết trong phần 2. Việc thiết kế và thực hiện lọc FIR có những thuận lợi sau đây: Đáp ứng pha là tuyến tính. Dễ thiết kế do không gặp các vấn đề ổn định (lọc FIR luôn ổn định). Việc thực hiện rất hiệu quả. Có thể sử dụng DFT để thực hiện Đáp ứng pha là tuyến tính (linear phase response) mang lại những thuận lợi sau: Bài toán thiết kế chỉ gồm các phép tính số học thực chứ không cần phép tính số học phức. Bộ lọc pha tuyến tính không có méo trễ nhóm và chỉ bị trễ một khoảng không đổi. Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 6 Đối với bộ lọc có chiều dài M (hoặc bậc M-1) số phép toán có bậc M/2 như đã khảo sát trong thực hiện bộ lọc có pha tuyến tính. 1.2. Cấu trúc của bộ lọc FIR Một bộ lọc đáp ứng xung hữu hạn với hàm hệ thống có dạng: 1M 0n n n M1 1M 1 10 zbzbzbb)z(H (1.3) Như vậy đáp ứng xung h(n) là: else Mnb nh n 0 10 )( (1.4) Và phương trình sai phân là: )1Mn(xb)1n(xb)n(xb)n(y 1M10 (1.5) Đây chính là tích chập tuyến tính của các dãy hữu hạn. Bậc của bộ lọc là M-1, trong khi chiều dài của bộ lọc là M (bằng với số lượng các hệ số). Các cấu trúc bộ lọc FIR luôn luôn ổn định, và tương đối đơn giản hơn so với các cấu trúc bộ lọc IIR. Hơn thế nữa, các bộ lọc FIR có thể được thiết kế để có một đáp ứng pha tuyến tính và đó là điều cần thiết trong một số ứng dụng. Chúng ta sẽ xem xét lần lượt các cấu trúc của bộ lọc FIR sau đây: a. Cấu trúc dạng trực tiếp Phương trình sai phân được thực hiện bởi một dãy liên tiếp các bộ trễ do không có đường phản hồi: )1Mn(xb)1n(xb)n(xb)n(y 1M10 (1.6) Do mẫu thức bằng đơn vị nên ta chỉ có một cấu trúc dạng trực tiếp duy nhất. Cấu trúc dạng trực tiếp được cho trong hình 1.2 với M = 5: b 0 z - 1 b 1 z - 1 b 2 z - 1 b 3 z - 1 b 4 y(n) x(n) Hình 1.2 Cấu trúc lọc FIR dạng trực tiếp Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 7 b. Cấu trúc dạng ghép tầng: Hàm hệ thống H(z) được biến đổi thành các tích của các khâu bậc 2 với các hệ số thực. Các khâu này được thực hiện ở dạng trực tiếp và bộ lọc tổng thể có dạng ghép tầng của các khâu bậc 2. M1 0 1M 1 0 1 0 M1 1M 1 10 z b b z b b 1bzbzbb)z(H K 1k 2 2,k 1 1,k0 )zBzB1(b (1.7) trong đó 2 M K , B k,1 và B k,2 là các số thực đại diện cho các hệ số của các khâu bậc 2. Cấu trúc dạng ghép tầng được cho trong hình 1.3 với M = 7: c. Cấu trúc dạng pha tuyến tính: Đối với các bộ lọc chọn tần, người ta mong muốn có đáp ứng pha là hàm tuyến tính theo tần số, nghĩa là: )e(H j (1.8) trong đó 0 hoặc 2 và là một hằng số. Đối với bộ lọc FIR nhân quả có đáp ứng xung trong khoảng [0, M-1], thì các điều kiện tuyến tính là: 1Mn0,0);n1M(h)n(h (1.9) 1Mn0,2/);n1M(h)n(h (1.10) Xét phương trình sai phân được cho trong phương trình (1.5) với đáp ứng xung đối xứng trong phương trình (1.9), ta có: )1Mn(xb)2Mn(xb)1n(xb)n(xb)n(y 0110 )]2Mn(x)1n(x[b)]1Mn(x)n(x[b 10 B 1,1 z - 1 z - 1 z - 1 y(n) x(n) B 2,1 B 3,1 b 0 B 1,2 z - 1 z - 1 z - 1 B 2,2 B 3,2 Hình 1.3 Cấu trúc lọc FIR dạng ghép tầng Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 8 Sơ đồ khối thực hiện phương trình sai phân trên được mô tả trong hình 1.4 dưới đây đối với cả M lẻ và M chẵn: Đối với M lẻ: M = 7, còn đối với M chẵn: M = 6 Rõ ràng, với cùng một bậc của bộ lọc (cùng M) cấu trúc pha tuyến tính sẽ tiết kiệm được 50% các bộ nhân so với cấu trúc dạng trực tiếp. 1.3. Các đặc tính của bộ lọc FIR pha tuyến tính Trong phần này chúng ta sẽ thảo luận về hình dạng của đáp ứng xung, đáp ứng tần số trong hàm hệ thống của các bộ lọc FIR pha tuyến tính. Cho h(n), trong đó 0 n M – 1, là đáp ứng xung có chiều dài M thì hàm truyền hệ thống là: 1M 0n n1M)1M( 1M 0n n z)n(hzz)n(h)z(H (1.11) có (M-1) điểm cực ở gốc (trivial poles) và M-1 điểm không nằm ở vị trí bất kỳ trên mặt phẳng z. Đáp ứng tần số là: ,e)n(h)e(H 1M 0n njj (1.12) b 0 z - 1 b 1 z - 1 b 2 x(n) z - 1 z - 1 z - 1 y(n) b 0 z - 1 b 1 z - 1 b 2 b 3 y(n) x(n) z - 1 z - 1 z - 1 z - 1 M=7 M=6 Hình 1.4 Cấu trúc lọc FIR pha tuyến tính với các hệ số M chẵn và lẻ Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 9 a. Đáp ứng xung h(n): Chúng ta có thể đưa ra ràng buộc pha tuyến tính: ,)e(H j (1.13) trong đó: là một hằng số trễ pha. Ta đã biết rằng h(n) phải đối xứng, nghĩa là: 2 1M ,1Mn0),n1M(h)n(h (1.14) Do đó h(n) là đối xứng theo , là chỉ số đối xứng. Có hai kiểu đối xứng: M lẻ: Trong trường hợp này, 2 1M là một số nguyên. Đáp ứng xung được mô tả trong hình 1.5 dưới đây: M chẵn: Trong trường hợp này, 2 1M không phải là một số nguyên. Đáp ứng xung được mô tả bằng hình 1.6 dưới đây: Hình 1.5 Đáp ứng xung đối xứng, M lẻ Hình 1.6 Đáp ứng xung đối xứng, M chẵn Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 10 Ta cũng có bộ lọc FIR pha tuyến tính loại hai nếu ta yêu cầu đáp ứng pha j eH thoả mãn điều kiện: )e(H j với (1.15) Đáp ứng pha là đường thẳng nhưng không đi qua gốc. Trong trường hợp này không phải là hằng số trễ pha, nhưng: d )e(Hd j (1.16) là hằng số, chính là trễ nhóm ( là một hằng số trễ nhóm). Trong trường hợp này, các tần số được làm trễ với một tốc độ không đổi. Nhưng một số tần số có thể được làm trễ với tốc độ lớn hơn hoặc nhỏ hơn. Đối với kiểu pha tuyến tính này, có thể thấy rằng: 1Mn0),n1M(h)n(h và 2 , 2 1M (1.17) Điều này có nghĩa rằng đáp ứng xung h(n) là phản đối xứng (antisymmetric). Chỉ số đối xứng vẫn là 2 1M . Một lần nữa chúng ta lại có 2 kiểu, cho M lẻ và M chẵn. M lẻ: Trong trường hợp này, 2 1M là một số nguyên. Đáp ứng xung được mô tả bằng hình 1.7 dưới đây: Lưu ý rằng mẫu h() tại 2 1M phải bằng 0, nghĩa là, 0 2 1M h . Hình 1.7 Đáp ứng xung phản đối xứng, M lẻ [...]... số Trang 21 c) Kết quả Kết quả chương trình được thực hiện với ws=0.4pi, wp=0.65pi, Rp=0.4, As=45 Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 22 TÀI LIỆU THAM KHẢO - [1] Ngô Văn Sỹ, giảng Xử lý Tín hiệu số, Đại học Bách Khoa Đà Nẵng [2] Monson H Hayes (Sách dịch - Tống Văn On), Lý thuyết và Bài tập Xử lý Tín hiệu số, Nhà xuất bản Lao động xã hội [3] Quách Tuấn Ngọc, Xử lý Tín. .. lý Tín hiệu số, Đại học Bách Khoa Hà Nội [4] Nguyễn Quốc Trung (2002), Xử lý Tín hiệu và Lọc số, Nhà xuất bản Khoa học và Kỹ thuật [5] Nguyễn Phùng Quang, Matlab & Simulink dành cho kỹ sư điều khiển tự động, NXB Khoa Học Kỹ Thuật, 2004 [6] Vinay K Ingle & John G Proalis, Digital Signal Procesing Using MATLAB, Brooks/Cole Thomson Learning [7] Hồ Văn Sung, Thực hành Xử Lý Số Tín Hiệu Trên Máy Tính PC... kỹ thuật lấy mẫu tần số : 1 0 1 2 1 3 4 5 6 7 0 1 2 4 3 5 6 7 8 Từ hình trên ta nhận thấy : Lỗi xấp xỉ là hiệu của đáp ứng lý tưởng và đáp ứng thực tế bằng không tại các tần số được lấy mẫu Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 15 Lỗi xấp xỉ ở tất cả các tần số khác phụ thuộc vào hình dạng của đáp ứng tần số lý tưởng; nghĩa là, đáp ứng tần số lý tưởng càng sắc nét... viết và chạy trên nền Matlab, với việc sử dụng một số hàm hỗ trợ có sẵn của Matlab cho xử lý tín hiệu số, và một số hàm viết thêm được tham khảo từ tài liệu [1] (các hàm dưới dạng các file m) a Lưu đồ thuật toán: Thiết kế lọc FIR thông cao bằng phương pháp lấy mẫu tần số Trang 18 BEGIN Nhập các chỉ tiêu ws, wp As, Rp No Chỉ tiêu có hợp lệ không? Yes Tính số mẫu M Tìm gía trị của T1, T2 tốt nhất ( ứng với... (high pass filter) Điều kiện: ws< wp . ĐẠI HỌC QUY NHƠN oOo BÀI TẬP LỚN MÔN HỌC: XỬ LÝ SỐ TÍN HIỆU ĐỀ TÀI: THIẾT KẾ BỘ LỌC FIR THÔNG CAO BẰNG PHƯƠNG PHÁP LẤY MẪU TẦN SỐ Người. phương pháp lấy mẫu tần số Trang 2 LỜI NÓI ĐẦU Với xu hướng số hóa các hệ thống thông tin hiện nay,việc xử lý tín hiệu số ngày càng trở nên quan trọng với khả năng xử lý thông tin một các ưu. tín hiệu số và các phương pháp xử lý. Một trong những kiến thức quan trọng đó là thiết kế bộ lọc số - hệ thống có thể làm thay đổi tín hiệu để phù hợp với mục đích của con người. Trong xử