-1

`

function [LLF, LL, H  ] = garchlike(data)


p=1;
q=1;


 if isempty(q)
    m=p;
 else
     m  =  max(p,q);   
 end

%this are the parameters
psi=[0.01;0.01];
 gamma=[0;0.02;0.02;0];
 lambda=[0;0.03;0.03;0];
  A=[1;-0.04;-0.03;1];

 stdestimate = std(data,1);                      
data   =  [stdestimate ; data];
T =  size(data,1);

B=(reshape(A,2,2))^-1;
%I am squaring each element of the inverse of A
B2=((reshape(A,2,2))^-1).^2;
%getting the Bl matrix(squared interactions)

Bl=[B2(1,:);B(1,:).*B(2,:);B2(2,:)];

garchc=Bl*psi;
garcha=Bl*(reshape(gamma,2,2))*reshape(A.^2,2,2);
garchb= Bl*(reshape(lambda,2,2))*reshape(A.^2,2,2);

H(1,1).day=[0.001; 0; 0; 0.002; ];
%convert to matrix

  H(1,1).day=ivech(H(1,1).day);

for t = (m + 1):T


  H(1,t).day= ivech(garchc + garcha*(diag(H(1,t-1).day))+ garchb*(data((t-1),:).^2)');

end
% Getting the loglikelihood
LLF=zeros(1,T);
%the loklikelihood
for t = (m + 2):T

  LLF(t)= log(det(H(1,t).day))+ data(t,:)/(H(1,t).day)*data(t,:)';

end
t=(m+2):T;
LLF=-LLF(t);
LL=sum(LLF);
`

我有这个功能并尝试使用 fmincon 最小化：

`  [parameters,LLF,EXITFLAG, OUTPUT,HESSIAN] =  fmincon('garchlike',    [psi;gamma;lambda;A],sumA ,sumB ,[] , [] , LB , UB,[],options, data, p , q);`

我不知道为什么不工作并给我一个错误：“错误使用 garchlike 错误使用 garchlike 输入参数太多。fmincon 中的错误（第 631 行）initVals.f = feval(funfcn{3},X,varargin{: }); 原因：初始用户提供的目标函数评估失败。FMINCON 无法继续。”

4

1 回答 1

0

“错误使用 garchlike 输入参数太多。” 是你的问题。

fmincon正在调用garchlike([psi;gamma;lambda;A],data, p , q)，但您garchlike只需要 1 个参数。

专业提示：学习如何使用调试器。从设置开始dbstop if error。

编辑：好的，这里有更多细节。请仔细阅读。

您garchlike将 1 个参数作为输入，而您的 fmincon 使用 4 个参数调用它。要么（1）你需要从 fmincon 中删除这 3 个额外的参数（即 data、p、q），要么（2）让 garchlike 多带 3 个参数！

于 2013-10-09T15:07:16.730 回答