Matlab的fmincon一个问题

uelophim · 发表于 2020-8-3 14:12

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

求教大神：
function f=nsllh(x)
close=xlsread('HW_3_Part2Data','#3','D3

1257');
nvix=xlsread('HW_3_Part2Data','#3','H4:H1257');
alpha=zeros(1,1);beta=zeros(1,1);theta=zeros(1,1);gamma=zeros(1,1);
x=[alpha,beta,theta,gamma];
r=zeros(1254,1);
for i=1:1254
r(i+1,1)=log(close(i+1,1)/close(i,1));
end
sigma=var(r);
omega=sigma*(1-alpha*(1+theta^2)-beta);
npsigma=ones(1254,1);
npsigma(1,1)=sigma;
for i=1:1253
npsigma(i+1,1)=omega+alpha*(r(i,1)-theta*npsigma(i,1)^0.5)^2+beta*npsigma(i,1)+gamma*nvix(i,1)^2/252;
end
alpha=x(1);beta=x(2);theta=x(3);gamma=x(4);
for i=1:1254
nllh(i,1)=0.5*log(2*3.1415926)+0.5*log(npsigma(i,1))+0.5*(r(i,1)^2/npsigma(i,1));
end
f=sum(nllh);

另外nonlcon的code是：
function [c,ceq]=mycon(x)
alpha=zeros(1,1);beta=zeros(1,1);theta=zeros(1,1);
x=[alpha,beta,theta];
c=alpha*(1+theta^2)+beta-1;
ceq=[];

run出来显示：
x0=[0.04;0.5;2;0.07];
Optimization terminated: first-order optimality measure less
than options.TolFun and maximum constraint violation is less
than options.TolCon.
No active inequalities.

x =

0.0400
0.5000
2.0000
0.0700

fval =

-3.6963e+003

exitflag =

   1

output =

   iterations: 1
      funcCount: 10
      stepsize: 1
      algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
firstorderopt: 0
   cgiterations: []
      message: [1x144 char]

lamda =

      lower: [4x1 double]
      upper: [4x1 double]
      eqlin: [1x0 double]
   eqnonlin: [1x0 double]
   ineqlin: [1x0 double]
ineqnonlin: 0

为什么会这样呢？fval已经改变了但是x不会变？而且iteration只有1

大小的小 · 发表于 2020-8-3 14:34

帮你顶一下

shuddkk · 发表于 2020-8-5 13:17

1、没有你的xls数据文件，没有进行测试；
2、你在哪儿使用了fmincon，没加上啊；
3、调整初值，再测试；
4、设置fmincon的options进行测试；

小白的白 · 发表于 2020-8-6 13:21

来学习一下

uelophim · 发表于 2020-8-6 13:25

shuddkk 发表于 2020-8-5 13:17
9 O7 Y9 s6 f3 a2 I1、没有你的xls数据文件，没有进行测试；
% W' s8 D. c/ e2 h, ~- ~/ h( k# h2、你在哪儿使用了fmincon，没加上啊；8 H9 r$ [& K9 O/ E2 m. k6 n' C
3、调整初值，再测试 ...

好的，我先去试试，谢谢

帐号		自动登录	找回密码
密码			注册

Matlab的fmincon一个问题

EDA365欢迎您登录！

点评

浏览过的版块

推荐内容 /1