bekk matlab,有木有高手懂Bekk-garch模型的malab代码的？可以给注释一下不？

>> function [parameters, loglikelihood, Ht, likelihoods, stdresid, stderrors, A, B, scores] = diagonal_bekk_T_mvgarch(data,p,q,BEKKoptions);

% PURPOSE:

%      To Estimate a diagonal BEKK multivariate GARCH model with T-dist errors.  ****SEE WARNING AT END OF HELP FILE****

%

% USAGE:

%      [parameters, loglikelihood, Ht, likelihoods, stdresid, stderrors, A, B, scores]  = diagonal_bekk_T_mvgarch(data,p,q,options);

%

% INPUTS:

%      data          - A t by k matrix of zero mean residuals

%      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+q*k+1 vector of estimated parameteters.

%