%function?[YUPQdeltaSijaSSijSjisumdeltaS]=NRchaoliu（）
function?[Ua]=NRchaoliu（）
%電力系統潮流計算程序；
%輸出：U——節點電壓，P--節點有功，Q--節點無功，deltaSij--支路功率損耗，
%Sij--從節點i流向節點j的功率S--節點復功率sumdeltaS--網絡總損耗
%輸入參數：point為節點信息矩陣，zhilu為支路信息矩陣；
[x]=xlsread（‘NRinput.xls‘‘A2:A2‘）;%從exel中讀取節點數x
[y]=xlsread（‘NRinput.xls‘‘B2:B2‘）;%從exel中讀取支路數y
e=xlsread（‘NRinput.xls‘‘B4:B4‘）;%誤差要求
[point]=xlsread（‘NRinput.xls‘‘C3:H50‘）;%從exel中讀取節點信息矩陣?zys
[zhilu]=xlsread（‘NRinput.xls‘‘J3:Q50‘）;%從exel中讀取支路信息矩陣
sizeHL=size（zhilu）;
zhilu0=zeros（sizeHL（1）2）;%增加2列存儲首末新節點編號zys
zhilu=[zhilu?zhilu0];%zys
%-------節點優化編號定義start??zys-------------------
nSW?=?0;???????????????????%?nSW為平衡節點個數?
nPV?=?0;???????????????????%?nPV為PV節點個數?
nPQ?=?0;???????????????????%?nPQ為PQ節點個數?
for?lx=1:x
????type0=point（lx2）;
????if?type0?==?3??????
????????nSW?=?nSW?+?1;?????
????????SW（nSW:）=point（lx:）;%?計算并儲存平衡節點???
????elseif?type0?==?2???????
????????nPV?=?nPV?+1;??????
????????PV（nPV:）=point（lx:）;%?計算并儲存PV節點???
????else???????
????????nPQ?=?nPQ?+?1;?????
????????PQ（nPQ:）=point（lx:）;%?計算并儲存PQ節點???
????end?
end
point=[PQ;PV;SW];????????????%?對bus矩陣按PQ、PV、平衡節點的順序重新排序
nodenum=[[1:x]‘?point（:1）];%?生成新舊節點對照表
for?hh=1:y?????
????for?ll=1:2?????????
????????for?k=1:x?????????????
????????????if?zhilu（hhll）==nodenum（k2）
????????????????zhilu（hh8+ll）=nodenum（k1）;????????????????
????????????????break?????????????
????????????end?????????
????????end?????
????end?
end????
%-------節點優化編號定義end?zys-------------------


TYPE=zeros（x1）;%TYPE為節點類型矩陣
U=zeros（x1）;%U為節點電壓矩陣
a=zeros（x1）;%a為節點電壓相角矩陣
P=zeros（x1）;%P為節點有功功率
Q=zeros（x1）;%Q為節點無功功率
S=ones（x1）+j*ones（x1）;%S為節點復功率
I=zeros（y1）;%I為起始節點編號矩陣
J=zeros（y1）;%J為終止節點編號矩陣
Rij=zeros（y1）;%R為線路電阻
Xij=zeros（y1）;%X為線路電抗
Zij=Rij+j*Xij;??%Yij為線路阻抗
Y=zeros（x）;??????%Y為n階節點導納方陣
G=zeros（x）;??????%G為n階節點電導方陣
B=zeros（x）;??????%B為n階節點電納方陣
B0=zeros（y1）;?%B0為n*1階線路對地電納值
RT=zeros（y1）;%RT為ij支路y（?矩陣zhilu的行數）*1階變壓器電阻
XT=zeros（y1）;%XT為ij支路y*1階變壓器電抗
ZT=RT+j*XT;%求變壓器阻抗
KT=zeros（y1）;?%K為ij支路y*1階變壓器變比，若k=0表示無變壓器K=1則為標準變比，k?不等于1為非標準變比
%------------------------------矩陣賦初值：

%TYPE=point（:1）;%將point矩陣的第一列賦給TYPE，以下類似
TYPE=point（:2）;%將point矩陣的第一列賦給TYPE，以下類似
%U=point（:2）;a=point（:3）;P=point（:4）;Q=point（:5）;
U=point（:3）;
a=point（:4）;
P=point（:5）;
Q=point（:6）;
%I=zhilu（:1）;J=zhilu（:2）;
I=zhilu（:9）;
J=zhilu（:10）;
Rij=zhilu（:3）;
Xij=zhilu（:4）;
Zij=Rij+j*Xij;
B0=zhilu（:5）;
RT=zhilu（:6）;
XT=zhilu（:7）;
ZT=RT+j*XT;
KT=zhilu（:8）;
tic??%程序運行時間開始計時
%------------------------------求節點導納矩陣Y
for?m=1:y??????%求Y中非對角元元素Yij
????if?KT（m）==0%若無變壓器則Yij直接為線路阻抗分之一取負值．
????????Y（I（m）J（m））=-1/Zij（m）;
????????Y（J（m）I（m））=-1/Zij（m）;
????else?%有變壓器時，Ｙij為線路阻抗乘以ＫＴ后分之一再取負值
??????????Y（I（m）J（m））=-1/（KT（m）*ZT（m））;
??????????Y（J（m）I（m））=-1/（KT（m）*ZT（m））;
???end
end
for?m=1:x????%求Y中的Ｙii
????for?n=1:y
???????if?KT（n）==0%無變壓器時Ｙii為Ｙij加上線路對地電導乘j
?????????if（I（n）==m|J（n）==m）
????????????Y（mm）=Y（

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????15553??2017-09-19?14:22??牛-拉法潮流\NRchaoliu.m
?????文件???????32256??2017-09-19?14:34??牛-拉法潮流\NRinput.xls
?????文件???????21671??2017-09-19?15:00??牛-拉法潮流\outputNR.txt
?????目錄???????????0??2019-03-04?20:43??牛-拉法潮流\