(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets(Luận văn thạc sĩ) Nhận dạng và phân loại các tín hiệu quá độ dựa vào mạng Neuron kết hợp với phân tích wavelets
LUẬN VĂN TỐT NGHIỆP i LỜI CẢM ƠN Lời đầu tiên, người thực xin bày tỏ lòng biết ơn Thầy NGUYỄN HỮU PHÚC, người tận tình hướng dẫn trình thực đề tài Người thực chân thành cảm ơn Thầy Quyền Uy nh, chủ nhiệm ngành giúp đỡ tạo nhiều điều kiện để đề tài hoàn thành thời hạn Người thực xin gởi lới cám ơn đến quý Thầy Cô Trường ĐH Sư Phạm Kỹ Thuật TPHCM, Trường ĐH Bách Khoa TPHCM giảng dạy giúp đỡ trình học thực đề tài Người thực xin gởi đến quý Thầy Cô Anh Chị Phòng QLKH -–QHQT - Sau ĐH Trường ĐH Sư Phạm Kỹ Thuật TPHCM có nhiều giúp đỡ trình học thực đề tài Xin gởi đến bạn học viên khóa lớp CH Điện Trường ĐH Sư Phạm Kỹ Thuật TPHCM có nhiều ý kiến quý báu cho người thực trình thực đề tài Sau cùng, người thực xin gởi đến người thân gia đình lời cám ơn động viên giúp đỡ suốt thời gian thực đề tài TPHCM, Ngày 30 tháng 08 năm 2005 Người Thực Hiện NGUYỄN TẤN ĐỜI TP HỒ CHÍ MINH 30 – 08 – 2005 LUẬN VĂN TỐT NGHIỆP ii TÓM TẮT Nội dung luận văn trình bày theo chương, bao gồm chương Mở đầu Kết luận, với chương nội dung đề tài đưa Chương 1: Giới thiệu đề tài đưa hướng giải nêu lên phạm vi giới hạn đề tài Chương 2: Trình bày chất lượng điện với đặc điểm tượng độ xảy hệ thống điện Chương 3: Giới thiệu lý thuyết wavelets phân tích wavelets Chương 4: Giới thiệu mạng nơron đặc điểm loại mạng nơron, qui tắc học hoạt động mạng Chương 5: Giới thiệu Simulink Matlab phần mềm ATP–EMTP Mô Matlab ATP–EMTP mạch tạo dạng sóng độ: Sóng Sin chuẩn, Võng điện áp, Tăng điện áp, Méo họa tần, Nhấp nháy điện áp, Mất điện tạm thời Đóng cắt tụ điện Chương 6: Trình bày việc áp dụng phân tích wavelets đa phân giải MRA mạng nơron xác suất PNN Đưa nội dung tính toán cách thực việc phân tích phân loại tượng độ phân tích wavelets mạng nơron Chương 7: Trình bày kết mô nội dung tính toán chương 6, với thực giao diện cho người sử dụng Chương 8: Rút kết luận đưa hướngphát triển cho đề tài tương lai TP HỒ CHÍ MINH 30 – 08 – 2005 LUẬN VĂN TỐT NGHIỆP iii MỤC LỤC LỜI CẢM ƠN TÓM TẮT MUÏC LUÏC CHƯƠNG 1: MỞ ĐẦU 1.1 ĐẶT VẤN ĐỀ 1.2 GIẢI QUYẾT VẤN ĐỀ 1.3 GIỚI HẠN ĐỀ TÀI 1.4 BỐ CỤC CỦA LUẬN VĂN CHƯƠNG 2: GIỚI THIỆU CÁC HIỆN TƯNG QUÁ ĐỘ TRÊN HỆ THỐNG ĐIỆN 2.1 KHÁI NIỆM 2.1.1 Định nghóa chất lượng điện 2.1.2 Tầm quan trọng chất lượng điện 2.2 CÁC HIỆN TƯNG QUÁ ĐỘ TRÊN HỆ THỐNG ĐIỆN 2.2.1 Dao động độ 2.2.2 Thay đổi áp thời gian daøi 2.2.3 Thay đổi áp thời gian ngắn 2.2.4 Mất cân điện áp 2.2.5 Méo họa tần 2.2.6 Nhấp nháy điện áp 2.3 ĐÁNH GIÁ CHẤT LƯNG ĐIỆN NĂNG 2.3.1 Giới thiệu 2.3.2 Đánh giá chất lượng điện 2.3.3 Thiết bị đo chất lượng điện CHƯƠNG 3: WAVELETS PHÂN TÍCH WAVELETS 3.1 GIỚI THIỆU VỀ WAVELETS 3.1.1 Định nghóa Wavelets 3.1.2 Phân tích Wavlets 3.2 BIẾN ĐỔI WAVELETS 3.2.1 Biến đổi Wavelets liên tục 3.2.2 Biến đổi Wavelets rời rạc CHƯƠNG 4: MẠNG NƠRON 4.1 GIỚI THIỆU VỀ MẠNG NƠRON 4.1.1 Định nghóa mạng nơron 4.1.2 Hoạt động mạng nơron 4.1.3 Các thành phần mạng nơron 4.2 HUẤN LUYỆN MẠNG NƠRON 4.2.1 Học có giám sát/ Học có thầy TP HỒ CHÍ MINH 30 – 08 – 2005 i ii iii 1 3 3 5 9 10 13 13 14 15 15 18 21 21 22 23 24 24 LUẬN VĂN TỐT NGHIỆP iv 4.2.2 Học không giám sát/ Học thầy 4.2.3 Tốc độ học 4.2.4 Các luật học 4.3 GIỚI THIỆU CÁC LOẠI MẠNG NƠRON 4.3.1 Mạng Lan truyền ngược 4.3.2 mạng LVQ 4.3.3 Mạng Xác suất 4.3.4 Maïng Hopffield 4.3.5 Mạng Tự xếp CHƯƠNG 5: MÔ PHỎNG CÁC HIỆN TƯNG QUÁ ĐỘ 5.1 MÔ PHỎNG BẰNG PHẦN MỀM MATLAB 5.1.1 Giới thiệu Power System Blockset Toolbox 5.1.2 Mạch mô phoûng 5.2 MÔ PHỎNG BẰNG PHẦN MỀM ATP 5.2.1 Giới thiệu ATP 5.2.2 Mạch kết mô CHƯƠNG 6: ỨNG DỤNG PHÂN TÍCH WAVELETS MẠNG NƠRON TRONG PHÂN TÍCH PHÂN LOẠI QUÁ ĐỘ 6.1 ỨNG DỤNG PHÂN TÍCH WAVELETS 6.1.1 Phân tích đa phân giải 6.1.2 Định lý Parseval 6.1.3 Phân tích tín hiệu độ Wavelets 6.2 PHÂN LOẠI TỰ ĐỘNG BẰNG MẠNG NƠRON 6.2.1 Mô hình mạng nơron Xác suất 6.2.2 Nhận dạng tượng độ mạng nơron Xác suất CHƯƠNG 7: KẾT QUẢ PHÂN TÍCH VÀ NHẬN DẠNG 7.1 XÂY DỰNG GIẢI THUẬT 7.2 KẾT QUẢ PHÂN TÍCH VÀ NHẬN DẠNG 7.2.1 Nhận dạng qua đặc trưng lượng 7.2.2 Mẫu huấn luyện nhận dạng tự động 7.3 CHƯƠNG TRÌNH GIAO DIỆN 7.3.1 Giới thiệu GUI 7.3.2 Phaân tích nhận dạng tín hiệu giao diện CHƯƠNG 8: KẾT LUẬN 8.1 KẾT LUẬN 8.2 HƯỚNG PHÁT TRIỂN ĐỀ TÀI TÀI LIỆU THAM KHẢO PHỤ LỤC TP HỒ CHÍ MINH 30 – 08 – 2005 25 25 25 25 25 27 28 29 30 32 32 33 34 34 38 41 41 42 42 47 47 49 50 52 52 61 63 63 65 68 68 69 LUẬN VĂN TỐT NGHIỆP Chương MỞ ĐẦU 1.1 ĐẶT VẤN ĐỀ : Các vấn đề chất lượng điện quan tâm nhiều nước giới có nhiều tài liệu đề cập đến đề tài Việt Nam nay, vấn đề chất lượng điện chưa quan tâm nhiều tương lai tiêu chí quan trọng để người tiêu dùng lựa chọn nhà cung cấp điện Do đó, tìm hiểu chất lượng điện vấn đề liên quan công việc cần thiết thời gian Hiện có nhiều kỹ thuật để phân tích, xử lý phân loại, nhận dạng tín hiệu Đặc biệt kỹ thuật phân tích wavelets hiệu phân tích tín hiệu động, không mang tính chu kỳ, mạng nơron thích hợp cho việc phân loại tự động nhiều tín hiệu khác Nhận thấy tầm quan trọng vấn đề trên, người thực chọn đề tài Nhận Dạng Phân Loại Các Tín Hiệu Quá Độ Dựa Vào Mạng Nơron Kết Hợp Với Phân Tích Wavelets để thực luận văn tốt nghiệp Đề tài dùng để nhận dạng phân loại nhanh tượng độ hệ thống điện như: Võng điện áp, Tăng điện áp, Méo họa tần, Nhấp nháy điện áp, Mất điện áp Đóng cắt tụ điện đường dây Khi nhận biết tượng giúp người vận hành hệ thống có thao tác xử lý thích hợp để nâng cao chất lượng điện 1.2 GIẢI QUYẾT VẤN ĐỀ: Được phân công Phòng quản lý sau đại học Ban chủ nhiệm ngành với đồng ý giáo viên hướng dẫn, người thực giải đề tài theo hướng kết hợp nghiên cứu lý thuyết mô máy tính 1.3 GIỚI HẠN ĐỀ TÀI: Do thời gian thực đề tài 06 tháng bị hạn chế thiết bị nên người thực phải thu hẹp phạm vi nghiên cứu Nội dung đề tài giới hạn việc mô kết phân loại số tín hiệu có sẵn máy tính Chấp nhận liệu tín hiệu độ tạo máy tính, chưa thử nghiệm việc xử lý thời gian thực Ngoài đề tài giới hạn việc nhận biết tượng độ khác nhau, không quan tâm đến việc nhận biết khác biệt loại độ TP HỒ CHÍ MINH 30 – 08 – 2005 LUẬN VĂN TỐT NGHIỆP 1.4 BỐ CỤC CỦA LUẬN VĂN: Vì nội dung đề tài nói nhận dạng tượng độ hệ thống điện nên luận văn trình bày theo hướng tìm hiểu, phân tích vấn đề lý thuyết liên quan làm tảng cho việc lập trình mô máy tính Cụ thể luận văn xếp theo nội dung sau: - Chương trình bày vấn đề chất lượng điện Từ nêu lên tượng độ hệ thống ảnh hưởng đến chất lượng điện Đây đối tượng nhận dạng phần sau - Chương trình bày khái niệm wavelets biến đổi wavelets Lý thuyết làm tảng cho việc nghiên cứu ứng dụng wavelets - Chương giới thiệu mạng nơron, đặc điểm hoạt động loại mạng nơron Lý thuyết giúp cho việc chọn lựa tối ưu ứng dụng loại mạng - Chương giới thiệu hai phần mềm: ATP – EMTP Matlab, đồng thời nêu lên cách thực mô dạng sóng tín hiệu độ Những tín hiệu làm liệu cho chương trình phân tích phân loại - Chương nêu cụ thể cách ứng dụng phân tích wavelets đa phân giải - - phân tích tín hiệu độ Trình bày ứng dụng định lý Parseval vào biến đổi wavelets để trích đặc trưng lượng tín hiệu độ Điều làm giảm bớt số lượng mẫu nhận dạng Ngoài ra, chương trình bày cách thực mô hình mạng nơron xác suất để nhận dạng tín hiệu, cách lựa chọn liệu huấn luyện mạng Chương nêu lên toàn kết mô đạt được, bao gồm kết nhận dạng qua đặc trưng lượng tín hiệu kết huấn luyện thử nghiệm nhận dạng mạng nơron Các đặc trưng lượng tín hiệu biểu diễn qua đồ thị phân bố lượng 13 cấp phân tích wavelets db4 Cuối cùng, chương nêu lên nhận xét, kết luận người thực đề tài, với định hướng phát triển đề tài tương lai Toàn nội dung chương trình mô phỏng, với số liệu liên quan trình bày phần phụ lục luận văn TP HỒ CHÍ MINH 30 – 08 – 2005 LUẬN VĂN TỐT NGHIỆP Chương GIỚI THIỆU CÁC HIỆN TƯNG QUÁ ĐỘ TRÊN HỆ THỐNG ĐIỆN 2.1 KHÁI NIỆM: 2.1.1 Định nghóa chất lượng điện năng: Bất kỳ sai lệch dòng điện điện áp so với dạng sóng sin chuẩn nguồn cung cấp gây hư hỏng tác động sai thiết bị Người sử dụng yêu cầu phải cung cấp nguồn điện có chất lượng, dạng sóng sin chuẩn có biên độ ổn định đảm bảo độ tin cậy cho thiết bị điện tử có độ nhạy cao Có thể định nghóa khác chất lượng điện - Các công ty cung cấp điện định nghóa chất lượng điện độ tin cậy cung cấp điện - Nhà máy sản xuất thiết bị điện định nghóa chất lượng điện đặc tính nguồn cho phép thiết bị hoạt động phù hợp, đặc tính khác loại thiết bị khác tùy vào nhà máy - Tuy nhiên cuối chất lượng điện khách hàng đưa Do chất lượng điện định nghóa vấn đề hệ thống điện liên quan đến điện áp, dòng điện tần số nhằm đảm bảo cho thiết bị hoạt động cách hiệu Thực tế, chất lượng điện chất lượng điện áp hệ thống điện Các tiêu chuẩn lãnh vực chất lượng điện dành cho việc trì điện áp với giới hạn xác định 2.1.2 Tầm quan trọng chất lượng điện năng: - Sự gia tăng thiết bị kỹ thuật số có độ nhạy cao đòi hỏi nghiêm ngặt chất lượng nguồn cung cấp Các lãnh vực liên quan công nghệ bán dẫn, máy tính điện tử, người máy, thiết bị điều khiển có lập trình Ngoài thiết bị điện tử lại làm tăng ảnh hưởng đến chất lượng điện Thị trường điện tạo tính cạnh tranh làm người ta quan tâm nhiều đến chất lượng điện - Chất lượng nguồn điện ảnh hưởng trực tiếp lên nhiều nhà tiêu thụ công nghiệp, người sử dụng thiết bị tự động đại - Lý cuối để ta quan tâm đến chất lượng điện giá trị kinh tế, ảnh hưởng kinh tế lên công ty điện, khách hàng, nhà cung cấp thiết bị 2.2 CÁC HIỆN TƯNG QUÁ ĐỘ TRÊN HỆ THỐNG ĐIỆN: 2.2.1 Dao động độ: TP HỒ CHÍ MINH 30 – 08 – 2005 LUẬN VĂN TỐT NGHIỆP Khái niệm độ dùng để kiện ta không mong muốn xảy thoáng qua hệ thống điện Quá độ thường liên quan đến tượng dao động tắt dần mạch R, L C Quá độ chia làm loại: xung độ dao động qua độ, tùy vào dạng sóng điện áp dòng điện độ Xung độ: Hình 2.1: Xung độ Xung độ thay đổi đột ngột điện áp dòng điện trạng thái xác lập, mang tính đơn cực Xung độ thường mô tả thời gian xuất thới gian suy giảm Ví dụ xung độ 1.2x50s 2000V có biên độ tăng từ đến max 2000V thời gian 1.2s, sau suy giảm biên độ thời gian 50s Nguyên nhân gây xung độ sét đánh Do có chứa thành phần tần số cao nên hình dạng xung độ thay đổi tùy theo thành phần mạch, khác đáng kể vị trí khác hệ thống Xung độ kích thích tần số gây dao động độ Dao động độ: Dao động độ thay đổi đột ngột điện áp dòng điện trạng thái xác lập hệ thống điện, mang tính lưỡng cực dương âm Dao động độ gồm giá trị tức thời điện áp dòng điện thay đổi cực tính liên tục Nó mô tả thành phần chứa phổ tần, gồm tần số, thời gian biên độ Dựa vào tần số người ta chia dao động độ thành loại: tần số cao, tần số trung bình tần số thấp Hình 2.2 minh họa dạng sóng dao động độ đóng cắt tụ Hình 2.2: Dao động độ TP HỒ CHÍ MINH 30 – 08 – 2005 LUẬN VĂN TỐT NGHIỆP Trong hình vẽ mục sau này, trục hoành biểu diễn thời gian số điểm lấy mẫu tín hiệu Tốc độ lấy mẫu 256điểm/chu kỳ - Các dao động độ có tần số >500KHz, thời gian tính s, xem dao động độ tần số cao, chúng thường đáp ứng hệ thống xung độ - Các dao động có tần số từ – 500KHz với thời gian xảy hàng chục s, xem dao động độ tần số trung bình, thường đóng cắt tụ đóng cắt đường dây gây - Các dao động độ có tần số < 5KHz, thời gian khoảng 0.3 – 50ms, xem dao động độ tần số thấp Các dao động thường xuất hệ thống truyền dẫn phân phối, thường nhiều kiện khác gây 2.2.2 Thay đổi điện áp thời gian dài: Các thay đổi làm lệch điện áp dụng thời gian dài 01 phút gọi thay đổi điện áp thời gian dài Chúng bao gồm điện áp thấp điện áp - Quá áp – Overvoltage: Một tượng gọi điện áp giá trị hiệu dụng tăng lớn 110% định mức thời gian dài 01 phút Quá điện áp thường xảy đóng cắt tải có giá trị lớn, nạp điện tụ nhánh Ngoài việc cài đặt sai đầu dây máy biến áp gây áp hệ thống - Thấp áp – Undervoltage: Một tượng xem thấp áp giá trị hiệu dụng giảm nhỏ 90% định mức thời gian dài 01 phút Thấp điện áp xảy kiện ngược với điện áp Việc đóng tải có giá trị lớn cắt tụ nhánh gây thấp điện áp Thấp áp áp gây cố hệ thống làm ảnh hưởng đến thiết bị điện tử công suất, máy tính, thiết bị xử lý - Mất áp kéo dài – Sustained Interruption: Hiện tượng điện áp hệ thống giảm zero kéo dài thời gian 01 phút gọi điện áp kéo dài Hiện tượng thường xảy cần phải có can thiệp người để phục hồi lại hệ thống 2.2.3 Thay đổi điện áp thời gian ngắn: Thay đổi điện áp thời gian ngắn chia làm nhiều loại tùy vào thời gian xảy ra, gồm: Instantaneous (3– 60s), Momentary (0.6– 3s), Temporary (0.1– 0.6s) Caùc thay đổi xuất điều kiện cố, dòng khởi động tải lớn Tùy vào vị trí xảy cố điều kiện hệ thống mà cố gây giảm áp tạm thời (sag), tăng áp tạm thời (swell), điện tạm thời (interruption) - Mất điện tạm thời – Voltage Interruption: Hiện tượng xảy điện áp nguồn dòng điện tải giảm đến giá trị nhỏ 0.1pu (power unit) thời gian 01 phút Như hình vẽ 2.3 TP HỒ CHÍ MINH 30 – 08 – 2005 LUẬN VĂN TỐT NGHIỆP Hiện tượng xảy cố hệ thống, cố thiết bị, điều khiển sai thiết bị Mất điện áp thường cố ngắn mạch hệ thống điện gây Sự điện đo thời gian xảy hệ thống, từ lúc biên độ điện áp nhỏ 10% định mức Thời gian điện xác định thời gian hoạt động thiết bị bảo vệ, thời gian không Hình 2.3: Mất điện tạm thời - Võng điện áp – Voltage Sag: Là giảm điện áp dòng điện hiệu dụng lại khoảng 0.1 – 0.9 pu thời gian từ 0.5 đến 30 chu kỳ Như hình 2.4 Hình 2.4: Võng điện áp Võng điện áp thường kèm theo cố hệ thống đóng điện tải nặng trình khởi động động lớn Các công ty điện, người sử dụng, nhà sản xuất thiết bị áp dụng nhiều biện pháp để giảm bớt cố gây võng điện áp giảm bớt ảnh hưởng lên thiết bị - Tăng điện áp – Voltage Swell: Là tăng dòng điện điện áp hiệu dụng đến giá trị khoảng 1.1 – 1.8 pu thời gian từ 0.5 đến 30 chu kỳ Như hình 2.5 Hình 2.5: Tăng điện áp TP HỒ CHÍ MINH 30 – 08 – 2005 PHỤ LỤC D MÃ NGUỒN CHƯƠNG TRÌNH: I CHƯƠNG TRÌNH PHÂN TÍCH WAVELETS DWT db4: Chương trình phân tích DWT 13 cấp, vẽ dạng sóng theo thời gian, vẽ phân bố lượng tính thời gian xảy cho loại độ Đây chương trình tổng quát cho loại tín hiệu: Pure Sin, Voltage Sag, Voltage Swell, Harmonic, Voltage Flicker, Voltage Interruption Capacitor Switching, biến liệu lệnh Load đặt X Muốn phân tích cho loại thay biến X tên tập tin chứa liệu Ví dụ: Pure Sin: X=sin01; Voltage Sag: X=SA01; Voltage Swell: X=SW01; Harmonic: X=HA10; Flicker: X=FL10; Voltage Interruption: X=IN10; Capacitor Switching: X=CL00 % nhap du lieu qua load X; s=vU; % phan tich wavelet db4 [C,L] = wavedec(s,13,'db4'); for i=1:13 P{i} = wrcoef('d',C,L,'db4',i);% dinh ly Parseval E(i)=sqrt((sum(P{i}.^2))/length(P{i}));% trich dac trung nang luong End % ve hinh cac thong so phan tich figure(1); subplot(4,1,1);plot(s);axis([1 5001 min(s) max(s)]); title('PHAN TICH WAVELET db4'); subplot(4,1,2);plot(abs(P{1}));axis([1 5001 max(abs(P{1}))]); ylabel('P1'); subplot(4,1,3);plot(abs(P{2}));axis([1 5001 max(abs(P{2}))]); ylabel('P2'); subplot(4,1,4);plot(abs(P{3}));axis([1 5001 max(abs(P{3}))]); ylabel('P3'); figure(2); bar(E);axis([0 14 ceil(max(E*10))/10]); xlabel('13 MUC NANG LUONG CUA X'); title('CAC MUC NANG LUONG'); % tinh thoi gian xay qua k=length(s)\5001; maxchuky2=0; temp=abs(P{3}); temp(1:10)=0; temp(k*4991:k*5001)=0; TP HOÀ CHÍ MINH 30 – 08 – 2005 PHỤ LỤC totaltime=k*(5001-ceil(5001/18/4)); begintime=ceil(k*5001/18/8); deltatime=ceil(k*5001/18/4); runtime=begintime; i=1; while (runtime2)|(length(pointmax)==3 & ( temp1>2 |temp2>2 ) ) | ( length(pointmax)==4 & temp1>2) tbegin=(begintime+(pointmax(1))*deltatime+find(temp(begintime+deltatime* (pointmax(1)1):begintime+deltatime*pointmax(1))==point(pointmax(1))))*0.3/5000; tend=(begintime+(pointmax(2)-1)*deltatime+find(temp(begintime+deltatime* (pointmax(2)-1) :begintime+deltatime*pointmax(2)) ==point(pointmax(2))))*0.3/5000; t=tend-tbegin; else pointmax=find(point>1/10*max(point)); pointm=find(point>0.8*max(point)); if length(pointmax)2 ) ) | ( length(pointmax)==4 & temp1>2 ) tbegin = (begintime+(pointmax(1)-1)*deltatime+find(temp(begintime+deltatime* (pointmax(1)-1):begintime+deltatime*pointmax(1))==point(pointmax(1))))* 0.3/5000; tend = (begintime+(pointmax(2)1)*deltatime+find(temp(begintime+deltatime*(pointmax(2)-1):begintime+ deltatime*pointmax(2))==point(pointmax(2))))*0.3/5000; t=tend-tbegin;%thoi gian xay qua else pointmax=find(point>1/10*max(point)); pointm=find(point>0.8*max(point)); if length(pointmax)