ieee118节点潮流计算程序及节点数据

包含ieee118节点以及线路详细数据 以及 MATLAB潮流计算程序

%%      运用牛顿-拉夫逊法解IEEE118 节点潮流

%%%%%%   首先把原始文件分割成两个独立的纯数据文件分别保存为Voltage.txt和Resis.txt
%%%%%%   计算n节点时，当Voltage.txt和Resis.txt准备完毕后，只需将14改为n，28改为2*n，20改为线路条数即可
%%  读取Voltage中的数据
clear 
load Voltage.txt
Bus=Voltage（:1） ;                 %设定节点数
Voltbase=Voltage（:2） ;            %设定电压标准值
VType=Voltage（:5）;                %设定节点类型
V=Voltage（:6）*100;                    %节点电压标幺值
Angle=Voltage（:7）;                %电压角
Pload=Voltage（:8）;                %负载有功
Qload=Voltage（:9）;                %负载无功
Pgen=Voltage（:10）;                %发电机发出有功
Qgen=Voltage（:11）;                %发电机发出无功
                                   %电压设定点
Vset=Voltage（:13）;                %节点所接并联电容器的电纳
Qgmax=Voltage（:14）;               %Q的最大值
Qgmin=Voltage（:15）;               %Q的最小值
Qsh=Voltage（:17）*100;                 %电容
baseMVA = 100;                     %功率标准值
%%  读取Resis中的数据
load Resis.txt
II=Resis（:1）;                  
JJ=Resis（:2）;                  %两端点
RType=Resis（:5）;               %线路类型
R=Resis（:6）;                   %两点间电阻
X=Resis（:7）;                   %两点间电抗
Btwo=Resis（:8）/2;              %线路两端电纳
K=Resis（:14）;                  % 变压器变比
Kmax=Resis（:16）;
Kmin=Resis（:17）; 
Kset=Resis（:18）;               %变比要求
tic;
t=clock;
%%   开始求导纳 Y
y0=zeros（118） ;                  %两点间的电阻
yt=zeros（118） ;                  %变压器之间的导
Qsh=i*Qsh;                      %电容电抗标幺值
K=1./K ;
for x=1:179 
    switch RType（x）             %比较线路类型  
        case 2                  %含变压器的
            y0（II（x）JJ（x））=1./（K（x）.*（R（x）+i*X（x））） ;                   
            y0（JJ（x）II（x））=y0（II（x）JJ（x））    ;               
            yt（II（x）JJ（x））=（1-K（x））./（（R（x）+i*X（x））.*K（x）^2）;  %II端的电压 
            yt（JJ（x）II（x））=（K（x）-1）./（（R（x）+i*X（x））.*K（x））;    %JJ端的电压 
        case 1                 %不含变压器的
            y0（II（x）JJ（x））=1/（R（x）+i*X（x））;                 
            y0（JJ（x）II（x））=1/（R（x）+i*X（x））;                  
    end
end
Y=zeros（118）;                  %导纳矩阵
clear x  
for x=1:118
    Y（xx）=sum（y0（x:）+yt（x:））+i*Qsh（x）/baseMVA ;     %求出对角导纳（不含B/2）
end
Y=Y-y0;
clear x
for x=1:179
    Y（II（x）II（x））=Y（II（x）II（x））+i*Btwo（x） ;
    Y（JJ（x）JJ（x））=Y（JJ（x）JJ（x））+i*Btwo（x） ;          %加B
end
clear x                           % Y导纳矩阵结束
%%     
U=zeros（1118）;                            %节点电压
for x=1:118
    U（x）=V（x）*exp（i*Angle（x）*pi/180）  ;   %极坐标电压
end
clear x
e=real（U）;
f=imag（U）;           %定义实部e跟虚部f 瓜分U
G=real（Y）;
B=imag（Y）;           %定义实部G跟虚部B 瓜分Y
D=ones（2361）;       %dQ dP dU^2 的矩阵
Ps=zeros（1118）;     
Qs=zeros（1118）;                 
N=0;
for m=1:118
    Ps（m）=（Pgen（m）-Pload（m））/baseMVA;
    Qs（m）=（Qgen（m）-Qload（m））/baseMVA;
end
clear m                       
while max（abs（D））>0.00001                        % 开始迭代过程求解
    for m=1:118
        switch VType（m）                     %按照节点类型分情况计算
            case 1
                D（2*m-1）=Ps（m）-e（m）*sum（G（m:）.*e-B（m:）.*f）-f（m

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       19411  2018-01-30 17:33  118节点电压数据.xlsx
     文件       27284  2018-01-21 10:20  118节点电阻数据.xlsx
     文件        7271  2018-01-31 09:51  ieee118.m