用对数极大似然函数估计写出类似GARCH模型的程序代码，
该模型中与GARCH模型不同的是在波动方程中新增了一个虚拟的解释变量（分别研究该解释变量在正和负对波动的影响）
不知道如何编写程序啊