使用python实现画AR模型时序图

概述使用python实现画AR模型时序图 背景: 用python画AR模型的时序图. 结果: 代码: import numpy as np import matplotlib.pyplot as plt """ AR(1)的时序图:x[t]=a*x[t-1]+e """ num = 2000 e = np.random.rand(num) x = np.empty(num) """ 平稳AR(1) """ a = -0 ...

背景：

用python画AR模型的时序图。

结果：

代码：

import numpy as npimport matplotlib.pyplot as plt"""AR(1)的时序图：x[t]=a*x[t-1]+e"""num = 2000e = np.random.rand(num)x = np.empty(num)"""平稳AR(1)"""a = -0.5x[0] = 2for i in range(1,num): x[i] = a*x[i-1]+e[i]plt.subplot(321,Title = "AR({0}):x[t]={1}*x[t-1]+e".format(1,a))plt.plot(x,"or")"""非平稳AR(1)"""a = -1.01x[0] = 2for i in range(1,num): x[i] = a*x[i-1]+e[i]plt.subplot(322,".b")"""平稳AR(2)"""a = -0.2b = 0.7x[0] = 2for i in range(2,num): x[i] = a*x[i-1]+b*x[i-2]+e[i]plt.subplot(323,Title = "AR({0}):x[t]={1}*x[t-1]+{2}*x[t-2]+e".format(2,a,b))plt.plot(x,"og")"""非平稳AR(2)"""a = -0.3b = 0.8x[0] = 2for i in range(2,num): x[i] = a*x[i-1]+b*x[i-2]+e[i]plt.subplot(324,".y")"""非平稳AR(2)"""a = -0.2b = 0.8x[0] = 2for i in range(2,num): x[i] = a*x[i-1]+b*x[i-2]+e[i]plt.subplot(313,"+",color="purple")plt.show()

以上这篇使用python实现画AR模型时序图就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持我们。

总结

以上是内存溢出为你收集整理的使用python实现画AR模型时序图全部内容，希望文章能够帮你解决使用python实现画AR模型时序图所遇到的程序开发问题。

如果觉得内存溢出网站内容还不错，欢迎将内存溢出网站推荐给程序员好友。