最优潮流内点法

在matlab上面编写的，基于原对偶内点法的（跟踪中心路径参数内点法）最优潮流，具体编程思想请参见王锡凡《现代电力系统分析》

function [A dh_dx dg_dx H_ d2h_dx_y d2g_dx_c d2f_dx temp L_Z U_W] = Coeff（num_node num_PQ Y B X A2）
[num_branch ~]  = size（B）;
len_x = 2*（num_node-num_PQ+num_node）; % 14
num_gen = num_node-num_PQ; % 2
num_equa = 2*num_node;
num_inequa = 2*num_gen+num_node+num_branch;

z = X（1:num_inequa）; l = X（num_inequa+1:2*num_inequa）; w = X（2*num_inequa+1:3*num_inequa）;
u = X（3*num_inequa+1:4*num_inequa）; x = X（4*num_inequa+1:4*num_inequa+len_x）; y = X（4*num_inequa+len_x+1:4*num_inequa+len_x+num_equa）;
Xtilde = x（2*num_gen+1:len_x）;
%% 计算等式约束Jacobian矩阵
dh_dPG = zeros（22*num_node）;
dh_dQR = zeros（22*num_node）;
for k = num_PQ+1:num_node
   dh_dPG（k-num_PQk*2-1） = 1;
   dh_dQR（k-num_PQk*2） = 1;
end

dh_dXtilde = zeros（2*num_node2*num_node）;
for ii = 1:num_node
    f

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

     文件       8814  2017-10-16 09:20  内点法\Coeff.m

     文件       1207  2017-10-15 21:54  内点法\dX.m

     文件       8554  2019-02-28 15:59  内点法\main.m

     文件        855  2017-10-07 10:56  内点法\makeY.m

     文件       1097  2017-10-27 14:41  内点法\solution.txt

     目录          0  2019-02-28 15:56  内点法

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

                20527                    6