(Đề tài NCKH) tối ưu hóa vị trí của thiết bị FACTS để nâng cao khả năng làm việc của hệ thống điện

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ PHAN THÁI TỒN TỐI ƯU HĨA VỊ TRÍ CỦA THIẾT BỊ FACTS ĐỂ NÂNG CAO KHẢ NĂNG LÀM VIỆC CỦA HỆ THỐNG ĐIỆN NGÀNH: KỸ THUẬT ĐIỆN - 60520202 SKC004337 Tp Hồ Chí Minh, tháng 10/2014 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ PHAN THÁI TỒN TỐI ƯU HĨA VỊ TRÍ CỦA THIẾT BỊ FACTS ĐỂ NÂNG CAO KHẢ NĂNG LÀM VIỆC CỦA HỆ THỐNG ĐIỆN NGÀNH: KỸ THUẬT ĐIỆN - 60520202 Tp Hồ Chí Minh, tháng 10/2014 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ PHAN THÁI TỒN TỐI ƯU HĨA VỊ TRÍ CỦA THIẾT BỊ FACTS ĐỂ NÂNG CAO KHẢ NĂNG LÀM VIỆC CỦA HỆ THỐNG ĐIỆN NGÀNH: KỸ THUẬT ĐIỆN - 60520202 Hướng dẫn khoa học: PGS.TS NGUYỄN HỒNG VIỆT Tp Hồ Chí Minh, tháng 10/2014 LÝ LỊCH KHOA HỌC I LÝ LỊCH SƠ LƯỢC: Họ & tên: Phan Thái Tồn Giới tính: Nam Ngày, tháng, năm sinh: 03/03/1988 Nơi sinh: Tiền Giang Quê quán: Tiền giang Dân tộc: Kinh Chỗ riêng địa liên lạc: Q7, TPHCM E-mail: phanthaitonntt@yahoo.com II QUÁ TRÌNH ĐÀO TẠO: Cao Đẳng: Hệ đào tạo: Chính Quy Thời gian đào tạo từ 10/2006 đến 10/2009 Nơi học (trường, thành phố): Q4, Thành Phố Hồ Chí Minh Ngành học: Công nghệ Kỹ Thuật Điện-Điện Tử Đại học: Hệ đào tạo: Chính Quy Thời gian đào tạo từ 10/2009 đến 10/2011 Nơi học (trường, thành phố): Quận Bình Thạnh, Thành Phố Hồ Chí Minh Ngành học: Cơng nghệ Kỹ Thuật Điện-Điện Tử III Q TRÌNH CƠNG TÁC CHUN MƠN KỂ TỪ KHI TỐT NGHIỆP ĐẠI HỌC: Thời gian 11/200903/2011 03/2011- Nay LỜI CAM ĐOAN Tôi cam đoan cơng trình nghiên cứu tơi 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 … năm 201… (Ký tên ghi rõ họ tên) LỜI CẢM TẠ Trước tiên em xin bày tỏ lòng biết ơn sâu sắc đến cố PGS.TS Nguyễn Hoàng Việt dù điều kiện sức khỏe có phần sa sút thầy cố gắng tận tình hướng dẫn em hồn thành luận văn tốt nghiệp Xin chân thành cảm ơn đến tất quý thầy cô Trường Đại Học Sư Phạm Kỹ Thuật trường Đại Học Bách Khoa TPHCM truyền đạt cho chúng em kiến thức vô quý báo làm hành trang cho em vững vàng công việc Xin cảm ơn tập thể lớp 2012b giúp đỡ, chia sẻ với khó khăn hoc tập, tạo cho em niềm tin nỗ lực phấn đấu hồn thành khóa hoc Chân Thành Cảm Ơn! TĨM TẮT Luận văn trình bày phương pháp tìm vị trí thiết bị FACTS tối ưu hệ thống điện sử dụng giải thuật di truyền GA (Genetic Alglorithms), để từ tiến hành viết giao diện nhỏ ngơn ngữ lập trình Matlab tính tốn tổn thất cơng suất phản kháng chi phí đầu tư hệ thống điện trước sau lắp thiết bị FACTS lưới điện IEEE 30 bus Thiết bị FACTS cụ thể luận văn SVC TCSC Kết tính tốn chương trình xác ứng dụng cho nhiều lưới điện khác Trong chương giới thiệu dạng ổn định hệ thống điện, nguyên tắt điều khiển hệ thống điện Chương giới thiệu tổng quan thiết bị FACTS điều khiển thiết bị FACTS Chương giới thiệu giải thuật di truyền phân bố trào lưu cơng suất truyền thống làm sở cho lập trình chương trình tính tốn Chương trình bày kết chương trình tính tốn thay đổi trọng số toán tối ưu Chương đánh giá kết tính tốn đề suất hướng phát triển để xây dựng nên cơng trình nghiên cứu tìm vị trí tối ưu thiết bị FACTS hệ thống điện hoàn chỉnh ABSTRACT This thesis presents methods to find the optimal location FACTS devices in power systems using genetic algorithm GA (Genetic Alglorithms), then proceed to write a small interface on the programming language Matlab computational costs reactive power losses and investment costs of the electricity system before and after installation of FACTS devices IEEE 30-bus grid FACTS devices specifically in this thesis is the SVC and TCSC Calculation results of this program is very accurate and can be applied to many different grids In Chapter introduces the types of power system stability, the material off of the power system control Chapters and present an overview of device controllers FACTS FACTS devices Chapter introduces genetic algorithms and distribution power flow as the basis for traditional programming calculation program Chapter presents the results of the calculation program to change the weights of the optimization problem Chapter reviews the results of calculations and proposal development continued to build up a study to find the optimal location of FACTS devices in power system complete MỤC LỤC TRANG Trang tựa i Quyết định giao đề tài ii Lý lịch cá nhân iii Lời cam đoan iv Cảm tạ v Tóm tắt vi Mục lục vii Danh sách chữ viết tắt viii Danh sách hình ix Danh sách bảng x PHẦN A: GIỚI THIỆU LUẬN VĂN Đặt vấn đề Mục tiêu nhiệm vụ luận văn Phạm vi nghiên cứu Phương pháp giải toán Điểm luận văn Giá trị thực tiễn Bố cục luận văn PHẦN B: NỘI DUNG CHƯƠNG 1: ĐIỀU KHIỂN VÀ CÁC GIỚI HẠN ỔN ĐỊNH CỦA HỆ THỐNG ĐIỆN 1.1 Đặc điểm lưới điện truyền tải 1.2 Cơ sở kiến thức 1.3 Điều khiển phân bố công suất 1.4 Giới hạn ổn định hệ thống 1.4.1Giới hạn nhiệt 1.4.2 Giới hạn điện áp 1.4.3 Giới hạn ổn định 1.4.3.1Ổn định độ 10 1.4.3.2 Ổn định dao động bé 10 1.4.4 Ổn định điện áp 11 1.5 Nguyên lý bù đường dây truyền tải 11 1.5.1 Bù song song 12 1.5.2 Bù nối tiếp 13 1.6 Tầm quan trọng việc giảm tổn hao công suất phản kháng 13 CHƯƠNG 2: TỔNG QUÁT VỀ THIẾT BỊ FACTS 2.1 Khái niệm thiết bị FACTS 15 2.1.1 Đặt vấn đề 15 2.1.2 Khái niệm 16 2.2 Các điều khiển thiết bị FACTS 17 2.2.1 Ưu điểm thiết bị FACTS 17 2.2.2 Phân loại thiết bị FACTS 17 2.2.2.1 Bộ bù ngang tĩnh SVC ( Static Var Compensator) 18 2.2.2.2 Bộ bù ngang tĩnh STATCOM (Static Synchronous Compensator) 22 2.2.2.3 Bộ bù dọc tĩnh TCSC( Thyristor Controlled Series Capacitor) 23 2.2.2.4 Bộ tích hợp điều khiển luồng công suất UPFC(Unified Power Flow Controlled) 27 CHƯƠNG 3: MƠ HÌNH CÁC BỘ ĐIỀU KHIỂN THIẾT BỊ FACTS 3.1 Giới Thiệu 30 3.2 Các điều khiển dựa Thyristor truyền thống 31 3.2.1Cuộn kháng điều khiển Thyristor (TCR) 31 3.2.2 Bộ bù tĩnh SVC ( Static Var Compensator) 36 function edit_SVC_MIN_Callback(hObject, % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_SVC_MIN as text % str2double(get(hObject,'String')) returns contents of edit_SVC_MIN as a double % - Executes during object creation, after setting all properties function edit_SVC_MIN_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_SVC_MAX_Callback(hObject, % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_SVC_MAX as text % str2double(get(hObject,'String')) returns contents of edit_SVC_MAX as a double % - Executes during object creation, after setting all properties function edit_SVC_MAX_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_NUMBER_OF_TCSC_Callback(hObject, eventdata, % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_NUMBER_OF_TCSC as text % str2double(get(hObject,'String')) returns contents of edit_NUMBER_OF_TCSC as a double % - Executes during object creation, after setting all properties function edit_NUMBER_OF_TCSC_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_TCSC_MIN_Callback(hObject, % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_TCSC_MIN as text % str2double(get(hObject,'String')) returns contents of edit_TCSC_MIN as a double % - Executes during object creation, after setting all properties function edit_TCSC_MIN_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_TCSC_MAX_Callback(hObject, % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_TCSC_MAX as text % str2double(get(hObject,'String')) returns contents of edit_TCSC_MAX as a double % - Executes during object creation, after setting all properties function edit_TCSC_MAX_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_ANPHA_1_Callback(hObject, % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_ANPHA_1 as text % str2double(get(hObject,'String')) returns contents of edit_ANPHA_1 as a double % - Executes during object creation, after setting all properties function edit_ANPHA_1_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_ANPHA_2_Callback(hObject, % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_ANPHA_2 as text % str2double(get(hObject,'String')) returns contents of edit_ANPHA_2 as a double % - Executes during object creation, after setting all properties function edit_ANPHA_2_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_LOCATION_OF_SVC_Callback(hObject, eventdata, handles) % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_LOCATION_OF_SVC as text % str2double(get(hObject,'String')) returns contents of edit_LOCATION_OF_SVC as a double % - Executes during object creation, after setting all properties function edit_LOCATION_OF_SVC_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_RATE_OF_SVC_Callback(hObject, eventdata, % hObject % eventdata % handles % text Hints: get(hObject,'String') returns contents of edit_RATE_OF_SVC as % str2double(get(hObject,'String')) returns contents of edit_RATE_OF_SVC as a double % - Executes during object creation, after setting all properties function edit_RATE_OF_SVC_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) end set(hObject,'BackgroundColor','white'); function edit_LOCATION_OF_TCSC_Callback(hObject, eventdata, handles) % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_LOCATION_OF_TCSC as text % str2double(get(hObject,'String')) returns contents of edit_LOCATION_OF_TCSC as a double % - Executes during object creation, after setting all properties function edit_LOCATION_OF_TCSC_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) end set(hObject,'BackgroundColor','white'); function edit_RATE_OF_TCSC_Callback(hObject, % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_RATE_OF_TCSC as text % str2double(get(hObject,'String')) returns contents of edit_RATE_OF_TCSC as a double % - Executes during object creation, after setting all properties function edit_RATE_OF_TCSC_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_X_OF_TCSC_Callback(hObject, eventdata, handles) % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_X_OF_TCSC as text % str2double(get(hObject,'String')) returns contents of edit_X_OF_TCSC as a double % - Executes during object creation, after setting all properties function edit_X_OF_TCSC_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) end set(hObject,'BackgroundColor','white'); function edit_QLOSS_Callback(hObject, % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_QLOSS as text % str2double(get(hObject,'String')) returns contents of edit_QLOSS as a double % - Executes during object creation, after setting all properties function edit_QLOSS_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_COST_Callback(hObject, eventdata, handles) % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_COST as text % as str2double(get(hObject,'String')) returns contents of edit_COST a double % - Executes during object creation, after setting all properties function edit_COST_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) end set(hObject,'BackgroundColor','white'); function edit_VMIN_Callback(hObject, eventdata, handles) % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_VMIN as text % as str2double(get(hObject,'String')) returns contents of edit_VMIN a double % - Executes during object creation, after setting all properties function edit_VMIN_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_QLOSS_INITIAL_Callback(hObject, eventdata, handles) % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_QLOSS_INITIAL as text % str2double(get(hObject,'String')) returns contents of edit_QLOSS_INITIAL as a double % - Executes during object creation, after setting all properties function edit_QLOSS_INITIAL_CreateFcn(hObject, eventdata, handles) % hObject % % eventdata handles % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_VMIN_INITIAL_Callback(hObject, % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_VMIN_INITIAL as text % str2double(get(hObject,'String')) returns contents of edit_VMIN_INITIAL as a double % - Executes during object creation, after setting all properties function edit_VMIN_INITIAL_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) end set(hObject,'BackgroundColor','white'); % - Executes on button press in pushbutton_START function pushbutton_START_Callback(hObject, eventdata, handles) % hObject % eventdata % handles tcscsvc(handles); % - Executes on button press in pushbutton_RESET function pushbutton_RESET_Callback(hObject, eventdata, handles) % hObject % eventdata % handles set(handles.edit_QLOSS,'string','MVAr'); set(handles.edit_COST,'string','$'); set(handles.edit_RATE_OF_TCSC,'string','MVAr'); set(handles.edit_LOCATION_OF_TCSC,'string','Branch'); set(handles.edit_X_OF_TCSC,'string','K'); set(handles.edit_LOCATION_OF_SVC,'string','Bus' ); set(handles.edit_RATE_OF_SVC,'string','MVar'); set(handles.edit_QLOSS_INITIAL,'string','Initial Qloss'); set(handles.edit_VMIN_INITIAL,'string','Initial Vmin'); set(handles.edit_VMIN,'string','New Vmin'); set(handles.edit_OPTIMAL_INDEX,'string','Index'); function edit_OPTIMAL_INDEX_Callback(hObject, eventdata, handles) % hObject % eventdata % handles % Hints: get(hObject,'String') returns contents of edit_OPTIMAL_INDEX as text % str2double(get(hObject,'String')) returns contents of edit_OPTIMAL_INDEX as a double % - Executes during object creation, after setting all properties function edit_OPTIMAL_INDEX_CreateFcn(hObject, eventdata, handles) % hObject % % % eventdata handles Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function [Qloss,rateTCSC]=main_tcscsvc(XY,busdata2,linedata2,old) busdata=busdata2; linedata=linedata2; lfybus lfnewton lineflow Ploss=deltaP;Qloss=deltaQ; new=thien; nhanh=[linedata2(branch_k,1),linedata2(branch_k,2)]; [m,n]=size(new); for i=1:m tachnhanh=[new(i,1),new(i,2)]; if nhanh==tachnhanh rateTCSC=abs(old(i,5))- abs(new(i,5)); break; end end function tcscsvc(handles) clear all; clc %doc du lieu ban dau busdata0; linedata0; %goi chuong trinh tinh toan popsize=get(handles.edit_POPSIZE,'string'); popsize=str2num(popsize); maxit=get(handles.edit_ITERATION,'string'); maxit=str2num(maxit); mutrate=get(handles.edit_MUTATION_RATE,'string'); mutrate=str2num(mutrate); selection=get(handles.edit_SELECTION,'string'); selection=str2num(selection); anpha1=get(handles.edit_ANPHA_1,'string'); anpha1=str2num(anpha1); anpha2=get(handles.edit_ANPHA_2,'string'); anpha2=str2num(anpha2); npar=4; varlo=[1 -0.8 -100]; %[min position TCSC, rate TCSC, position_SVC, varhi=[41 0.2 30 100]; busdata=busdata0; linedata=linedata0; lfnewton lfybus lineflow busout_plotV; vmin_Initial = vmin; old_Voltage = thuan(:,2); Old_Qloss=deltaQ; old=thien; [bestcost,bestsolution,globalbestcost,rate_TCSC,new_Qloss,cost] = GA(npar,varhi, varlo,popsize,maxit,busdata0,linedata0,old,mutrate,selection,anpha1, anpha2); figure(123); plot(globalbestcost,'b'); grid on title('\fontsize{12}\bf Optimal location of FACTS devices'); xlabel('\fontsize{12}\bf Iteration');ylabel('\fontsize{12}\bf Optimal Index'); legend('\fontsize{10}\bf Genetic Algorithim',1); set(handles.edit_OPTIMAL_INDEX,'string',num2str(bestcost)); set(handles.edit_QLOSS,'string',num2str(new_Qloss)); set(handles.edit_COST,'string',num2str(cost)); set(handles.edit_RATE_OF_TCSC,'string',num2str(-rate_TCSC)); set(handles.edit_LOCATION_OF_TCSC,'string',num2str(bestsolution(1,1))); set(handles.edit_X_OF_TCSC,'string',num2str(bestsolution(1,2))); set(handles.edit_LOCATION_OF_SVC,'string',num2str(bestsolution(1,3))); set(handles.edit_RATE_OF_SVC,'string',num2str(bestsolution(1,4))); set(handles.edit_QLOSS_INITIAL,'string',num2str(Old_Qloss)); set(handles.edit_VMIN_INITIAL,'string',num2str(vmin_Initial)); [V_plot,Vmin]=take_new_vmin(bestsolution,busdata0,linedata0);%chuong trinh lay vmin va ve V sau uu set(handles.edit_VMIN,'string',num2str(Vmin)); figure(456); plot(V_plot,'-r'); title('\fontsize{12}\bf Compare voltages after optimization'); xlabel('\fontsize{12}\bf Bus');ylabel( '\fontsize{12}\bf P.U.'); legend('\fontsize{10}\bf Optimal',1); hold on; plot(old_Voltage,'-b'); grid on; function [bestcost,bestsolution,globalbestcost,rate_TCSC,new_Qloss,money ]=GA(npar,varhi, varlo, popsize,maxit) % _ ff='testfunction'; % objective function Nt=npar; % keep=floor(selection*popsize); % # (lam tron so nguyen nho nhat) nmut = ceil((popsize-1)*Nt*mutrate); % (lam tron ve so nguyen lon hon) M=ceil((popsize-keep)/2); % so luong sinh san % _ iga=0; % Nvar = size(FieldD,2); aux = rand(popsize,npar); m=[-1 1]*FieldD; ublb=ones(popsize,1)*m; lb=ones(popsize,1)*FieldD(1,:); par=ublb.*aux+lb; [cost,TCSC]=feval(ff,par,busdata1,linedata1,old,anpha1, anpha2); % [cost,ind]=sort(cost); % cost in element sap xep theo thu tu tang dan, % _ while iga