如何在Python中将ARIMA模型的方差设置为GARCH

Question

在Matlab中，我可以将ARIMA模型的方差拟合为GARCH模型，例如：

rng(1);
y=randn(100,1)
y=.01+cumsum(y)*.0025;
mdl=arima(3,0,3);
mdl.Variance=garch(1,1)

如何在python中以这种方式拟合arima-garch模型？ 谢谢！