【优化电价】内点法求解实时电价最优问题【含Matlab源码 1161期】

⛄一、内点法简介

法。而且无论是面对LP还是QP，内点法都显示出了相当的极好的性能，例如多项式的算法复杂度。

⛄二、部分源代码

clear;
clc;
errArr=[];
%%
%初始化！！！
initial;
% Start clock
t1 = clock;
%%
ROU=sl’MU_MIN+su’MU_MAX;
MUt=SIGMAROU/(2length(sl));%初始对偶因子与惩罚因子计算%
ik=0;%计迭代次数！！！
%迭代循环过程！！
while(abs(ROU)>=err)
%%
%Calcute h,g matrix
ROU=sl’MU_MIN+su’MU_MAX;
errArr=[errArr;ROU;];
SIGMA=0;
MU=SIGMAROU/(2length(sl)); %中心参数置零%
for i=1:30
temp=0;
for j=1:30
temp=temp-V(j)*aY(i,j)*cos(Vth(i)-Vth(j)-Yth(i,j));
end
if (i>6)
tPg=0;
else
tPg=