作者jajagood (蓝海白沙)
站内MATLAB
标题[问题]有关fmincon的初始值设定
时间Tue Feb 25 01:44:07 2014
最近在写有关fmincon的语法，遇到初始值的问题
以下是方程式
rpha=[0.01];
betaa=[0.05];
pAQL=[100];
pLQL=[2000];
initial_n=10;
x0=[initial_n 1 1.5];
options=optimset('Display','iter','Algorithm','sqp');
[x,f_xStr,exitflag]=fmincon(@QSSobjfuntion,x0,[],[],[],[],[],[],@QSSfuntion,options);
可以解出x=[53 0.9594 1.1856]
如果解不出来，将初始值initial_n一直往上叠加也可以找到答案
但是在某些数据标准情况下
EX
rpha=[0.01];
betaa=[0.05];
pAQL=[1000];
pLQL=[2000];
initial_n=400;
x0=[initial_n 1 1.5];
options=optimset('Display','iter','Algorithm','sqp');
[x,f_xStr,exitflag]=fmincon(@QSSobjfuntion,x0,[],[],[],[],[],[],@QSSfuntion,options)
已知答案应该是x[719 0.9594 1.0155]
但是初始值已经叠加到400能找不出此解，跑出以下问题，我想与初始值这块有关
但是不太懂这意思。
Warning: Infinite or Not-a-Number value encountered.
> In quadgk>midpArea at 319
In quadgk at 216
In QSSfuntion at 18
In fmincon at 722
In QSS at 21
Error using sqpLineSearch (line 25)
Nonlinear constraint function is undefined at initial point. Fmincon cannot
continue.
Error in fmincon (line 910)
[X,FVAL,EXITFLAG,OUTPUT,LAMBDA,GRAD,HESSIAN] =
sqpLineSearch(funfcn,X,full(A),full(B),full(Aeq),full(Beq), ...
Error in QSS (line 21)
[x,
f_xStr,exitflag]=fmincon(@QSSobjfuntion,x0,[],[],[],[],[],[],@QSSfuntion,options)
请版上的高手大大帮我找出问题所在。谢谢您!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.43.244.176
1F：→ LiamIssac:function长怎样? convex? 02/25 02:26
※ 编辑: jajagood 来自: 114.43.244.176 (02/25 02:56)
2F：推 letibe:3变数而已 先画图看看吧 03/02 16:37
3F：→ jajagood:可是数值解不出来无法画图吧? 我不知道PAQL的值对解 03/04 17:43
4F：→ jajagood:有甚麽样的问题 03/04 17:43