%?IEEE-9標準試驗系統潮流求解，直角坐標系求解
%?1號節點為平衡節點
%?基準功率100MW，基準電壓230KV

clear

%?time=zeros（1001）;
%?for?t=1:100
%?tic

%%?讀Data1中數據
load?ieee9_node.txt
Bus=ieee9_node（:1）;????????????????????%節點號
Vtype=ieee9_node（:3）;??????????????????%節點類型
Pload=ieee9_node（:4）;??????????????????%負載有功功率
Qload=ieee9_node（:5）;??????????????????%負載無功功率
Pgen=ieee9_node（:6）;???????????????????%發電機發出的有用功率
Qgen=ieee9_node（:7）;???????????????????%發電機發出的無功功率
Vset=ieee9_node（:9）;???????????????????%電壓設定值
Qsh=ieee9_node（:12）;???????????????????%并聯電容器的電納（標幺值）

swing=find（Vtype==3）;???????????????????%平衡節點編號
PV=find（Vtype==2）;??????????????????????%PV節點編號
PQ=find（Vtype==1）;??????????????????????%PQ節點編號
PVPQ=find（Vtype<3）;?????????????????????%PV/PQ節點編號
%%?讀取Data2中的數據
load?ieee9_branch.txt
II=ieee9_branch（:1）;??????????????????????%支路首端節點
JJ=ieee9_branch（:2）;??????????????????????%支路末端節點
Ltype=ieee9_branch（:3）;???????????????????%支路類型
R=ieee9_branch（:4）;???????????????????????%支路電阻
X=ieee9_branch（:5）;???????????????????????%支路電抗
B_=ieee9_branch（:6）;??????????????????????%1/2線路電納
K=ieee9_branch（:7）;???????????????????????%變壓器非標準變比
%%?節點導納矩陣
n=length（Bus）;??????????????????????%節點數量
y1=zeros（n）;????????????????????????%支路/變壓器導納
y3=zeros（n）;????????????????????????%支路對地導納

y1（（JJ-1）*n+II）=1./（（R+1i*X）.*K）;
y1（（II-1）*n+JJ）=y1（（JJ-1）*n+II）;
y3（（JJ-1）*n+II）=（1-K）./（（R+1i*X）.*K.^2）+1i*B_;
y3（（II-1）*n+JJ）=（K-1）./（（R+1i*X）.*K）+1i*B_;
Y=diag（sum（y12）+sum（y32）+1i*Qsh）;???????%節點導納矩陣
Y=Y-y1;
G=real（Y）;??????????????????????????????????????????%節點電導
B=imag（Y）;??????????????????????????????????????????%節點電納
%%?設電壓初始值
U=Vset;
e=real（U）;
f=imag（U）;
%%?牛頓-拉夫遜法迭代求解
Ps=Pgen-Pload;??????????????????????????????????????%節點注入有功功率
Qs=Qgen-Qload;??????????????????????????????????????%節點注入無功功率
D=ones（2*（n-1）1）;??????????????????????????????????%節點功率不平衡量
N=0;????????????????????????????????????????????????%迭代次數
J=zeros（2*（n-1））;???????????????????????????????????%雅克比矩陣
tolerance=1e-5;?????????????????????????????????????%誤差精度
while?max（abs（D））>tolerance
????%計算節點功率和節點電壓平方的不平衡量（不包括平衡節點）
????
????D（2*PVPQ-3）=Ps（PVPQ）-e（PVPQ）.*（G（PVPQ:）*e-B（PVPQ:）*f）-f（PVPQ）.*（G（PVPQ:）*f+B（PVPQ:）*e）;?????%有功不平衡
????D（2*PQ-2）=Qs（PQ）-f（PQ）.*（G（PQ:）*e-B（PQ:）*f）+e（PQ）.*（G（PQ:）*f+B（PQ:）*e）;?????????????????????%PQ節點無功不平衡
????D（2*PV-2）=?Vset（PV）.^2-（e（PV）.^2+f（PV）.^2）;?????????????????????????????????????????????????????%PV節點電壓不平衡
????%計算雅克比矩陣
????for?a=1:n-1
????????j=PVPQ（a）;
????????for?b=1:n-1
????????????k=PVPQ（b）;
????????????if?j~=k
????????????????J（2*j-32*k-3）=-B（jk）*e（j）+G（jk）*f（j）;
????????????????J（2*j-32*k-2）=G（jk）*e（j）+B（jk）*f（j）;
????????????????if?Vtype（j）==1?????????????????????%PQ節點
????????????????????J（2*j-22*k-3）=-J（2*j-32*k-2）;
????????????????????J（2*j-22*k-2）=J（2*j-32*k-3）;
????????????????end
????????????else
????????????????J（2*j-32*k-3）=-B（jk）*e（j）+G（jk）*f（j）+G（j:）*f+B（j:）*e;
????????????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????4982??2018-11-15?11:00??ieee9_1.m
?????文件????????5883??2018-11-15?10:50??ieee9_12.m
?????文件????????7040??2018-11-15?10:50??ieee9_2.m
?????文件????????9298??2018-11-15?10:50??ieee9_22.m
?????文件????????1932??2018-11-11?21:36??ieee9_branch.txt
?????文件????????2080??2018-11-14?21:25??ieee9_node.txt