program.rar

通过matlab计算综合能源系统（IEEE39节点天然气系统和比利时20节点天然气系统），并通过电力系统指标衡量天然气系统对综合能源系统的影响

function [V converged i lam] = cpf_corrector（Ybus Sbusb V_hat ref pv pq ...
                lam_hat Sbust Vprv lamprv z step parameterization mpopt）
%CPF_CORRECTOR  Solves the corrector step of a continuation power flow
%解决连续潮流的校正步长
%   [V CONVERGED I LAM] = CPF_CORRECTOR（YBUS SBUSB V_HAT REF PV PQ ...
%                                       LAM_HAT SBUST VPRV LPRV Z ...
%                                       STEP PARAMETERIZATION MPOPT）
%
%   Computes the corrector step of a continuation power flow using a
%   full Newton method with selected parameterization scheme.
%
%   Inputs:
%       YBUS : complex bus admittance matrix
%       SBUSB : handle of function returning nb x 1 vector of complex
%               base case injections in p.u. and derivatives w.r.t. |V|
%       V_HAT :  predicted complex bus voltage vector
%       REF : vector of indices for REF buses
%       PV : vector of indices of PV buses
%       PQ : vector of indices of PQ buses
%       LAM_HAT : predicted scalar lambda
%       SBUST : handle of function returning nb x 1 vector of complex
%               target case injections in p.u. and derivatives w.r.t. |V|
%       VPRV : complex bus voltage vector at previous solution
%       LAMPRV : scalar lambda value at previous solution
%       STEP : continuation step length
%       Z : normalized tangent prediction vector
%       STEP : continuation step size
%       PARAMETERIZATION : Value of cpf.parameterization option.
%       MPOPT : Options struct
%
%   Outputs:
%       V : complex bus voltage solution vector
%       CONVERGED : Newton iteration count
%       I : Newton iteration count
%       LAM : lambda continuation parameter
%
%   See also RUNCPF.

%   MATPOWER
%   Copyright （c） 1996-2016 Power Systems Engineering Research Center （PSERC）
%   by Ray Zimmerman PSERC Cornell
%   Shrirang Abhyankar Argonne National Laboratory
%   and Alexander Flueck IIT
%
%   Modified by Alexander J. Flueck Illinois Institute of Technology
%   2001.02.22 - corrector.m （ver 1.0） based on newtonpf.m （MATPOWER 2.0）
%
%   Modified by Shrirang Abhyankar Argonne National Laboratory
%   （Updated to be compatible with MATPOWER version 4.1）
%
%   This file is part of MATPOWER.
%   Covered by the 3-clause BSD License （see LICENSE file for details）.
%   See https://matpower.org for more info.

%% default arguments
if nargin < 14
    mpopt = mpoption;
end

%% options
tol     = mpopt.pf.tol;
max_it  = mpopt.pf.nr.max_it;

%% initialize
converged = 0;
i = 0;
V = V_hat;          %% initialize V with predicted V
Va = angle（V）;
Vm = abs（V）;
lam = lam_hat;      %% initialize lam with predicted lam

%% set up indexing for updating V
npv = length（pv）;
npq = length（pq）;
nb = length（V）;         %% number of buses
j1 = 1;         j2 = npv;           %% j1:j2 - V angle of pv buses
j3 = j2 + 1;    j4 = j2 + npq;      %% j3:j4 - V angle of pq buses
j5 = j4 + 1;    j6 = j4 + npq;      %% j5:j6 - V mag of pq buses
j7 = j6 + 1;   

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

    .CA....      5805  2020-03-19 18:49  program\01  cpf\cpf_corrector.m

    .CA....      1439  2020-03-19 18:53  program\01  cpf\cpf_predictor.m

    .CA....       509  2020-03-06 21:15  program\01  cpf\exemple.m

    .CA....        84  2020-03-05 18:32  program\01  cpf\foring.m

    .CA....      4597  2019-06-21 09:41  program\02  case\case14.m

    .CA....     12797  2020-04-08 16:47  program\02  case\case20.m

    .CA....      9628  2019-06-21 09:41  program\02  case\case39.m

    .CA....      2235  2019-06-21 09:41  program\02  case\case9.m

    .CA....      1909  2019-06-21 09:41  program\02  case\case9Q.m

    .CA....      1884  2019-06-21 09:41  program\02  case\case9target.m

    .CA....      3368  2020-04-30 10:59  program\03  IE\case20.m

    .CA....      9628  2019-06-21 09:41  program\03  IE\case39.m

    .CA....      2235  2019-06-21 09:41  program\03  IE\case9.m

    .CA....      2620  2020-05-07 16:04  program\03  IE\IEcf.m

    .CA....      1607  2020-05-07 11:10  program\03  IE\IEf.m

    .CA....      1202  2020-04-20 19:02  program\03  IE\IEof.m

    .CA....       912  2020-04-19 22:38  program\03  IE\Integrated.m

    .CA....      3056  2020-05-05 15:10  program\03  IE\rungf.m

    .CA....      1552  2020-05-07 16:26  program\04  Quota\quota.m

    .CA....       944  2020-05-08 10:30  program\04  Quota\Quota_mar.m

    .CA....       801  2020-05-08 10:33  program\04  Quota\Quota_vol.m

    .CA....       830  2020-05-08 10:56  program\04  Quota\Quota_vrp.m

    .C.D...         0  2020-07-23 22:37  program\01  cpf

    .C.D...         0  2020-07-23 22:37  program\02  case

    .C.D...         0  2020-07-23 22:38  program\03  IE

    .C.D...         0  2020-07-23 22:38  program\04  Quota

    .C.D...         0  2020-07-23 22:38  program

----------- ---------  ---------- -----  ----

                69642                    27



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