%data支路信息：
%第一列存儲存儲支路序號
%第二列存儲支路得一個端點
%第三例存儲支路得另一個端點
%第四利存儲支路電阻
%第五列存儲支路電抗
%第六列存儲支路對地電導
%第七列存儲支路對地電納
%第八列存儲變壓器變比

%B2pie節點信息：
%第一列為節點編號
%第二列為發電機有功功率
%第三列為發電機無功功率
%第四列為負荷有功功率
%第五列為負荷無功功率
%第六列為節點電壓初始值
%第七列為節點電壓相角初始值
%第八例為該點得類型：1為PQ節點，2為PV節點，3為平衡節點負荷節點排在前面，PV節點排在后面，平衡節點最后
%B3=input（‘\n請輸入對地導納信息：B3=‘）;
%B4=input（‘\n請輸入阻抗信息：B4=‘）;
%k為迭代次數；
clc
clear
%?close?all
?[dataY1]=xlsread（‘G:\電力系統課程潮流資料\inputdata2.xls‘）;
%?Filepath=‘e:\PQmethoddata\inputdata.txt‘;
%?data=load（Filepath）;
nbus=data（11）;??????????????%節點數
nline=data（12）;??????????%支路數
ngen=data（13）;??????????????%發電機節點數
nload=data（14）;???????????????%?負荷節點數
basemva=data（15）;??????????%?容量基準值?
jingdu=data（16）;%計算精度
detaY=ones（nbus+nload-11）;
B1pie=data（2:nline+11:8）
for?i=1:nline
????B1（i1:3）=B1pie（i1:3）;
????B1（i4）=B1pie（i4）+B1pie（i5）*sqrt（-1）;
????B1（i5）=B1pie（i6）+B1pie（i7）*sqrt（-1）;
????B1（i6）=B1pie（i8）;
end
Y=zeros（nbus）;
for?a=1:nline
????p=B1（a2）;
????q=B1（a3）;
????Y（pp）=Y（pp）+1./B1（a4）+B1（a5）./2;
????Y（pq）=Y（pq）-1./（B1（a4）*B1（a6））;
????Y（qp）=Y（pq）;
????Y（qq）=Y（qq）+1./（B1（a4）*B1（a6）^2）+B1（a5）./2;
end
Y?????????%節點導納陣中線路對地導納是正嗎？
G=real（Y）;
B=imag（Y）;
B2pie=data（nline+2:nline+1+nbus1:8）;
for?i=1:nbus
B2（i1）=（B2pie（i2）-B2pie（i4）+（B2pie（i3）-B2pie（i5））*sqrt（-1））/100;
B2（i2）=B2pie（i6）*cos（B2pie（i7））+（B2pie（i6）*sin（B2pie（i7）））*sqrt（-1）;
B2（i3）=B2pie（i8）;
end
for?i=1:nbus
????V（i）=abs（B2（i2））;???????%各節點電壓幅值
????phase（i）=angle（B2（i2））;?%各節點電壓相角
????P0（i）=real（B2（i1））????%各節點的注入有功功率
????Q0（i）=imag（B2（i1））????%各節點的注入無功功率
end
%根據初始值計算各節點注入功率及功率不平衡量
ndiedai=0;
while?max（abs（detaY））>jingdu
????ndiedai=ndiedai+1;
?Q=zeros（nbus1）;
?P=zeros（nbus1）;
?for?i=1:nbus
???for??j=1:nbus
???P（i）=V（i）*V（j）*（G（ij）*cos（phase（i）-phase（j））+B（ij）*sin（phase（i）-phase（j）））+P（i）;
???Q（i）=V（i）*V（j）*（G（ij）*sin（phase（i）-phase（j））-B（ij）*cos（phase（i）-phase（j）））+?Q（i）;
???end
?end
for?i=1:nbus
?????
????switch?B2（i3）
????????case?1
???????????
??????????????detaP（i）=P0（i）-P（i）;
??????????????detaQ（i）=Q0（i）-Q（i）;
????????case?2
?????????????
?????????????detaP（i）=P0（i）-P（i）;
????????otherwise
????????pinghengjd=i;
????end
end
detaP=remove（detaPpinghengjd）;%?在1行，1列向量或方陣中去除平衡節點
%?牛頓法的因變量不平衡量組成的向量????
?for?i=1:nbus-1???
detaY（i）=detaP（i）;????
?end
?for?i=1:（nbus-ngen）
?????detaY（i+nbus-1）=detaQ（i）
?end
?
??????????
%接著生成雅可比矩陣
for?i=1:（nbus-1）
????for?j=1:（nbus-1）
??????if?i~=j
??????????H（ij）=-V（i）*V（j）*（G（ij）*sin（phase（i）-phase（j））-B（ij）*cos（phase（i）-phase（j）））;
??????else
??????????H（ii）=V（i）^2*B（ii）+Q（i）;
??????end
????end
end

