how to resolve the issue with fmincon

8 vues (au cours des 30 derniers jours)
dav
dav le 11 Mai 2014
Commenté : dav le 11 Mai 2014
Hi,
I need to use the following code in my program. however, it gives an error and the warning. It happens because I use SIZE() when assigning the initial values. but i must use it that way.
Can someone please help me overcome this?
code:
theta0 = ones(size(svhat),1);
Av = ones(size(svhat),1)';
Av(1,1) = 0;
lbb = zeros(size(svhat),1);
lbb(1,1) = 1;
lbb = eps*lbb;
bv = 0.99999;
opts = optimset('Display','off','Algorithm','sqp');
%function
[theta, opt] = fmincon(@(theta)alach(theta,epshat),theta0,Av,bv,[],[],lbb,[],[],opts);
error:
Warning: Input arguments must be scalar.
> In test at 146
Warning: Input arguments must be scalar.
> In test at 148
Warning: Input arguments must be scalar.
> In test at 151
Error using sqpLineSearch (line 22)
Objective function is undefined at initial point. Fmincon cannot continue.
Error in fmincon (line 851)
[X,FVAL,EXITFLAG,OUTPUT,LAMBDA,GRAD,HESSIAN] = sqpLineSearch(funfcn,X,full(A),full(B),full(Aeq),full(Beq), ...
Error in test (line 161)
[theta, opt] = fmincon(@(theta)alach(theta,epshat),theta0,Av,bv,[],[],lbb,[],[],opts);

Connectez-vous pour répondre à cette question.

Réponse acceptée

Matt J
Matt J le 11 Mai 2014
Modifié(e) : Matt J le 11 Mai 2014
I think you might really want
theta0 = ones(numel(svhat),1);
or
theta0=ones(size(svhat));
As for the fmincon error, your function @(theta)alach(theta,epshat) is not returning a finite value at theta0. You should check what it is actually returning.
  6 commentaires
Afficher 4 commentaires plus anciens
dav
dav le 11 Mai 2014
Modifié(e) : dav le 11 Mai 2014
This is what i got. thanks
theta0 =
1 1 1 1 1 1 1
ftheta0 =
NaN
Av_dims =
7 1
Furthermore, Av should be a row vector while lbb,theta0 are column vectors.
dav
dav le 11 Mai 2014
Dear Matt,
The following code worked just fine.
Thank you very much!
code:
theta0 = ones(numel(svhat),1);
Av = ones(numel(svhat),1)';
Av(1,1) = 0;
%Av = Av';
lbb = ones(numel(svhat),1);
lbb(1,1) = 1;
lbb = eps*lbb;
%lbb=lbb';
bv = 0.99999;
opts = optimset('Display','off','Algorithm','sqp');
%function
[theta, opt] = fmincon(@(theta)alach(theta,epshat),theta0,Av,bv,[],[],lbb,[],[],opts);

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Systems of Nonlinear Equations dans Help Center et File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by