algorithm - Matlab 中的障碍函数

Question

我在哪里可以找到 Matlab 中障碍函数的实现？我试图查看算法 interior-point 是如何实现的，这是我在 fmincon.m

的末尾找到的

elseif strcmpi(OUTPUT.algorithm,interiorPoint)
    defaultopt.MaxIter = 1000; defaultopt.MaxFunEvals = 3000; defaultopt.TolX = 1e-10;
    defaultopt.Hessian = 'bfgs';
    mEq = lin_eq + sizes.mNonlinEq + nnz(xIndices.fixed); % number of equalities
    % Interior-point-specific options. Default values for lbfgs memory is 10, and 
    % ldl pivot threshold is 0.01
    options = getIpOptions(options,sizes.nVar,mEq,flags.constr,defaultopt,10,0.01); 

    [X,FVAL,EXITFLAG,OUTPUT,LAMBDA,GRAD,HESSIAN] = barrier(funfcn,X,A,B,Aeq,Beq,l,u,confcn,options.HessFcn, ...
        initVals.f,initVals.g,initVals.ncineq,initVals.nceq,initVals.gnc,initVals.gnceq,HESSIAN, ...
        xIndices,options,optionFeedback,finDiffFlags,varargin{:});

所以我想看看障碍中有什么但失败了。

edit barrier.m

我得到:

Answer 1

障碍函数在 p-file 中定义(精确位于 MATLABROOT/toolbox/optim/optim/barrier.p)。

不幸的是，p 文件的关键在于它们被混淆了，即您无法阅读源代码。这是一个关于 SO 的反复出现的问题，请参阅 this thread例如。

恐怕您无法阅读 barrier 中的内容。也许如果您好心询问 Mathworks，他们可以为您提供有关内容的一些信息。

最佳