matlab牛顿拉夫逊计算潮流

%initialing the parameters
clear;
%input（‘Please input the parameters of the Line and the Transformers‘）;
Nline=input（‘Please input the number of the bus\nNline=‘）;
Ntrans=input（‘Please input the number of the transformer\nNtrans=‘）;
Line=zeros（Nline5）;
Transformer=zeros（Ntrans6）;
for m=1:Nline
    Line（m:）=input（‘Please input the parameters of the bus\nLine= ‘）;
end
for m=1:Ntrans
    Transformer（m:）=input（‘Please input the parameters of the transformer\nTransformer= ‘）;
end
clc;
PQ0=input（‘Please input the value of the PQ0\nPQ0=‘）;
PQ0=PQ0‘;
EF=input（‘Please input the value of the EF\nEF=‘）;
EF=EF‘;
%initial some variables that will be used later
Nnode=Nline+Ntrans;
k=0;
error=1;
Sslack=0;
%form the Matrix of Y
Y=zeros（Nnode）;
if Nline>=1               %compute by adding the nodes of the lines
    for k=1:Nline
        t1=Line（k2）;
        t2=Line（k3）;
        Y1=1/（Line（k4）+j*Line（k5））;
        Y（t1t1）=Y（t1t1）+Y1;
        Y（t1t2）=Y（t1t2）-Y1