Ứng dụng lọc thích nghivới mô hình laguerre trong trệt tiếng ồn tích cực

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA -0O0 - LÊ HOÀNG MINH ỨNG DỤNG LỌC THÍCH NGHI VỚI MÔ HÌNH LAGUERRE TRONG TRIỆT TIẾNG ỒN TÍCH CỰC CHUYÊN NGÀNH: KỸ THUẬT VÔ TUYẾN-ĐIỆN TỬ MÃ SỐ NGÀNH: 2.07.01 LUẬN VĂN THẠC SĨ TP Hồ Chí Minh, tháng 10 năm 2003 BỘ GIÁO DỤC VÀ ĐÀO TẠO ĐẠI HỌC QUỐC GIA TP.HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA - CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM ĐỘC LẬP –TỰ DO – HẠNH PHÚC o0o NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên : LÊ HOÀNG MINH Ngày, tháng, năm sinh : 23 - 10 - 1974 Chuyên ngành: KỸ THUẬT VÔ TUYẾN-ĐIỆN TỬ Khóa :12 I Phái : Nữ Nơi sinh : Nam Định TÊN ĐỀ TÀI : ỨNG DỤNG LỌC THÍCH NGHI VỚI MÔ HÌNH LAGUERRE TRONG TRIỆT TIẾNG ỒN TÍCH CỰC II NHIỆM VỤ VÀ NỘI DUNG: • Tìm hiểu giải thuật lọc thích nghi • Xây dựng giải thuật lọc thích nghi dùng mô hình FIR, Laguerre • Mô hệ thống triệt tiếng ồn truyền thẳng thích nghi dùng mô hình FIR mô hình Laguerre III NGÀY GIAO NHIỆM VU Ï: 20/4/2003 IV NGÀY HOÀN THÀNH NHIỆM VỤ : 25/10/2003 V HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN : TS DƯƠNG HOÀI NGHĨA VI HỌ VÀ TÊN CÁN BỘ NHẬN XÉT 1: TS VII HỌ VÀ TÊN CÁN BỘ NHẬN XÉT 2: TS CÁN BỘ HƯỚNG DẪN CÁN BỘ NHẬN XÉT CÁN BỘ NHẬN XÉT TS DƯƠNG HOÀI NGHĨA TS TS Nội dụng đề cương luận văn thạc só thông qua Hội Đồng Chuyên Ngành Ngày …… tháng …… năm 2003 PHÒNG QLKH - SAU ĐẠI HỌC CHỦ NHIỆM NGÀNH NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN NHẬN XÉT CỦA GIÁO VIÊN PHẢN BIỆN NHẬN XÉT CỦA GIÁO VIÊN PHẢN BIỆN LỜI CẢM ƠN Em xin gởi lời biết ơn sâu sắc đến tất Thầy Cô Trường Đại Học Bách Khoa TP Hồ Chí Minh tận tình giảng dạy em thời gian học tập thực luận văn Cảm ơn Ban Giám Hiệu, Khoa Điện, Bộ Môn Điện tử Trường Đại Học Sư Phạm Kỹ Thuật tạo điều kiện thuận lợi để em học tập suốt thời gian qua Đặc biệt, em xin chân thành cảm ơn Thầy, Tiến Só Dương Hoài Nghóa tận tình hướng dẫn giúp đỡ em trình thực luận văn Xin chân thành cảm ơn Bạn bè , đồng nghiệp, người thân gia đình nguồn động viên , hỗ trợ suốt khóa học Lê Hoaøng Minh ABSTRACT This paper presents an adaptive filter algorithm using Laguerre model and its application to the active noise control problem Simulations show that, with respect to the algorithm using FIR model, the one using Laguerre model converges more quickly, has smaller error, and yields better performance of noise cancellation MỤC LỤC Trang Mở đầu CHƯƠNG 1: TRIỆT TIẾNG ỒN TÍCH CỰC 1.1 Giới thiệu 1.2 Các hệ thống triệt tiếng ồn 1.2.1 Hệ thống truyền thẳng 1.2.2 Hệ thống hồi tiếp 1.2.3 Hệ thống truyền thẳng hồi tiếp 1.3 Xây dựng hệ thống triệt tiếng ồn tích cực 10 CHƯƠNG 2: LỌC THÍCH NGHI 10 Giới thiệu 11 2.2 Lọc tối ưu Wiener 12 2.3 Giải thuật lọc thích nghi giảm nhanh 14 2.4 Giải thuật lọc thích nghi LMS 16 2.5 Giải thuật LMS lọc thích nghi FIR bậc p 17 CHƯƠNG 3: MÔ HÌNH LAGUERRE 18 3.1 Giới thiệu 19 3.2 Mô hình Laguerre 20 3.2.1 Hệ thống tuyến tính liên tục 21 3.2.2 Hệ thống tuyến tính rời rạc 22 3.23 Những biến thể mô hình Laguerre 23 CHƯƠNG 4: ỨNG DỤNG LỌC THÍCH NGHI VỚI MÔ HÌNH LAGUERRE TRONG TRIỆT TIẾNG ỒN TÍCH CỰC 24 4.1 Nguyên lý làm việc 25 41.1 Lọc FIR 26 4.1.2 Lọc Laguerre 27 4.2 Mô hệ thống triệt tiếng ồn thích nghi 4.2.1 Hệ thống triệt tiếng ồn thích nghi mô 28 28 Matlab công cụ Simulink 4.2.2 Mô hình nhận dạng LMS nguồn âm thứ cấp ứng với lọc FIR 29 4.2.3 Mô hình nhận dạng LMS nguồn âm thứ cấp ứng với lọc Laguerre 30 4.3 Kết mô 31 4.3.1 Thông số mô 31 4.3.2 Dữ liệu mô 31 4.4 Nhận dạng dùng mô hình FIR kết mô 32 4.4.1 Hiệu triệt tiếng ồn dùng lọc thích nghi với mô hình FIR 32 4.4.2 Đánh giá sai số mô hình nhận dạng theo bậc mô hình 35 4.4.3 Biểu đồ Bode lọc FIR với bậc khác 36 4.4.4 Biến thiên thông số mô hình FIR với bậc khác trình nhận dạng 4.5 Nhận dạng dùng mô hình Laguerre kết mô 50 50 4.5.1 Hiệu triệt tiếng ồn dùng lọc thích nghi với mô hình Laguerre 51 4.5.2 Đánh giá sai số mô hình nhận dạng theo bậc mô hình 52 4.5.3 Biểu đồ Bode lọc Laguerre với bậc khác 53 4.5.4 Biến thiên thông số mô hình Laguerre với bậc khác trình nhận dạng 56 4.6 nh hưởng hệ số μ 62 4.7 nh hưởng hệ số a 63 KẾT LUẬN 65 PHỤ LỤC 66 TÀI LIỆU THAM KHẢO 67 LUẬN VĂN CAO HỌC th1 = XX(27,:); th2 = XX(28,:); th3 = XX(29,:); th4 = XX(30,:); th5 = XX(31,:); th6 = XX(32,:); th7 = XX(33,:); th8 = XX(34,:); th9 = XX(35,:); th10 = XX(36,:); th11 = XX(37,:); th12 = XX(38,:); th13 = XX(39,:); th14 = XX(40,:); th15 = XX(41,:); th16 = XX(42,:); th17 = XX(43,:); th18 = XX(44,:); th19 = XX(45,:); th20 = XX(46,:); th21 = XX(47,:); myB20(1) = th1h(size(th1h,2)); myB20(2) = th2h(size(th2h,2)); myB20(3) = th3h(size(th3h,2)); myB20(4) = th4h(size(th4h,2)); myB20(5) = th5h(size(th5h,2)); myB20(6) = th6h(size(th6h,2)); myB20(7) = th7h(size(th7h,2)); myB20(8) = th8h(size(th8h,2)); myB20(9) = th9h(size(th9h,2)); myB20(10) = th10h(size(th10h,2)); myB20(11) = th11h(size(th11h,2)); myB20(12) = th12h(size(th12h,2)); myB20(13) = th13h(size(th13h,2)); myB20(14) = th14h(size(th14h,2)); myB20(15) = th15h(size(th15h,2)); myB20(16) = th16h(size(th16h,2)); myB20(17) = th17h(size(th17h,2)); myB20(18) = th18h(size(th18h,2)); Phụ Lục LUẬN VĂN CAO HỌC myB20(19) = th19h(size(th19h,2)); myB20(20) = th20h(size(th20h,2)); myB20(21) = th21h(size(th21h,2)); myA20 = [1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; SYS20 = TF(myB20,myA20,Tsamp); figure(1), bode(SYS20,{0.01,100000},'k'), hold on figure(2) subplot(2,2,1), plot(time,th1,'k'), hold on, plot(time,th1h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th2,'k'), hold on, plot(time,th2h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,3), plot(time,th3,'k'), hold on, plot(time,th3h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,4), plot(time,th4,'k'), hold on, plot(time,th4h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(3) subplot(2,2,1), plot(time,th5,'k'), hold on, plot(time,th5h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th6,'k'), hold on, plot(time,th6h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) subplot(2,2,3), plot(time,th7,'k'), hold on, plot(time,th7h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,4), plot(time,th8,'k'), hold on, plot(time,th8h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(4) subplot(2,2,1), plot(time,th9,'k'), hold on, plot(time,th9h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th10,'k'), hold on, plot(time,th10h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,3), plot(time,th11,'k'), hold on, plot(time,th11h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,4), plot(time,th12,'k'), hold on, plot(time,th12h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(5) subplot(2,2,1), plot(time,th13,'k'), hold on, plot(time,th13h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th14,'k'), hold on, plot(time,th14h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) Phuï Lục LUẬN VĂN CAO HỌC subplot(2,2,3), plot(time,th15,'k'), hold on, plot(time,th15h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,4), plot(time,th16,'k'), hold on, plot(time,th16h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(6) subplot(2,2,1), plot(time,th17,'k'), hold on, plot(time,th17h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th18,'k'), hold on, plot(time,th18h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,3), plot(time,th19,'k'), hold on, plot(time,th19h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,4), plot(time,th20,'k'), hold on, plot(time,th20h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(7) subplot(2,2,1), plot(time,th21,'k'), hold on, plot(time,th21h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) figure(8) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.5,0.5]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load FIR25.mat XX time = XX(1,:); x = XX(2,:); e = XX(3,:); d = XX(4,:); dh = XX(5,:); th1h = XX(6,:); th2h = XX(7,:); th3h = XX(8,:); th4h = XX(9,:); th5h = XX(10,:); th6h = XX(11,:); th7h = XX(12,:); th8h = XX(13,:); th9h = XX(14,:); th10h = XX(15,:); th11h = XX(16,:); th12h = XX(17,:); th13h = XX(18,:); th14h = XX(19,:); Phụ Lục LUẬN VĂN CAO HỌC th15h = XX(20,:); th16h = XX(21,:); th17h = XX(22,:); th18h = XX(23,:); th19h = XX(24,:); th20h = XX(25,:); th21h = XX(26,:); th22h = XX(27,:); th23h = XX(28,:); th24h = XX(29,:); th25h = XX(30,:); th26h = XX(31,:); th1 = XX(32,:); th2 = XX(33,:); th3 = XX(34,:); th4 = XX(35,:); th5 = XX(36,:); th6 = XX(37,:); th7 = XX(38,:); th8 = XX(39,:); th9 = XX(40,:); th10 = XX(41,:); th11 = XX(42,:); th12 = XX(43,:); th13 = XX(44,:); th14 = XX(45,:); th15 = XX(46,:); th16 = XX(47,:); th17 = XX(48,:); th18 = XX(49,:); th19 = XX(50,:); th20 = XX(51,:); th21 = XX(52,:); th22 = XX(53,:); th23 = XX(54,:); th24 = XX(55,:); th25 = XX(56,:); th26 = XX(57,:); myB25(1) = th1h(size(th1h,2)); Phụ Lục LUẬN VĂN CAO HỌC myB25(2) = th2h(size(th2h,2)); myB25(3) = th3h(size(th3h,2)); myB25(4) = th4h(size(th4h,2)); myB25(5) = th5h(size(th5h,2)); myB25(6) = th6h(size(th6h,2)); myB25(7) = th7h(size(th7h,2)); myB25(8) = th8h(size(th8h,2)); myB25(9) = th9h(size(th9h,2)); myB25(10) = th10h(size(th10h,2)); myB25(11) = th11h(size(th11h,2)); myB25(12) = th12h(size(th12h,2)); myB25(13) = th13h(size(th13h,2)); myB25(14) = th14h(size(th14h,2)); myB25(15) = th15h(size(th15h,2)); myB25(16) = th16h(size(th16h,2)); myB25(17) = th17h(size(th17h,2)); myB25(18) = th18h(size(th18h,2)); myB25(19) = th19h(size(th19h,2)); myB25(20) = th20h(size(th20h,2)); myB25(21) = th21h(size(th21h,2)); myB25(22) = th22h(size(th22h,2)); myB25(23) = th23h(size(th23h,2)); myB25(24) = th24h(size(th24h,2)); myB25(25) = th25h(size(th25h,2)); myB25(26) = th26h(size(th26h,2)); myA25 = [1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; SYS25 = TF(myB25,myA25,Tsamp); figure(1), bode(SYS25,{0.01,100000},'k'), hold on figure(2), clf subplot(211), plot(time,e,'k'), hold on, ylabel('e'),hold on, xlabel('t'),axis([0,Tsim,-0.2,0.2]) figure(2) subplot(2,2,1), plot(time,th1,'k'), hold on, plot(time,th1h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th2,'k'), hold on, plot(time,th2h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,3), plot(time,th3,'k'), hold on, plot(time,th3h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) Phụ Lục LUẬN VĂN CAO HOÏC subplot(2,2,4), plot(time,th4,'k'), hold on, plot(time,th4h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(3) subplot(2,2,1), plot(time,th5,'k'), hold on, plot(time,th5h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th6,'k'), hold on, plot(time,th6h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) subplot(2,2,3), plot(time,th7,'k'), hold on, plot(time,th7h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,4), plot(time,th8,'k'), hold on, plot(time,th8h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(4) subplot(2,2,1), plot(time,th9,'k'), hold on, plot(time,th9h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th10,'k'), hold on, plot(time,th10h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,3), plot(time,th11,'k'), hold on, plot(time,th11h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,4), plot(time,th12,'k'), hold on, plot(time,th12h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(5) subplot(2,2,1), plot(time,th13,'k'), hold on, plot(time,th13h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th14,'k'), hold on, plot(time,th14h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,3), plot(time,th15,'k'), hold on, plot(time,th15h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,4), plot(time,th16,'k'), hold on, plot(time,th16h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(6) subplot(2,2,1), plot(time,th17,'k'), hold on, plot(time,th17h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th18,'k'), hold on, plot(time,th18h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,3), plot(time,th19,'k'), hold on, plot(time,th19h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,4), plot(time,th20,'k'), hold on, plot(time,th20h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(7) Phuï Luïc LUẬN VĂN CAO HỌC subplot(2,2,1), plot(time,th21,'k'), hold on, plot(time,th21h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th22,'k'), hold on, plot(time,th22h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,3), plot(time,th23,'k'), hold on, plot(time,th23h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) subplot(2,2,4), plot(time,th24,'k'), hold on, plot(time,th24h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(8) subplot(2,2,1), plot(time,th25,'k'), hold on, plot(time,th25h,'k'), hold on,axis([0,Tsim,-0.02,0.1]) subplot(2,2,2), plot(time,th26,'k'), hold on, plot(time,th26h,'k'), hold on, axis([0,Tsim,-0.02,0.1]) figure(9) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.5,0.5]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% SYS0 = TF(myB,myA,Tsamp) figure(1), bode(SYS0,{0.01,100000},'k'), hold on, ,% bode cua ham truyen moi truong %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Order = [3 ,5 ,9 ,12 ,15 ,20 ,25 ]; JJ = [0.01573 , 0.01332 ,0.004658 ,0.002122 , 0.001434 , 0.001244, 0.001221 ]; figure(3), plot(Order,JJ,'k'), axis([0,25,0,0.02]), hold on,% bac bo loc FIR %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load LAGUERRE3.mat XX time = XX(1,:); x = XX(2,:); e = XX(3,:); d = XX(4,:); dh = XX(5,:); th1h = XX(6,:); th2h = XX(7,:); th3h = XX(8,:); myB3(1) = th1h(size(th1h,2)); myB3(2) = th2h(size(th2h,2)); myB3(3) = th3h(size(th3h,2)); myA3 = [1,0,0]; SYS3 = TF(myB3,myA3,Tsamp); Phụ Lục LUẬN VĂN CAO HỌC figure(1), bode(SYS3,{0.01,100000},'k'), hold on figure(2) subplot(2,2,1), plot(time,th1h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,th2h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,3), plot(time,th3h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) figure(3) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.5,0.5]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load LAGUERRE4.mat XX time = XX(1,:); x = XX(2,:); e = XX(3,:); d = XX(4,:); dh = XX(5,:); th1h = XX(6,:); th2h = XX(7,:); th3h = XX(8,:); th4h = XX(9,:); myB4(1) = th1h(size(th1h,2)); myB4(2) = th2h(size(th2h,2)); myB4(3) = th3h(size(th3h,2)); myB4(4) = th4h(size(th4h,2)); myA4 = [1,0,0,0]; SYS4 = TF(myB4,myA4,Tsamp); figure(1), bode(SYS4,{0.01,100000},'k'), hold on figure(2) subplot(2,2,1), plot(time,th1h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,th2h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,3), plot(time,th3h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,4), plot(time,th4h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) figure(3) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.5,0.5]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load LAGUERRE5.mat XX time = XX(1,:); x = XX(2,:); e = XX(3,:); d = XX(4,:); dh = XX(5,:); Phụ Lục LUẬN VĂN CAO HOÏC th1h = XX(6,:); th2h = XX(7,:); th3h = XX(8,:); th4h = XX(9,:); th5h = XX(10,:); myB5(1) = th1h(size(th1h,2)); myB5(2) = th2h(size(th2h,2)); myB5(3) = th3h(size(th3h,2)); myB5(4) = th4h(size(th4h,2)); myB5(5) = th1h(size(th5h,2)); myA5 = [1,0,0,0,0]; SYS5 = TF(myB5,myA5,Tsamp); figure(1), bode(SYS5,{0.01,100000},'k'), hold on figure(2) subplot(2,2,1), plot(time,th1h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,th2h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,3), plot(time,th3h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,4), plot(time,th4h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) figure(3) subplot(2,2,1), plot(time,th5h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) figure(4) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.5,0.5]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load LAGUERRE8.mat XX time = XX(1,:); x = XX(2,:); e = XX(3,:); d = XX(4,:); dh = XX(5,:); th1h = XX(6,:); th2h = XX(7,:); th3h = XX(8,:); th4h = XX(9,:); th5h = XX(10,:); th6h = XX(11,:); th7h = XX(12,:); th8h = XX(13,:); myB8(1) = th1h(size(th1h,2)); myB8(2) = th2h(size(th2h,2)); Phụ Lục LUẬN VĂN CAO HỌC myB8(3) = th3h(size(th3h,2)); myB8(4) = th4h(size(th4h,2)); myB8(5) = th1h(size(th5h,2)); myB8(6) = th2h(size(th6h,2)); myB8(7) = th3h(size(th7h,2)); myB8(8) = th4h(size(th8h,2)); myA8 = [1,0,0,0,0,0,0,0]; SYS8 = TF(myB8,myA8,Tsamp); figure(1), bode(SYS8,{0.01,100000},'k'), hold on figure(2) subplot(2,2,1), plot(time,th1h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,th2h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,3), plot(time,th3h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,4), plot(time,th4h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) figure(3) subplot(2,2,1), plot(time,th5h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,th6h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,3),plot(time,th7h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,4), plot(time,th8h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) figure(4) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.5,0.5]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load LAGUERRE10.mat XX time = XX(1,:); x = XX(2,:); e = XX(3,:); d = XX(4,:); dh = XX(5,:); th1h = XX(6,:); th2h = XX(7,:); th3h = XX(8,:); th4h = XX(9,:); th5h = XX(10,:); th6h = XX(11,:); th7h = XX(12,:); th8h = XX(13,:); th9h = XX(14,:); th10h = XX(15,:); myB10(1) = th1h(size(th1h,2)); Phụ Lục LUẬN VĂN CAO HOÏC myB10(2) = th2h(size(th2h,2)); myB10(3) = th3h(size(th3h,2)); myB10(4) = th4h(size(th4h,2)); myB10(5) = th1h(size(th5h,2)); myB10(6) = th2h(size(th6h,2)); myB10(7) = th3h(size(th7h,2)); myB10(8) = th4h(size(th8h,2)); myB10(9) = th1h(size(th9h,2)); myB10(10) = th2h(size(th10h,2)); myA10 = [1,0,0,0,0,0,0,0,0,0]; SYS10 = TF(myB10,myA10,Tsamp); figure(1), bode(SYS10,{0.01,100000},'k'), hold on figure(2) subplot(2,2,1), plot(time,th1h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,th2h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,3), plot(time,th3h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,4), plot(time,th4h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) figure(3) subplot(2,2,1), plot(time,th5h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,th6h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,3), plot(time,th7h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,4), plot(time,th8h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) figure(4) subplot(2,2,1), plot(time,th9h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,th10h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) figure(5) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.5,0.5]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load LAGUERRE5.mat XX time = XX(1,:); x = XX(2,:); e = XX(3,:); d = XX(4,:); dh = XX(5,:); th1h = XX(6,:); th2h = XX(7,:); th3h = XX(8,:); th4h = XX(9,:); th5h = XX(10,:); Phụ Lục LUẬN VĂN CAO HỌC th1 = XX(12,:); th2 = XX(13,:); th3 = XX(14,:); th4 = XX(15,:); th5 = XX(16,:); myB5(1) = th1h(size(th1h,2)); myB5(2) = th2h(size(th2h,2)); myB5(3) = th3h(size(th3h,2)); myB5(4) = th4h(size(th4h,2)); myB5(5) = th5h(size(th5h,2)); myA5 = [1,0,0,0,0]; SYS5 = TF(myB5,myA5,Tsamp); figure(1), bode(SYS5,{0.01,100000},'k'), hold on figure(2) subplot(2,2,1), plot(time,0.0257,'k'), hold on, plot(time,th1h,'k'),hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,0.1065,'k'), hold on, plot(time,th2h,'k'), hold on,axis([0,Tsim,-0.05,0.2]) subplot(2,2,3), plot(time,0.1355,'k'), hold on, plot(time,th3h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,4), plot(time,0.0434,'k'), hold on, plot(time,th4h,'k'), hold on,axis([0,Tsim,-0.05,0.2]) figure(3) subplot(2,2,1), plot(time,0.0011,'k'), hold on, plot(time,th5h,'k'),hold on, axis([0,Tsim,-0.05,0.2]) figure(4) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.5,0.5]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load LAGUERRE4.mat XX time = XX(1,:); x = XX(2,:); e = XX(3,:); d = XX(4,:); dh = XX(5,:); th1h = XX(6,:); th2h = XX(7,:); th3h = XX(8,:); th4h = XX(9,:); th1 = XX(10,:); Phụ Lục LUẬN VĂN CAO HỌC th2 = XX(11,:); th3 = XX(12,:); th4 = XX(13,:); myB4(1) = th1h(size(th1h,2)); myB4(2) = th2h(size(th2h,2)); myB4(3) = th3h(size(th3h,2)); myB4(4) = th4h(size(th4h,2)); myA4 = [1,0,0,0]; SYS4 = TF(myB4,myA4,Tsamp); figure(1), bode(SYS4,{0.01,100000},'k'), hold on figure(2) subplot(2,2,1), plot(time,0.0257,'k'), hold on, plot(time,th1h,'k'),hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,0.1065,'k'), hold on, plot(time,th2h,'k'), hold on,axis([0,Tsim,-0.05,0.2]) subplot(2,2,3), plot(time,0.1355,'k'), hold on, plot(time,th3h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,4), plot(time,0.0434,'k'), hold on, plot(time,th4h,'k'), hold on,axis([0,Tsim,-0.05,0.2]) figure(3) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.5,0.5]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load LAGUERRE3.mat XX time = XX(1,:); x = XX(2,:); e = XX(3,:); d = XX(4,:); dh = XX(5,:); th1h = XX(6,:); th2h = XX(7,:); th3h = XX(8,:); th1 = XX(9,:); th2 = XX(10,:); th3 = XX(11,:); myB3(1) = th1h(size(th1h,2)); myB3(2) = th2h(size(th2h,2)); myB3(3) = th3h(size(th3h,2)); myA3 = [1,0,0]; SYS3 = TF(myB3,myA3,Tsamp); Phụ Lục LUẬN VĂN CAO HỌC figure(1), bode(SYS3,{0.01,100000},'k'), hold on figure(2) subplot(2,2,1), plot(time,0.0257,'k'), hold on, plot(time,th1h,'k'),hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,0.1065,'k'), hold on, plot(time,th2h,'k'), hold on,axis([0,Tsim,-0.05,0.2]) subplot(2,2,3), plot(time,0.1355,'k'), hold on, plot(time,th3h,'k'), hold on, axis([0,Tsim,-0.05,0.2]) figure(3) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.2,0.2]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load LAGUERRE2.mat XX time = XX(1,:); x = XX(2,:); e = XX(3,:); d = XX(4,:); dh = XX(5,:); th1h = XX(6,:); th2h = XX(7,:); th1 = XX(8,:); th2 = XX(9,:); myB2(1) = th1h(size(th1h,2)); myB2(2) = th2h(size(th2h,2)); myA2 = [1,0]; SYS2 = TF(myB2,myA2,Tsamp); figure(1), bode(SYS2,{0.01,100000},'k'), hold on figure(2) subplot(2,2,1), plot(time,0.0257,'k'), hold on, plot(time,th1h,'k'),hold on, axis([0,Tsim,-0.05,0.2]) subplot(2,2,2), plot(time,0.1065,'k'), hold on, plot(time,th2h,'k'), hold on,axis([0,Tsim,-0.05,0.2]) figure(3) subplot(211), plot(time,e,'k'),hold on, axis([0,Tsim,-0.5,0.5]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% SYS0 = TF(myB,myA,Tsamp) figure(1), bode(SYS0,{0.01,100000},'k'), hold on, ,% bode cua ham truyen moi truong %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Order = [1 ,2 ,3 ,4 ,5 ]; Phụ Lục LUẬN VĂN CAO HOÏC JJ = [0.05349 ,0.005351 , 0.00005678 , 0.0002069 ,0.000195 ]; figure(1), plot(Order,JJ,'k'), hold on, axis([0,5,0,0.02]), ,% bac bo loc LAGUERRE a=0.7 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Order = [1 ,2 ,3 ,4 ,5 ]; JJ = [0.06177 ,0.02513 , 0.0000271 , 0.0001077, 0.0001192 ]; figure(1), plot(Order,JJ,'k'), hold on, axis([0,5,0,0.02]), ,% bac bo loc LAGUERRE a=0.6 Phuï Luïc ... giải thuật lọc thích nghi dẫn đến giải thuật lọc thích nghi LMS Chương : Lọc Laguerre Tìm hiểu mô hình lọc Laguerre Chương : Ứng dụng lọc thích nghi với mô hình Laguerre triệt tiếng ồn tích cực... thống Chương 3: Mô hình Laguerre 22 LUẬN VĂN CAO HỌC CHƯƠNG IV: ỨNG DỤNG LỌC THÍCH NGHI VỚI MÔ HÌNH LAGUERRE TRONG TRIỆT TIẾNG ỒN TÍCH CỰC Chương 4: Ứng dụng lọc thích nghi triệt tiếng ồn 24 LUẬN... [8],[9] Hình 4.1: Hệ thống triệt tiếng ồn tích cực Chương 4: Ứng dụng lọc thích nghi triệt tiếng ồn 25 LUẬN VĂN CAO HỌC Để thực lọc Fˆ (z-1) ta sử dụng lọc FIR lọc Laguerre 41.1 Lọc FIR Hình 4.2: Mô