(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện

(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện(Luận văn thạc sĩ) Ứng dụng toán lai GA HS cho bài toán phân bố công suất trong hệ thống điện

LỜI CAM ĐOAN Tôi xin cam đoan đề tài nghiên cứu riêng Các số liệu, kết nêu luận văn trung thực chƣa đƣợc công bố đề tài khác Tôi xin cam đoan giúp đỡ cho việc thực luận văn đƣợc cảm ơn thơng tin trích dẫn luận văn đƣợc rõ nguồn gốc Học viên thực luận văn LÊ VŨ TRỌNG BẢO iii LỜI CÁM ƠN Tơi xin bày tỏ lịng biết ơn sâu sắc đến cá nhân giúp tơi hồn thành tốt thời hạn luận văn Thạc sĩ, cụ thể: - Thầy TS NGƠ CAO CƢỜNG – Trƣởng Phịng Tổ Chức Trƣờng Đại học Kỹ Thuật Công Nghệ Tp.HCM, ngƣời trực tiếp hƣớng dẫn luận văn - Thầy Th.s Lê Đình Lƣơng Thầy Th.s Nguyễn Vinh Quan tận tình giúp đỡ tơi thời gian thực Luận văn - Các thầy cô giảng dạy, truyền đạt kiến thức cho tơi q trình học tập nghiên cứu trình học cao học trƣờng Đại Học Sƣ Phạm Kỹ Thuật Thành Phố Hồ Chí Minh - Ban giám hiệu, Phịng Đào Tạo - Sau đại học Khoa Điện – Điện Tử Trƣờng Đại Học Sƣ Phạm Kỹ Thuật Thành Phố Hồ Chí Minh giúp đỡ, cung cấp thông tin tạo điều kiện thuận lợi cho tơi q trình học tập - Các Anh/Chị em, bạn bè đồng hành suốt trời gian thực luận văn TP.Hồ Chí Minh, tháng 04 năm 2013 Ngƣời thực Lê Vũ Trọng Bảo iv TÓM TẮT LUẬN VĂN Trong Luận văn này, phƣơng pháp lai đƣợc phát triển cách sử dụng thuật toán di truyền (GA) thuật tốn tìm kiếm hịa hợp (HS), đƣợc gọi tên GAHS Thuật toán dùng để giải toán điều độ kinh tế có xét hiệu ứng van điểm Thuật toán HS thuật toán meta-heuristic phát triển gần đây, thành công loạt vấn đề tối ƣu hóa Cách tìm kiếm giải pháp thuật toán HS dễ rơi vào giải pháp tối ƣu địa phƣơng giải vấn đề phức tạp Tuy nhiên, hoạt động lựa chọn, lai tạo đột biến thuật toán GA tạo giải pháp đa dạng, tránh đƣợc khả cực tiểu địa phƣơng Xác suất chiến lƣợc cập nhật dựa thuật toán HS cho nhớ hài hòa cao giai đoạn ban đầu để đẩy nhanh tốc độ tìm kiếm khả chiến lƣợc cập nhật dựa GA cao giai đoạn cuối để thoát khỏi giải pháp tối ƣu địa phƣơng Phƣơng pháp GA-HS đƣợc ứng dụng vào để giải toán Điều độ kinh tế với hiệu ứng điểm van Kết tính tốn cho thấy thuật tốn lai đề xuất có hiệu hội tụ tốt v ABSTRACT In this paper, a hybrid heuristic method is developed using the Genetic Algorithm (GA) and Harmony Search (HS), called GA/HS In this algorithm to solved the economic power dispatch with valve-point effects The HS algorithm is a recently developed meta-heuristic algorithm, and has been very successful in a wide variety of optimization problems Refresh strategy of solution vectors is so simple that HS algorithm is easy to fall into local optimum solution when complex problem is solved However, selection, crossover and mutation operations in GA can generate diverse solutions, so premature stagnation behavior can be avoided in some extent The probability of update strategy based on HS algorithm for harmony memory is higher during the initial stage in order to accelerate the searching speed, and the probability of update strategy based on GA is higher during the final stage in order to escape from local optimum solution Calculation results show that the proposed hybrid algorithm is effective and it converges well vi MỤC LỤC Trang tựa TRANG Quyết định giao đề tài LÝ LỊCH KHOA HỌC i LỜI CAM ĐOAN iii LỜI CÁM ƠN iv TÓM TẮT LUẬN VĂN v ABSTRACT vi MỤC LỤC vii DANH SÁCH CÁC CHỮ VIẾT TẮT x DANH SÁCH CÁC HÌNH xi DANH SÁCH CÁC BẢNG xii Chƣơng TỔNG QUAN 1.1TỔNG QUAN CHUNG 1.2MỤC TIÊU CỦA ĐỀ TÀI 1.3NHIỆM VỤ ĐỀ TÀI VÀ GIỚI HẠN ĐỀ TÀI 1.4PHƢƠNG PHÁP NGHIÊN CỨU Chƣơng2 BÀI TỐN ĐIỀU PHỐI CƠNG SUẤT ED 2.1GIỚI THIỆU 2.2BÀI TOÁN ĐIỀU PHỐI KINH TẾ CỔ ĐIỂN 2.2.1Hàm mục tiêu 2.2.2Ràng buộc đẳng thức 2.2.3Ràng buộc bất đẳng thức 2.3BÀI TOÁN ĐIỀU PHỐI KINH TẾ VỚI HÀM CHI PHÍ NHIÊN LIỆU KHƠNG TRƠN 2.3.1.Đặc điểm toán điều phối kinh tế với điểm van công suất 2.3.2.Biểu thức điều phối kinh tế với điểm van công suất Chƣơng THUẬT TOÁN DI TRUYỀN 10 3.1GIỚI THIỆU CHUNG 10 vii 3.2TƢƠNG QUAN GIỮA THUẬT TỐN GEN VÀ Q TRÌNH TIẾN HĨA CỦA SINH VẬT 10 3.3THUẬT TOÁN DI TRUYỀN 11 3.3.1.Cấu trúc tổng quát thuật toán di truyền 12 3.3.2.Các thành phần thuật toán di truyền 13 3.3.3.Các phƣơng pháp mã hóa 13 3.3.4.Xác định độ phù hợp (Fitness) 15 3.3.5.Điều chỉnh độ phù hợp (Fitness scaling) 15 3.3.5.1.Phƣơng pháp ánh xạ tuyến tính 15 3.3.5.2.Phƣơng pháp xác định độ phù hợp theo thứ tự giá trị hàm mục tiêu (Ranking) 16 3.3.6.Tiêu chuẩn ngừng lặp 17 3.4CÁC TOÁN TỬ DI TRUYỀN 17 3.4.1.Chọn lọc cá thể (Selection) 17 3.4.1.1.Phƣơng pháp bánh xe Roulette 17 3.4.1.2.Phƣơng pháp lấy mẫu ngẫu nhiên 19 3.4.2 Lai tạo 20 3.4.2.1.Toán tử lai cho phƣơng pháp mã hóa nhị phân thập phân 20 3.4.2.2.Tốn tử lai cho biến điều khiển biểu diễn số thực 20 3.4.3.Đột biến 21 Chƣơng THUẬT TỐN TÌM KIẾM SỰ HÀI HỊA (HS) 23 4.1GIỚI THIỆU CHUNG 23 4.2GIẢI THUẬT CỦA THUẬT TOÁN HS 25 4.3ƢU VÀ NHƢỢC ĐIỂM THUẬT TOÁN HS 29 4.3.1 Ƣu điểm HS 29 4.3.2Nhƣợc điểm HS 29 4.4 PHÂN TÍCH VÍ DỤ 29 Chƣơng THUẬT TOÁN LAI GA-HS VÀ ỨNG DỤNG GIẢI BÀI TOÁN ED VỚI HIỆU ỨNG ĐIỂM VAN 31 viii 5.1.THUẬT TOÁN LAI GA-HS 31 5.1.1Giới thiệu 31 5.1.2 Các bƣớc phƣơng pháp lai thuật toán GA – HS 32 5.2LƢU ĐỒ GIẢI THUẬT 34 5.2.1Lƣu đồ HS cho toán ED với hiệu ứng điểm van 34 5.2.2Lƣu đồ GA-HS cho toán ED với hiệu ứng điểm van 36 5.3ỨNG DỤNG THUẬT TOÁN LAI GA-HS GIẢI BÀI TOÁN ED VỚI HIỆU ỨNG ĐIỂM VAN 39 5.3.1Hệ thống 13 nút 39 5.3.2Hệ thống 40 nút 42 Chƣơng KẾT LUẬN VÀ HƢỚNG NGHIÊN CỨU PHÁT TRIỂN 47 6.1KẾT LUẬN 47 6.2HƢỚNG NGHIÊN CỨU VÀ PHÁT TRIỂN 47 6.3LỜI KẾT 48 TÀI LIỆU THAM KHẢO 49 PHỤ LỤC 52 ix DANH SÁCH CÁC CHỮ VIẾT TẮT DE Differential Evolution ED Economic Dispatch fw Fret Width GA Genetic algorithm HM Harmony Memory HMCR Harmony Memory Considering Rate HMS Harmony Memory Size HS Harmony Search IFEP Improved Fast Evolutionary Program IHS Improved Harmony Search LP Linear programming MPPSO Mean personal-best oriented PSO NLP Nonlinear programming PAR Pitch Adjusting Rate PPSO Personal best-oriented PSO PSO Particle Swarm Optimization x DANH SÁCH CÁC HÌNH HÌNH TRANG Hình 2.1: Đƣờng cong chi phí phổ biến nhà máy nhiệt điện Hình 2.2: Hàm chi phí nhiên liệu nhà máy nhiệt điện với van nạp Hình 3.1: Cơ chế hoạt động thuật tốn di truyền 11 Hình 3.2: Lƣu đồ giải thuật thuật toán di truyền 12 Hình 3.3: Các khơng gian làm việc thuật tốn di truyền 14 Hình 4.1: Harmony Search mối tƣơng đồng hài hòa âm nhạc tối ƣu kỹ thuật 23 Hình 4.2: Sơ đồ khối đơn giản thuật toán HS 26 Hình 4.3: Sơ đồ chi tiết thuật toán HS 28 Hình 4.4: Đồ thị hàm số Himmelblau 29 Hình 4.5: Đồ thị kết hàm Himmelblau 30 Hình 5.1: Lƣu đồ giải thuật HS cho tốn Valve point 36 Hình 5.2: Lƣu đồ giải thuật GA - HS cho toán Valve point 37 Hình 5.3: Đặc tính hội tụ mạng điện 13 máy phát 41 Hình 5.4: Đặc tính hội tụ mạng điện 40 máy phát 45 xi DANH SÁCH CÁC BẢNG BẢNG TRANG Bảng 5.1: Thông số hệ thống 13 nút với hiệu ứng điểm van 39 Bảng 5.2: Công suất phát hệ thống 13 máy phát 40 Bảng 5.3: Kết so sánh hệ thống 13 máy phát có xét ảnh hƣởng điểm van công suất 41 Bảng 5.4: Thông số hệ thống 40 nút với hiệu ứng điểm van 42 Bảng 5.5: Công suất phát hệ 40 máy phát 43 Bảng 5.6: Kết so sánh hệ thống 40 máy phát có xét ảnh hƣởng điểm van công suất 45 xii Chƣơng KẾT LUẬN VÀ HƢỚNG NGHIÊN CỨU PHÁT TRIỂN 6.1 KẾT LUẬN Phƣơng pháp tối ƣu GA-HS dự sở sử dụng hai phƣơng pháp tối ƣu HS GA, gọi phƣơng pháp lai GA-HS Phƣơng pháp đƣợc tạo thành dự sở thuật toán tìm kiếm hịa hợp (HS)[3] sử dụng thuật tốn GA để tạo giá trị cá thể Mục đích việc đƣa thuật tốn GA[2] vào bƣớc thuật toán HS nhằm tăng cƣờng khả thăm dò để đẩy nhanh tốc độ hội tụ thuật toán HS Thuật toán GA tạo giải pháp đa dạng tránh đƣợc điểm cực tiểu địa phƣơng để tiến đến giải pháp tối ƣu toàn cục nhanh Từ kết đạt đƣợc, cho thấy thuật toán lai GA- HSđã giải tốt toán điều độ kinh tế cơng suất có ảnh hƣởng điểm van Trongphƣơng pháp cải tiếnmới, thuật toánGA-HS đƣợc sử dụng với ƣu điểm hai thuật toán GA HS phƣơng pháp làm tăng tốc độhộitụvới lời giảitồn cụcmột cách nhanh chóngkhơng phụ thuộc vàohình dạngcủa hàm chi phí Tuy cịn hạn chế mặt thời gian tính tốn so với phƣơng pháp tối ƣu khác Nhƣng kết tính tốnđã rarằng thuật tốnGA-HSđãthể ƣu hơnnhiều phƣơng pháp kháctrong tài liệu tham khảo việc giải quyếtcácbài toán điều độ kinh tế(ED) có ảnh hƣởng điểm van cơng suấtcho tất cáctrƣờng hợp thử nghiệmtrongđiều kiện củatổng giá trịhàm chi phí tính tốn 6.2 HƢỚNG NGHIÊN CỨU VÀ PHÁT TRIỂN Thuật toán lai GA-HS chứng tỏ đƣợc khả tốn tối ƣu hóa ED hệ thống điện Từ kết đạt đƣợc trên, xin đƣa hƣớng phát triển đề tài: 47  Giải tốn phân bố cơng suất tối ƣu mạng điện có số lƣợng nút lớn  Cải thiện thời gian tính tốn thuật tốn  Nghiên cứu thuật toán khác để đƣa đến ý tƣởng với kinh nghiệm từ phƣơng pháp lai GA-HS để tạo phƣơng pháp tối ƣu khác  Áp dụng thuật toán lai GA-HS vào giải toán khác nhƣ: toán ED với hàm chi phí có xét vùng cấm tổn hao, toán quy hoạch phát triển hệ thống điện, tốn phân bố cơng suất tối ƣu thị trƣờng điện cạnh tranh, tốn phân bố cơng suất tối ƣu hệ thống Thủy điện- Nhiệt điện, 6.3 LỜI KẾT Luận văn tốt nghiệp xem cơng trình khoa học quan trọng tơi Trong q trình tìm hiểu nghiên cứu thuật tốn tối ƣu HS GA, tơi có đƣợc đam mê Tơi cố gắng hồn thành tốt luận văn Tuy nhiên tránh đƣợc sai sót, nhầm lẫn mong chia đóng góp Hội đồng nhƣ ngƣời đọc để tơi hồn chỉnh tƣơng lai Xin chân thành cảm ơn! 48 TÀI LIỆU THAM KHẢO TIẾNG VIỆT [1] PGS.TS Quyền Huy Ánh, Vận Hành Tối Ưu Hệ Thống Điện Trƣờng Đại Học Sƣ Phạm Kỹ Thuật Tp.Hồ Chí Minh, 01/2012 [2] Phạm Việt Cƣờng,Ứng dụng thuật tốn di truyền phân bố cơng suất tối ưu hệ thống điện Đại học Bách Khoa Tp Hồ Chí Minh, Luận văn Thạc sĩ, 07/2003, 700998 TIẾNG NƢỚC NGOÀI [3] Zong Woo Geem, “Optimal Scheduling of Multiple Dam System Using Harmony Search Algorithm” Johns Hopkins University, Environmental Planning and Management Program, 729 Fallsgrove Drive #6133, Rockville, Maryland 20850, USA [4] Kang Seok Lee, Zong Woo Geem, “A new meta-heuristic algorithm for continuous engineering optimization: harmony search theory and practice” Comput Methods Appl Mech Engrg 194 (2005) 3902–393 [5] A Vasebi, M Fesanghary, S.M.T Bathaee, “Combined heat and power economic dispatch by harmony search algorithm” Electrical Power and Energy Systems 29 (2007) 713–719 [6] Leandro dos Santos Coelho, Viviana Cocco Mariani, “An improved harmony search algorithm for power economic load dispatch” Energy Conversion and Management 50 (2009) 2522–2526 [7] Nekooei, M.M Farsangi, H Nezamabadi-pour, “An Improved Harmony Search Approach To Economic Dispatch” K ISSN 2077-3528, IJTPE Journal September 2011, Issue 8, Volume 3, Number 3, Pages 25-31 [8] Tarek Bouktir, Linda Slimani, M Belkacemi, “A Genetic Algorithm for Solving the Optimal Power Flow Problem”, Leonardo Journal of Sciences, Issue 4, January-June 2004, p 44-58 49 [9] Chao-Lung Chiang, “Improved Genetic Algorithm for Power Economic Dispatch of Units With Valve-Point Effects and Multiple Fuels” IEEE Transactions On Power Systems, Vol 20, No 4, November 2005 [10] PhanTu Vu, DinhLuong Le, NgocDieu Vo and Josef Tlust, “A Novel Weight-Improved Particle Swarm Optimization Algorithm for Optimal Power Flow and Economic Load Dispatch Problems” [11]Wei -Wei Shi, Wei Han, Wei -Chao Si, “A Hybrid Genetic Algorithm Based on Harmony Search and its Improving” Informatics and Management Science ILecture Notes in Electrical Engineering Volume 204, 2013, pp 101-109 [12] Li Liu, Hongtao Yu, Lixia Li, “Distribution Network Reconfiguration Based on Harmony Search/Genetic Hybrid Algorithm” Shenyang Institute of Engineering,2012 China International Conference on Electricity Distribution (CICED 2012) Shanghai, 5-6 Sep 2012 [13] Khang Nguyen, Phuc Nguyen and Nuong Tran “A hybrid algorithm of Harmony Search and Bees Algorithm for a University Course Timetabling Problem” International Journal of Computer Science Issues (IJCSI), Vol 9, Issue 1, No 1, January 2012 ISSN (Online): 1694-0814 [14] A Kaveh and A Nasrollahi “Engineering Design Optimization Using A Hybrid PSO And HS Algorithm” Center of Excellence for Fundamental Studies in Structural Engineering, Iran University of Science and Technology, Narmak, Tehran, Iran Sian Journad Of Civil Engineering (BHRC) Vol 14, No (2013) Pages 201-223 [15] N Sinha, R Chakrabarti, and P K Chattopadhyay, “Improved Fast Evolutionary Program for Economic Load Dispatch with Non-smooth Cost Curves”, IE (I) Journal-EL [16] C H Chen, and S N Yeh, “Particle Swarm Optimization for Economic Power Dispatch With Valve-Point Effects” Senior Member, IEEE 2006 IEEE PES Tranmission and Distribution Conference and Eaposition Latin America, Venezuela 50 [17]M Mahdavi, M Fesanghary, E Damangir, “An improved harmony search algorithm for solving optimization problems”, Applied Mathematics and Computation 188 (2007) 1567–1579 [18] N Sinha, R Chakrabarti, and P K Chattopadhyay, “Evolutinary programming techniques for economic load dispatch”, IEEE Trans on Evolutionary Compution, Vol 7, No.1, Feb., pp.83-94, 2003 51 PHỤ LỤC I THÔNG SỐ CÀI ĐẶT CHO LẬP TRÌNH THUẬT TỐN PD = 10500; Ploss = PD = 1800; Ploss = Hệ thống 13 máy phát HS GA-HS 20 20 500.000 50.000 Thông số cài đặt HMS MaxImp Hệ thống 40 máy phát HS GA-HS 40 40 800.000 500.000 HCMR 0.85 0.85 0.95 0.95 PAR b Kích thƣớc quần thể Số hệ 0.45 (pU-pL)./1000 0.45 (pU-pL)./1000 0.55 (pU-pL)./10000 0.45 (pU-pL)./10000 15 40 30 30 II KẾT QUẢ TÍNH TỐN BẰNG PHƢƠNG PHÁP HS Thuật tốn HS đƣợc áp dụng vào tính tốn mạng điện khác bao gồm: hệ thống 13 máy phát với ảnh hƣởng điểm van, hệ thống 40 máy phát với ảnh hƣởng điểm van (thông số hệ thống đƣợc trình bày Bảng 5.1 Bảng 5.40 Đối với hệ thống, thuật toán đề xuất đƣợc tiến hành chạy thử 50 lần độc lập (chƣơng trình MatLab R2012a (7.14.0.739) máy tính có cấu hình CORE i5 1.7Ghz) Hệ thống 13 nút xét hiệu ứng điểm van: Máy phát Pi,min Pi,max Công suất (MW) 680 610.1237 360 148.8922 360 223.4919 60 180 109.8640 60 180 60.0000 60 180 109.8362 60 180 109.8226 52 60 180 109.7852 60 180 128.1458 10 40 120 40.0000 11 40 120 40.0000 12 55 120 55.0385 13 55 120 55.0000 Tổng công suất phát (MW) 10500.0000 Chi phí cực tiểu ($/h) 17964.7989 Chi phí cực đại ($/h) 18065,5947 Chi phí trung bình ($/h) 17993.6875 Độ lệch tiêu chuẩn ($/h) 11,07665 Thời gian tính tốn (s) 15 1.9 x 10 1.88 Total cost($) 1.86 1.84 1.82 1.8 1.78 0.5 1.5 2.5 Iterations 3.5 4.5 x 10 Hình 1: Đặc tính hội tụ mạng điện 13 máy phát (HS) 53 Hệ thống 40 nút xét hiệu ứng điểm van: Máy phát Máy phát Pi,min Pi,max 21 254 550 Công suất (MW) 523.2745 110.9917 22 254 550 523.2749 120 97.4078 23 254 550 523.2707 80 190 179.7427 24 254 550 523.2832 47 97 96.9153 25 254 550 523.2866 68 140 139.9902 26 254 550 523.3347 110 300 259.6243 27 10 150 10.0040 135 300 284.6256 28 10 150 10.0000 135 300 284.6290 29 10 150 10.0000 10 130 300 204.7997 30 47 97 93.6602 11 94 375 168.8010 31 60 190 179.4145 12 94 375 284.8912 32 60 190 95.7411 13 125 500 304.4778 33 60 190 190.0000 14 125 500 304.5208 34 90 200 164.8952 15 125 500 394.2747 35 90 200 164.8691 16 125 500 304.5606 36 90 200 91.9741 17 220 500 489.2791 37 25 110 109.9950 18 220 500 489.3013 38 25 110 109.9986 19 242 550 511.3030 39 25 110 110.0000 20 242 550 511.2755 40 242 550 511.2593 Pi,min Pi,max 36 114 36 114 60 Công suất (MW) 87.0531 Tổng công suất phát (MW) Chi phí cực tiểu ($/h) 10,500 121,424.37 Chi phí cực đại ($/h) 121,550.35 Chi phí trung bình ($/h) 121,520.08 Độ lệch tiêu chuẩn ($/h) 34,07665 Thời gian tính tốn (s) 54 54 x 10 2.8 2.6 Total cost($) 2.4 2.2 1.8 1.6 1.4 1.2 0.5 1.5 2.5 Iterations 3.5 4.5 5 x 10 Hình 2: Đặc tính hội tụ mạng điện 40 máy phát (HS) III CODE TỐI ƢU HÀM HIMMELBLAU Gồm có file sau: Script (runable) files: hs.m Thể thông số thuật tốn HS tối ƣu hóa hàm hs_plot.m Thể thơng số thuật toán HS với kết đồ MHB thị hội tụ Function (subroutine) files: harmony.m Tối ƣu hóa HS cổ điển harmony_plot.m Tối ƣu hóa HS cổ điển với kết là: lịch sử hội tụ, giá trị tốt hàm mục tiêu MHB.m (Modified Himmelblau Himmelblau Code file MHB.m function f = MHB(x) % Định nghĩa hàm Himmelblau 55 function) : Định nghĩa hàm f = (x(1)^2 + x(2) - 11)^2 + (x(1) + x(2)^2 - 7)^2 + 0.1*( (x(1)-3)^2 + (x(2)-2)^2 ); Code file hs.m % HS: Harmony Search tối ưu % Use the Matlab startup random number sequence: rand('twister',5489); % Commment out to start rand from current position % Xác định hàm mục tiêu giới hạn f = @MHB; xL = [-6 -6]; xU = [6 6]; % thiết lập thơng số thuật tốn HSparams.HMS = 30; HSparams.MaxImp = 100000; HSparams.HMCR = 0.95; HSparams.PAR = 0.3; HSparams.b = (xU-xL)/100000; % Thực tối ưu [xbest,fbest] = harmony(f,xL,xU,HSparams); fprintf('Best solution found:\n') disp(xbest) fprintf('Function value = %e', fbest) Code file harmony.m function [xbest,fbest] = harmony(f,xL,xU,HSparams) % HARMONY tối ưu dùng HS với biến liên tục % thành phần vector tham số HMS = HSparams.HMS; MaxImp = HSparams.MaxImp; HMCR = HSparams.HMCR; PAR = HSparams.PAR; b = HSparams.b; % kích N = HM = F = xnew = thước mảng length(xL); % Số lượng biến định zeros(HMS,N); zeros(HMS,1); zeros(1,N); % Ngẫu nhiên khởi tạo ma trận HM,đảm bảo để biến giới hạn % Đánh giá hàm mục tiêu tương ứng for j = 1:HMS HM(j,:) = xL + (xU-xL).*rand(1,N); 56 F(j) = f(HM(j,:)); end % Vòng ngẫu hứng qua MaxImp for j = 1:MaxImp % Ngẫu hứng tạo HM mới: vòng lặp qua biến for i = 1:N % ngẫu nhiên thực hoạt động với vectorxnew(i) HM if rand < HMCR % Bộ nhớ xem xét lựa chọn ngẫu nhiên vector lưu trữ HM xnew(i) = HM(ceil(rand*HMS),i); if rand < PAR %Hiệu chỉnh cường độ: điều chỉnh ngẫu nhiên cường độ nhẹ % within +/- b(i), đảm bảo giới hạn thỏa mãn xnew(i) = xnew(i) + (2*rand-1)*b(i); xnew(i) = min(max(xnew(i),xL(i)),xU(i)); end else % Random playing: ngẫu nhiên lựa chọn điều chỉnh giới hạn xnew(i) = xL(i) + rand*(xU(i)-xL(i)); end end % kết thúc trinh ngẫu hứng % HM update: kiểm tra hài hòa vector có tốt cũ tệ % hài hịa có HM fnew = f(xnew); [fworst,idxworst] = max(F); if fnew < fworst HM(idxworst,:) = xnew; F(idxworst) = fnew; end end % số lượng tối đa ngẫu hứng đạt % trở lại ngẫu hứng tốt đạt [fbest,idxbest] = min(F); xbest = HM(idxbest,:); Code file hs_plot.m % HS_PLOT: Harmony Search minimization with functionplots % Use the Matlab startup random number sequence: 57 rand('twister',5489); % Commment out to start rand from current position % thiết lập hàm mục tiêu giới hạn f = @MHB; xL = [-6 -6]; xU = [6 6]; % Thiết lập thơng số thuật tốn HSparams.HMS = 10; HSparams.MaxImp = 10000; HSparams.HMCR = 0.8; HSparams.PAR = 0.4; HSparams.b = (xU-xL)/1000; % Tối ưu hóa [xbest,fbest,xnew_history,fbest_history] = harmony_plot(f,xL,xU,HSparams); % kết đồ thị % đồ thị có dánh dấu vị trí với lần lặp x1 = linspace(xL(1),xU(1)); x2 = linspace(xL(2),xU(2)); [X1,X2] = meshgrid(x1,x2); F = zeros(length(x1),length(x2)); for i=1:length(x1) for j=1:length(x2) F(i,j) = f([X1(i,j),X2(i,j)]); end end figure contour(X1,X2,F,75) xlabel('{\itx}_1') ylabel('{\itx}_2') hold on % Show iterations ImpStep = 100; for i = 1:ImpStep:HSparams.MaxImp c = [0.7 0.7 0.7]*(1 - i/HSparams.MaxImp); % Vary color gray to black plot(xnew_history(i:(i+ImpStep1),1),xnew_history(i:(i+ImpStep-1),2), '+','MarkerEdgeColor',c,'markersize',3) pause(0.01) end % đồ thị hội tụ figure 58 semilogy(fbest_history) xlabel('{\itImp}') ylabel('{\itf}({\bfx_{best}})') set(gca,'YGrid','on') Code file harmony_plot.m: function [xbest,fbest,xnew_history,fbest_history] = harmony_plot(f,xL,xU,HSparams) % HARMONY_PLOT basic Harmony Search minimization for continuous variables, % ghi lại vector giá trị tốt lần lặp % thành phần vector tham số HMS = HSparams.HMS; MaxImp = HSparams.MaxImp; HMCR = HSparams.HMCR; PAR = HSparams.PAR; b = HSparams.b; % Kích thước mảng N = length(xL); % Số lượng biến định HM = zeros(HMS,N); F = zeros(HMS,1); xnew = zeros(1,N); xnew_history = zeros(MaxImp,N); fbest_history = zeros(MaxImp,1); % Ngẫu nhiên khởi tạo ma trận HM,đảm bảo để biến giới hạn % Đánh giá hàm mục tiêu tương ứng for j = 1:HMS HM(j,:) = xL + (xU-xL).*rand(1,N); F(j) = f(HM(j,:)); end % Vòng ngẫu hứng qua MaxImp for j = 1:MaxImp % Ngẫu hứng tạo HM mới: vòng lặp qua biến for i = 1:N % ngẫu nhiên thực hoạt động với vector xnew(i) HM if rand < HMCR % Bộ nhớ xem xét lựa chọn ngẫu nhiên vector lưu trữ HM xnew(i) = HM(ceil(rand*HMS),i); 59 if rand < PAR % Hiệu chỉnh cường độ: điều chỉnh ngẫu nhiên cường độ nhẹ % within +/- b(i), đảm bảo giới hạn thỏa mãn xnew(i) = xnew(i) + (2*rand-1)*b(i); xnew(i) = min(max(xnew(i),xL(i)),xU(i)); end else % Random playing: ngẫu nhiên lựa chọn điều chỉnh giới hạn xnew(i) = xL(i) + rand*(xU(i)-xL(i)); end end % Kết thúc ngẫu hứng hài hòa % Ghi lại hài hòa xnew_history(j,:) = xnew; % HM update: kiểm tra hài hịa vector có tốt cũ tệ % hài hòa fnew = f(xnew); [fworst,idxworst] = max(F); if fnew < fworst HM(idxworst,:) = xnew; F(idxworst) = fnew; end % Ghi lại giá trị hàm số tốt fbest_history(j) = min(F); end % Maximum number of improvisations reached % Trở lại hài hòa tốt [fbest,idxbest] = min(F); xbest = HM(idxbest,:); 60 ... 34 5.2.1Lƣu đồ HS cho toán ED với hiệu ứng điểm van 34 5.2.2Lƣu đồ GA- HS cho toán ED với hiệu ứng điểm van 36 5. 3ỨNG DỤNG THUẬT TOÁN LAI GA- HS GIẢI BÀI TOÁN ED VỚI HIỆU ỨNG ĐIỂM VAN ... 5.3 ỨNG DỤNG THUẬT TOÁN LAI GA- HS GIẢI BÀI TOÁN ED VỚI HIỆU ỨNG ĐIỂM VAN Bài toán ED liên quan đến việc tối ƣu kế hoạch phát công suất máy phát hệ thống điện để tối thiểu hóa chi phí đáp ứng. .. cơng suất mức công suất phát Điều phối công suất kỹ thuật để điều chỉnh biến điều khiển phân phối công suất cho hệ thống điện hoạt động tối ƣu Điều phối công suất có hai cách: điều phối cơng suất