![](https://csdnimg.cn/release/download_crawler_static/35455464/bg1.jpg)
多元 GARCH 模型预测的 Matlab 程序
function [parameters,loglikelihood, Ht, likelihoods, stdresid, stderrors, A, B, scores]
= full_bekk_mvgarch(data,p,q, BEKKoptions);
% PURPOSE:
% To Estimatea full BEKK multivariate GARCH model. %
%
% USAGE:
% [parameters,loglikelihood, Ht, likelihoods, stdresid,stderrors,A, B, scores]
= full_bekk_mvgarch(data,p,q,options);
%
%
% INPUTS:
% data - A t by k matrix of zeromeanresiduals
% p - The lag length of the innovation process
% q - The lag length of the AR process
% options - (optional) Options for the optimization(fminunc)
%
% OUTPUTS:
% parameters - A (k*(k+1))/2+p*k^2+q*k^2 vector of estimated
parameteters.F
% or any k^2 setof Innovation or AR parametersX,
% reshape(X,k,k) will give the correct matrix
% To recoverC, useivech(parmaeters(1:(k*(k+1))/2)
% loglikelihood - The loglikelihood of the function at the optimum
% Ht - A k x k x t 3 dimension matrix of conditional covariances
% likelihoods - A t by 1 vector of individual likelihoods
% stdresid - A t by k matrix of multivariate standardizedresiduals
% stderrors - A numParams^2 square matrix of robust Standad
Errors(A^(-1)*B*A^(-1)*t^(-1))
% A - The estimatedinverseof the non-robust Standarderrors
% B - The estimatedcovarianceof teh scores
% scores - A t by numParamsmatrix of individual scores
% needto try and get somesmartstartgin values