Thanks Prof Ray...
RegardsSabhan K
Pada Kamis, 19 September 2019 20.33.54 WIB, Ray Zimmerman
<[email protected]> menulis:
You can always use the diary command to save all of the output of a MATLAB
session to a file. It sounds like that may be what you’re looking for.
You could also keep a copy of the results struct for the best solution so far,
replacing it each time you find an improved solution. Then at the end, use
printpf() to re-print the output of the best solution.
Ray
On Sep 18, 2019, at 10:11 PM, Sabhan Kanata <[email protected]> wrote:
Dear Prof Ray Zimmerman,
I did the optimization process using metaheuristics with the calculation of
power flow using matpower. When I use the number of iterations above 5, the
runpf output value is overwritten in the windows command. As a result, the
power flow that shows the minimum PLOSS and VD optimization results is erased
from history. How to withdraw all previous power flow.
Ex: PLOSS and voltage deviation (VD) with the minimum value, namely = 120.5285
MW and 0,18971.
PLOSS #1: Min = 120.5285
Max = 170.736
Range = 50.2075
St.D. = 22.1753
Mean = 131.1021
VD #2:
Min = 0.18971
Max = 0.38363
Range = 0.19392
St.D. = 0.091242
Mean = 0.30149
Example of a matpower command output that does not bring up the minimum PLOSS
and VD values because it is overwritten from the previous iteration output
Newton's method power flow converged in 4 iterations.
Converged in 0.01 seconds
================================================================================
| Bus Data |
================================================================================
Bus Voltage Generation Load
# Mag(pu) Ang(deg) P (MW) Q (MVAr) P (MW) Q (MVAr)
----- ------- -------- -------- -------- -------- --------
1 1.011 1.164 2411.84 7434.88 148.43 29.76
2 1.009 0.848 - - 295.80 225.46
3 0.996 -3.946 - - 449.98 -43.59
4 0.996 -3.501 - - 865.82 -172.55
5 0.998 -3.556 - - 586.13 195.76
6 1.013 -4.729 - - 552.61 45.05
7 1.007 -4.707 - - 459.79 78.89
8 1.037 -3.176 1253.00 2431.39 - -
9 1.000 -3.711 - - 709.19 422.36
10 0.978 -3.237 533.00 -212.17 626.00 290.00
11 0.973 -2.941 447.00 183.51 - -
12 0.949 -2.750 - - 660.94 268.66
13 0.924 5.009 - - 300.89 218.78
14 0.936 10.087 - - 583.82 315.34
15 0.976 19.288 2432.00 802.54 497.00 165.00
16 1.001 15.164 - - 1303.75 545.75
17 1.007 15.382 552.00 499.06 215.20 23.00
18 0.996 -3.473 - - 372.62 -73.09
19 0.965 3.015 - - 301.18 89.10
20 0.925 6.761 - - 1288.57 475.95
21 0.978 16.599 - - 544.44 199.85
22 1.044 24.386 4239.00 787.61 1039.30 209.00
23 1.044 19.334 418.00 1345.14 576.30 338.37
24 0.994 -1.529 - - 827.95 177.43
25 0.981 13.619 - - 220.82 110.72
26 0.936 -1.164 - - 302.54 113.09
27 0.966 7.614 566.00 214.87 - -
28 0.966 7.766 565.00 10.46 - -
29 0.961 6.745 - - 122.80 52.00
30 0.927 7.002 - - - -
31 0.974 1.168 575.00 -6773.36 - -
-------- -------- -------- --------
Total: 13991.84 6723.93 13851.87 4300.09
================================================================================
| Branch Data |
================================================================================
Brnch From To From Bus Injection To Bus Injection Loss (I^2 * Z)
# Bus Bus P (MW) Q (MVAr) P (MW) Q (MVAr) P (MW) Q (MVAr)
----- ----- ----- -------- -------- -------- -------- -------- --------
1 1 2 826.93 270.06 -826.47 -264.87 0.463 5.18
2 1 24 1376.53 378.49 -1369.20 -308.09 7.327 70.40
3 1 31 59.95 6756.58 -57.27 -6506.64 2.678 249.94
4 31 24 632.27 -266.72 -631.91 302.34 0.358 35.62
5 2 5 530.67 39.41 -527.00 -19.65 3.665 41.00
6 4 3 450.29 -40.09 -449.98 43.59 0.312 3.50
7 18 4 147.19 -0.19 -147.18 0.26 0.008 0.07
8 7 5 -445.46 260.20 446.62 -249.00 1.167 11.21
9 8 5 184.51 666.30 -181.75 -639.78 2.760 26.52
10 11 5 179.26 -544.98 -177.89 534.53 1.370 15.32
11 6 7 14.56 341.25 -14.33 -339.09 0.224 2.16
12 8 6 569.75 411.10 -567.17 -386.30 2.581 24.80
13 8 9 498.74 1353.99 -493.28 -1301.52 5.462 52.47
14 10 9 218.16 -857.61 -215.91 879.16 2.243 21.55
15 11 10 311.50 -352.13 -311.16 355.44 0.344 3.31
16 12 11 46.18 -1053.58 -43.76 1080.62 2.417 27.04
17 13 12 85.96 -28.72 -84.68 24.15 1.280 12.30
18 13 26 938.11 -173.31 -927.32 264.85 10.783 103.60
19 26 12 624.78 -377.94 -622.44 400.47 2.344 22.52
20 30 13 1329.07 63.07 -1324.96 -16.75 4.117 46.31
21 14 30 1335.66 135.73 -1329.07 -63.07 6.587 72.66
22 15 14 1935.00 637.54 -1932.03 -345.53 2.967 331.93
23 16 14 310.96 186.41 -306.80 -184.99 4.152 39.89
24 14 20 586.05 71.83 -582.45 -37.26 3.598 34.56
25 16 17 -336.33 -471.57 336.80 476.06 0.469 4.49
26 23 16 1785.11 917.14 -1770.37 -752.30 14.735 164.84
27 18 5 146.15 -177.67 -146.11 178.14 0.049 0.47
28 19 18 673.47 -253.74 -665.97 250.94 7.504 83.96
29 29 19 980.39 -132.86 -974.65 164.64 5.743 64.26
30 27 29 1130.87 223.83 -1129.30 -205.77 1.568 18.06
31 28 27 565.00 10.46 -564.87 -8.96 0.130 1.50
32 20 29 -25.04 -331.16 26.11 286.63 1.070 11.97
33 21 20 692.42 204.53 -681.08 -107.52 11.341 126.88
34 22 21 1254.76 530.59 -1236.86 -404.38 17.898 200.24
35 22 23 1944.94 48.02 -1943.41 89.62 1.533 171.53
36 24 4 1173.16 -171.68 -1168.93 212.39 4.237 40.70
37 25 14 269.23 162.05 -266.69 -167.50 2.537 24.38
38 16 25 492.00 291.50 -490.05 -272.77 1.949 18.73
-------- --------
Total: 139.969 2185.88
PLOSS =
139.9686
sum_of_V_deviations =
0.3173
RegardsSabhan K
Pada Rabu, 18 September 2019 20.21.56 WIB, Ray Zimmerman <[email protected]>
menulis:
Use the printpf() function to print the results of a previously run solution
that you’ve saved in a results struct. For example …
results = runpf('case118');...printpf(results);
Hope this helps,
Ray
On Sep 18, 2019, at 5:50 AM, Sabhan Kanata <[email protected]> wrote:
Dear Prof Ray Zimmerman,
How to display the entire matpower power flow in the matlab command from the
previous iteration that was overwritten?
Thanks
RegardsSabhan K
