Nghiên cứu giải pháp điều khiển bị động kết cấu với hệ cản thay đổi khối lượng mtmd (multiple tuned mass damper)

Đại Học Quốc Gia Tp Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA PHẠM HỮU NGHĨA NGHIÊN CỨU GIẢI PHÁP ĐIỀU KHIỂN BỊ ĐỘNG KẾT CẤU VỚI HỆ CẢN THAY ĐỔI KHỐI LƯNG MTMD (MULTIPLE TUNED MASS DAMPER) Chuyên ngành: XÂY DỰNG DÂN DỤNG & CÔNG NGHIỆP Mã số ngành: 23 04 10 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, tháng 11 năm 2007 CÔNG TRÌNH ĐƯC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học: PGS TS BÙI CÔNG THÀNH Cán chấm nhận xeùt 1: Cán chấm nhận xét 2: Luận văn thạc só bảo vệ HỘI ĐỒNG CHẤM LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA, ngày tháng năm CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐỘC LẬP – TỰ DO – HẠNH PHÚC - - - oOo - - - - - - oOo - - Tp HCM, ngày 05 tháng 11 năm 2007 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: PHẠM HỮU NGHĨA Giới tính: Nam Ngày, tháng, năm sinh: 12 – 08 - 1977 Nơi sinh: Lâm Đồng Chuyên ngành: XÂY DỰNG DÂN DỤNG VÀ CÔNG NGHIỆP Khố (Năm trúng tuyển) : 2004 I TÊN ĐỀ TÀI : NGHIÊN CỨU GIẢI PHÁP ĐIỀU KHIỂN BỊ ĐỘNG KẾT CẤU VỚI HỆ CẢN THAY ĐỔI KHỐI LƯNG MTMD (MULTIPLE TUNED MASS DAMPER) II NHIỆM VỤ LUẬN VĂN: • Nghiên cứu làm việc nhiều hệ cản MTMD điều khiển dao động kết cấu • Xét thông số hệ cản TMD hệ MTMD, mô hình hệ cản MTMD vị trí bố trí hệ cản TMD hệ MTMD • Sử dụng phương pháp Newmark để phân tích dao động kết cấu • Xây dựng chương trình tính toán với ngôn ngữ kỹ thuật Matlab • Ví dụ kiểm chứng minh hoạ • Kết luận kiến nghị III NGÀY GIAO NHIỆM VỤ: 05 /02/ 2007 IV NGÀY HOÀN THÀNH NHIỆM VỤ: 05 /11 / 2007 V CÁN BỘ HƯỚNG DẪN: PGS TS BÙI CÔNG THÀNH Nội dung đề cương luận văn thạc só Hội đồng chuyên ngành thông qua CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM BỘ MÔN QUẢN LÝ CHUYÊN NGÀNH PGS TS BÙI CÔNG THÀNH LỜI CẢM ƠN Sau ngày tháng cố gắng học hỏi nghiên cứu, cuối hoàn thành luận văn kết thúc khóa học thạc só trường Đại Học Bách Khoa TP.Hồ Chí Minh Tôi xin chân thành cảm ơn Thầy Bùi Công Thành người thầy bảo hướng dẫn thực luận văn này.Với hướng dẫn tận tình, vốn kiến thức sâu rộng khoa học kinh nghiệm sẵn có thầy hình thành cho ý tưởng ban đầu đường nghiên cứu, suốt trình nghiên cứu luận văn thầy thường xuyên theo dõi nhắc nhở hướng dẫn giải vấn đề khó khăn mà gặp phải Chính lời khuyên bảo vô qúy báu thầy động lực để hoàn thành luận văn Tôi xin gửi lòng biết ơn đến thầy Đỗ Kiến Quốc người thầy truyền đạt cho kiến thức động lực học làm sở nghiên cứu luận văn xin chân thành cảm ơn thầy cô truyền đạt kiến thức học kinh nghiệm qúy báu cho hai năm qua Tôi tin học vô may mắn có để làm tảng cho bước vững tương lai Xin chân thành cảm ơn Ban Giám Hiệu, Phòng Đào Tạo Sau Đại Học trường ĐH Bách Khoa TP.HCM tạo điều kiện thuận lợi để học tập Tôi xin bày tỏ lòng biết ơn sâu sắc đến Cha Mẹ người cho tất để có kiến thức ngày hôm nay, cuối xin cảm ơn đến anh em, bạn bè, đồng nghiệp hỗ trợ , giúp đỡ tạo điều kiện cho suốt thời gian học tập thực luận văn Chân thành cảm ơn ! Tác giả Phạm Hữu Nghóa TÓM TẮT LUẬN VĂN o0o -Trong giới đại ngày , với phát triển khoa học kỹ thuật ngày cao, nhu cầu sử dụng thẩm mỹ người sử dụng ngày cao tạo điều kiện cho công trình cao vật liệu nhẹ, vật liệu cường độ cao phát triển Sự phát triển yếu tố làm cho độ cứng công trình cao vốn thấp lại thấp, nên công trình nhạy với tác động từ bên trận động đất gió lớn thường xuyên xảy Nên vấn đề dao động công trình quan tâm ngành xây dựng Giải pháp điều khiển kết cấu xem làø ý tưởng tiến “hay nói bước tiến công nghệ nhà cao tầng nói riêng ngày xây dựng nói chung” giới đại ngày Trên cở sở đó, tác giả nghiên cứu làm việc hệ cản MTMD(Multiple tuned mass dampers) với mục tiêu khắc phục số nhược điểm hệ TMD, hệ MTMD điều khiển dải rộng tần số khác nhau, TMD MTMD điều khiển tần số khác nhau, nên hệ phản ứng tốt với tải ngẫu nhiên có tần số thay đổi trước tải thường gặp thực tế Trong luận văn tác giả khảo sát hiệu làm việc MTMD chung với kết cấu, bắt đầu việc thiết lập ma trận độ cứng tương ứng với vị trí bố trí TMD, sau phân tích tìm thông số tối ưu β , fopt, ξ opt thông qua phương pháp số dựa việc cực tiểu hoá chuyển vị gia tốc Để đạt hiệu toán điều khiển, luận văn tác giả xác định mode cần điều khiển bố trí hệ điều khiển để điều khiển tần số mode dao động cho có lợi hiệu cao tìm vị trí thông số hệ MTMD (khối lượng, độ cứng , hệ số cản, β , fopt, ξ opt ) yếu tố quan trọng cần thiết toán điều khiển MTMD Dựa kết phân tích, luận văn đưa nhận xét kết luận phạm vi sử dụng hiệu hệ cản MTMD đề số biện pháp nhằm cải tiến nâng cao hiệu hướng phát triển đề tài Nghiên cứu giải pháp điều khiển bị động với hệ cản MTMD MỤC LỤC TỜ GIAO NHIỆM VỤ TỜ NHẬN XÉT LỜI CẢM ƠN TÓM TẮT LUẬN VĂN MỤC LỤC Chương 1: TỔNG QUAN …………………………………………………… … 1.1 Giới thiệu …… …………………………….………………………… 1.2 Tổng quan điều khiển dao động …………………………………… 1.3 Điều khiển bị động ……………………………………………………… 1.3.1 Hệ cô lập ……………………………………………………… 1.3.2 Hệ tiêu tán lượng ……………………………………… 1.4 Điều khiển chủ động - bán chủ động – kết hợp ………………………… 1.4.1 Hệ chủ động …………………………………………………… 1.4.2 Hệ bán chủ động ……………………………………………… 10 1.4.3 Hệ kết hợp ……………………………………………………… 11 1.5 Hệ cản MTMD …………………………………………………………… 13 1.5.1 Cấu tạo, chế hoạt động hệ MTMD ………………………… 13 1.5.2 Lịch sử phát triển tình hình nghiên cứu hệ cản MTMD …… 15 1.5.3 Một số công trình sử dụng MTMD …………………………… 18 1.6 Nhiệm vụ luận văn ……………………………………………………… 21 Chương 2: CƠ SỞ LÝ THUYẾT ……………………………………………………………… …………………… 22 2.1 Sự tiêu tán lượng hệ cản bị động ……………………………… 22 2.2 Thiết lập phương trình dao động ………………………………………… 26 2.2.1 Thiết lập phương trình hệ BTD ……………………………… 26 2.2.2 Thiết lập phương trình hệ BTD + 1TMD ……………………… 28 2.2.3 Thiết lập phương trình hệ BTD + nTMD ……………………… 29 2.2.4 Thiết lập phương trình hệ n BTD + nTMD tầng ………… 33 2.2.5 Thiết lập phương trình hệ n BTD + nTMD nhiều tầng …… 36 Chương 3: XÁC ĐỊNH THÔNG SỐ TỐI ƯU CHO HỆ MTMD ………….……… 38 3.1 Mô hình tính …………………………………………………………… 38 3.2 Định nghóa thông số MTMD ……………………………………… 40 3.3 Phương pháp xác định thông số tối ưu MTMD …………………… 41 3.4 Bố trí hệ MTMD điều khiển dao động ………………………………… 42 3.5 Các bước thiết kế MTMD ……………………………………………… 43 Chương 4: PHÂN TÍCH ĐỘNG …………………………………………………………………… ………… 44 4.1 Các giả thiết tính toán …………………………………………………… 44 4.2 Xác định ma trận tính chất …………………………………………… 45 4.3 Các phương pháp phân tích động ……………………………………… 47 4.4 Thuật giải Newmark …………………………………………………… 48 Chương 5: KIỂM CHỨNG CHƯƠNG TRÌNH ……………………………………………………… 50 5.1 Kiểm tra tính suy biến MTMD ………………………………………… 50 5.2 Kiểm kết với nghiên cứu trước …………………………………… 55 5.3 Kiểm tra Sap2000 ………………………………………………………… 58 Chương 6: TÍNH TOÁN BẰNG SỐ ……………………………………………………………….…………… 62 6.1 Bài toán điều khiển 1BTD + nTMD n=3 ……………………………… 62 6.2 Bài toán điều khiển 1BTD + nTMD n=3.5 với mô hình MTMD …… 76 6.3 Bài toán điều khiển n BTD + nTMD n=3.5.7.11 ……………………… 89 6.4 Bài toán điều khiển n BTD+nTMD n=3 xét vị trí bố trí hệ cản TMD … 111 Chương 7: KẾT LUẬN, KIẾN NGHỊ – HƯỚNG PHÁT TRIỂN ĐỀ TÀI … 121 7.1 Kết luận kiến nghị ………………………………………………… 121 7.2 Hướng phát triển đề tài ……………………………………………… 123 TÀI LIỆU THAM KHẢO ……………………………………………………….……………………………………… 126 PHỤ LỤC TÍNH TOÁN ………………………………………………………….……………………………………… 130 Chương Tổng quan Chương 1: TỔNG QUAN 1.1 Giới Thiệu Ngày thảm họa trận động đất mối đe dọa người nói chung trận động đất thường xuyên xảy mối nguy hiểm công nghiệp xây dựng nhà cao tầng Các trận động đất gây thiệt hại lớn người mà kinh tế Theo số liệu thống kê từ trận động đất[34] Bảng số liệu Trận động Địa điểm Ngày tháng Độ lớn đất Thiệt hại Người ShaanXi Trung Quốc 23.1.1956 830.000 Tangshang Trung Quốc 27.1.1976 7.5 250.000 Izmit Thổ Nhó Kỳ 17.8.1999 7.6 17.118 Spitak Armania 7.12.1988 6.8 25.000 Gujarat Ấn Độ 26.1.2001 6.6 20.085 Sumatra An Độ Dương 26.12.2004 230.000 Kashmir Pakistan 8.10.2005 7.8 100.000 Hình ảnh tàn phá số trận động đất Động đất Kobe Động đất Nhật 26.7.2007 Hình 1.1 Nhà sụp đổ nhật trận động đất Kobe 1995 & động đất 7.2007 Nghiên cứu giải pháp điều khiển bị động Với hệ cản MTMD Chương Tổng quan Động đất Northridge Động đất Kobe Động đất Marina Động đất Pakistan 8.10.05 Hình 1.2 Các nhà bị sụp đổ trận động đất Trước phát triển khoa học kỹ thuật công nghệ tiên tiến tiền đề cho phát triển công trình cao tầng mà đặc điểm bật công trình sử dụng vật liệu cường độ cao sắt, kết cấu liên hợp … với liên kết cứng nửa cứùng nên phản ứng nhạy với tác động bên nên nguy hiểm trước trận động đất ta giải pháp khắc phục thích hợp Trên sở quan điểm điều khiển kết cấu đời xem quan điểm tiến công nghiệp nhà cao tầng yếu tố thiết kế quan trọng công trình xây dựng Nghiên cứu giải pháp điều khiển bị động Với hệ cản MTMD Chương Tổng quan 1.2 Tổng quan điều khiển dao động Khi chịu tác động bên động đất phần lớn lượng trận động đất truyền nên công trình gây biến dạng chuyển vị làm hư hỏng công trình Để khắc phục yếu tố ta có giải pháp làm tiêu tán lượng tác nhân bên gây giải pháp điều khiển kết cấu Hệ thống giải pháp điều khiển kết cấu mô sau Sơ đồ giải pháp điều khiển kết cấu (Structure Control) Điều khiển kết cấu Structure Control ĐK Bị động Passive control ĐK chủ động +Bán chủ động Active + Semi- active ĐK kết hợp (Hybrid) Active + Passive Hệ cô lập Seimic Isolate + Gối đàn hồi + Gối ma sát Hệ tiêu tán Energy dissipation + Hệ cản kim loại ( Metallic damper) + Hệ cản masát ( Friction damper) + Hệ cản chất lỏng ( Visco – viscoelestic damper) Hệ cản khối lượng Mass damper + Hệ cản khối lượng ( TMD ) + Hệ cản chất lỏng ( TLCD) 1.3 Hệ thống điều khiển bị động( Passive control) Hệ điều khiển bị động hệ thống làm giảm phản ứng kết cấu hấp thụ phần lượng kết cấu, tác nhân bên (gió , động Nghiên cứu giải pháp điều khiển bị động Với hệ cản MTMD Phụ Lục FF=length(f) for f1=1:FF FT(f1,1)=f(1,f1); end % -Ssi= Si_s: delta_Si: Si_e; % Si khoang ti so can SS=length(Ssi) for s1=1:SS SI(s1,1)=Ssi(1,s1); end % -% - Bat dau vong lap tinh toan cac thong so cua TMD lai -% -for i2=1:be % Beta khoang tan so beta2(i2)=BE(i2,1); for i4=1:FF % Ft khoang ti so tan so ft(i4)=FT(i4,1); for i3=1:SS % Si khoang ti so can si(i3)=SI(i3,1); mega0=ft(i4)*mega; mega1=mega0*(1-0.5*beta2(i2)); mega2=mega0*(1+0.5*beta2(i2)); k11=mm3*mega0*mega0; k22=mm3*mega1*mega1; k33=mm3*mega2*mega2; c11=2*si(i3)*sqrt(mm3*k11); c22=2*si(i3)*sqrt(mm3*k22); c33=2*si(i3)*sqrt(mm3*k33); 3.4 Giải phương trình trị riêng tìm tần số mode dao động FILE NAME: eigen.m % -% -Giai Phuong trinh tri rieng tim tan so va mode dao dong -% -[V,D]=eig(KK0,MM0); [lambda,k]=sort(diag(D));V=V(:,k); Factor=diag(V'*MM0*V); Phi=V*inv(sqrt(diag(Factor))); Phi1=Phi; Omega=diag(sqrt(V'*KK0*V)) 3.5 Xác định ma trận cản theo Rayleigh FILE NAME: damping.m 14 Mã Nguồn Matlab Phụ Luïc % -% - Xac Dinh Ma Tran He so Can % -si=0.05; % Ti so can a0=2*si*Omega(1)*Omega(2)/(Omega(1)+Omega(2)); a1=2*si/(Omega(1)+Omega(2)); CC0=a0*MM0+a1*KK0; % -3.6 Tính thông số ban đầu thuật giải Newmark FILE NAME: firstNewmark.m %= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = % firstNewmark;% Tinh Gia Tri Ban Dau %= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = xt0=zeros(nel3,1); vt0=zeros(nel3,1); FFt0=zeros(nel3,1); at0=inv(MM)*(FFt0-CC*vt0-KK*xt0); % step 1.1 gama=0.56; beta=0.25*(0.5+gama)^2; a1=1/beta/dtat^2; a2=1/beta/dtat; a3=1/beta/2-1; a4=gama/beta/dtat; a5=gama/beta-1; a6=dtat*(gama/2/beta-1); xti=xt0;vti=vt0;ati=at0; k=KK+CC*a4+MM*a1; % step 1.3 an=MM*a2+CC*gama/beta; bn=MM*1/2/beta+CC*a6; 3.7 Gán Tải Trọng Cho Kết Cấu FILE NAME: load.m %= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = % load;% Tinh Tai Trong %= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = for s=1:1:nel3 if i8 10') 21 Mã Nguồn Matlab Phụ Lục DỮ LIỆU PHỔ GIA TỐC ĐỘNG ĐẤT ELCENTRO o0o Data of Elcentro May 18, 1940 North South Component 2829 point at equal spacing of 0.02 sec The units are (g) Elcen =9.81*[ 0.0011 0.0012 0.0013 0.0011 0.0013 0.0011 0.0012 0.0013 0.0012 0.0012 0.0012 0.0012 0.0012 0.0011 -0.0003 -0.0118 -0.0151 -0.0188 -0.0136 -0.0153 0.0357 0.0360 0.0413 -0.0238 0.1258 0.0864 -0.2004 -0.1646 0.0811 0.3272 -0.1887 0.1217 -0.0584 -0.1244 -0.0083 -0.0108 0.0019 0.0266 -0.1278 -0.1088 -0.0328 0.0254 0.1216 -0.0023 0.1144 0.1324 -0.2430 0.0784 0.0967 -0.1088 -0.1353 -0.1168 -0.0924 0.0011 0.0016 0.0011 0.0012 0.0010 0.0011 0.0011 0.0012 0.0012 0.0012 0.0012 0.0013 0.0011 0.0011 -0.0104 -0.0099 -0.0136 -0.0188 -0.0196 -0.0154 0.0488 0.0433 0.0419 -0.0044 0.1578 0.0933 -0.2058 -0.1354 0.1199 0.3495 -0.1102 0.1796 0.0241 -0.1013 -0.0424 0.0191 0.0353 0.0696 -0.0705 -0.1169 -0.0336 -0.0005 0.1497 -0.0126 0.1464 -0.1041 -0.2068 0.1304 0.1419 -0.0929 -0.1143 -0.1153 -0.0964 0.0013 0.0012 0.0012 0.0016 0.0011 0.0013 0.0011 0.0012 0.0012 0.0011 0.0011 0.0012 0.0011 0.0013 -0.0092 -0.0072 -0.0096 -0.0050 -0.0255 -0.0054 0.0513 0.0554 0.0068 0.0154 0.1486 0.1027 -0.1831 -0.1090 0.1639 0.2884 -0.0760 0.0602 -0.0666 -0.0765 -0.0828 0.0431 0.0574 0.0709 -0.0554 -0.0781 -0.0105 -0.0224 0.1764 0.0143 0.1656 -0.0553 -0.1874 0.1815 0.1857 -0.0506 -0.0809 -0.0759 -0.0449 0.0014 0.0013 0.0012 0.0015 0.0013 0.0003 0.0009 0.0013 0.0011 0.0012 0.0012 0.0014 0.0011 0.0012 -0.0078 -0.0075 -0.0071 0.0047 -0.0321 0.0041 0.0437 0.0666 -0.0525 0.0328 0.1183 0.1249 -0.1741 -0.0781 0.2008 0.2383 -0.0168 -0.2639 -0.1987 -0.0526 -0.0946 0.0688 0.0899 0.1332 0.0061 -0.0569 0.0020 0.0077 0.0510 0.0214 0.1970 -0.0391 -0.1367 0.0095 0.2458 -0.1198 -0.0538 -0.0569 -0.0175 22 0.0013 0.0013 0.0013 0.0013 0.0012 0.0014 0.0013 0.0016 0.0013 0.0012 0.0009 0.0013 0.0012 0.0012 -0.0085 -0.0123 -0.0029 0.0161 -0.0299 0.0172 0.0374 0.0759 -0.0793 0.0524 0.0963 0.0332 -0.1772 -0.0422 0.2468 -0.1208 0.0121 -0.1548 -0.1652 -0.0267 -0.0719 -0.0073 0.1151 0.1381 0.0672 -0.0224 0.0300 0.0217 0.0077 0.0502 0.1904 -0.0316 -0.0999 -0.0127 0.1785 -0.2049 -0.0542 0.0011 0.0286 0.0012 0.0011 0.0011 0.0012 0.0014 0.0013 0.0010 0.0013 0.0013 0.0011 0.0011 0.0012 0.0016 0.0014 -0.0111 -0.0171 -0.0058 -0.0048 -0.0156 0.0257 0.0286 0.0676 -0.0594 0.0741 0.0918 -0.1516 -0.1773 -0.0003 0.2790 -0.2382 0.0552 -0.1747 -0.1703 -0.0040 -0.0606 -0.0355 0.1388 0.2068 -0.1016 -0.0124 0.0496 0.0569 0.0269 0.0655 0.2021 -0.1070 -0.0371 0.0271 -0.1117 -0.1663 -0.1146 -0.0724 0.0624 0.0012 0.0013 0.0013 0.0011 0.0013 0.0013 0.0016 0.0013 0.0011 0.0016 0.0014 0.0012 0.0011 0.0012 -0.0133 -0.0187 -0.0125 -0.0126 -0.0189 0.0271 0.0248 0.0619 -0.0476 0.1029 0.0957 -0.2089 -0.1826 0.0382 0.3106 -0.1644 0.0920 -0.1024 -0.1498 0.0196 -0.0336 -0.0034 0.0264 -0.0841 -0.1472 -0.0655 0.0621 0.0836 0.0306 0.0967 0.1828 -0.1631 0.0123 0.0448 -0.0992 -0.1707 -0.1206 -0.1554 0.0871 Mã Nguồn Matlab Phụ Lục 0.0531 0.1346 0.0383 -0.0237 0.0014 0.0076 0.0004 0.0090 -0.0129 -0.0101 0.0784 -0.0163 0.0417 -0.0059 0.0182 0.0051 0.0427 -0.0224 0.0134 -0.0313 -0.0102 0.0323 -0.1215 -0.0495 0.0152 0.1274 -0.0618 0.1144 -0.1060 0.1509 0.0735 0.0306 0.0078 0.0014 -0.0149 -0.0153 -0.0054 -0.0268 -0.0406 -0.0175 -0.0310 -0.0047 -0.0344 -0.0818 0.0360 0.1256 0.0127 -0.1107 -0.0874 0.0461 0.0288 0.0172 -0.0054 0.0496 0.0472 0.0216 -0.0075 0.0223 -0.0173 0.0048 -0.0668 0.0096 0.0248 0.1613 0.0501 -0.0156 -0.0366 0.0204 0.0078 -0.0025 -0.0024 -0.0087 0.0293 -0.0414 0.0186 0.0047 0.0404 0.0098 0.0551 -0.0366 -0.0129 -0.0269 0.0056 0.0320 -0.1341 -0.0304 0.0089 -0.0529 -0.0654 0.0589 -0.1149 0.0244 0.0084 0.0526 -0.0261 0.0032 0.0112 -0.0273 0.0203 -0.0253 -0.0118 0.0029 -0.0376 -0.0338 -0.0279 -0.0656 0.0559 0.1334 -0.0438 -0.0982 -0.0873 0.0426 0.0221 0.0099 -0.0229 0.0504 0.0563 -0.0023 0.0271 -0.0344 -0.0283 0.0214 -0.0872 0.0346 0.0014 0.1871 0.0278 -0.0071 -0.0153 0.0372 0.0053 -0.0259 0.0184 -0.0009 -0.0169 -0.0424 -0.0097 0.0131 0.0497 0.0239 0.0739 -0.0403 -0.0420 -0.0258 0.0249 0.0170 -0.1361 -0.0101 0.0229 -0.0483 -0.0577 0.0639 -0.0900 0.0013 -0.0448 0.0719 -0.0418 -0.0064 0.0431 -0.0226 0.0438 -0.0329 0.0041 0.0286 -0.0276 -0.0579 -0.0437 -0.0516 0.0726 0.1528 -0.0574 -0.0821 -0.0882 0.0290 0.0391 -0.0041 -0.0258 0.0496 0.0601 -0.0234 0.0068 -0.0722 -0.0285 0.0333 -0.0779 0.0234 0.0275 0.1379 -0.0031 0.0134 0.0075 0.0509 -0.0068 -0.0396 -0.0055 0.0072 -0.0125 -0.0243 -0.0179 0.0256 0.0263 0.0359 0.0459 -0.0453 -0.0508 -0.0259 0.0315 0.0211 -0.1239 -0.0078 0.0517 -0.0432 -0.0156 0.0284 -0.0490 -0.0019 -0.0266 0.0806 -0.0178 -0.0223 0.0768 -0.0666 0.0418 -0.0487 0.0006 0.0544 -0.0111 -0.0658 -0.0583 -0.0384 0.0891 0.1754 -0.0640 -0.0792 -0.0681 0.0406 0.0413 -0.0089 -0.0109 0.0334 0.0744 -0.0374 -0.0286 -0.0492 -0.0280 0.0689 -0.0639 -0.0234 23 0.0442 0.0748 -0.0144 0.0362 0.0223 0.0217 -0.0029 -0.0258 -0.0080 0.0223 -0.0047 -0.0068 -0.0088 0.0249 0.0131 0.0513 0.0429 -0.0431 -0.0390 -0.0323 0.0363 0.0166 -0.1089 -0.0151 0.0812 -0.0282 0.0274 -0.0006 -0.0085 0.0012 -0.0177 0.0674 0.0094 -0.0378 0.1074 -0.0925 0.0191 -0.0659 -0.0003 0.0663 0.0076 -0.0613 -0.0804 -0.0247 0.1034 0.1986 -0.0871 -0.0838 -0.0187 0.0560 0.0332 -0.0033 0.0098 0.0237 0.0828 -0.0412 -0.0204 -0.0201 -0.0237 0.0932 -0.0431 -0.0397 0.0770 0.0281 -0.0106 0.0570 -0.0074 0.0012 -0.0006 -0.0232 -0.0031 0.0350 0.0139 0.0135 -0.0177 0.0119 0.0016 0.0553 0.0171 -0.0448 -0.0395 -0.0308 0.0363 0.0348 -0.0888 -0.0144 0.1139 -0.0786 0.0907 -0.0393 0.0409 0.0122 -0.0063 0.0409 0.0459 -0.0367 0.0876 -0.0706 0.0009 -0.0771 -0.0103 0.0409 0.0301 -0.0561 -0.1033 -0.0058 0.1133 0.1574 -0.1143 -0.0859 0.0144 0.0568 0.0387 0.0048 0.0330 0.0273 0.0911 -0.0332 0.0040 -0.0014 -0.0122 0.0541 -0.0228 -0.0528 0.1031 0.0328 -0.0207 0.0308 -0.0041 -0.0084 0.0041 -0.0177 -0.0087 0.0673 0.0076 0.0305 -0.0138 0.0022 -0.0007 0.0467 -0.0124 -0.0129 -0.0321 -0.0226 0.0372 -0.0578 -0.0699 0.0039 0.1602 -0.0715 0.1532 -0.0734 0.1016 0.0441 0.0049 0.0131 0.0238 -0.0278 0.0423 -0.0424 -0.0146 -0.0609 -0.0318 0.0073 0.0298 -0.0454 -0.0959 0.0146 0.1195 0.0799 -0.1204 -0.0865 0.0239 0.0445 0.0333 0.0067 0.0444 0.0359 0.0633 -0.0270 0.0315 -0.0024 -0.0046 -0.0069 -0.0042 -0.0587 Mã Nguồn Matlab Phụ Lục -0.0444 0.1007 0.0478 -0.0149 -0.0066 -0.0523 -0.0111 0.0198 0.0119 0.0061 0.0416 -0.0572 -0.0406 -0.0071 0.0541 -0.0344 0.0306 -0.0326 -0.0314 -0.0081 0.0845 -0.0033 -0.0575 -0.0218 -0.0479 0.0146 0.0082 0.0208 0.0104 -0.0161 0.0064 0.0330 -0.0076 -0.0535 -0.0132 0.0163 -0.0301 0.0204 -0.0287 0.0125 -0.0420 -0.0021 0.0447 0.0207 -0.0216 0.0102 0.0530 -0.0563 0.0383 -0.0005 0.0578 -0.0379 -0.0170 -0.0028 -0.0338 0.0044 0.0396 -0.0190 0.0154 0.0103 -0.0128 0.0144 -0.0134 0.0526 0.0244 -0.0242 0.0071 -0.0399 -0.0116 -0.0073 -0.0098 -0.0021 0.0440 -0.0538 -0.0424 -0.0194 0.0606 -0.0111 0.0252 -0.0253 -0.0326 -0.0144 0.0693 0.0072 -0.0650 -0.0362 -0.0450 0.0584 -0.0021 0.0174 -0.0044 -0.0291 0.0082 0.0375 -0.0056 -0.0309 -0.0268 0.0208 -0.0284 0.0395 -0.0259 -0.0025 -0.0480 0.0037 0.0685 -0.0086 0.0024 -0.0093 0.0637 -0.0753 0.0609 -0.0174 0.0503 -0.0467 -0.0142 -0.0118 -0.0416 0.0291 0.0362 -0.0018 0.0043 0.0064 -0.0258 0.0364 0.0367 0.0101 0.0008 -0.0397 0.0128 -0.0162 0.0114 -0.0202 -0.0316 0.0118 0.0586 -0.0319 -0.0284 -0.0320 0.0495 0.0166 0.0223 -0.0094 -0.0228 -0.0033 0.0413 0.0171 -0.0536 -0.0427 -0.0394 0.0791 -0.0074 0.0048 -0.0154 -0.0233 0.0119 0.0435 0.0053 0.0036 -0.0436 0.0289 -0.0286 0.0374 -0.0051 -0.0160 -0.0461 -0.0021 0.0750 -0.0297 0.0290 -0.0298 0.0544 -0.0913 0.0584 -0.0119 0.0304 -0.0516 -0.0093 -0.0219 -0.0496 0.0472 0.0306 0.0208 -0.0003 0.0027 -0.0418 0.0504 0.0741 -0.0134 -0.0117 -0.0605 0.0191 0.0081 0.0364 -0.0223 -0.0406 0.0303 0.0687 -0.0133 -0.0125 -0.0229 0.0186 0.0330 0.0204 0.0059 -0.0197 0.0194 0.0149 0.0155 -0.0405 -0.0373 -0.0331 0.0864 0.0009 -0.0045 -0.0113 -0.0080 0.0160 0.0466 -0.0048 0.0279 -0.0424 0.0358 -0.0309 0.0245 0.0171 -0.0277 -0.0404 -0.0086 0.0780 -0.0408 0.0481 -0.0355 0.0382 -0.0813 0.0402 0.0019 0.0143 -0.0459 -0.0033 -0.0329 -0.0574 0.0649 -0.0024 0.0401 0.0035 -0.0018 -0.0556 0.0504 24 0.1081 -0.0062 -0.0088 -0.0644 0.0338 0.0330 0.0638 -0.0154 -0.0281 0.0491 0.0473 0.0039 0.0056 -0.0021 -0.0344 0.0239 0.0126 0.0071 -0.0171 0.0491 -0.0044 -0.0014 -0.0246 -0.0330 -0.0268 0.0809 0.0066 -0.0010 0.0053 -0.0030 0.0202 0.0364 -0.0334 0.0199 -0.0148 0.0339 -0.0304 0.0120 0.0412 -0.0334 -0.0386 -0.0144 0.0734 -0.0468 0.0583 -0.0129 0.0151 -0.0579 0.0361 0.0149 -0.0042 -0.0357 0.0013 -0.0367 -0.0624 0.0581 -0.0281 0.0498 0.0065 -0.0014 -0.0633 0.0378 0.1351 0.0131 -0.0426 0.0281 0.0310 0.0720 -0.0028 -0.0112 0.0574 0.0053 0.0017 0.0119 0.0181 -0.0679 0.0234 -0.0024 -0.0143 -0.0112 0.0791 -0.0157 -0.0191 -0.0090 -0.0354 -0.0209 0.0673 0.0113 0.0088 0.0145 0.0018 0.0246 0.0219 -0.0571 -0.0012 0.0048 0.0096 -0.0194 -0.0056 0.0546 -0.0361 -0.0303 -0.0044 0.0621 -0.0421 0.0433 0.0068 -0.0091 -0.0273 0.0315 0.0324 -0.0193 -0.0269 0.0011 -0.0326 -0.0534 0.0478 -0.0451 0.0375 0.0118 0.0016 -0.0277 0.0354 0.1321 0.0374 -0.0047 -0.0229 -0.0197 0.0104 0.0513 0.0146 0.0073 0.0503 -0.0304 -0.0207 0.0026 0.0384 -0.0540 0.0340 -0.0209 -0.0250 -0.0029 0.0906 -0.0112 -0.0388 -0.0067 -0.0448 -0.0109 0.0343 0.0148 0.0132 0.0045 0.0044 0.0289 0.0071 -0.0614 -0.0061 0.0109 -0.0171 -0.0013 -0.0178 0.0369 -0.0389 -0.0150 0.0194 0.0452 -0.0377 0.0268 0.0317 -0.0346 0.0139 0.0213 0.0456 -0.0287 -0.0189 0.0011 -0.0287 -0.0360 0.0398 -0.0419 0.0276 0.0140 -0.0033 -0.0019 0.0322 Mã Nguồn Matlab Phụ Lục 0.0300 -0.0287 -0.0284 0.0123 -0.0429 0.0167 0.0278 0.0213 -0.0413 -0.0186 0.0513 -0.0896 0.0333 0.0545 -0.0082 -0.0728 0.0497 0.0374 -0.0458 0.0771 0.0888 -0.0951 0.0553 0.0694 0.0504 -0.1200 -0.0161 -0.0024 0.0483 -0.0025 0.0262 -0.0153 -0.0224 -0.0451 -0.0083 -0.0452 0.0308 -0.0243 0.0311 -0.0204 0.0110 -0.0149 -0.0375 0.0418 0.0009 -0.0137 0.0044 0.0175 -0.0011 -0.0093 0.0013 0.0122 -0.0023 -0.0462 0.0194 0.0045 -0.0159 0.0080 0.0330 0.0302 -0.0330 -0.0117 0.0264 -0.0263 -0.0153 0.0121 -0.0406 0.0173 0.0258 0.0348 -0.0218 -0.0401 0.0844 -0.1191 0.0246 0.0584 -0.0229 -0.0793 0.0363 0.0211 -0.0529 0.0754 0.0961 -0.1336 0.0290 0.1058 0.0184 -0.0989 -0.0126 -0.0047 0.0601 -0.0061 0.0284 -0.0103 -0.0164 -0.0421 -0.0084 -0.0389 0.0383 -0.0216 0.0383 -0.0100 0.0033 -0.0137 -0.0357 0.0426 -0.0063 -0.0203 0.0123 0.0216 0.0032 -0.0100 0.0071 0.0143 -0.0041 -0.0383 0.0213 -0.0085 -0.0108 0.0070 0.0313 0.0318 -0.0203 -0.0133 0.0157 -0.0234 0.0054 0.0104 -0.0280 0.0074 0.0081 0.0499 -0.0049 -0.0591 0.1115 -0.0841 0.0123 0.0620 -0.0301 -0.0819 0.0256 0.0051 -0.0764 0.0571 0.1001 -0.0826 -0.0034 0.1155 -0.0227 -0.0743 -0.0083 -0.0019 0.0575 -0.0098 0.0273 -0.0056 -0.0118 -0.0319 -0.0140 -0.0338 0.0366 -0.0113 0.0319 0.0034 -0.0002 -0.0062 -0.0234 0.0400 -0.0124 -0.0271 0.0223 0.0255 0.0118 -0.0105 0.0147 0.0168 -0.0079 -0.0274 0.0276 -0.0188 -0.0116 0.0062 0.0264 0.0208 -0.0129 -0.0118 -0.0025 -0.0265 0.0206 0.0093 -0.0163 -0.0064 -0.0053 0.0438 0.0235 -0.0761 0.0926 -0.0569 0.0255 0.0639 -0.0327 -0.0528 0.0196 0.0026 -0.0585 0.0481 0.0948 -0.0410 -0.0374 0.1126 -0.0608 -0.0538 -0.0034 0.0007 0.0469 -0.0197 0.0159 -0.0122 -0.0183 -0.0191 -0.0160 -0.0246 0.0273 -0.0052 0.0184 0.0179 -0.0011 0.0015 -0.0148 0.0338 -0.0118 -0.0348 0.0273 0.0146 0.0106 -0.0123 0.0218 0.0148 -0.0179 -0.0086 0.0354 -0.0226 -0.0115 0.0079 0.0229 0.0110 0.0009 -0.0087 25 -0.0225 -0.0313 0.0171 0.0007 -0.0019 -0.0208 -0.0198 0.0153 0.0258 -0.0573 0.0466 -0.0178 0.0375 0.0600 -0.0373 -0.0164 0.0418 -0.0076 -0.0142 0.0557 0.0788 0.0143 -0.0313 0.0958 -0.1058 -0.0478 0.0038 0.0099 0.0301 -0.0168 -0.0016 -0.0157 -0.0262 -0.0051 -0.0238 -0.0152 0.0101 0.0038 0.0022 0.0317 -0.0029 -0.0016 0.0013 0.0190 -0.0102 -0.0351 0.0208 0.0060 0.0048 -0.0162 0.0216 0.0130 -0.0326 0.0076 0.0347 -0.0304 -0.0067 0.0145 0.0174 -0.0041 0.0024 -0.0056 -0.0349 -0.0359 0.0146 -0.0138 0.0120 -0.0180 -0.0050 -0.0091 0.0121 -0.0214 0.0074 0.0201 0.0536 0.0452 -0.0504 0.0168 0.0619 -0.0290 0.0238 0.0668 0.0231 0.0642 0.0010 0.0816 -0.1339 -0.0391 0.0093 0.0229 0.0112 0.0051 -0.0097 -0.0190 -0.0320 0.0043 -0.0327 0.0004 -0.0043 0.0114 -0.0161 0.0253 -0.0036 -0.0146 0.0181 0.0132 -0.0093 -0.0287 0.0172 -0.0061 0.0009 -0.0187 0.0182 0.0108 -0.0421 0.0151 0.0263 -0.0268 -0.0009 0.0199 0.0208 -0.0164 -0.0048 -0.0002 -0.0300 -0.0343 0.0148 -0.0261 0.0179 0.0098 0.0062 -0.0395 0.0107 -0.0436 0.0479 0.0605 0.0145 -0.0668 0.0411 0.0628 -0.0382 0.0555 0.0766 -0.0342 0.0829 0.0236 0.0619 -0.1234 -0.0274 0.0043 0.0356 0.0006 0.0157 -0.0208 -0.0236 -0.0339 0.0009 -0.0421 0.0184 -0.0124 0.0222 -0.0314 0.0142 -0.0082 -0.0226 0.0328 0.0077 -0.0074 -0.0208 0.0118 -0.0089 -0.0061 -0.0089 0.0154 0.0069 -0.0519 0.0176 0.0180 -0.0218 0.0056 0.0263 0.0252 -0.0325 -0.0095 0.0049 Mã Nguồn Matlab Phụ Lục 0.0103 0.0073 -0.0103 0.0088 0.0275 -0.0026 0.0174 -0.0092 0.0005 0.0104 -0.0034 0.0081 -0.0003 0.0007 -0.0038 -0.0101 0.0102 -0.0030 -0.0071 -0.0099 0.0129 0.0015 -0.0269 0.0073 0.0240 -0.0255 0.0009 0.0153 -0.0133 -0.0029 -0.0007 0.0066 -0.0081 0.0045 0.0089 0.0102 0.0151 -0.0061 -0.0266 -0.0209 -0.0055 0.0060 0.0078 0.0182 0.0031 0.0159 0.0148 -0.0018 0.0107 -0.0014 0.0085 0.0148 -0.0106 -0.0129 -0.0048 0.0004 0.0020 0.0047 0.0146 -0.0025 -0.0015 0.0031 0.0178 0.0016 -0.0039 0.0100 0.0339 -0.0045 0.0226 -0.0065 0.0002 0.0145 -0.0093 0.0067 0.0008 0.0026 -0.0068 -0.0056 0.0133 -0.0036 -0.0065 -0.0097 0.0170 -0.0064 -0.0266 0.0118 0.0176 -0.0246 0.0063 0.0134 -0.0189 -0.0036 -0.0001 0.0026 -0.0075 0.0042 0.0082 0.0085 0.0118 -0.0079 -0.0273 -0.0159 -0.0046 0.0054 0.0102 0.0171 0.0040 0.0158 0.0116 0.0071 0.0023 0.0106 0.0154 -0.0155 -0.0116 -0.0004 0.0026 0.0009 0.0063 0.0194 -0.0059 0.0012 0.0021 0.0263 -0.0007 0.0008 0.0116 0.0301 -0.0074 0.0154 -0.0041 0.0004 0.0172 -0.0086 0.0056 0.0006 0.0054 -0.0051 -0.0023 0.0080 -0.0056 -0.0104 -0.0085 0.0214 -0.0101 -0.0204 0.0162 0.0136 -0.0238 0.0108 0.0111 -0.0186 -0.0040 -0.0003 -0.0024 -0.0039 0.0054 0.0074 0.0063 0.0088 -0.0115 -0.0284 -0.0112 -0.0030 0.0054 0.0126 0.0159 0.0054 0.0138 0.0088 0.0034 0.0024 0.0026 0.0132 0.0111 -0.0193 -0.0115 0.0017 0.0043 -0.0016 0.0070 0.0155 -0.0089 0.0051 0.0013 0.0345 -0.0039 0.0010 0.0143 0.0240 -0.0038 0.0078 -0.0017 0.0009 0.0143 -0.0032 0.0028 -0.0005 0.0081 -0.0054 -0.0013 0.0046 -0.0116 -0.0129 -0.0088 0.0182 -0.0112 -0.0166 0.0203 0.0059 -0.0218 0.0110 0.0081 -0.0125 -0.0027 0.0046 -0.0073 -0.0009 0.0069 0.0086 0.0073 0.0051 -0.0156 -0.0270 -0.0096 -0.0011 0.0053 0.0151 0.0118 0.0071 0.0134 0.0048 0.0060 -0.0023 0.0038 0.0151 0.0034 -0.0179 -0.0111 0.0007 0.0068 -0.0009 0.0079 0.0111 -0.0079 0.0066 0.0001 26 0.0362 -0.0069 0.0012 0.0162 0.0108 0.0009 0.0028 0.0014 0.0114 0.0011 0.0009 -0.0008 0.0074 -0.0101 0.0014 -0.0001 -0.0107 -0.0159 -0.0049 0.0153 -0.0132 -0.0100 0.0251 -0.0049 -0.0165 0.0114 0.0026 -0.0085 -0.0006 0.0087 -0.0084 0.0030 0.0080 0.0101 0.0093 0.0020 -0.0203 -0.0260 -0.0084 0.0036 0.0048 0.0177 0.0076 0.0085 0.0138 0.0029 0.0095 -0.0081 0.0038 0.0156 -0.0029 -0.0173 -0.0106 0.0010 0.0066 0.0006 0.0089 0.0076 -0.0071 0.0056 -0.0002 0.0333 -0.0115 0.0027 0.0189 0.0006 0.0054 -0.0013 0.0013 0.0030 0.0067 0.0048 -0.0022 -0.0026 0.0063 -0.0155 0.0017 -0.0016 -0.0099 -0.0134 0.0015 0.0111 -0.0161 -0.0044 0.0283 -0.0151 -0.0121 0.0113 -0.0023 -0.0035 -0.0009 0.0136 -0.0084 0.0056 0.0099 0.0116 0.0117 -0.0016 -0.0233 -0.0249 -0.0084 0.0048 0.0044 0.0198 0.0041 0.0110 0.0137 0.0016 0.0119 -0.0080 0.0046 0.0152 -0.0053 -0.0161 -0.0106 -0.0004 0.0053 0.0017 0.0095 0.0042 -0.0061 0.0052 0.0016 0.0218 -0.0149 0.0061 0.0216 -0.0001 0.0114 -0.0066 0.0013 0.0068 0.0014 0.0059 -0.0018 -0.0021 0.0013 -0.0165 0.0051 -0.0022 -0.0087 -0.0121 0.0078 0.0073 -0.0213 0.0018 0.0269 -0.0210 -0.0066 0.0118 -0.0078 -0.0031 -0.0001 0.0107 -0.0084 0.0045 0.0099 0.0126 0.0133 -0.0037 -0.0252 -0.0243 -0.0081 0.0057 0.0060 0.0186 0.0018 0.0133 0.0147 -0.0001 0.0148 -0.0037 0.0061 0.0152 -0.0078 -0.0145 -0.0094 -0.0006 0.0040 0.0034 0.0105 0.0013 -0.0036 0.0042 0.0029 Mã Nguồn Matlab Phụ Lục 0.0046 0.0056 -0.0097 -0.0028 0.0036 -0.0123 -0.0206 0.0073 0.0043 0.0048 0.0086 0.0066 0.0155 0.0068 0.0037 0.0111 -0.0009 -0.0108 0.0004 0.0048 -0.0039 -0.0053 0.0066 -0.0044 0.0012 0.0020 -0.0016 -0.0068 0.0011 0.0056 0.0081 0.0107 0.0073 0.0044 0.0009 -0.0038 -0.0133 -0.0046 -0.0021 0.0069 0.0065 0.0011 -0.0091 -0.0063 0.0027 0.0022 -0.0040 0.0020 0.0076 0.0123 0.0075 0.0028 0.0035 -0.0038 0.0035 0.0083 -0.0035 -0.0089 -0.0079 -0.0060 -0.0053 0.0012 0.0034 0.0030 -0.0082 -0.0021 0.0011 -0.0146 -0.0188 0.0099 0.0012 0.0063 0.0088 0.0060 0.0146 0.0049 0.0056 0.0111 -0.0015 -0.0120 0.0003 0.0056 -0.0045 -0.0010 0.0061 -0.0031 0.0009 0.0010 -0.0016 -0.0071 0.0035 0.0062 0.0093 0.0102 0.0059 0.0028 0.0001 -0.0055 -0.0109 -0.0041 -0.0018 0.0068 0.0042 -0.0002 -0.0094 -0.0063 0.0058 -0.0006 -0.0041 0.0030 0.0094 0.0118 0.0067 0.0024 0.0021 -0.0033 0.0049 0.0058 -0.0061 -0.0093 -0.0079 -0.0058 -0.0046 0.0024 0.0028 0.0005 -0.0065 -0.0014 -0.0011 -0.0158 -0.0161 0.0104 -0.0017 0.0068 0.0091 0.0072 0.0135 0.0038 0.0069 0.0106 -0.0024 -0.0091 0.0008 0.0058 -0.0054 0.0022 0.0054 -0.0024 0.0009 0.0006 -0.0014 -0.0073 0.0030 0.0058 0.0107 0.0102 0.0046 0.0024 -0.0011 -0.0068 -0.0091 -0.0038 0.0001 0.0076 0.0038 -0.0094 -0.0055 0.0088 -0.0026 -0.0044 0.0033 0.0111 0.0111 0.0056 0.0018 0.0007 -0.0026 0.0076 0.0040 -0.0079 -0.0091 -0.0084 -0.0060 -0.0037 0.0029 0.0028 -0.0021 -0.0054 0.0001 -0.0030 -0.0168 -0.0120 0.0118 -0.0051 0.0071 0.0084 0.0094 0.0127 0.0016 0.0075 0.0101 -0.0039 -0.0076 0.0008 0.0055 -0.0054 0.0064 0.0053 -0.0012 0.0005 -0.0006 -0.0007 -0.0072 0.0034 0.0062 0.0113 0.0098 0.0049 0.0027 -0.0016 -0.0084 -0.0066 -0.0038 0.0024 0.0099 0.0043 -0.0002 -0.0096 -0.0054 0.0114 -0.0027 -0.0048 0.0038 0.0111 0.0106 0.0051 0.0014 -0.0006 -0.0012 0.0098 0.0014 -0.0077 -0.0088 -0.0086 -0.0058 -0.0033 0.0021 27 0.0040 -0.0054 -0.0053 0.0006 -0.0057 -0.0178 -0.0044 0.0106 -0.0033 0.0076 0.0081 0.0113 0.0110 0.0016 0.0086 0.0030 -0.0058 -0.0053 0.0019 0.0023 -0.0064 0.0079 0.0040 -0.0002 0.0008 -0.0021 -0.0016 -0.0064 0.0033 0.0061 0.0109 0.0091 0.0045 0.0025 -0.0020 -0.0098 -0.0055 -0.0031 0.0048 0.0116 0.0040 -0.0004 -0.0091 -0.0054 0.0097 -0.0033 -0.0034 0.0038 0.0116 0.0096 0.0043 0.0019 -0.0021 0.0009 0.0127 -0.0002 -0.0083 -0.0085 -0.0088 -0.0057 -0.0024 0.0014 0.0044 -0.0079 -0.0045 0.0018 -0.0079 -0.0188 0.0004 0.0091 0.0002 0.0078 0.0075 0.0130 0.0097 0.0021 0.0093 0.0009 -0.0070 -0.0031 0.0029 0.0002 -0.0071 0.0073 0.0009 0.0013 0.0012 -0.0018 -0.0034 -0.0037 0.0031 0.0061 0.0109 0.0092 0.0048 0.0024 -0.0024 -0.0113 -0.0054 -0.0030 0.0061 0.0093 0.0039 -0.0022 -0.0076 -0.0026 0.0076 -0.0036 -0.0016 0.0049 0.0118 0.0091 0.0039 0.0026 -0.0037 0.0023 0.0121 -0.0016 -0.0084 -0.0080 -0.0080 -0.0054 -0.0012 0.0009 0.0064 -0.0109 -0.0039 0.0031 -0.0108 -0.0199 0.0038 0.0069 0.0024 0.0082 0.0067 0.0141 0.0083 0.0022 0.0103 0.0008 -0.0091 -0.0006 0.0038 -0.0029 -0.0064 0.0071 -0.0019 0.0013 0.0014 -0.0016 -0.0050 -0.0015 0.0046 0.0064 0.0106 0.0089 0.0049 0.0023 -0.0028 -0.0131 -0.0047 -0.0028 0.0061 0.0082 0.0021 -0.0060 -0.0068 0.0003 0.0051 -0.0038 0.0064 0.0122 0.0085 0.0032 0.0028 -0.0044 0.0029 0.0096 -0.0022 -0.0087 -0.0077 -0.0066 -0.0055 -0.0001 -0.0003 Mã Nguồn Matlab Phụ Lục 0.0004 0.0013 0.0011 -0.0012 -0.0031 -0.0011 0.0082 0.0074 0.0047 -0.0003 0.0019 0.0070 0.0006 0.0010 0.0106 -0.0139 -0.0099 0.0106 -0.0022 0.0069 0.0048 0.0024 -0.0050 -0.0092 0.0030 0.0071 -0.0042 0.0079 0.0096 0.0003 0.0004 -0.0037 -0.0114 -0.0041 0.0011 0.0013 0.0013 0.0011 0.0012 0.0012 0.0008 0.0018 0.0013 -0.0001 -0.0013 -0.0041 -0.0007 0.0086 0.0059 0.0036 0.0011 -0.0004 0.0020 0.0084 0.0011 -0.0051 0.0104 0.0004 -0.0004 0.0098 -0.0055 -0.0037 0.0024 0.0015 -0.0021 -0.0067 0.0071 0.0051 -0.0043 0.0108 0.0074 -0.0016 -0.0010 -0.0053 -0.0089 0.0006 0.0012 0.0012 0.0014 0.0012 0.0012 0.0005 0.0029 0.0026 -0.0012 -0.0018 -0.0036 -0.0013 0.0093 0.0042 0.0018 0.0013 0.0001 0.0022 0.0092 0.0042 -0.0017 0.0093 0.0121 0.0082 0.0064 -0.0043 -0.0108 0.0010 0.0002 0.0009 -0.0012 0.0086 0.0035 -0.0042 0.0100 0.0050 -0.0038 -0.0021 -0.0078 -0.0066 0.0011 0.0012 0.0013 0.0013 0.0010 0.0010 0.0012 0.0023 0.0026 -0.0013 -0.0016 -0.0032 -0.0001 0.0095 0.0041 0.0004 0.0016 0.0003 0.0022 0.0074 0.0063 0.0074 0.0104 0.0111 0.0151 0.0036 -0.0009 -0.0086 0.0004 -0.0016 0.0046 0.0043 0.0084 0.0010 -0.0048 0.0093 0.0023 -0.0056 -0.0041 -0.0103 -0.0036 0.0010 0.0010 0.0013 0.0013 0.0012 0.0012 0.0015 28 0.0026 0.0036 -0.0006 -0.0011 -0.0028 0.0025 0.0094 0.0042 -0.0020 0.0012 0.0008 0.0034 0.0063 0.0096 0.0144 0.0027 0.0061 0.0078 0.0028 0.0025 -0.0050 0.0019 -0.0045 0.0048 0.0053 0.0088 -0.0023 -0.0024 0.0081 0.0013 -0.0054 -0.0050 -0.0128 -0.0021 0.0012 0.0013 0.0011 0.0012 0.0012 0.0013 0.0011 0.0020 0.0033 -0.0004 -0.0009 -0.0019 0.0055 0.0097 0.0044 -0.0023 0.0007 0.0011 0.0044 0.0046 0.0096 0.0239 -0.0074 0.0015 0.0033 0.0024 0.0059 -0.0014 0.0021 -0.0074 -0.0004 0.0039 0.0086 -0.0061 0.0008 0.0088 0.0015 -0.0034 -0.0046 -0.0131 -0.0030 0.0014 0.0011 0.0009 0.0015 0.0014 0.0011 0.0011 0.0016 0.0019 -0.0006 -0.0021 -0.0013 0.0079 0.0090 0.0047 -0.0008 0.0006 0.0018 0.0063 0.0029 0.0044 0.0198 -0.0194 -0.0059 0.0046 0.0002 0.0110 0.0043 0.0036 -0.0082 -0.0035 0.0024 0.0085 -0.0058 0.0044 0.0089 0.0018 -0.0018 -0.0045 -0.0130 -0.0038 0.0014 0.0012 0.0005 0.0013 0.0011 0.0012 0.0011 0] Mã Nguồn Matlab ... : NGHIÊN CỨU GIẢI PHÁP ĐIỀU KHIỂN BỊ ĐỘNG KẾT CẤU VỚI HỆ CẢN THAY ĐỔI KHỐI LƯNG MTMD (MULTIPLE TUNED MASS DAMPER) II NHIỆM VỤ LUẬN VĂN: • Nghiên cứu làm việc nhiều hệ cản MTMD điều khiển dao động. .. hệ điều khiển bán chủ động 10 Nghiên cứu giải pháp điều khiển bị động Với hệ cản MTMD Chương Tổng quan 1.4.3 Hệ điều khiển kết hợp (Hybrid control) [25] Hệ kết nối từ hệ điều khiển chủ động hệ. .. dụng hệ điều khiển kết hợp 11 Nghiên cứu giải pháp điều khiển bị động Với hệ cản MTMD Chương Tổng quan Kết luận : Tuy hệ cản bị động có nhược điểm so với hệ cản chủ động khả thích nghi với biến đổi