求助大神用matlab编写万有引力搜索算法（GSA）求解电力系统无功优化的程序，感谢！！俺的邮箱是[e-mail]

[Danielxue]

求助大神用matlab编写万有引力搜索算法（GSA）求解电力系统无功优化的程序，感谢！！俺的邮箱是[e-mail]

万有引力算法 Gravitational Search algorithm
源程序见：http://download.csdn.net/detail/xionghui0501/4239995

%测试系统为21节点系统，相关原始数据格式说明如下：
%n——节点个数；n1——支路条数；isb——平衡节点号；H——PQ节点个数（为后面形成PVU存储PV节点初始电压用）；pr——误差精度。
%B1——支路参数矩阵，其中第一列和第二列是起始节点编号和终点节点编号,第三列、第四列、第五列、第六列分别为：支路电阻、电抗、变压器变比、电纳。（不考虑电导）
%B2——节点参数矩阵，其中第一列和第二列为节点编号和节点类型；第三列到第六列分别为：注入有功、注入无功、电压幅值、电压相位。
%节点类型分类如下：“0”为平衡节点,“1”为PQ，“2”为PV节点；“3”为PQ（V）节点，“4”为PI节点。
n=21 ;
n1=20;
isb=1;
H=20;
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;
2 11 0.0164 0.01565i 1 0;
11 12 0.15042 0.13554i 1 0;
12 13 0.04095 0.04784i 1 0;
13 14 0.07089 0.09373i 1 0;
3 15 0.04512 0.03083i 1 0;
15 16 0.08980 0.07091i 1 0;
16 17 0.08960 0.07011i 1 0;
6 18 0.0203 0.01034i 1 0;
18 19 0.02842 0.01447i 1 0;
19 20 0.1059 0.09337i 1 0;
20 21 0.08042 0.07006i 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 3 -0.006 -0.0035 1 0;
11 1 -0.009 -0.004 1 0;
12 1 -0.009 -0.004 1 0;
13 1 -0.009 -0.004 1 0;
14 3 -0.009 -0.004 1 0;
15 1 -0.009 -0.005 1 0;
16 1 -0.042 -0.02 1 0;
17 4 -0.042 -0.02 1 0;
18 1 -0.006 -0.0025 1 0;
19 1 -0.006 -0.0025 1 0;
20 1 -0.006 -0.002 1 0;
21 4 -0.012 -0.007 1 0];