做GARCH(1,1)模型设定如下：
yt = mu + et, et ~ N(0, ht)
ht = w + a e_{t-1}^2 + b h_{t-1}
时间序列yt=常数值+扰动项，扰动项满足garch(1,1)模型，服从正态分布。要求用matlab程序实现联合估计参数mu、w、a、b~
我自己的思路是先写出极大似然函数likelihood，然后用fmincon得到上述参数，但是我写的程序总是出现各种各样的问题，这几天快要把我憋疯了！！！这个程序应该很简单的，不知为何总是出错，求大神们帮忙编个程序，我的邮箱464728396@qq.com。
programe是我编的主程序，likelihood是我编的极大似然估计函数，yt是原始数据。在我的程序基础上修改一下应该很快吧，楼主很捉急，这几天死磕就是磕不出来。。。。。要求用极大似然估计的方法联合估计参数，很急，谢谢！！！