如题，我想自己写一个最简单的GARCH（1，1）的程序，包括先产生数据
e(t)=sqrt(h(t))*epsil(t), h(t)=w+alpha*e^2(t-1)+beta*h(t-1)，然后用极大似然估计法去估计参数（w,alpha,beta）。我用的是newton法，提供gradient,和hessian,然后用迭代去找负的对数似然函数的最小值。但是郁闷的是，不管是用matlab自己提供的优化程序（fminunc,fmincon）,还是自己写的，都不能得到理想结果，最主要的是，程序停止的点gradient并不接近于0.我郁闷了很久。虽然matlab里面有现成的garchfit，可是因为我要做的不是基本的GARCH模型，所以一定要自己编的过关了才行。郁闷！！！
谁有估计GARCH(1,1)的程序，能发给我吗?快疯了！万分感谢先。邮箱limuyi1981@gmail.com