statsmodels.sandbox.tsa.fftarma.ArmaFft.invpowerspd

ArmaFft.invpowerspd(n)[source]

从谱密度计算自协方差

缩放是正确的,但为了数值精度,n需要很大 可能在fft中用零填充会更快 不进行切片时,返回带有fftshift的双边自协方差

>>> ArmaFft([1, -0.5], [1., 0.4], 40).invpowerspd(2**8)[:10]
array([ 2.08    ,  1.44    ,  0.72    ,  0.36    ,  0.18    ,  0.09    ,
        0.045   ,  0.0225  ,  0.01125 ,  0.005625])
>>> ArmaFft([1, -0.5], [1., 0.4], 40).acovf(10)
array([ 2.08    ,  1.44    ,  0.72    ,  0.36    ,  0.18    ,  0.09    ,
        0.045   ,  0.0225  ,  0.01125 ,  0.005625])

Last update: Oct 16, 2024