使用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使用os模块的os.walk遍历文件夹示例
Jan 27 Python
Python time模块详解(常用函数实例讲解,非常好)
Apr 24 Python
python中string模块各属性以及函数的用法介绍
May 30 Python
Python获取暗黑破坏神3战网前1000命位玩家的英雄技能统计
Jul 04 Python
Python编程实现删除VC临时文件及Debug目录的方法
Mar 22 Python
利用python对Excel中的特定数据提取并写入新表的方法
Jun 14 Python
python opencv实现图片旋转矩形分割
Jul 26 Python
python模拟登陆,用session维持回话的实例
Dec 27 Python
python 爬虫百度地图的信息界面的实现方法
Oct 27 Python
基于python实现简单网页服务器代码实例
Sep 14 Python
Python使用openpyxl复制整张sheet
Mar 24 Python
Python实现滑雪小游戏
Sep 25 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
php.ini中date.timezone设置分析
2011/07/29 PHP
解析如何修改phpmyadmin中的默认登陆超时时间
2013/06/25 PHP
Javascript 错误处理的几种方法
2009/06/13 Javascript
IE8下String的Trim()方法失效的解决方法
2013/11/08 Javascript
jquery validate在ie8下的bug解决方法
2013/11/13 Javascript
javascript检查浏览器是否支持flash的实现代码
2014/08/14 Javascript
Javascript表单验证要注意的事项
2014/09/29 Javascript
jQuery点击输入框显示验证码图片
2016/05/19 Javascript
AngularJS入门教程之Scope(作用域)
2016/07/27 Javascript
JQuery遍历元素的后代和同胞实现方法
2016/09/18 Javascript
微信小程序 教程之事件
2016/10/18 Javascript
JavaScript中关于iframe滚动条的去除和保留
2016/11/17 Javascript
详解本地Node.js服务器作为api服务器的解决办法
2017/02/28 Javascript
Laravel整合Bootstrap 4的完整方案(推荐)
2018/01/25 Javascript
详解创建自定义的Angular Schematics
2018/06/06 Javascript
详解vue使用vue-layer-mobile组件实现toast,loading效果
2018/08/31 Javascript
35个最好用的Vue开源库(史上最全)
2019/01/03 Javascript
详解vue 不同环境配置不同的打包命令
2019/04/07 Javascript
小程序云开发实现数据库异步操作同步化
2019/05/18 Javascript
JavaScript中0、空字符串、'0'是true还是false的知识点分享
2019/09/16 Javascript
微信小程序实现禁止分享代码实例
2019/10/19 Javascript
深入浅析vue全局环境变量和模式
2020/04/28 Javascript
python批量提交沙箱问题实例
2014/10/08 Python
Python中的引用和拷贝浅析
2014/11/22 Python
详解字典树Trie结构及其Python代码实现
2016/06/03 Python
使用PyV8在Python爬虫中执行js代码
2017/02/16 Python
用TensorFlow实现戴明回归算法的示例
2018/05/02 Python
查看python下OpenCV版本的方法
2018/08/03 Python
Python批量修改xml的坐标值全部转为整数的实例代码
2020/11/26 Python
大学生自荐信范文
2015/03/05 职场文书
毕业生自荐信范文
2015/03/05 职场文书
导盲犬小Q观后感
2015/06/11 职场文书
MySQL七种JOIN类型小结
2021/10/24 MySQL
Python识别花卉种类鉴定网络热门植物并自动整理分类
2022/04/08 Python
SQL Server使用PIVOT与unPIVOT实现行列转换
2022/05/25 SQL Server
numpy array找出符合条件的数并赋值的示例代码
2022/06/01 Python