使用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中的命令行参数解析工具之docopt详解
Mar 27 Python
Python制作刷网页流量工具
Apr 23 Python
python使用生成器实现可迭代对象
Mar 20 Python
Python向MySQL批量插数据的实例讲解
Mar 31 Python
Python利用公共键如何对字典列表进行排序详解
May 19 Python
详解pandas安装若干异常及解决方案总结
Jan 10 Python
python+pyqt5编写md5生成器
Mar 18 Python
python利用Opencv实现人脸识别功能
Apr 25 Python
python模拟斗地主发牌
Apr 22 Python
tensorflow 20:搭网络,导出模型,运行模型的实例
May 26 Python
python集合能干吗
Jul 19 Python
python如何遍历指定路径下所有文件(按按照时间区间检索)
Sep 14 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
Android ProgressBar进度条和ProgressDialog进度框的展示DEMO
2013/06/19 PHP
thinkphp中html:list标签传递多个参数实例
2014/10/30 PHP
基于PHP实现微信小程序客服消息功能
2019/08/12 PHP
laravel5.6 框架邮件队列database驱动简单demo示例
2020/01/26 PHP
JavaScript 特殊字符
2007/04/05 Javascript
JS 的应用开发初探(mootools)
2009/12/19 Javascript
JavaScript获取鼠标移动时的坐标(兼容IE8、chome谷歌、Firefox)
2014/09/13 Javascript
jQuery如何获取同一个类标签的所有值(默认无法获取)
2014/09/25 Javascript
javaScript知识点总结(必看篇)
2016/06/10 Javascript
多种方式实现js图片预览
2016/12/12 Javascript
基于JavaScript实现的折半查找算法示例
2017/04/14 Javascript
es7学习教程之Decorators(修饰器)详解
2017/07/21 Javascript
Vue resource中的GET与POST请求的实例代码
2017/07/21 Javascript
Angular @HostBinding()和@HostListener()用法
2018/03/05 Javascript
深入浅出理解JavaScript闭包的功能与用法
2018/08/01 Javascript
微信小程序 image组件遇到的问题
2019/05/28 Javascript
vue中watch和computed的区别与使用方法
2020/08/23 Javascript
windows下安装python paramiko模块的代码
2013/02/10 Python
python解析html开发库pyquery使用方法
2014/02/07 Python
基于Python闭包及其作用域详解
2017/08/28 Python
python使用__slots__让你的代码更加节省内存
2018/09/05 Python
Python2.7实现多进程下开发多线程示例
2019/05/31 Python
Gauss-Seidel迭代算法的Python实现详解
2019/06/29 Python
Python3显示当前时间、计算时间差及时间加减法示例代码
2019/09/07 Python
Python jieba库用法及实例解析
2019/11/04 Python
Python中BeautifulSoup通过查找Id获取元素信息
2020/12/07 Python
Python plt 利用subplot 实现在一张画布同时画多张图
2021/02/26 Python
CSS的pointer-events属性详细介绍(作用和注意事项)
2014/04/23 HTML / CSS
Canvas与图片压缩的示例代码
2017/11/28 HTML / CSS
SIXPAD智能健身仪英国官网:革命性的训练装备品牌
2018/09/27 全球购物
爱尔兰橄榄球店:Irish Rugby Store
2019/12/05 全球购物
swtich是否能作用在byte上,是否能作用在long上,是否能作用在String上?
2013/03/30 面试题
师范生免费教育协议书范本
2014/10/09 职场文书
幼儿园大班教学反思
2016/03/02 职场文书
Redis入门教程详解
2021/08/30 Redis
python 详解turtle画爱心代码
2022/02/15 Python