tic???????????????????%%%%%%%%%%%%%%%
clc;clear?all;
%相關原始數據格式說明如下：
%n——節點個數；n1——支路條數；isb——平衡節點號；H——PQ節點個數（為后面形成PVU存儲PV節點初始電壓用）；pr——誤差精度。
%B1——支路參數矩陣，其中第一列和第二列是起始節點編號和終點節點編號第三列、第四列、第五列、第六列分別為：支路電阻、電抗、變壓器變比、電納。（不考慮電導）
%B2——節點參數矩陣，其中第一列和第二列為節點編號和節點類型；第三列到第六列分別為：注入有功、注入無功、電壓幅值、電壓相位。
%節點類型分類如下：“0”為平衡節點“1”為PQ，“2”為PV節點；“3”為PQ（V）節點，“4”為PI節點。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%數據輸入：基準電壓10KV，基準功率10MW?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n=33?;???????
n1=32;
isb=1;
H=32;????????????????%%%%%%%%%%%%%18節點加DG????PQV處理
pr=0.0001;??????????
B1=[1?2?0.00922?0.0047i?1?0;
????2?3?0.00493?0.02511i?1?0;
????3?4?0.0366?0.01864i?1?0;
????4?5?0.03811?0.01941i?1?0;
????5?6?0.0819?0.0707i?1?0;
????6?7?0.01872?0.06188i?1?0;
????7?8?0.07114?0.02351i?1?0;
????8?9?0.103?0.074i?1?0;
????9?10?0.1044?0.074i?1?0;
????10?11?0.01966?0.0065i?1?0;
????11?12?0.03744?0.01238i?1?0;
????12?13?0.1468?0.1155i?1?0;
????13?14?0.05416?0.07129i?1?0;
????14?15?0.05910?0.0526i?1?0;
????15?16?0.07463?0.05450i?1?0;
????16?17?0.1289?0.1721i?1?0;
????17?18?0.0732?0.0574i?1?0;
????2?19?0.0164?0.01565i?1?0;
????19?20?0.15042?0.13554i?1?0;
????20?21?0.04095?0.04784i?1?0;
????21?22?0.07089?0.09373i?1?0;
????3?23?0.04512?0.03083i?1?0;
????23?24?0.08980?0.07091i?1?0;
????24?25?0.08960?0.07011i?1?0;
????6?26?0.0203?0.01034i?1?0;
????26?27?0.02842?0.01447i?1?0;
????27?28?0.1059?0.09337i?1?0;
????28?29?0.08042?0.07006i?1?0;
????29?30?0.05075?0.02585i?1?0;
????30?31?0.09744?0.0963i?1?0;
????31?32?0.03105?0.03619i?1?0;
????32?33?0.03410?0.05302i?1?0];
B2=[1?0?0?0?1.05?0;
????2?1?-0.01?-0.006?1?0;
????3?1?-0.009?-0.004?1?0;
????4?1?-0.012?-0.008?1?0;
????5?1?-0.006?-0.003?1?0;
????6?1?-0.006?-0.002?1?0;
????7?1?-0.02?-0.01?1?0;
????8?1?-0.02?-0.01?1?0;
????9?1?-0.006?-0.002?1?0;
????10?1?-0.006?-0.0035?1?0;
????11?1?-0.0045?-0.003?1?0;
????12?1?-0.006?-0.0035?1?0;
????13?1?-0.006?-0.0035?1?0;
????14?1?-0.012?-0.008?1?0;
????15?1?-0.006?-0.001?1?0;
????16?1?-0.006?-0.002?1?0;
????17?1?-0.006?-0.002?1?0;
????18?3?-0.009+0.1?-0.004?1?0;
????19?1?-0.009?-0.004?1?0;
????20?1?-0.009?-0.004?1?0;
????21?1?-0.009?-0.004?1?0;
????22?1?-0.009?-0.004?1?0;
????23?1?-0.009?-0.005?1?0;
????24?1?-0.042?-0.02?1?0;
????25?1?-0.042?-0.02?1?0;
????26?1?-0.006?-0.0025?1?0;
????27?1?-0.006?-0.0025?1?0;
????28?1?-0.006?-0.002?1?0;
????29?1?-0.012?-0.007?1?0;
????30?1?-0.02?-0.06?1?0;
????31?1?-0.015?-0.007?1?0;
????32?1?-0.021?-0.01?1?0;
????33?4?-0.006?-0.004?1?0];
T=60;%仿真代數???????????????遺傳算法
M=40;%?群體規模
pm=0.1;pc=0.8;%交叉變異概率
Ugmax=1.06;Ugmin=1.0;%參數取值范圍
Timax=1.1;Timin=0.9;
Qcmax=0.05;Qcmin=-0.05;
Vimax=1.07;Vimin=0.93;
Z1=1;%Z2=1;%%懲罰因子%%%%%%%%%%%%%%%&&&&&&&&&&&&&&&&&&&&&&&*********************具體數值待定
%QGmax=;QGmin=;
Long=5;%單個參數字串長度，總編碼長度3L
bestv=-inf;
bval=round（rand（M3*Long））;%初始種群?round?四舍五入取整函數??編碼后的初始種群
%%%%%%%%%%%%%bestv=-inf%最優適應度初值%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%