使用python实现画AR模型时序图
Python  /  管理员 发布于 7年前   163
背景:
用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,title = "AR({0}):x[t]={1}*x[t-1]+e".format(1,a))plt.plot(x,".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,title = "AR({0}):x[t]={1}*x[t-1]+{2}*x[t-2]+e".format(2,a,b))plt.plot(x,".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,title = "AR({0}):x[t]={1}*x[t-1]+{2}*x[t-2]+e".format(2,a,b))plt.plot(x,"+",color="purple") plt.show()
以上这篇使用python实现画AR模型时序图就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持。
122 在
学历:一种延缓就业设计,生活需求下的权衡之选中评论 工作几年后,报名考研了,到现在还没认真学习备考,迷茫中。作为一名北漂互联网打工人..123 在
Clash for Windows作者删库跑路了,github已404中评论 按理说只要你在国内,所有的流量进出都在监控范围内,不管你怎么隐藏也没用,想搞你分..原梓番博客 在
在Laravel框架中使用模型Model分表最简单的方法中评论 好久好久都没看友情链接申请了,今天刚看,已经添加。..博主 在
佛跳墙vpn软件不会用?上不了网?佛跳墙vpn常见问题以及解决办法中评论 @1111老铁这个不行了,可以看看近期评论的其他文章..1111 在
佛跳墙vpn软件不会用?上不了网?佛跳墙vpn常见问题以及解决办法中评论 网站不能打开,博主百忙中能否发个APP下载链接,佛跳墙或极光..
Copyright·© 2019 侯体宗版权所有·
粤ICP备20027696号