`
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无法继续。“
答案 0 :(得分:0)
“使用garchlike太多输入参数时出错。”是你的问题。
fmincon
正在调用garchlike([psi;gamma;lambda;A],data, p , q)
,但您的garchlike
只需要1个参数。
专业提示:了解如何使用调试器。首先设置dbstop if error
。
编辑:好的,这里有更多细节。请仔细阅读。
你的garchlike
接受1个参数作为输入,而你的fmincon用4个参数调用它。
要么(1)你需要从fmincon中删除那3个额外的参数(即数据,p,q),或者(2)使garchlike再获取3个参数!!