电力系统IEEE标准测试系统

本资源包含十多个不同节点的IEEE标准测试系统。每个系统详细的给出了系统结构与系统数据。可供电力系统的研究人员使用。

%data支路信息：
%第一列存储存储支路序号
%第二列存储支路得一个端点
%第三例存储支路得另一个端点
%第四利存储支路电阻
%第五列存储支路电抗
%第六列存储支路对地电导
%第七列存储支路对地电纳
%第八列存储变压器变比

%B2pie节点信息：
%第一列为节点编号
%第二列为发电机有功功率
%第三列为发电机无功功率
%第四列为负荷有功功率
%第五列为负荷无功功率
%第六列为节点电压初始值
%第七列为节点电压相角初始值
%第八例为该点得类型：1为PQ节点，2为PV节点，3为平衡节点负荷节点排在前面，PV节点排在后面，平衡节点最后
%B3=input（‘\n请输入对地导纳信息：B3=‘）;
%B4=input（‘\n请输入阻抗信息：B4=‘）;
%k为迭代次数；
clc
clear
% close all
 [dataY1]=xlsread（‘G:\电力系统课程潮流资料\inputdata2.xls‘）;
% Filepath=‘e:\PQmethoddata\inputdata.txt‘;
% data=load（Filepath）;
nbus=data（11）;              %节点数
nline=data（12）;          %支路数
ngen=data（13）;              %发电机节点数
nload=data（14）;               % 负荷节点数
basemva=data（15）;          % 容量基准值 
jingdu=data（16）;%计算精度
detaY=ones（nbus+nload-11）;
B1pie=data（2:nline+11:8）
for i=1:nline
    B1（i1:3）=B1pie（i1:3）;
    B1（i4）=B1pie（i4）+B1pie（i5）*sqrt（-1）;
    B1（i5）=B1pie（i6）+B1pie（i7）*sqrt（-1）;
    B1（i6）=B1pie（i8）;
end
Y=zeros（nbus）;
for a=1:nline
    p=B1（a2）;
    q=B1（a3）;
    Y（pp）=Y（pp）+1./B1（a4）+B1（a5）./2;
    Y（pq）=Y（pq）-1./（B1（a4）*B1（a6））;
    Y（qp）=Y（pq）;
    Y（qq）=Y（qq）+1./（B1（a4）*B1（a6）^2）+B1（a5）./2;
end
Y         %节点导纳阵中线路对地导纳是正吗？
G=real（Y）;
B=imag（Y）;
B2pie=data（nline+2:nline+1+nbus1:8）;
for i=1:nbus
B2（i1）=（B2pie（i2）-B2pie（i4）+（B2pie（i3）-B2pie（i5））*sqrt（-1））/100;
B2（i2）=B2pie（i6）*cos（B2pie（i7））+（B2pie（i6）*sin（B2pie（i7）））*sqrt（-1）;
B2（i3）=B2pie（i8）;
end
for i=1:nbus
    V（i）=abs（B2（i2））;       %各节点电压幅值
    phase（i）=angle（B2（i2））; %各节点电压相角
    P0（i）=real（B2（i1））    %各节点的注入有功功率
    Q0（i）=imag（B2（i1））    %各节点的注入无功功率
end
%根据初始值计算各节点注入功率及功率不平衡量
ndiedai=0;
while max（abs（detaY））>jingdu
    ndiedai=ndiedai+1;
 Q=zeros（nbus1）;
 P=zeros（nbus1）;
 for i=1:nbus
   for  j=1:nbus
   P（i）=V（i）*V（j）*（G（ij）*cos（phase（i）-phase（j））+B（ij）*sin（phase（i）-phase（j）））+P（i）;
   Q（i）=V（i）*V（j）*（G（ij）*sin（phase（i）-phase（j））-B（ij）*cos（phase（i）-phase（j）））+ Q（i）;
   end
 end
for i=1:nbus
     
    switch B2（i3）
        case 1
           
              detaP（i）=P0（i）-P（i）;
              detaQ（i）=Q0（i）-Q（i）;
        case 2
             
             detaP（i）=P0（i）-P（i）;
        otherwise
        pinghengjd=i;
    end
end
detaP=remove（detaPpinghengjd）;% 在1行，1列向量或方阵中去除平衡节点
% 牛顿法的因变量不平衡量组成的向量    
 for i=1:nbus-1   
