如何用程序修改matpower中算例参数并保存

[33nodes.rar] - 配电网潮流计算源程序，可直接运行的，修改输入数据即可 [hhu_pfcp.rar] - 基于MATPOWER修改而成的潮流计算软件，可以采用直流潮流法、牛顿拉夫逊法、PQ分解法与高斯－塞代尔法解算IEEE9-300节点的潮流。 [33node.rar] - 33节点的配电网潮流计算程序，程序自带33节点数据 [model14.rar] - 这是一个基于matlab应用于电力潮流计算的ieee14节点模型 [disflowtry.rar] - 配电网潮流计算，以IEEE33节点为例,运用前退回代法计算 [IEEE33jdmlabbc.rar] - 本程序是利用牛顿法编制的潮流计算程序员源码，33节点的数据，符合ieee33节点的要求。 [powerflow33.rar] - 33节点配电网前退回代法潮流计算程序。可以计算功率损耗，是至今为止比较完善的程序 [fzjd.rar] - 改程序模拟了，对前推回代方法计算配电网潮流计算，迭代次数少 [NRflowCal.rar] - 本程序用牛顿-拉夫逊法计算电力系统潮流，在本程序中支路中含有变压器可以生成导纳矩阵，程序中的算例使用了《电力系统分析》(下册 华中科技大学) P66页的数据，已通过运行，数据正确。 [10.rar] - 极坐标潮流程序，附有4节点，33节点，1047节点的数据，已编译通过。

% 配电网前推回代潮流计算程序 

% 使用IEEE 33节点配电系统作为算例，可实现弱环网情况下的潮流计算

countnum=0

BranchData = [

    1   2   0.0922  0.0470

    2   3   0.4930  0.2511

    3   4   0.3660  0.1864

    4   5   0.3811  0.1941

    5   6   0.8190  0.7070

    6   7   0.1872  0.6188

    7   8   0.7114  0.2351

    8   9   1.0300  0.7400

    9   10  1.0440  0.7400

    10  11  0.1966  0.0650

    11  12  0.3744  0.1238

    12  13  1.4680  1.1550

    13  14  0.5416  0.7129

    14  15  0.5910  0.5260

    15  16  0.7463  0.5450

    16  17  1.2890  1.7210

    17  18  0.7320  0.5740

    2   19  0.1640  0.1565

    19  20  1.5042  1.3554

    20  21  0.4095  0.4784

    21  22  0.7089  0.9373

    3   23  0.4512  0.3083

    23  24  0.8980  0.7091

    24  25  0.8960  0.7011

    6   26  0.2030  0.1034

    26  27  0.2842  0.1447

    27  28  1.0590  0.9337

    28  29  0.8042  0.7006

    29  30  0.5075  0.2585

    30  31  0.9744  0.9630

    31  32  0.3105  0.3619

    32  33  0.3410  0.5302

]                                     % 支路，阻抗

NodeData = [

    2   100.00  60.00

    3   90.00   40.00

    4   120.00  80.00

    5   60.00   30.00

    6   60.00   20.00

    7   200.00  100.00

    8   200.00  100.00

    9   60.00   20.00

    10  60.00   20.00

    11  45.00   30.00

    12  60.00   35.00

    13  60.00   35.00

    14  120.00  80.00

    15  60.00   10.00

    16  60.00   20.00

    17  60.00   20.00

    18  90.00   40.00

    19  90.00   40.00

    20  90.00   40.00

    21  90.00   40.00

    22  90.00   40.00

    23  90.00   50.00

    24  420.00  200.00

    25  420.00  200.00

    26  60.00   25.00

    27  60.00   25.00

    28  60.00   20.00

    29  120.00  70.00

    30  200.00  600.00

    31  150.00  70.00

    32  210.00  100.00

    33  60.00   40.00

]                                      % 节点，负荷

UB = 12.66                             % 电压基准 kV

SB = 10                                % 功率基准 MVA

ZB = UB^2/SB                           % 阻抗基准 ohm

BranchData(:,[3,4]) = BranchData(:,[3,4]) / ZB     % 阻抗标幺化

NodeData(:,[2,3]) = NodeData(:,[2,3]) / SB / 1000% 功率标幺化

NN = 33                                % 节点数

A0 = zeros(NN)

for n = 1:NN-1

    A0(BranchData(n,1),BranchData(n,2)) = 1

end                                     % 形成 A0

AssociatedMatrix=0

for n=2:NN-1

    AssociatedMatrix(n,n)=1

    temp=BranchData(n-1,1)

    AssociatedMatrix(n,1:n-1)=AssociatedMatrix(temp,1:n-1)

end

A0T = A0'                              % 形成 A0 的转置

S = [0-NodeData(:,2) - i*NodeData(:,3)]       % 形成 S

ZL = [0BranchData(:,3) + i*BranchData(:,4)]       % 形成 ZL

V = ones(NN,1)

V(1) = 1                            % 各个节点电压赋初值

IL(NN,1) = -conj(S(NN) / V(NN))          % 最末支路电流赋初值

Delta = 1                              % 收敛判据赋初值

TempV = V                              % 赋初值，用于记忆上次迭代结果

while Delta > 1e-8

    countnum=countnum+1

    IN = conj(S ./ V)                    % 节点注入电流

    

    for n = 1:NN-1

        IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n)

    end                                 % 电流回代过程

    

    for n = 2:NN

        V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n)

    end                                 % 电压前推过程

    

    Delta = max(abs(V-TempV))          % 更新收敛判据

    TempV = V                          % 记忆迭代结果

end

Vangle(:,1)=abs(V)

Vangle(:,2)=angle(V)/3.1415*180

for i=1:NN-1

    st=BranchData(i,1)

    en=BranchData(i,2)

    Sij(i,1)=V(st)*conj((V(st)-V(en))/ZL(i+1))

    Sji(i,1)=V(en)*conj((V(en)-V(st))/ZL(i+1))

end

希望你能用到

针对配电网的特点和常用前推回代潮流算法存在的不足,提出了一种改进的配电网牛顿潮流算法。该算法从基于基尔霍夫电流定律的Ybus算法出发,将迭代得到的结果作为电流注入模型的牛顿法的初值进行潮流分析。算法通过修改雅可比迭代矩阵元素对PV节点进行处理,通过变换系统节点电流方程来解决多平衡节点系统问题,既可用于放射形配网,又能用于含有环网的配电系统。该算法物理意义清晰,容易实现,收敛性好,给出的算例验证了算法的可行性。

---