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


Posted in Python onNovember 20, 2019

背景:

用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.5
x[0] = 2
for 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.01
x[0] = 2
for 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.2
b = 0.7
x[0] = 2
for 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.3
b = 0.8
x[0] = 2
for 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.2
b = 0.8
x[0] = 2
for 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模型时序图就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
Python中无限元素列表的实现方法
Aug 18 Python
Python中用PIL库批量给图片加上序号的教程
May 06 Python
Python Web框架Tornado运行和部署
Oct 19 Python
机器学习经典算法-logistic回归代码详解
Dec 22 Python
Python读取英文文件并记录每个单词出现次数后降序输出示例
Jun 28 Python
set在python里的含义和用法
Jun 24 Python
Django Rest framework权限的详细用法
Jul 25 Python
对pytorch中的梯度更新方法详解
Aug 20 Python
python+OpenCV实现车牌号码识别
Nov 08 Python
如何将你的应用迁移到Python3的三个步骤
Dec 22 Python
关于tensorflow的几种参数初始化方法小结
Jan 04 Python
Python闭包与装饰器原理及实例解析
Apr 30 Python
Python笔记之工厂模式
Nov 20 #Python
Python常用模块logging——日志输出功能(示例代码)
Nov 20 #Python
将python2.7添加进64位系统的注册表方式
Nov 20 #Python
10个Python面试常问的问题(小结)
Nov 20 #Python
python使用pip安装SciPy、SymPy、matplotlib教程
Nov 20 #Python
Python笔记之facade模式
Nov 20 #Python
将python安装信息加入注册表的示例
Nov 20 #Python
You might like
咖啡店都有些什么常规豆子呢?有什么风味在里面
2021/03/04 咖啡文化
简单的移动设备检测PHP脚本代码
2011/02/19 PHP
php学习之简单计算器实现代码
2011/06/09 PHP
PHP获取数组中重复最多的元素的实现方法
2014/11/11 PHP
PHP简单判断手机设备的方法
2016/08/23 PHP
YII2框架中查询生成器Query()的使用方法示例
2020/03/18 PHP
laravel数据库查询结果自动转数组修改实例
2021/02/27 PHP
一个简单的JS鼠标悬停特效具体方法
2013/06/17 Javascript
Jquery原生态实现表格header头随滚动条滚动而滚动
2014/03/18 Javascript
简单实现异步编程promise模式
2015/07/31 Javascript
jQuery实现的产品自动360度旋转展示特效源码分享
2015/08/21 Javascript
Node.js实用代码段之正确拼接Buffer
2016/03/17 Javascript
浅谈JavaScript 执行环境、作用域及垃圾回收
2016/05/31 Javascript
基于node.js依赖express解析post请求四种数据格式
2017/02/13 Javascript
JavaScript实现二分查找实例代码
2017/02/22 Javascript
微信小程序 wx.login解密出现乱码的问题解决办法
2017/03/10 Javascript
React Native基础入门之初步使用Flexbox布局
2018/07/02 Javascript
基于vue中keep-alive缓存问题的解决方法
2018/09/21 Javascript
webpack4.0 入门实践教程
2018/10/08 Javascript
微信实现自动跳转到用其他浏览器打开指定APP下载
2019/02/15 Javascript
微信小程序开发之转发分享功能
2019/10/22 Javascript
JS操作字符串转数字的常见方法示例
2019/10/29 Javascript
Vue+Vuex实现自动登录的知识点详解
2020/03/04 Javascript
Python Numpy库常见用法入门教程
2020/01/16 Python
django 模版关闭转义方式
2020/05/14 Python
一文读懂python Scrapy爬虫框架
2021/02/24 Python
美国宠物商店:Wag.com
2016/10/25 全球购物
意大利顶级奢侈品电商:LUISAVIAROMA(支持中文)
2020/05/26 全球购物
数控技术专科生自我评价
2014/01/08 职场文书
小孩百日宴答谢词
2014/01/15 职场文书
评析教师个人的自我评价
2014/02/19 职场文书
人力管理专业毕业生求职信
2014/02/27 职场文书
高等教育专业自荐信范文
2014/03/26 职场文书
党的群众路线教育实践活动学习计划
2014/11/03 职场文书
优秀教师先进事迹材料
2014/12/15 职场文书
优秀教师先进材料
2014/12/16 职场文书