资源简介

基于IEEE33节点所做的matlab潮流计算,具有通用性,每一行代码都带注释,推荐下载。

资源截图

代码片段和文件信息

%*********************%潮流计算*****************************************
function flow()
global X Y Z  Ploss Num V DG
format long
Branch_Data=dlmread(‘branchdata.txt‘);   %读取线路参数矩阵
Node_Data=dlmread(‘nodedata.txt‘);       %读取节点参数矩阵
UB = 12.66;                             % 电压基准 kV
SB = 10;                                % 功率基准 MVA
ZB = UB^2/SB;                           % 阻抗基准 ohm
Branch_Data(:[34]) = Branch_Data(:[34]) / ZB;     % 阻抗标幺值
Node_Data(:[23]) = Node_Data(:[23]) / SB / 1000;  % 功率标幺值
for j=1:DG
Node_Data(X(j)-12) = Node_Data(X(j)-12) - Y(j);        %节点编号从1开始
Node_Data(X(j)-13) = Node_Data(X(j)-13) - Z(j);
end
A = zeros(Num);
for n = 1:Num-1
    A(Branch_Data(n1)Branch_Data(n2)) = 1;
end                                     % 形成邻接矩阵A
AT = A‘;                                        % 形成A的转置矩阵
S = [0;Node_Data(:2) + i*Node_Data(:3)];           % 形成节点负荷矩阵,根节点赋值0
ZL = [0;Branch_Data(:3) + i*Branch_Data(:4)];       % 形成线路阻抗矩阵
 
V = ones(Num1);                         % 各个节点电压赋初值
Delta = 1;                              % 收敛判据赋初值
TempV = V;         

评论

共有 条评论