Failure in initial objective function evaluation. FMINCON cannot continue.
이전 댓글 표시
Attaching my code, I am just not able to understand why this is not working. Gives the error, Failure in initial objective function evaluation. FMINCON cannot continue.
function new_hyp = max_likelihood(model)
problem.f = @(x) log_likelihood([x(1:model.nv), model.hyp(end)]',model); % calling the log likelihodd function, we maximize log likelihood the find the
% problem.f = @(x) log_likelihood( model.hyp(end)',model);
A = [];
b = [];
Aeq = [];
beq = [];
% model.hyper_bound = log([0.01, 50]); % bounds of hyperparamter
model.hyper_bound = log([0.35,2.73]); % bounds of hyperparamter
lb= model.hyper_bound(1)*ones(1,model.nv);
ub= model.hyper_bound(2)*ones(1,model.nv);
% lb = [-4.60517018598809,-4.60517018598809];
% ub = [3.91202300542815,3.91202300542815];
x0 = zeros(1,model.nv);%[0,0]; %Starting point
options = optimset('Display', 'off') ;
nonlcon=[];
new_hyp = fmincon(problem.f,x0,A,b,Aeq,beq,lb,ub,nonlcon, options)'; % compuatation of hyperparameter that maximizes the log likelihood
% new_hyp = fminsearch(problem.f,x0)';
end
function [lml] = log_likelihood(hyper_param, model)
K = model.cov_model(hyper_param, model.X(1:model.n,:), model.X(1:model.n,:)) + eye(model.n)*model.noise; % Covariance matrix using squared exponential kernel for the input points
Kchol = chol(K);
ALPHA = (Kchol\(Kchol'\model.f));
lml = model.f'*ALPHA + 2*sum(log(diag(Kchol))); % computation of log likelihood ,
end
댓글 수: 3
Matt J
2022년 2월 11일
Before running fmincon, did you test your objective function at x0 to see if it works? (You should always do so).
Joshua Scott
2022년 11월 14일
Hey Matt,
What is the best way to test the function at x0 in matlab? (my apologies if this is a dumb question.)
Alan Weiss
2022년 11월 15일
Define your parameters hyper_param and model in your workspace. Then call
lml = log_likelihood(hyper_param, model)
Or set up your problem and x0 parameter, and call
problem.f(x0)
Alan Weiss
MATLAB mathematical toolbox documentation
답변 (0개)
카테고리
도움말 센터 및 File Exchange에서 Multibody Modeling에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!