plotvmpoint.m

這是螞蟻系統演算法裡面更新費洛蒙以及新的搜尋機制的程式碼

clear
clc
Vmpoint=[0.99801
0.998
0.9815
0.96228
0.9815
0.99037
0.97953
0.97274
0.96673
0.97087
0.95915
0.95204
0.94672
0.94977
0.95028
0.95133
0.95375
0.94543
0.9658
0.94715
0.98742
0.99329
0.96679
0.95487
];
VmpointAS=[1.0659
1.0489
1.0467
1.0324
1.0631
1.0244
1.0155
1.0036
1.0022
1.0118
1.0151
0.99749
0.99797
1.0038
1.0128
1.0136
0.99995
1.0057
1.006
0.98812
1.0149
1.0465
0.99483
0.98326
];
VmpointACS=[1.0604
1.0491
1.0436
1.0305
0.99754
1.0029
1.0064
0.9944
0.99301
0.9974
0.99567
0.98379
0.98161
0.98611
0.9952
0.9974
0.99499
1.0064
1.0057
0.98783
1.0139
1.042
0.99384
0.98225
];
inter=[3
4
6
7
9
10
12
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
];
plot(1:1:length(Vmpoint),Vmpoint,'bo');
hold on
plot(1:1:length(VmpointAS),VmpointAS,'g*');
hold on
plot(1:1:length(VmpointACS),VmpointACS,'rsquare');
set(gca, 'xtick', [1:length(Vmpoint)]);
set(gca, 'xticklabel', [inter]);
axis([-inf,inf,0.94,1.08]);
xlabel('P-Q Bus NO.');
ylabel('V(p.u.)');
legend('程ㄎて玡','AS','ACS');