for?i=1:（nbus-1）
????for?j=1:（nbus-ngen）
????????if?i~=j
?????????N（ij）=-V（i）*（G（ij）*cos（phase（i）-phase（j））+B（ij）*sin（phase（i）-phase（j）））;
????????else
?????????N（ii）=（-V（i）^2*G（ii）-P（i））/V（i）;
????????end
????end
end

for?i=1:（nbus-ngen）
????for?j=1:（nbus-1）
????????if?i~=j
?????????J（ij）=V（i）*V（j）*（G（ij）*cos（pha

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????139776??2012-01-01?19:32??IEEE標準測試系統\003_bus\003readme.doc

?????文件???????2267??2012-01-01?19:32??IEEE標準測試系統\003_bus\bpa\003bpa.DAT

?????文件????????495??2012-01-01?19:32??IEEE標準測試系統\003_bus\bpa\003bpaswi.DAT

?????文件????????828??2012-01-01?19:32??IEEE標準測試系統\003_bus\eurostag\003edfech.DAT

?????文件???????2346??2012-01-01?19:32??IEEE標準測試系統\003_bus\ieee\003ieee.dat

?????文件????????513??2012-01-01?19:32??IEEE標準測試系統\003_bus\th\003thdat.dat

?????文件????????221??2012-01-01?19:32??IEEE標準測試系統\003_bus\th\003thmac.dat

?????文件????????313??2012-01-01?19:32??IEEE標準測試系統\003_bus\th\003thres.dat

?????文件???????2086??2012-01-01?19:32??IEEE標準測試系統\003_bus.png

?????文件??????53248??2012-01-01?19:32??IEEE標準測試系統\004_bus\004readme.doc

?????文件???????2649??2012-01-01?19:32??IEEE標準測試系統\004_bus\bpa\004bpa.dat

?????文件??????45814??2012-01-01?19:32??IEEE標準測試系統\004_bus\bpa\004bpa.pfo

?????文件??????????1??2012-01-01?19:32??IEEE標準測試系統\004_bus\bpa\PWRFLO.DIS

?????文件????????927??2012-01-01?19:32??IEEE標準測試系統\004_bus\eurostag\004edfech.dat

?????文件???????2647??2012-01-01?19:32??IEEE標準測試系統\004_bus\ieee\004ieee.dat

?????文件????????381??2012-01-01?19:32??IEEE標準測試系統\004_bus\th\004thdat.dat

?????文件????????360??2012-01-01?19:32??IEEE標準測試系統\004_bus\th\004thres.dat

?????文件???????3272??2012-01-01?19:32??IEEE標準測試系統\004_bus.png

?????文件??????74752??2012-01-01?19:32??IEEE標準測試系統\005_bus\005readme.doc

?????文件???????3471??2012-01-01?19:32??IEEE標準測試系統\005_bus\bpa\005bpa.DAT

?????文件???????1208??2012-01-01?19:32??IEEE標準測試系統\005_bus\eurostag\005edfech.DAT

?????文件???????3550??2012-01-01?19:32??IEEE標準測試系統\005_bus\ieee\005ieee.dat

?????文件????????516??2012-01-01?19:32??IEEE標準測試系統\005_bus\th\005thdat.dat

?????文件????????410??2012-01-01?19:32??IEEE標準測試系統\005_bus\th\005thres.dat

?????文件???????4590??2012-01-01?19:32??IEEE標準測試系統\005_bus.png

?????文件?????143360??2012-01-01?19:32??IEEE標準測試系統\009_bus\009readme.doc

?????文件???????5907??2012-01-01?19:32??IEEE標準測試系統\009_bus\bpa\009bpa.DAT

?????文件????????624??2012-01-01?19:32??IEEE標準測試系統\009_bus\bpa\009bpaswi.DAT

?????文件???????1267??2012-01-01?19:32??IEEE標準測試系統\009_bus\eurostag\009edfdta.DAT

?????文件???????2051??2012-01-01?19:32??IEEE標準測試系統\009_bus\eurostag\009edfech.DAT

............此處省略241個文件信息