Fast Decouple
Post on 20-Feb-2016
20 Views
Preview:
DESCRIPTION
Transcript
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 1/21
clc Psp=[0;0.2;‐0.45;‐0.4;‐0.60] Qsp=[0;0.2;‐0.15;‐0.5;‐0.1] v=[1.06;1;1;1;1] V=[1.06;1;1;1;1] tol=0.001 ybus=[6.25‐18.695i ‐5+15i ‐1.25+3.75i 0 0 ‐5+15i 10.83334‐32.415i ‐1.6667+5i ‐1.6667+5i ‐2.5+7.5i ‐1.25+3.75i ‐1.6667+5i 12.91667‐38.695i ‐10+30i 0 0 ‐1.6667+5i ‐10+30i 12.91667‐38.695i ‐1.25+3.75i 0 ‐2.5+7.5i 0 ‐1.25+3.75i 3.75‐11.21i] B=‐imag(ybus(2:1:5,2:1:5)) delP=1 delQ=1 iter=0 delt=zeros(5,1) while abs(max(delP))>=tol||abs(max(delQ))>=tol iter=iter+1 I=ybus*V Scal=(conj(V)).*I delP=Psp(2:1:5)‐real(Scal(2:1:5)) delQ=Qsp(2:1:5)‐(‐imag(Scal(2:1:5))) if abs(max(delP))>=tol deldelt=inv(B)*delP./v(2:1:5) delt(2:1:5)=delt(2:1:5)+deldelt end if abs(max(delQ))>=tol delv=inv(B)*delQ./v(2:1:5) v(2:1:5)=v(2:1:5)+delv end for t=1:5 V(t)=v(t)*exp(i*delt(t)) end
end
Psp =
0 0.2000 ‐0.4500 ‐0.4000 ‐0.6000
Qsp =
0 0.2000 ‐0.1500 ‐0.5000 ‐0.1000
v =
1.0600
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 2/21
1.0000 1.0000 1.0000 1.0000
V =
1.0600 1.0000 1.0000 1.0000 1.0000
tol =
1.0000e‐03
ybus =
Columns 1 through 4
6.2500 ‐18.6950i ‐5.0000 +15.0000i ‐1.2500 + 3.7500i 0.0000 + 0.0000i ‐5.0000 +15.0000i 10.8333 ‐32.4150i ‐1.6667 + 5.0000i ‐1.6667 + 5.0000i ‐1.2500 + 3.7500i ‐1.6667 + 5.0000i 12.9167 ‐38.6950i ‐10.0000 +30.0000i 0.0000 + 0.0000i ‐1.6667 + 5.0000i ‐10.0000 +30.0000i 12.9167 ‐38.6950i 0.0000 + 0.0000i ‐2.5000 + 7.5000i 0.0000 + 0.0000i ‐1.2500 + 3.7500i
Column 5
0.0000 + 0.0000i ‐2.5000 + 7.5000i 0.0000 + 0.0000i ‐1.2500 + 3.7500i 3.7500 ‐11.2100i
B =
32.4150 ‐5.0000 ‐5.0000 ‐7.5000 ‐5.0000 38.6950 ‐30.0000 0 ‐5.0000 ‐30.0000 38.6950 ‐3.7500 ‐7.5000 0 ‐3.7500 11.2100
delP =
1
delQ =
1
iter =
0
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 3/21
delt =
0 0 0 0 0
iter =
1
I =
0.3750 ‐ 1.0667i ‐0.3001 + 0.9850i ‐0.0750 + 0.2800i ‐0.0000 + 0.0550i 0.0000 + 0.0400i
Scal =
0.3975 ‐ 1.1307i ‐0.3001 + 0.9850i ‐0.0750 + 0.2800i ‐0.0000 + 0.0550i 0.0000 + 0.0400i
delP =
0.5001 ‐0.3750 ‐0.4000 ‐0.6000
delQ =
1.1850 0.1300 ‐0.4450 ‐0.0600
deldelt =
‐0.0378 ‐0.0868 ‐0.0932 ‐0.1100
delt =
0 ‐0.0378
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 4/21
‐0.0868 ‐0.0932 ‐0.1100
delv =
0.0494 0.0205 0.0139 0.0323
v =
1.0600 1.0494 1.0205 1.0139 1.0323
V =
1.0600 1.0000 1.0000 1.0000 1.0000
V =
1.0600 + 0.0000i 1.0486 ‐ 0.0396i 1.0000 + 0.0000i 1.0000 + 0.0000i 1.0000 + 0.0000i
V =
1.0600 + 0.0000i 1.0486 ‐ 0.0396i 1.0167 ‐ 0.0885i 1.0000 + 0.0000i 1.0000 + 0.0000i
V =
1.0600 + 0.0000i 1.0486 ‐ 0.0396i 1.0167 ‐ 0.0885i 1.0095 ‐ 0.0944i 1.0000 + 0.0000i
V =
1.0600 + 0.0000i
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 5/21
1.0486 ‐ 0.0396i 1.0167 ‐ 0.0885i 1.0095 ‐ 0.0944i 1.0261 ‐ 0.1133i
iter =
2
I =
1.0376 + 0.0339i 0.5966 ‐ 0.1052i ‐0.4313 + 0.0295i ‐0.5305 + 0.4017i ‐0.6546 ‐ 0.0601i
Scal =
1.0998 + 0.0359i 0.6298 ‐ 0.0866i ‐0.4411 ‐ 0.0081i ‐0.5734 + 0.3555i ‐0.6649 ‐ 0.1358i
delP =
‐0.4298 ‐0.0089 0.1734 0.0649
delQ =
0.1134 ‐0.1581 ‐0.1445 ‐0.2358
deldelt =
‐0.0128 ‐0.0000 0.0025 ‐0.0023
delt =
0 ‐0.0506 ‐0.0869 ‐0.0907 ‐0.1123
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 6/21
delv =
‐0.0185 ‐0.0369 ‐0.0395 ‐0.0459
v =
1.0600 1.0308 0.9837 0.9745 0.9864
V =
1.0600 + 0.0000i 1.0486 ‐ 0.0396i 1.0167 ‐ 0.0885i 1.0095 ‐ 0.0944i 1.0261 ‐ 0.1133i
V =
1.0600 + 0.0000i 1.0295 ‐ 0.0521i 1.0167 ‐ 0.0885i 1.0095 ‐ 0.0944i 1.0261 ‐ 0.1133i
V =
1.0600 + 0.0000i 1.0295 ‐ 0.0521i 0.9800 ‐ 0.0853i 1.0095 ‐ 0.0944i 1.0261 ‐ 0.1133i
V =
1.0600 + 0.0000i 1.0295 ‐ 0.0521i 0.9800 ‐ 0.0853i 0.9705 ‐ 0.0883i 1.0261 ‐ 0.1133i
V =
1.0600 + 0.0000i 1.0295 ‐ 0.0521i 0.9800 ‐ 0.0853i 0.9705 ‐ 0.0883i 0.9802 ‐ 0.1105i
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 7/21
iter =
3
I =
1.3542 ‐ 0.3320i 0.1593 ‐ 0.3667i ‐0.4800 + 0.1842i ‐0.3873 + 0.6079i ‐0.6279 + 0.1990i
Scal =
1.4355 ‐ 0.3520i 0.1831 ‐ 0.3693i ‐0.4861 + 0.1395i ‐0.4295 + 0.5557i ‐0.6375 + 0.1257i
delP =
0.0169 0.0361 0.0295 0.0375
delQ =
‐0.1693 ‐0.0105 0.0557 0.0257
deldelt =
0.0057 0.0093 0.0099 0.0106
delt =
0 ‐0.0449 ‐0.0775 ‐0.0808 ‐0.1016
delv =
‐0.0060
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 8/21
‐0.0019 ‐0.0010 ‐0.0022
v =
1.0600 1.0249 0.9818 0.9735 0.9842
V =
1.0600 + 0.0000i 1.0295 ‐ 0.0521i 0.9800 ‐ 0.0853i 0.9705 ‐ 0.0883i 0.9802 ‐ 0.1105i
V =
1.0600 + 0.0000i 1.0238 ‐ 0.0460i 0.9800 ‐ 0.0853i 0.9705 ‐ 0.0883i 0.9802 ‐ 0.1105i
V =
1.0600 + 0.0000i 1.0238 ‐ 0.0460i 0.9789 ‐ 0.0761i 0.9705 ‐ 0.0883i 0.9802 ‐ 0.1105i
V =
1.0600 + 0.0000i 1.0238 ‐ 0.0460i 0.9789 ‐ 0.0761i 0.9703 ‐ 0.0786i 0.9802 ‐ 0.1105i
V =
1.0600 + 0.0000i 1.0238 ‐ 0.0460i 0.9789 ‐ 0.0761i 0.9703 ‐ 0.0786i 0.9791 ‐ 0.0999i
iter =
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 9/21
4
I =
1.2577 ‐ 0.4635i 0.1258 ‐ 0.1888i ‐0.4452 + 0.2057i ‐0.3418 + 0.5582i ‐0.5803 + 0.1801i
Scal =
1.3332 ‐ 0.4913i 0.1375 ‐ 0.1875i ‐0.4514 + 0.1675i ‐0.3755 + 0.5147i ‐0.5862 + 0.1184i
delP =
0.0625 0.0014 ‐0.0245 ‐0.0138
delQ =
0.0125 0.0175 0.0147 0.0184
deldelt =
0.0017 ‐0.0002 ‐0.0006 ‐0.0003
delt =
0 ‐0.0432 ‐0.0777 ‐0.0814 ‐0.1019
delv =
0.0030 0.0048 0.0051 0.0054
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 10/21
v =
1.0600 1.0279 0.9866 0.9785 0.9897
V =
1.0600 + 0.0000i 1.0238 ‐ 0.0460i 0.9789 ‐ 0.0761i 0.9703 ‐ 0.0786i 0.9791 ‐ 0.0999i
V =
1.0600 + 0.0000i 1.0269 ‐ 0.0444i 0.9789 ‐ 0.0761i 0.9703 ‐ 0.0786i 0.9791 ‐ 0.0999i
V =
1.0600 + 0.0000i 1.0269 ‐ 0.0444i 0.9836 ‐ 0.0766i 0.9703 ‐ 0.0786i 0.9791 ‐ 0.0999i
V =
1.0600 + 0.0000i 1.0269 ‐ 0.0444i 0.9836 ‐ 0.0766i 0.9753 ‐ 0.0796i 0.9791 ‐ 0.0999i
V =
1.0600 + 0.0000i 1.0269 ‐ 0.0444i 0.9836 ‐ 0.0766i 0.9753 ‐ 0.0796i 0.9845 ‐ 0.1007i
iter =
5
I =
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 11/21
1.2140 ‐ 0.4066i 0.1960 ‐ 0.1783i ‐0.4391 + 0.1881i ‐0.3630 + 0.5337i ‐0.5916 + 0.1558i
Scal =
1.2868 ‐ 0.4310i 0.2092 ‐ 0.1744i ‐0.4463 + 0.1514i ‐0.3965 + 0.4916i ‐0.5981 + 0.0938i
delP =
‐0.0092 ‐0.0037 ‐0.0035 ‐0.0019
delQ =
0.0256 0.0014 ‐0.0084 ‐0.0062
deldelt =
‐0.0009 ‐0.0012 ‐0.0013 ‐0.0012
delt =
0 ‐0.0441 ‐0.0789 ‐0.0827 ‐0.1031
delv =
1.0e‐03 *
0.7809 0.1272 ‐0.0182 ‐0.0222
v =
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 12/21
1.0600 1.0287 0.9867 0.9785 0.9896
V =
1.0600 + 0.0000i 1.0269 ‐ 0.0444i 0.9836 ‐ 0.0766i 0.9753 ‐ 0.0796i 0.9845 ‐ 0.1007i
V =
1.0600 + 0.0000i 1.0277 ‐ 0.0454i 0.9836 ‐ 0.0766i 0.9753 ‐ 0.0796i 0.9845 ‐ 0.1007i
V =
1.0600 + 0.0000i 1.0277 ‐ 0.0454i 0.9836 ‐ 0.0778i 0.9753 ‐ 0.0796i 0.9845 ‐ 0.1007i
V =
1.0600 + 0.0000i 1.0277 ‐ 0.0454i 0.9836 ‐ 0.0778i 0.9752 ‐ 0.0808i 0.9845 ‐ 0.1007i
V =
1.0600 + 0.0000i 1.0277 ‐ 0.0454i 0.9836 ‐ 0.0778i 0.9752 ‐ 0.0808i 0.9844 ‐ 0.1019i
iter =
6
I =
1.2293 ‐ 0.3890i
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 13/21
0.1942 ‐ 0.2072i ‐0.4426 + 0.1855i ‐0.3689 + 0.5415i ‐0.5954 + 0.1620i
Scal =
1.3030 ‐ 0.4124i 0.2089 ‐ 0.2041i ‐0.4498 + 0.1481i ‐0.4035 + 0.4983i ‐0.6026 + 0.0988i
delP =
‐0.0089 ‐0.0002 0.0035 0.0026
delQ =
‐0.0041 ‐0.0019 ‐0.0017 ‐0.0012
deldelt =
1.0e‐03 *
‐0.2122 0.0632 0.1255 0.1275
delt =
0 ‐0.0443 ‐0.0789 ‐0.0826 ‐0.1030
delv =
1.0e‐03 *
‐0.4497 ‐0.6055 ‐0.6428 ‐0.6361
v =
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 14/21
1.0600 1.0282 0.9861 0.9779 0.9890
V =
1.0600 + 0.0000i 1.0277 ‐ 0.0454i 0.9836 ‐ 0.0778i 0.9752 ‐ 0.0808i 0.9844 ‐ 0.1019i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0456i 0.9836 ‐ 0.0778i 0.9752 ‐ 0.0808i 0.9844 ‐ 0.1019i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0456i 0.9830 ‐ 0.0777i 0.9752 ‐ 0.0808i 0.9844 ‐ 0.1019i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0456i 0.9830 ‐ 0.0777i 0.9745 ‐ 0.0807i 0.9844 ‐ 0.1019i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0456i 0.9830 ‐ 0.0777i 0.9745 ‐ 0.0807i 0.9838 ‐ 0.1017i
iter =
7
I =
1.2349 ‐ 0.3973i
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 15/21
0.1835 ‐ 0.2062i ‐0.4433 + 0.1875i ‐0.3657 + 0.5446i ‐0.5928 + 0.1641i
Scal =
1.3090 ‐ 0.4211i 0.1979 ‐ 0.2035i ‐0.4504 + 0.1499i ‐0.4004 + 0.5012i ‐0.5999 + 0.1011i
delP =
0.0021 0.0004 0.0004 ‐0.0001
delQ =
‐0.0035 ‐0.0001 0.0012 0.0011
deldelt =
1.0e‐03 *
0.1391 0.1509 0.1597 0.1379
delt =
0 ‐0.0442 ‐0.0787 ‐0.0824 ‐0.1029
delv =
1.0e‐04 *
‐0.8391 0.1473 0.3836 0.5767
v =
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 16/21
1.0600 1.0281 0.9861 0.9779 0.9891
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0456i 0.9830 ‐ 0.0777i 0.9745 ‐ 0.0807i 0.9838 ‐ 0.1017i
V =
1.0600 + 0.0000i 1.0271 ‐ 0.0454i 0.9830 ‐ 0.0777i 0.9745 ‐ 0.0807i 0.9838 ‐ 0.1017i
V =
1.0600 + 0.0000i 1.0271 ‐ 0.0454i 0.9831 ‐ 0.0775i 0.9745 ‐ 0.0807i 0.9838 ‐ 0.1017i
V =
1.0600 + 0.0000i 1.0271 ‐ 0.0454i 0.9831 ‐ 0.0775i 0.9746 ‐ 0.0805i 0.9838 ‐ 0.1017i
V =
1.0600 + 0.0000i 1.0271 ‐ 0.0454i 0.9831 ‐ 0.0775i 0.9746 ‐ 0.0805i 0.9838 ‐ 0.1016i
iter =
8
I =
1.2325 ‐ 0.3993i
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 17/21
0.1847 ‐ 0.2020i ‐0.4430 + 0.1877i ‐0.3650 + 0.5434i ‐0.5926 + 0.1628i
Scal =
1.3064 ‐ 0.4232i 0.1988 ‐ 0.1991i ‐0.4500 + 0.1502i ‐0.3995 + 0.5002i ‐0.5996 + 0.1000i
delP =
0.0012 0.0000 ‐0.0005 ‐0.0004
delQ =
1.0e‐03 *
0.8693 0.2175 0.1896 0.0128
deldelt =
1.0e‐04 *
0.2158 ‐0.1669 ‐0.2619 ‐0.3101
delt =
0 ‐0.0442 ‐0.0787 ‐0.0824 ‐0.1029
V =
1.0600 + 0.0000i 1.0271 ‐ 0.0454i 0.9831 ‐ 0.0775i 0.9746 ‐ 0.0805i 0.9838 ‐ 0.1016i
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 18/21
V =
1.0600 + 0.0000i 1.0271 ‐ 0.0454i 0.9831 ‐ 0.0775i 0.9746 ‐ 0.0805i 0.9838 ‐ 0.1016i
V =
1.0600 + 0.0000i 1.0271 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9746 ‐ 0.0805i 0.9838 ‐ 0.1016i
V =
1.0600 + 0.0000i 1.0271 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9746 ‐ 0.0805i 0.9838 ‐ 0.1016i
V =
1.0600 + 0.0000i 1.0271 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9746 ‐ 0.0805i 0.9838 ‐ 0.1016i
iter =
9
I =
1.2322 ‐ 0.3994i 0.1858 ‐ 0.2017i ‐0.4429 + 0.1877i ‐0.3655 + 0.5433i ‐0.5931 + 0.1627i
Scal =
1.3061 ‐ 0.4233i 0.2000 ‐ 0.1988i ‐0.4500 + 0.1502i ‐0.4000 + 0.5000i ‐0.6000 + 0.0998i
delP =
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 19/21
1.0e‐04 *
‐0.3505 0.0053 0.0060 ‐0.0427
delQ =
0.0012 0.0002 0.0000 ‐0.0002
delv =
1.0e‐04 *
0.7005 0.6736 0.6818 0.5770
v =
1.0600 1.0282 0.9862 0.9780 0.9891
V =
1.0600 + 0.0000i 1.0271 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9746 ‐ 0.0805i 0.9838 ‐ 0.1016i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9746 ‐ 0.0805i 0.9838 ‐ 0.1016i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9746 ‐ 0.0805i 0.9838 ‐ 0.1016i
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 20/21
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9747 ‐ 0.0805i 0.9838 ‐ 0.1016i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9747 ‐ 0.0805i 0.9839 ‐ 0.1016i
iter =
10
I =
1.2318 ‐ 0.3980i 0.1862 ‐ 0.2029i ‐0.4429 + 0.1875i ‐0.3655 + 0.5432i ‐0.5931 + 0.1629i
Scal =
1.3057 ‐ 0.4219i 0.2005 ‐ 0.1999i ‐0.4500 + 0.1500i ‐0.4000 + 0.5000i ‐0.6001 + 0.1000i
delP =
1.0e‐03 *
‐0.5003 ‐0.0348 0.0267 0.1096
delQ =
1.0e‐04 *
0.5049 0.1063 ‐0.0742 ‐0.2469
7/15/2015 fast_decouple
file:///C:/Users/Siddhu/Documents/MATLAB/html/fast_decouple.html 21/21
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9747 ‐ 0.0805i 0.9839 ‐ 0.1016i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9747 ‐ 0.0805i 0.9839 ‐ 0.1016i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9747 ‐ 0.0805i 0.9839 ‐ 0.1016i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9747 ‐ 0.0805i 0.9839 ‐ 0.1016i
V =
1.0600 + 0.0000i 1.0272 ‐ 0.0454i 0.9831 ‐ 0.0776i 0.9747 ‐ 0.0805i 0.9839 ‐ 0.1016i
Published with MATLAB® R2014a
top related