Nghiên cứu phương pháp mô phỏng và đánh giá ổn định động hệ thống điện

Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện MỤC LỤC MỤC LỤC LỜI CẢM ƠN LỜI CAM ĐOAN TỪ VIẾT TẮT DANH MỤC CÁC BẢNG BIỂU .7 DANH MỤC CÁC HÌNH VẼ MỞ ĐẦU 11 M.1 Lý chọn đề tài 11 M.2 Lịch sử nghiên cứu 11 M.3 Mục đích, đối tượng, phạm vi nghiên cứu .12 M.4 Ý nghĩa khoa học thực tiễn đề tài .13 M.5 Phương pháp nghiên cứu .13 M.6 Nội dung luận văn 13 Chƣơng – TỔNG QUAN VỀ CÁC PHƢƠNG PHÁP NGHIÊN CỨU ỔN ĐỊNH ĐỘNG CỦA HỆ THỐNG ĐIỆN 14 1.1 Khái niệm chung 14 1.2 Mơ hình phần tử hệ thống điện dùng đánh giá ổn định động 16 1.2.1 Cấu trúc chung mơ hình hệ thống điện 16 1.2.2 Mơ hình máy phát, hệ thống kích từ tự động điều chỉnh kích từ, hệ thống điều chỉnh tốc độ quay tuabin .17 1.2.3 Mơ hình lưới điện, phụ tải thiết bị điều khiển hệ thống truyền tải24 1.2.4 Mơ hình tổng thể hệ thống điện .27 1.3 Các phương pháp tính tốn ổn định động 27 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện 1.3.1 Phương pháp tích phân số 28 1.3.2 Phương pháp cân diện tích 29 1.3.3 Phương pháp trực tiếp 34 1.4 Nghiên cứu ổn định động hệ thống nhiều máy phát .36 Chƣơng – NGHIÊN CỨU PHƢƠNG PHÁP TÌM THỜI GIAN CẮT TỚI HẠN VÀ ỨNG DỤNG ĐÁNH GIÁ ỔN ĐỊNH ĐỘNG CHO HTĐ MẪU 37 2.1 Một số phương pháp tính thời gian cắt tới hạn CCT .37 2.2 Đánh giá ổn định động dựa kết phân tích 40 2.3 Ứng dụng chương trình tính CCT đánh giá ổn định động cho HTĐ mẫu .43 2.3.1 Giới thiệu hệ thống điện IEEE nút 39 nút 43 2.3.2 Các bước thực 44 2.3.3 Các chương trình sử dụng luận văn 45 2.3.4 Kết tính tốn mơ chi tiết 45 2.3.5 Đánh giá ảnh hưởng thông số chế độ xác lập đến thời gian cắt tới hạn 57 2.4 Nhận xét 59 Chƣơng - ỨNG DỤNG MẠNG NƠRON NHÂN TẠO ĐÁNH GIÁ ỔN ĐỊNH QUÁ ĐỘ HỆ THỐNG ĐIỆN 60 3.1 Tổng quan mạng nơron nhân tạo mạng lan truyền thẳng MLP 60 3.1.1 Cấu trúc mạng MLP 61 3.1.2 Quá trình học mạng MLP 63 3.2 Ứng dụng mạng nơron đánh giá ổn định động HTĐ .66 3.3 Kết tính tốn với lưới máy nút 68 3.3.1 Đầu vào Pg 69 3.3.2 Đầu vào Pg-Qg 70 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện 3.3.3 Đầu vào U θ 71 3.4 Kết tính tốn với lưới 10 máy 39 nút 73 3.4.1 Đầu vào Pg .73 3.4.2 Đầu vào Pg Qg 76 3.4.3 Đầu vào U θ 79 KẾT LUẬN 83 PHỤ LỤC 84 PHỤ LỤC 98 TÀI LIỆU THAM KHẢO 103 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện LỜI CẢM ƠN  Đầu tiên cho gửi lời cám ơn đến tồn thể thầy giáo môn Hệ thống điện – Đại học ách khoa Hà Nội đ tạo điều kiện cho tơi hồn thành luận văn thạc s này, hội tốt thực hành k đ học giảng đường c ng giúp ngày tự tin vào thân  Tơi xin gửi lời cám ơn chân thành tới TS Ngu n Đức Hu suốt thời gian qua đ nhiệt tình dạy, giúp đ tơi hồn thành tốt luận văn thạc s  Tôi c ng xin gửi lời cám ơn sâu sắc tới toàn thể bạn b , người thân, nh ng người đ bên cạnh tôi, ủng hộ suốt thời gian qua  Do thời gian có hạn, chắn luận văn khơng tránh khỏi nh ng thiếu sót Em kính mong thầy bảo, đóng góp ý kiến để em hồn thiện, tiếp tục nghiên cứu phát triển đề tài Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện LỜI CAM ĐOAN Tôi xin cam đoan nh ng vấn đề trình bày luận văn nh ng nghiên cứu riêng cá nhân Các số liệu thống kê, báo cáo, tài liệu khoa học luận văn sử dụng công trình khác đ nghiên cứu thích đủ, quy định Hà Nội, ngày tháng năm 2015 Tác giả luận văn Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện TỪ VIẾT TẮT  ANN: Mạng nơron nhân tạo (Artifical Neural Networks)  AVR: Tự động điều chỉnh điện áp (Automatic Voltage Regulator)  CCT: Thời gian cắt tới hạn (Critical Clearing Time)  COA: Tâm góc (Center Of Angle)  COI: Tâm quán tính Center Of Inertia)  E C: Phương pháp cân diện tích Equal Area Criterion)  FACTS: Hệ thống truyền tải điện linh hoạt (Flexible Alternating Current Transmission System)  HVDC: Truyền tải điện cao áp chiều (High Voltage Direct Current)  HTĐ: Hệ thống điện  MAPE: Sai số trung bình phần trăm tuyệt đối (Mean Absolute Percent Error)  PSS: Bộ ổn định hệ thống nguồn (Power System Stabilizer)  SIME: Mơ hình đẳng trị máy phát (Single Machine Equivalent)  SVC: Tụ bù tĩnh Static VAR Compensator)  TĐK: Tự động điều chỉnh kích từ  UPFC: Hệ thống điều chỉnh dịng cơng suất (Unified Power Flow Controller) Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện DANH MỤC CÁC BẢNG BIỂU ảng 2.1 – Thông số d liệu nút hệ thống điện IEEE nút ảng 2.2 – Thông số d liệu nút hệ thống điện IEEE 39 nút ảng 2.3 – CCT máy phát cố với lưới IEEE nút ảng 2.4 – CCT máy phát cố với lưới IEEE 39 nút ảng 2.5 – Khảo sát CCT máy phát thay đổi thông số chế độ xác lập ảng 3.1 – Thông số chi tiết số liệu huấn luyện cho mạng nơ ron ảng 3.2 – Sai số huấn luyện máy cố, lưới máy nút ảng 3.3 – Sai số huấn luyện máy cố, lưới 10 máy 39 nút Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện DANH MỤC CÁC HÌNH VẼ Hình 1.1 – Đánh giá ổn định động hệ thống điện Hình 1.2 – Dao động góc lệch tuyệt đối a), tương đối (b) Hình 1.3 – Chuyển động rotor máy phát Hình 1.4 – Mơ hình máy phát đẳng trị hệ tọa độ vng góc Hình 1.5 – Mơ hình hệ thống kích từ TĐK Hình 1.6 – Mơ hình máy phát nối với lưới vơ lớn Hình 1.7 – Đồ thị cơng suất - góc Hình 1.8 – Đồ thị cơng suất - góc Hình 1.9 – Ngun lý ổn định bóng trượt khơng ma sát Hình 2.1 – Sơ đồ mơ tả qui trình tìm CCT Hình 2.2 – Sơ đồ khối tính CCT Hình 2.3 – Giải thuật phát ổn định động Hình 2.4 – Quan hệ gi a góc rotor máy phát thời gian máy cố Hình 2.5 – Quan hệ gi a tốc độ rotor máy phát thời gian máy cố Hình 2.6 – Quan hệ gi a điện áp máy phát thời gian máy cố Hình 2.7 – Quan hệ gi a góc rotor máy phát thời gian máy cố Hình 2.8 – Quan hệ gi a tốc độ rotor máy phát thời gian máy cố Hình 2.9 – Quan hệ gi a điện áp máy phát thời gian máy cố Hình 2.10 – Quan hệ gi a góc rotor máy phát thời gian máy cố Hình 2.11 – Quan hệ gi a tốc độ rotor máy phát thời gian máy 1sự cố Hình 2.12 – Quan hệ gi a điện áp máy phát thời gian máy cố Hình 2.13 – Quan hệ gi a góc rotor máy phát thời gian máy cố Hình 2.14 – Quan hệ gi a tốc độ rotor máy phát thời gian máy cố Hình 2.15 – Quan hệ gi a điện áp máy phát thời gian máy cố Hình 2.16 – Quan hệ gi a góc rotor máy phát thời gian máy cố Hình 2.17 – Quan hệ gi a tốc độ rotor máy phát thời gian máy cố Hình 2.18 – Quan hệ gi a điện áp máy phát thời gian máy cố Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện Hình 2.19 – Quan hệ gi a góc rotor máy phát thời gian máy cố Hình 2.20 – Quan hệ gi a tốc độ rotor máy phát thời gian máy cố Hình 2.21 – Quan hệ gi a điện áp máy phát thời gian máy cố Hình 2.22 – Quan hệ gi a góc rotor máy phát thời gian máy cố Hình 2.23 – Quan hệ gi a tốc độ rotor máy phát thời gian máy cố Hình 2.24 – Quan hệ gi a điện áp máy phát thời gian máy cố Hình 3.1 – Mơ hình mạng MLP lớp ẩn Hình 3.2 – CCT xấp xỉ mơ phỏng, máy cố đầu vào Pg, lưới IEEE nút Hình 3.3 – Phân bố sai số xấp xỉ %), máy cố đầu vào Pg, lưới IEEE nút Hình 3.4 – CCT xấp xỉ mơ phỏng, máy cố đầu vào Pg-Qg, lưới IEEE nút Hình 3.5 – Phân bố sai số xấp xỉ %), máy cố đầu vào Pg-Qg, lưới IEEE nút Hình 3.6 – CCT xấp xỉ mơ phỏng, máy cố đầu vào U-θ, lưới IEEE nút Hình 3.7 – Phân bố sai số xấp xỉ %), máy cố đầu vào U-θ, lưới IEEE nút Hình 3.8 – CCT xấp xỉ mô phỏng, máy cố, đầu vào Pg lưới IEEE 39 nút Hình 3.9 – Phân bố sai số xấp xỉ %), máy cố, đầu vào Pg lưới IEEE 39 nút Hình 3.10 – CCT xấp xỉ mô phỏng, máy cố, đầu vào Pg lưới IEEE 39 nút Hình 3.11 – Phân bố sai số xấp xỉ %), máy cố, đầu vào Pg lưới IEEE 39 nút Hình 3.12 – CCT xấp xỉ mô phỏng, máy cố, đầu vào Pg lưới IEEE 39 nút Hình 3.13 – Phân bố sai số xấp xỉ %), máy cố, đầu vào Pg lưới IEEE 39 nút Hình 3.14 – CCT xấp xỉ mô phỏng, máy cố, đầu vào Pg-Qg lưới IEEE 39 nút Hình 3.15 – Phân bố sai số xấp xỉ %), máy cố, đầu vào Pg-Qg lưới IEEE 39 nút Hình 3.16 – CCT xấp xỉ mô phỏng, máy cố, đầu vào Pg-Qg lưới IEEE 39 nút Hình 3.17 – Phân bố sai số xấp xỉ %), máy cố, đầu vào Pg-Qg lưới IEEE 39 nút Hình 3.18 – CCT xấp xỉ mô phỏng, máy cố, đầu vào Pg-Qg lưới IEEE 39 nút Hình 3.19 – Phân bố sai số xấp xỉ %), máy cố, đầu vào Pg-Qg lưới IEEE 39 nút Hình 3.20 – CCT xấp xỉ mô phỏng, máy cố, đầu vào U-θ lưới IEEE 39 nút Hình 3.21 – Phân bố sai số xấp xỉ %), máy cố, đầu vào U-θ lưới IEEE 39 nút Hình 3.22 – CCT xấp xỉ mơ phỏng, máy cố, đầu vào U-θ lưới IEEE 39 nút Hình 3.23 – Phân bố sai số xấp xỉ %), máy cố, đầu vào U-θ lưới IEEE 39 nút Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện Hình 3.24 – CCT xấp xỉ mô phỏng, máy cố, đầu vào U-θ lưới IEEE 39 nút Hình 3.25 – Phân bố sai số xấp xỉ %), máy cố, đầu vào U-θ lưới IEEE 39 nút 10 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện if output; fprintf('> Allocate memory '); end chunk = 500; Time = zeros(chunk,1); Time(1,:) = t; Errest = zeros(chunk,1); Errest(1,:) = errest; Stepsize = zeros(chunk,1); Stepsize(1,:) = stepsize; % System variables Voltages = zeros(chunk, length(U0)); Voltages(1,:) = U0.'; % Generator Angles = zeros(chunk,ngen); Angles(1,:) = Xgen0(:,1).*180./pi; Speeds = zeros(chunk,ngen); Speeds(1,:) = Xgen0(:,2)./(2.*pi.*freq); Eq_tr = zeros(chunk,ngen); Eq_tr(1,:) = Xgen0(:,3); Ed_tr = zeros(chunk,ngen); Ed_tr(1,:) = Xgen0(:,4); Pe = zeros(chunk,ngen); Pe(1,:)=Pe0; % Exciter and governor Efd = zeros(chunk,ngen); Efd(1,:) = Efd0(:,1); PM = zeros(chunk,ngen); PM(1,:) = Pm0(:,1); %Stability Checker Sc=1; STOP = 0; %% Main stability loop while (t < stoptime + stepsize) if i>3 CO(i) = Angles(i,:)*Pgen0(:,4)/sum(Pgen0(:,4)); if i > 50 if abs(max(Angles(i-1,:))-CO(i-1)-max(Angles(i-46,:))+CO(i46)) > 100 STOP = 1; end end; end %% Output i=i+1; if mod(i,45)==0 && output fprintf('\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b> %6.2f%% completed', t/stoptime*100) end %% Numerical Method switch method case 90 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện [Xgen0, Pgen0, Vgen0, Xexc0, Pexc0, Vexc0, Xgov0, Pgov0, Vgov0, U0, t, newstepsize] = ModifiedEuler(t, Xgen0, Pgen0, Vgen0, Xexc0, Pexc0, Vexc0, Xgov0, Pgov0, Vgov0, Ly, Uy, Py, gbus, genmodel, excmodel, govmodel, stepsize); case [Xgen0, Pgen0, Vgen0, Xexc0, Pexc0, Vexc0, Xgov0, Pgov0, Vgov0, U0, t, newstepsize] = RungeKutta(t, Xgen0, Pgen0, Vgen0, Xexc0, Pexc0, Vexc0, Xgov0, Pgov0, Vgov0, Ly, Uy, Py, gbus, genmodel, excmodel, govmodel, stepsize); case [Xgen0, Pgen0, Vgen0, Xexc0, Pexc0, Vexc0, Xgov0, Pgov0, Vgov0, U0, errest, failed, t, newstepsize] = RungeKuttaFehlberg(t, Xgen0, Pgen0, Vgen0, Xexc0, Pexc0, Vexc0, Xgov0, Pgov0, Vgov0, U0, Ly, Uy, Py, gbus, genmodel, excmodel, govmodel, tol, maxstepsize, stepsize); case [Xgen0, Pgen0, Vgen0, Xexc0, Pexc0, Vexc0, Xgov0, Pgov0, Vgov0, U0, errest, failed, t, newstepsize] = RungeKuttaHighamHall(t, Xgen0, Pgen0, Vgen0, Xexc0, Pexc0, Vexc0, Xgov0, Pgov0, Vgov0, U0, Ly, Uy, Py, gbus, genmodel, excmodel, govmodel, tol, maxstepsize, stepsize); case [Xgen0, Pgen0, Vgen0, Xexc0, Pexc0, Vexc0, Xgov0, Pgov0, Vgov0, U0, t, eulerfailed, newstepsize] = ModifiedEuler2(t, Xgen0, Pgen0, Vgen0, Xexc0, Pexc0, Vexc0, Xgov0, Pgov0, Vgov0, Ly, Uy, Py, gbus, genmodel, excmodel, govmodel, stepsize); end if eulerfailed fprintf('\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b> Error: No solution found Try lowering tolerance or increasing maximum number of iterations in ModifiedEuler2 Exiting \n') return; end if failed t = t-stepsize; end % End exactly at stop time if t + newstepsize > stoptime newstepsize = stoptime - t; elseif stepsize < minstepsize 91 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện fprintf('\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b> Error: No solution found with minimum step size Exiting \n') return; end; %% Allocate new memory chunk if matrices are full if i>size(Time,1) Stepsize = [Stepsize; zeros(chunk,1)];Errest = [Errest; zeros(chunk,1)];Time = [Time; zeros(chunk,1)]; Voltages = [Voltages; zeros(chunk,length(U0))];Efd = [Efd; zeros(chunk,ngen)];PM = [PM; zeros(chunk,ngen)]; Angles=[Angles;zeros(chunk,ngen)];Speeds=[Speeds;zeros(chunk,ngen)];Eq_tr =[Eq_tr;zeros(chunk,ngen)];Ed_tr=[Ed_tr;zeros(chunk,ngen)]; end %% Save values Stepsize(i,:) = stepsize.'; Errest(i,:) = errest.'; Time(i,:) = t; Voltages(i,:) = U0.'; % exc Efd(i,:) = Xexc0(:,1).*(genmodel>1); % Set Efd to zero when using classical generator model % gov PM(i,:) = Xgov0(:,1); % gen Angles(i,:) = Xgen0(:,1).*180./pi; Speeds(i,:) = Xgen0(:,2)./(2.*pi.*freq); Eq_tr(i,:) = Xgen0(:,3); Ed_tr(i,:) = Xgen0(:,4); Pe(i,:) = Vgen0(:,3).'; %% Adapt step size if event will occur in next step if ~isempty(event) && ev = event(ev,1) if event(ev,1) - t < newstepsize newstepsize = event(ev,1) - t; end end end 92 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện %% Check for events if ~isempty(event) && ev 10*eps || ev > size(event,1) % that happen on time t break; else eventhappened = true; end switch event(ev,2) case bus(buschange(ev,2),buschange(ev,3)) = buschange(ev,4); case branch(linechange(ev,2),linechange(ev,3)) = linechange(ev,4); end ev=ev+1; end if eventhappened % Refactorise [Ly, Uy, Py] = AugYbus(baseMVA, bus, branch, xd_tr, gbus, bus(:,PD)./baseMVA, bus(:,QD)./baseMVA, U00); U0 = SolveNetwork(Xgen0, Pgen0, Ly, Uy, Py, gbus, genmodel); [Id0,Iq0,Pe0] = MachineCurrents(Xgen0, Pgen0, U0(gbus), genmodel); Vgen0 = [Id0,Iq0,Pe0]; Vexc0 = abs(U0(gbus)); % decrease stepsize after event occured if method==3 || method==4 newstepsize = minstepsize; end i=i+1; % if event occurs, save values at t- and t+ %% Save values Stepsize(i,:) = stepsize.'; Errest(i,:) = errest.'; Time(i,:) = t; Voltages(i,:) = U0.'; % exc 93 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện Efd(i,:) = Xexc0(:,1).*(genmodel>1); % Set Efd to zero when using classical generator model % gov PM(i,:) = Xgov0(:,1); % gen Angles(i,:) = Xgen0(:,1).*180./pi; Speeds(i,:) = Xgen0(:,2)./(2.*pi.*freq); Eq_tr(i,:) = Xgen0(:,3); Ed_tr(i,:) = Xgen0(:,4); Pe(i,:)= Vgen0(:,3).'; eventhappened = false; end end %% Advance time stepsize = newstepsize; t = t + stepsize; if STOP && (t > tc + 0.2) stoptime = t-stepsize * 2; STOP = 0; end; end % end of main stability loop %% Output if output fprintf('\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b> 100%% completed') else fprintf('\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b') end Angles = Angles(1:i,:); Speeds = Speeds(1:i,:); Eq_tr = Eq_tr(1:i,:); Ed_tr = Ed_tr(1:i,:); Efd = Efd(1:i,:); PM = PM(1:i,:); Voltages = Voltages(1:i,:); Pe=Pe(1:i,:); Stepsize = Stepsize(1:i,:); Errest = Errest(1:i,:); Time = Time(1:i,:); if plots 94 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện figure xlabel('Time [s]') ylabel('Angle [deg]') hold on CO(i) = Angles(i,:)*Pgen0(:,4)/sum(Pgen0(:,4)); plot(Time,Angles); hold on;plot(Time,CO','LineWidth',2,'Color','Red'); hold on;plot(Time,Angles(:,gen_fault),'LineWidth',2,'Color','Black'); axis([0 Time(end) -1 1]) axis 'auto y' figure xlabel('Time [s]') ylabel('Speed [pu]') hold on plot(Time,Speeds) axis([0 Time(end) -1 1]) axis 'auto y' figure xlabel('Time [s]') ylabel('Voltage [pu]') hold on plot(Time,abs(Voltages)) axis([0 Time(end) -1 1]) axis 'auto y' end ns_ = numel(CO); if abs(max(Angles(ns_-1,:))-CO(ns_-1)-max(Angles(ns_-46,:))+CO(ns_-46)) > 60 Sc = -1; [~,spidx] = max(Speeds(end,:)); disp(['Max speed at = ' num2str(spidx) ' machine']) end return; 95 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện Chƣơng trình tính CCT function [CCT] = CCT_finder_v3(tc,delta0,tol,casefile_pf, casefile_dyn, casefile_ev,gen_fault) %% tc Clearing time % delta0 leap % tol Tolerence % tic t=tc; t0 = casefile_ev.event(1,1); delta=delta0; mpopt=mpoption; mpopt(31)=0; mpopt(32)=0; casefile_pf = runpf(casefile_pf,mpopt); Sc=rundyn_ss_v3(casefile_pf, casefile_dyn, casefile_ev,t,gen_fault); if Sc > lastknownSTABLE = t - t0; else lastknownUNSTABLE = t - t0; end; i=0; t1=[]; Sc1=[]; flag=1; while delta>tol i=i+1; t=t+delta*Sc; t1=[t1 t]; Sc0=Sc; Sc=rundyn_ss_v3(casefile_pf, casefile_dyn, casefile_ev,t,gen_fault); Sc1=[Sc1 Sc]; if Sc > lastknownSTABLE = t - t0; else lastknownUNSTABLE = t - t0; end; if Sc*Sc040 display('Qua nhieu vong lap'); return end end [t1;Sc1] CCT = (lastknownUNSTABLE + lastknownSTABLE)/2; End 97 Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện PHỤ LỤC Các chƣơng trình tính tốn với mạng nơron nhân tạo  Tính tốn cho lƣới IEEE nút clear all; close all %% load case9_data_may02 CCTdat.mat; Macidx PGIDX = 1:3; QGIDX = 4:6; UBIDX = 7:15; THETAIDX = 16:24; PQLOADIDX = 25:42; CCTIDX = 43:44; = 2; % input = [data(:,PGIDX)]; % input = [data(:,PGIDX) data(:,QGIDX)]; input = [data(:,UBIDX) data(:,THETAIDX)]; output = data(:,CCTIDX(Macidx)); %% Remove zero inputs x = input(1,:); pridx = find(x ~= 0); input = input(:,pridx); %% output1=output; min_err=inf; desired_err = 2.6; % Target error while min_err>desired_err for i = 1:30 net=newff(minmax(input'),output1',[40],{'tansig' % Tao mang MLP lop an net.divideParam.trainRatio=99/100; net.divideParam.valRatio =15/100; net.divideParam.testRatio=15/100; net.trainParam.epochs=200; net.trainParam.show=50; net.trainParam.goal=1e-5; net.inputs{1}.processFcns = {}; net.outputs{2}.processFcns = {}; net =trainlm(net,input',output1'); 98 'purelin'}, 'trainlm'); Nghiên cứu phương pháp mô đánh giá ổn định động hệ thống điện %Du bao Y=sim(net,input'); % Kiem tra sai so huan luyen err = output1-Y'; errpct = abs(err)./output1*100; errpct1 = (err)./output1*100; MAE = mean(abs(err)); MAPE = nanmean(errpct(~isinf(errpct))); fprintf('Sai so huan luyen MAPE: %0.3f%%\n',MAPE); if MAPE