function?[PlossTotal0VmIL]=TideCal（Z）?????%潮流計算輸出系統總損耗值和各節點電壓幅值
%PlossTotal??為系統總損耗
%Vm為各節點電壓幅值
[Z]=OpDF_;
[ml]=size（Z）;
b=m;
n=m+1;
Vb=11;?????????????????%Vb=11kV?跟系統電壓等級有關
Sb=100000;?????????????%Sb=100MVA=100000kVA
Zb=（Vb^2/Sb）*1000;?????%Zb單位為Ω????電流單位為A
??%Z=[支路數?首節點?尾節點?尾節點類型?單位長度的支路阻抗?支路長度?尾節點的有功功率?尾節點的無功功率]???電源點：0??末梢點：1??T接點：2???負荷點：3
h=0;
Vchu=ones（n1）;
Vnode=Vchu;
Iin=zeros（b1）;%負荷電流（節點注入電流）
IL=zeros（b1）;???%支路電流
Sloss=zeros（b1）;%線路損耗
S=zeros（b1）;
Vold=zeros（n1）;%記錄上次的電壓幅值
QC=zeros（b2）;???%電容器
kk=100;
while??（kk>0.000000001）
????%求節點注入電流
????for?L=1:b
????????k=Z（L3）;
????????Ua=Vnode（k+11）;
????????S（k1）=Z（k7）+i*Z（k8）;?%對應尾節點有功和無功????????????%負荷的視在功率
????????Iin（k1）=conj（S（k1）/Ua）-（QC（k1）/Sb）/Ua-（QC（k2）/Sb）/Ua;%%%
????end
????Ii