%%??????運用牛頓-拉夫遜法解IEEE118?節點潮流

%%%%%%???首先把原始文件分割成兩個獨立的純數據文件分別保存為Voltage.txt和Resis.txt
%%%%%%???計算n節點時，當Voltage.txt和Resis.txt準備完畢后，只需將14改為n，28改為2*n，20改為線路條數即可
%%??讀取Voltage中的數據
clear?
load?Voltage.txt
Bus=Voltage（:1）?;?????????????????%設定節點數
Voltbase=Voltage（:2）?;????????????%設定電壓標準值
VType=Voltage（:5）;????????????????%設定節點類型
V=Voltage（:6）*100;????????????????????%節點電壓標幺值
Angle=Voltage（:7）;????????????????%電壓角
Pload=Voltage（:8）;????????????????%負載有功
Qload=Voltage（:9）;????????????????%負載無功
Pgen=Voltage（:10）;????????????????%發電機發出有功
Qgen=Voltage（:11）;????????????????%發電機發出無功
???????????????????????????????????%電壓設定點
Vset=Voltage（:13）;????????????????%節點所接并聯電容器的電納
Qgmax=Voltage（:14）;???????????????%Q的最大值
Qgmin=Voltage（:15）;???????????????%Q的最小值
Qsh=Voltage（:17）*100;?????????????????%電容
baseMVA?=?100;?????????????????????%功率標準值
%%??讀取Resis中的數據
load?Resis.txt
II=Resis（:1）;??????????????????
JJ=Resis（:2）;??????????????????%兩端點
RType=Resis（:5）;???????????????%線路類型
R=Resis（:6）;???????????????????%兩點間電阻
X=Resis（:7）;???????????????????%兩點間電抗
Btwo=Resis（:8）/2;??????????????%線路兩端電納
K=Resis（:14）;??????????????????%?變壓器變比
Kmax=Resis（:16）;
Kmin=Resis（:17）;?
Kset=Resis（:18）;???????????????%變比要求
tic;
t=clock;
%%???開始求導納?Y
y0=zeros（118）?;??????????????????%兩點間的電阻
yt=zeros（118）?;??????????????????%變壓器之間的導
Qsh=i*Qsh;??????????????????????%電容電抗標幺值
K=1./K?;
for?x=1:179?
????switch?RType（x）?????????????%比較線路類型??
????????case?2??????????????????%含變壓器的
????????????y0（II（x）JJ（x））=1./（K（x）.*（R（x）+i*X（x）））?;???????????????????
????????????y0（JJ（x）II（x））=y0（II（x）JJ（x））????;???????????????
????????????yt（II（x）JJ（x））=（1-K（x））./（（R（x）+i*X（x））.*K（x）^2）;??%II端的電壓?
????????????yt（JJ（x）II（x））=（K（x）-1）./（（R（x）+i*X（x））.*K（x））;????%JJ端的電壓?
????????case?1?????????????????%不含變壓器的
????????????y0（II（x）JJ（x））=1/（R（x）+i*X（x））;?????????????????
????????????y0（JJ（x）II（x））=1/（R（x）+i*X（x））;??????????????????
????end
end
Y=zeros（118）;??????????????????%導納矩陣
clear?x??
for?x=1:118
????Y（xx）=sum（y0（x:）+yt（x:））+i*Qsh（x）/baseMVA?;?????%求出對角導納（不含B/2）
end
Y=Y-y0;
clear?x
for?x=1:179
????Y（II（x）II（x））=Y（II（x）II（x））+i*Btwo（x）?;
????Y（JJ（x）JJ（x））=Y（JJ（x）JJ（x））+i*Btwo（x）?;??????????%加B
end
clear?x???????????????????????????%?Y導納矩陣結束
%%?????
U=zeros（1118）;????????????????????????????%節點電壓
for?x=1:118
????U（x）=V（x）*exp（i*Angle（x）*pi/180）??;???%極坐標電壓
end
clear?x
e=real（U）;
f=imag（U）;???????????%定義實部e跟虛部f?瓜分U
G=real（Y）;
B=imag（Y）;???????????%定義實部G跟虛部B?瓜分Y
D=ones（2361）;???????%dQ?dP?dU^2?的矩陣
Ps=zeros（1118）;?????
Qs=zeros（1118）;?????????????????
N=0;
for?m=1:118
????Ps（m）=（Pgen（m）-Pload（m））/baseMVA;
????Qs（m）=（Qgen（m）-Qload（m））/baseMVA;
end
clear?m???????????????????????
while?max（abs（D））>0.00001????????????????????????%?開始迭代過程求解
????for?m=1:118
????????switch?VType（m）?????????????????????%按照節點類型分情況計算
????????????case?1
????????????????D（2*m-1）=Ps（m）-e（m）*sum（G（m:）.*e-B（m:）.*f）-f（m

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????19411??2018-01-30?17:33??118節點電壓數據.xlsx
?????文件???????27284??2018-01-21?10:20??118節點電阻數據.xlsx
?????文件????????7271??2018-01-31?09:51??ieee118.m