我想估计汇率的 GARCH 模型。我有对数似然的公式，并且我有一个平方误差的向量。但是我不知道如何估算参数，因为我只有\sigma_1。
我要最大化：
-\frac{N}{2}log(2\pi) - \frac{1}{2} \sum_{t=1}^N log \sigma_t^2 - \frac{1}{2}\sum_{t=1}^N \frac{\epsilon_t^2}{\sigma_t^2}
我想用python编程，并想使用scipy.optimize，但是我不确定要使用哪种方法。你有提示吗？
但是我的主要问题是\sigma必须每次都更新，但是随后我需要尚未估计的参数。