detaY（i）=detaP（i）;    
 end
 for i=1:（nbus-ngen）
     detaY（i+nbus-1）=detaQ（i）
 end
 
          
%接着生成雅可比矩阵
for i=1:（nbus-1）
    for j=1:（nbus-1）
      if i~=j
          H（ij）=-V（i）*V（j）*（G（ij）*sin（phase（i）-phase（j））-B（ij）*cos（phase（i）-phase（j）））;
      else
          H（ii）=V（i）^2*B（ii）+Q（i）;
      end
    end
end

for i=1:（nbus-1）
    for j=1:（nbus-ngen）
        if i~=j
         N（ij）=-V（i）*（G（ij）*cos（phase（i）-phase（j））+B（ij）*sin（phase（i）-phase（j）））;
        else
         N（ii）=（-V（i）^2*G（ii）-P（i））/V（i）;
        end
    end
end

for i=1:（nbus-ngen）
    for j=1:（nbus-1）
        if i~=j
         J（ij）=V（i）*V（j）*（G（ij）*cos（pha

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     139776  2012-01-01 19:32  IEEE标准测试系统\003_bus\003readme.doc

     文件       2267  2012-01-01 19:32  IEEE标准测试系统\003_bus\bpa\003bpa.DAT

     文件        495  2012-01-01 19:32  IEEE标准测试系统\003_bus\bpa\003bpaswi.DAT

     文件        828  2012-01-01 19:32  IEEE标准测试系统\003_bus\eurostag\003edfech.DAT

     文件       2346  2012-01-01 19:32  IEEE标准测试系统\003_bus\ieee\003ieee.dat

     文件        513  2012-01-01 19:32  IEEE标准测试系统\003_bus\th\003thdat.dat

     文件        221  2012-01-01 19:32  IEEE标准测试系统\003_bus\th\003thmac.dat

     文件        313  2012-01-01 19:32  IEEE标准测试系统\003_bus\th\003thres.dat

     文件       2086  2012-01-01 19:32  IEEE标准测试系统\003_bus.png

     文件      53248  2012-01-01 19:32  IEEE标准测试系统\004_bus\004readme.doc

     文件       2649  2012-01-01 19:32  IEEE标准测试系统\004_bus\bpa\004bpa.dat

     文件      45814  2012-01-01 19:32  IEEE标准测试系统\004_bus\bpa\004bpa.pfo

     文件          1  2012-01-01 19:32  IEEE标准测试系统\004_bus\bpa\PWRFLO.DIS

     文件        927  2012-01-01 19:32  IEEE标准测试系统\004_bus\eurostag\004edfech.dat

     文件       2647  2012-01-01 19:32  IEEE标准测试系统\004_bus\ieee\004ieee.dat

     文件        381  2012-01-01 19:32  IEEE标准测试系统\004_bus\th\004thdat.dat

     文件        360  2012-01-01 19:32  IEEE标准测试系统\004_bus\th\004thres.dat

     文件       3272  2012-01-01 19:32  IEEE标准测试系统\004_bus.png

     文件      74752  2012-01-01 19:32  IEEE标准测试系统\005_bus\005readme.doc

     文件       3471  2012-01-01 19:32  IEEE标准测试系统\005_bus\bpa\005bpa.DAT

     文件       1208  2012-01-01 19:32  IEEE标准测试系统\005_bus\eurostag\005edfech.DAT

     文件       3550  2012-01-01 19:32  IEEE标准测试系统\005_bus\ieee\005ieee.dat

     文件        516  2012-01-01 19:32  IEEE标准测试系统\005_bus\th\005thdat.dat

     文件        410  2012-01-01 19:32  IEEE标准测试系统\005_bus\th\005thres.dat

     文件       4590  2012-01-01 19:32  IEEE标准测试系统\005_bus.png

     文件     143360  2012-01-01 19:32  IEEE标准测试系统\009_bus\009readme.doc

     文件       5907  2012-01-01 19:32  IEEE标准测试系统\009_bus\bpa\009bpa.DAT

     文件        624  2012-01-01 19:32  IEEE标准测试系统\009_bus\bpa\009bpaswi.DAT

     文件       1267  2012-01-01 19:32  IEEE标准测试系统\009_bus\eurostag\009edfdta.DAT

     文件       2051  2012-01-01 19:32  IEEE标准测试系统\009_bus\eurostag\009edfech.DAT

............此处省略241个文件信息