最优潮流主程序代码matlab

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