%*********************%潮流計算*****************************************
function?flow（）
global?X?Y?Z??Ploss?Num?V?DG
format?long
Branch_Data=dlmread（‘branchdata.txt‘）;???%讀取線路參數矩陣
Node_Data=dlmread（‘nodedata.txt‘）;???????%讀取節點參數矩陣
UB?=?12.66;?????????????????????????????%?電壓基準?kV
SB?=?10;????????????????????????????????%?功率基準?MVA
ZB?=?UB^2/SB;???????????????????????????%?阻抗基準?ohm
Branch_Data（:[34]）?=?Branch_Data（:[34]）?/?ZB;?????%?阻抗標幺值
Node_Data（:[23]）?=?Node_Data（:[23]）?/?SB?/?1000;??%?功率標幺值
for?j=1:DG
Node_Data（X（j）-12）?=?Node_Data（X（j）-12）?-?Y（j）;????????%節點編號從1開始
Node_Data（X（j）-13）?=?Node_Data（X（j）-13）?-?Z（j）;
end
A?=?zeros（Num）;
for?n?=?1:Num-1
????A（Branch_Data（n1）Branch_Data（n2））?=?1;
end?????????????????????????????????????%?形成鄰接矩陣A
AT?=?A‘;????????????????????????????????????????%?形成A的轉置矩陣
S?=?[0;Node_Data（:2）?+?i*Node_Data（:3）];???????????%?形成節點負荷矩陣，根節點賦值0
ZL?=?[0;Branch_Data（:3）?+?i*Branch_Data（:4）];???????%?形成線路阻抗矩陣
?
V?=?ones（Num1）;?????????????????????????%?各個節點電壓賦初值
Delta?=?1;??????????????????????????????%?收斂判據賦初值
TempV?=?V;?????????