(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện(Luận văn thạc sĩ) Ứng dụng giải thuật di truyền để tính toán tối ưu dung lượng bù cho hệ thống điện
LỜI CAM ĐOAN Tơi cam đoan cơng trình nghiên cứu Các số liệu, kết nêu luận văn trung thực chưa cơng bố cơng trình khác TP Hồ Chí Minh, ngày tháng Học viên (Ký tên ghi rõ họ tên) NGUYỄN TẤN HÙNG ii năm 2018 LỜI CẢM ƠN Qua trình học tập, nghiên cứu hồn thành luận văn, em kính gửi lời cảm ơn chân thành sâu sắc đến: PGS.TS TRƯƠNG ĐÌNH NHƠN thầy/cơ tận tình dạy, tạo điều kiện động viên em suốt q trình thực Q thầy, giáo tham gia công tác giảng dạy, hướng dẫn em thành viên lớp Cao học chuyên ngành Kỹ Thuật Điện KDD17B tồn khố học Q thầy, cô giáo giảng dạy khoa Điện - Điện Tử, Phòng Đào tạo – phận sau đại học Trường Đại học Sư phạm Kỹ thuật Tp Hồ Chí Minh giúp đỡ em thực thời gian học tập nghiên cứu trường Kính gửi lời cảm tạ tới Ban Giám Hiệu Trường Đại học Sư phạm Kỹ thuật Tp Hồ Chí Minh tạo điều kiện thuận lợi cho học viên trường học tập nghiên cứu Kính chúc Quý thầy, cô giáo thật nhiều sức khỏe Thành phố Hồ Chí Minh, ngày tháng Học viên thực NGUYỄN TẤN HÙNG iii năm 2018 TÓM TẮT Đề tài “ ỨNG DỤNG GIẢI THUẬT DI TRUYỀN ĐỂ TÍNH TỐN TỐI ƯU DUNG LƯỢNG BÙ CHO HỆ THỐNG ĐIỆN ” tiến hành khoảng thời gian năm Trường Đại Học Sư Phạm Kỹ Thuật TP HCM Sau thời gian nghiên cứu đề tài triển khai tập trung giải vấn đề sau: Tìm hiểu thuật tốn di truyền (GA) áp dụng tính toán tối ưu dung lượng bù cho hệ thống điện dựa nguyên tắc hệ thống mạch vòng yêu cầu vận hành hình tia, nhằm mục đích giảm tổn thất công suất, tổn thất điện áp, nâng cao chất lượng điện hệ thống điện Tìm hiểu thiết bị bù STATCOM ứng dụng hệ thống điện Mô hệ thống điện 14-bus ứng dụng thuật tốn di truyền (GA) bù cơng suất phản kháng Mô hệ thống điện 14-bus ứng dụng thuật tốn di truyền (GA) bù cơng suất phản kháng có kết hợp với thiết bị bù STATCOM Kiểm tra so sánh độ xác thuật toán lưới điện mẫu IEEE (hệ thống 14-bus) Học viên thực NGUYỄN TẤN HÙNG iv ABSTRACT Research “Applying genetic algorithm to optimize the amount of capacity compensated for the power system” was conducted at Ho Chi Minh University of Technology and Education within the period of year The research was subsequently deployed and applied to solve the following issues: Understanding genetic algorithm (GA) and its application in optimizing the amount of capacity compensated for the power system by using the principles of loop system and x-ray tube operation is aimed to reduce energy loss, voltage drop and improve the quality of electric power system Learning about the static synchronous compensator (STATCOM) and its application in the power system Simulating 14-bus power system which is applied genetic algorithm (GA) to compensate reactive power Simulating 14-bus power system which is applied genetic algorithm (GA) to compensate reactive power and combined with the static synchronous compensator (STATCOM) Checking and comparing the precision of the algorithm on the grid pattern of IEEE (14-bus power system) Author NGUYỄN TẤN HÙNG v MỤC LỤC TRANG Trang tựa Quyết định giao đề tài Lý lịch khoa học ……………………………………………………………… .i Lời cam đoan…………………………………………………………………… ii Lời cảm ơn……………………………………………………………………… iii Tóm tắt………………………………………………………………………….…iv Mục lục ………………………………………………………………………… vi Danh mục chữ viết tắt……………………………………………………… vii Danh sách bảng …………………………………………………….……… viii Danh sách hình………………………………………………………….…… ix Chương 1: TỔNG QUAN………………………………………………… … 1.1 Giới thiệu .1 1.2 Đặt vấn đề 1.3 Mục tiêu đề tài 1.4 Nhiệm vụ giới hạn đề tài 1.5 Phương pháp nghiên cứu 1.6 Nội dung đề tài .5 Chương 2: CƠ SỞ LÝ THUYẾT……………………………………… ………6 2.1 Mục tiêu lợi ích bù cơng suất phản kháng 2.2 Ổn định điện áp hệ thống điện………………………………………6 2.3 Các giới hạn ổn định hệ thống điện …………… ……………….7 2.3.1 Giới hạn điện áp ……………….………………………………… 2.3.2 Giới hạn nhiệt……………………………………………………… 2.3.3 Giới hạn ổn định………….…………………………… ………… Chương 3: ỨNG DỤNG GIẢI THUẬT DI TRUYỀN ĐỂ TÍNH TỐN TỐI ƯU DUNG LƯỢNG BÙ CHO HỆ THỐNG ĐIỆN……………………… … 14 3.1 Giới thiệu giải thuật di truyền (GA) …………………………………….14 vi 3.2 Giới thiệu hệ thống điện 14-bus 15 3.3 Lưu đồ giải thuật ……………………… ………………………….…….16 Chương 4: BỘ BÙ ĐỒNG BỘ TĨNH STATCOM VÀ ỨNG DỤNG TRONG HỆ THỐNG ĐIỆN ………………………………………… …….18 4.1 Giới thiệu STATCOM 18 4.2 Cấu trúc nguyên lí hoạt động STATCOM 20 4.2.1 Cấu trúc STATCOM…………………………………20 4.2.2 Nguyên lí hoạt động STATCOM………………………… 21 4.3 Kết mơ .24 Chương 5: KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN ………………….……36 TÀI LIỆU THAM KHẢO………………………………………………………37 PHỤ LỤC……………… ………………………………………………………38 vii DANH MỤC CÁC CHỮ VIẾT TẮT GO: Global Optimization GA: Genetic Algorithms FACTS: Flexible AC Tranmission System STATCOML: Static Synchronous Compensator SVC: Static Var Compensator VSC: Voltage Source Converter IGBT: Insulated gate bipolar transistor PWM: Pulse Width Modulated PI: Proportional integral viii DANH SÁCH CÁC BẢNG Bảng 3.1: Nhiệm vụ bước ……………………………………………… 17 ix DANH SÁCH CÁC HÌNH Hình 2.1: Phân loại ổn định hệ thống điện ………………………………… Hình 2.2: Hệ thống điện …………………………………………………… …….9 Hình 2.3: Đường cong cơng suất - góc 10 Hình 2.4: Sự thay đổi hệ thống ổn định độ ổn định 11 Hình 2.5: Độ thay đổi góc hệ thống ổn định dao động bé (a), hệ thống ổn định dao động (b) hệ thống ổn định (c) 12 Hình 2.6: Giới hạn vận hành đường dây theo mức điện áp … ……… …13 Hình 3.1: Sơ đồ đơn tuyến hệ thống 14 - bus 15 Hình 3.2: Lưu đồ giải thuật 16 Hình 4.1: Mạch điện tương đương STATCOM 20 Hình 4.2: Cấu trúc STATCOM 21 Hình 4.3: Ngun lí hoạt động STATCOM 21 Hình 4.4: Nguyên lí bù bù tích cực 22 Hình 4.5: Trạng thái phát công suất phản kháng bù 23 Hình 4.6: Sơ đồ hệ thống điện 14 nút xây dựng Matlab 24 Hình 4.7: Giá trị điện áp nút trước bù 25 Hình 4.8: Giá trị điện áp nút sau lắp đặt STATCOM nút 13 28 Hình 4.9: Mơ hình mơ phần tử khối STATCOM 28 Hình 4.10: Khối điều khiển STATCOM 29 Hình 4.11: Bên khối điều khiển STATCOM 30 Hình 4.12: Bộ đo giá trị dịng điện điện áp 31 Hình 4.13: Khâu phản hồi 31 Hình 4.14: Khâu tạo xung đồng 32 Hình 4.15: Khâu khuếch đại xung 32 Hình 4.16: Thơng số mơ trường hợp tải động 34 Hình 4.17: Thông số mô điện áp bus gần 35 x Chương TỔNG QUAN 1.1 Giới thiệu Hiện nay, có nhiều cơng trình nghiên cứu ứng dụng thuật tốn trí tuệ nhân để tính tốn tối ưu hóa cho hệ thống điện, đặc biệt sử dụng thiết bị phát nguồn công suất phản kháng cho hệ thống lưới điện nhằm đảm bảo ổn định điện áp cho hệ thống Tuy nhiên, việc đánh giá lựa chọn thiết bị phát công suất hợp lý, dung lượng bù tối ưu phân tích chế độ xác lập, độ chưa quan tâm sâu sắc Theo thực tế nay, hệ thống điện sử dụng hệ thống điện xoay chiều Hệ thống điện xoay chiều hệ thống điện phức tạp, gồm có máy phát đồng bộ, đường dây truyền tải, máy biến áp, thiết bị bù phụ tải…và chia thành ba khâu: sản xuất, truyền tải phân phối Một hệ thống điện xoay chiều hoạt động phải thỏa yêu cầu sau: - Các máy phát điện làm việc chế độ đồng - Điện áp vận hành nằm giới hạn cho phép theo qui định - Tần số vận hành nằm giới hạn cho phép theo qui định - Các phụ tải phải cung cấp nguồn điện đầy đủ - Các đường dây phải vận hành điều kiện bình thường khơng bị tải Trong hệ thống điện, công suất truyền tải đường dây phụ thuộc vào tổng trở đường dây, điện áp góc truyền tải điểm đầu điểm cuối đường dây, đại lượng giới hạn cơng suất truyền tải đường dây Vì vậy, khả truyền tải công suất đường dây cải thiện đáng kể việc tăng công suất phản kháng phía phụ tải, lắp cuộn kháng bù ngang (mắc song song), lắp tụ điện bù dọc (mắc nối tiếp) vào đường dây để điều khiển điện áp dọc theo chiều dài đường dây Hình 4.12: Bộ đo giá trị dòng điện điện áp Bộ đo giá trị dòng điện điện áp nằm bên khối điều khiển có tham số đầu vào điện áp Vabc, góc pha sin_cos dịng điện Iabc; tham số đầu dòng điện IdIq điện áp MagV Hình 4.13: Khâu phản hồi 31 Tham số đầu vào khâu phản hồi IdIq dòng điện tham chiếu IdIq_Ref đưa vào khối Sum qua khối hiệu chỉnh PID tham số ngõ điện áp VdVq Hình 4.14: Khâu tạo xung đồng Hình 4.15: Khâu khuếch đại xung Trường hợp 1: Giả sử cố ngắn mạch pha xảy phía bus đường dây nối bus bus thời gian chu kỳ bắt đầu thời điểm 0.2 giây với STATCOM gắn bus 13 Các kết mô trường hợp chưa lắp STATCOM (đường màu xanh) có STATCOM (đường màu đỏ) giá trị điện áp dòng điện số bus quan trọng trình bày hình 4.16 32 5.5 0.8 4.5 I bus (p.u.) Vbus (p.u.) 1.2 0.6 0.4 0.2 0 3.5 0.2 0.4 0.6 0.8 2.5 1.0 t (s) 0.2 0.4 0.6 0.8 1.0 0.8 1.0 t (s) a Điện áp bus b Dòng điện bus 1.2 0.25 0.2 0.8 I bus (p.u.) Vbus (p.u.) 0.6 0.4 0.1 0.2 0 0.15 0.2 0.4 0.6 0.8 1.0 0.05 t (s) 0.6 d Dòng điện bus 1.2 0.035 0.03 0.8 0.025 I bus 12 (p.u.) Vbus 12 (p.u.) 0.4 t (s) c Điện áp bus 0.6 0.4 0.2 0 0.2 0.02 0.015 0.01 0.2 0.4 0.6 0.8 0.005 1.0 t (s) 0.2 0.4 0.6 0.8 t (s) e Điện áp bus 12 f Dòng điện bus 12 33 1.0 1.2 0.25 0.2 0.8 I bus 13 (p.u.) Vbus 13 (p.u.) 0.6 0.4 0.1 0.05 0.2 0 0.15 0.2 0.4 0.6 0.8 1.0 t (s) 0 0.6 0.8 1.0 h Dòng điện bus 13 1.2 0.025 0.02 0.8 I bus 14 (p.u.) Vbus 14 (p.u.) 0.4 t (s) g Điện áp bus 13 0.6 0.4 0.015 0.01 0.005 0.2 0 0.2 0.2 0.4 0.6 0.8 0 1.0 0.2 t (s) 0.4 0.6 0.8 1.0 t (s) i Điện áp bus 14 j Dòng điện bus 14 Hình 4.16: Thơng số mơ trường hợp tải động Từ kết hình 4.16, ta kết thấy độ dao động độ hệ thống tăng cao đặc biệt điệ áp bus gần tải bus 12, bus 13 bus 14 Tuy nhiên, gắn STATCOM vào dao động bị loại trừ nhằm đảm bảo chất lượng điện Đây ưu điểm STATCOM Trường hợp 2: Giả sử trạng thái làm việc trường hợp cố ngắn mạch xảy bus 13 Các kết mô trình bày hình 4.17 Có thể thấy rõ từ kết mơ có cố ngắn mạch bus 13 nên điện áp giảm xuống p.u Đây trường hợp cố nghiêm trọng công suất STATCOM nhỏ nên 34 nâng mức điện áp lúc ngắn mạch lên Tuy nhiên, hệ thống khôi phục tốt dao 1.2 1 0.8 0.8 Vbus 13 (p.u.) 1.2 0.6 0.4 0.2 0 0.6 0.4 0.2 0.2 0.4 0.6 0.8 1.0 0 0.2 0.4 t (s) 0.6 0.8 t (s) a Điện áp bus 12 b Điện áp bus 13 1.2 Vbus 14 (p.u.) Vbus 12 (p.u.) động có giảm khơng đáng kể 0.8 0.6 0.4 0.2 0 0.2 0.4 0.6 0.8 1.0 t (s) c Điện áp bus 14 Hình 4.17: Thơng số mô điện áp bus gần 35 1.0 Chương KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN 5.1 Kế t luâ ̣n Luận văn đã nghiên cứu và giải quyế t đươ ̣c những vấ n đề cu ̣ thể sau: - Nghiên cứu yếu tố ảnh hưởng đến ổn định điện áp hệ thống điện - Nghiên cứu giải thuật di truyền GA ứng dụng giải thuật di truyền GA hệ thống điện để tìm vị trí dung lượng cơng suất phản kháng cần bù tối ưu Kết mô ứng dụng giải thuật di truyền GA cho hệ thống điện chuẩn 14 nút IEEE cho kết dung lượng bù tối ưu 0,706 (p.u) vị trí nút 13 - Nghiên cứu cấu tạo, nguyên lý làm việc đặt tính làm việc thiết bị bù STATCOM - Xây dựng mơ hình mơ hệ thống điện chuẩn 14 nút IEEE có ứng dụng STATCOM vị trí nút 13 với dung lượng bù tối ưu 0,706 (p.u) nâng cao chất lượng điện áp hệ thống lúc bình thường lúc xảy cố ngắn mạch ba pha cân bằng, có kết mơ điện áp nút hệ thống tốt thỏa điều kiện sai lệch điện áp cho phép ±5% so với trường hợp hệ thống chưa bù công suất phản kháng hay bù tụ điện GA đề xuất hệ thống điện chuẩn tìm giải pháp tối ưu vị trí giá trị công suất phản kháng cần bù cho hệ thống 5.2 Hướng phát triển Luận văn phát triển theo hướng sau: - Nghiên cứu áp dụng cho hệ thống thực Việt Nam có xét đến tính kinh tế - Nghiên cứu thuật tốn điều khiển tối ưu cho STATCOM việc nâng cao chất lượng điện cho hệ thống 36 TÀI LIỆU THAM KHẢO [1] Yong Hua Song and Allan T Johns - Flexible AC Transmission Systems (FACTS) [2] Hồ Văn Hiến “Hệ Thống Điện Truyền Tải Và Phân Phối” NXB ĐH Quốc Gia Tp HCM [3] TS Quyền Huy Ánh “Giáo trình giải tích mạng ” (Trường Đại học Sư phạm Kỹ thuật Tp.HCM) [4] Nguyễn Văn Nhờ “Giáo trình điện tử cơng suất 1” Nhà xuất Đại học Quốc gia TP Hồ Chí Minh 2002 [5] Holland JH Adaptation in natural and artificial systems Ann Arbor, MI: The University of Michigan Press; 1975 [6] Goldberg E Genetic algorithms in search, optimization & machine learning Reading, MA: Addison-Wesley; 1989 [7] Wenjuan Zhang, Fangxing Li, Leon M Tolbert, “Optimal allocation of Shunt Dynamic Var Source SVC and STATCOM: A survey”, Student Member, IEEE, Senior Member, IEEE, Senior Member, IEEE 2003 [8] Goldberg D Genetic algorithms in search, optimization & machine learning Reading, MA: Addison-Wesley; 1998 [9] Kalyanmoy D Multi-objective genetic algorithms: problem difficulties and construction of test problems Evol Computat 1999; 7(3):205–230 [10] Das DB, Patvardhan C New multi-objective stochastic search technique for economic load dispatch IEE Proc Generation Transmission Distribution 1998; 145(6):747–752 [11] How FACTS controllers benefits AC transmission systems: John J Paserba, Fellow IEEE [12] How FACTS improve the performance of electrical grid: Rolf Grunbaum, Ake Petersson, Bjom Thorvaldsson (ABB Review3/2002) 37 PHỤ LỤC Chương trình matlab sử dụng giải thuật GA cho hệ thống điện IEEE 14 bus clc; %% clean the screen command clear; close all; %% delete all current data tic; %% Add open Fourteen_bus_Copy.mdl; %% Mở mô hình P10 = 0.35; P11 = 0.35; P12 = 0.35; P13 = 0.35; P14 = 0.35; Q10 = 0; Q11 = 0; Q12 = 0; Q13 = 0; Q14 = 0; Q10temp = Q10; Q11temp = Q11; Q12temp = Q12; Q13temp = Q13; Q14temp = Q14; %% Add %% Algorithm Parameters SelMethod = 1; CrossMethod = 1; PopSize = 10; %% Number of Individuals MaxIteration = 25; %% Maximun Number of Iterations CrossPercent = 60; MutatPercent = 30; %% Percent of Mutation ElitPercent = 100 - CrossPercent - MutatPercent; CrossNum = round((CrossPercent/100)*PopSize); %% làm tròn số % if mod(CrossNum,2)~=0; Percent if mod(CrossNum,2)~=0; %% Chia lấy số dư 38 CrossNum = CrossNum - 1; end MutatNum = round((MutatPercent/100)*PopSize); % ElitNum = PopSize - CrossNum - MutatNum; ElitNum = PopSize - CrossNum - MutatNum; %% Problem Satement % VarMin = [1 2]; % VarMax = [3 5 3]; VarMin = [0 0 0]; %% Giá trị Min Q10-Q14 VarMax = [1 1 1]; %% Giá trị Max Q10-Q14 b=repmat(VarMin,PopSize,1); s=repmat(VarMax,PopSize,1); % DimNum = 6; DimNum = 5; % CostFuncName = @rastriginsfcn; %% Initial Population Pop = rand(PopSize,DimNum).* (s - b) +b; % Cost = feval(CostFuncName,Pop); Voltage_Value = Voltage_Fcn(Pop); %% Hàm tính tốn tổng giá trị điện áp tìm so với giá trị chuẩn (=1) % [Cost Indx] = sort(Cost); [Voltage_Value Indx] = sort(Voltage_Value); Pop = Pop(Indx,:); %% Main Loop MeanMat = []; MinMat = []; for Iter = 1:MaxIteration %% Elitism ElitPop = Pop(1:ElitNum,:); %% Cross Over CrossPop = []; % ParentIndexes = SelectParents_Fcn(Cost,CrossNum,SelMethod); ParentIndexes = SelectParents_Fcn(Voltage_Value,CrossNum,SelMethod); %% Chọn cá thể cha mẹ lai ghép lần for ii = 1:CrossNum/2 Par1Indx = ParentIndexes(ii*2-1); Par2Indx = ParentIndexes(ii*2); 39 Par1 = Pop(Par1Indx,:); Par2 = Pop(Par2Indx,:); [Off1 Off2] = MyCrossOver_Fcn(Par1,Par2,CrossMethod); CrossPop = [CrossPop ; Off1 ; Off2]; end %% Mutation % MutatPop = rand(MutatNum,DimNum).*(s(1:100,:) - b(1:100,:)) + b(1:100,:); MutatPop = rand(MutatNum,DimNum).*(s(1:MutatNum,:) - b(1:MutatNum,:)) + b(1:MutatNum,:); %% New Population Pop = [ElitPop ; CrossPop ; MutatPop]; % Cost = feval(CostFuncName,Pop); % [Cost Indx] = sort(Cost); Voltage_Value = Voltage_Fcn(Pop); [Voltage_Value Indx] = sort(Voltage_Value); Pop = Pop(Indx,:); %% Algorithm Progress disp(' ') BestP = Pop(1,:) % BestC = Cost(1) BestC = Voltage_Value(1); % MinMat(Iter) = Cost(1); MinMat(Iter) = Voltage_Value(1); % MeanMat(Iter) = mean(Cost); MeanMat(Iter) = mean(Voltage_Value); plot(MinMat,' r','linewidth',2); hold on % plot(MeanMat,' k','linewidth',2); % hold off pause(.5) end % ylim([0 5]) %% Results BestSolution = Pop(1,:) %% Add for i11 = 0:0 %% Repeat one time Q10 = Pop(1,1); %% Lấy vị trí hàng cột Pop(x,y): x hàng, y cột 40 Q11 = Q11temp; Q12 = Q12temp; Q13 = Q13temp; Q14 = Q14temp; sim( 'Fourteen_bus_Copy.mdl'); %% Mô B10temp = B10.signals(1,1).values(20001,1); %% Lấy giá trị điện áp Bus B11temp = B11.signals(1,1).values(20001,1); B12temp = B12.signals(1,1).values(20001,1); B13temp = B13.signals(1,1).values(20001,1); B14temp = B14.signals(1,1).values(20001,1); %% Tính tổng giá trị lệch so với chuẩn (=1) S11_temp1 = abs(1 - B10temp) + abs(1 - B11temp) + abs(1 - B12temp) + abs(1 B13temp) + abs(1 - B14temp); temp1 = 2; %% Cho gia tri tam =2 if B10temp >=0.95 && B11temp >=0.95 && B12temp >=0.95 && B13temp >=0.95 && B14temp >=0.95 if B10temp =0.95 if B10temp =0.95 if B10temp =0.95 if B10temp =0.95 if B10temp