⛄一、内点法简介
法。而且无论是面对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=