用matlab怎么算quasi maximum likelihood参数估计渐进分布的协方差

一个小问题~

最近在做QMLE，在求参数估计渐进分布的协方差时卡住了。协防差矩阵里面有个
\[T^-1\sum{h(\hat{\theta},y_t) \cdot h(\hat{\theta},y_t)}\]
这个不知道怎么用matlab实现这个公式。协方差的具体公式hamilton里有（P143）

\[h(\hat{\theta},y_t) \]是yt的loglikelihood对参数的一阶导数。

我用gradient函数不知道为什么整不出来，也不知道用这个函数对不对。

先谢过大家！