công suất có xét đến ràng buộc an ninh:
Dựa vào ý tưởng cải tiến phương pháp tối ưu bầy đàn bằng kỹ thuật gradient giả đã trình bày trong chương 2 và kỹ thuật áp dụng thuật toán tối ưu vào bài toán phân bố tối ưu công suất như trình bày trong phần trên, ta xây dựng thuật toán hoàn chỉnh như sau:
- Bƣớc 1: Khởi tạo các giá trị điều khiển cho thuật toán tối ưu bầy đàn cải tiến: Số lượng cá thể NP = 20, số vòng lặp tối đa: Maxiter = 250, các hằng số gia tốc c1 = c2 = 2.05, các số phạt Kp, Kq, Kv, Ks được gán bằng 106.
- Bƣớc 2: Khởi tạo vị trí và tốc độ ban đầu cho các thuật toán tối ưu bầy đàn theo công thức (4.3), (4.4).
- Bƣớc 3: Tính giá trị ban đầu theo hàm Fitness của từng phần tử. Ta sử dụng phương pháp Newton-Raphson để tính phân bố công suất của hệ thống theo từng cá thể. Từ đó xác định giá trị hàm Fitness dựa theo biểu thức (4.1). Giá trị ban đầu này của hàm Fitness được gán vào biến FTPbest. - Bƣớc 4: Dùng biến Pbest để lưu trữ vị trí ban đầu của các phần tử và biến
Gbest để lưu trữ vị trí của phần tử tốt nhất dựa vào giá trị của hàm Fitness min
Gbest Pbest
FT FT .
- Bƣớc 5: Gán giá trị ban đầu cho các gradient giả bằng 0 và khởi tạo biến đếm k = 1.
- Bƣớc 6: Tính toán giá trị mới của vị trí và vận tốc từng cá thể dựa vào biểu thức (3.18), (3.19) và (3.20).
- Bƣớc 7: Dùng phương pháp Newton-Raphson để tính phân bố công suất cho từng phần tử.
- Bƣớc 8: Tính lại giá trị hàm Fitness cho từng phần tử và so sánh với
FTPbest và FTGbest để cập nhật lại vị trí tốt nhất của từng phần tử và vị trí
- Bƣớc 9: Xác định giá trị gradient giả cho từng cá thể theo biểu thức (3.15), (3.16) và (3.17).
- Bƣớc 10: Kiểm tra nếu k < Maxiter thì k = k + 1 và quay lại bước 6, ngược lại thì dừng vòng lặp và xuất kết quả.
Để giải quyết cụ thể cho bài toán phân bố tối ưu công suất có xét đến ràng buộc an ninh, ta thực hiện theo trình tự theo 3 phần như sau:
Phần 1: Giải bài toán phân bố tối ưu công suất trong trường hợp bình