Minimize vector using Fmincon
조회 수: 21 (최근 30일)
이전 댓글 표시
I want to minimize a cost function using Fmincon. My goal is to minimize "obj" such that my model has the best fit to the experimental data.
C = importdata('Experiment6.txt')
S = C(:,1); % Column 84x1; experimental data
L = C(:,2); % Column 84x1; experimental data
F = C(:,3); % Column 84x1; experimental data
MEdp = 0; % Fixed
MEp = 196; % Fixed
E = ones(1,11); % Initial guess (Xo)
t = logspace(-4, 1,11); % Fixed vector
Eoptimized = fmincon(@objective, E, [],[])
function obj = objective(E, MEdp, t, MEp, F, L, S)
% Calculate model, returns 84x1 columns
for i = 1:length(t)
MEp = MEp + ( E(i) * t(i).^2 * F.^2) ./ ( 1 + ( t(i)^2 * F.^2) );
MEdp = MEdp + ( E(i) * t(i)* F) ./ ( 1 + ( t(i)^2 * F.^2) );
end
% Cost function
obj = sum ( ( ( (L - MEdp) -1 ).^2 ) + ( ( (S - MEp) -1 ) .^2 ) );
end
Objection returns a single value. I want to have E returned as an optimized vector. However, when I try to run the the code it states that it does not have enough input arguments.
% Without this code the function works propely:
Eoptimized = fmincon(@objective, E, [],[])
The following error message appears:
Not enough input arguments.
Error in DiscreteSpectrum>objective(line 264)
for i = 1:length(t)
Error in fmincon (line 552)
initVals.f = feval(funfcn{3},X,varargin{:});
Error in DiscreteSpectrum (line 244)
Eoptimized = fmincon(@objective, E, [], []);
Caused by:
Failure in initial objective function evaluation. FMINCON cannot continue.
댓글 수: 0
채택된 답변
추가 답변 (1개)
Alan Weiss
2020년 6월 8일
I think that you called the objective function incorrectly. Try
Eoptimized = fmincon(@(E)objective(E, MEdp, t, MEp, F, L, S), E)
That said, I think you would do better to use the correct solver for this problem, namely lsqnonlin. You will find that lsqnonlin is faster and converges better than fmincon for this problem.
Alan Weiss
MATLAB mathematical toolbox documentation
참고 항목
카테고리
Help Center 및 File Exchange에서 Nonlinear Optimization에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!