使用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实现简单的Markdown转换器
Jul 16 Python
win10 64bit下python NLTK安装教程
Sep 19 Python
python 内置模块详解
Jan 01 Python
Python常见的pandas用法demo示例
Mar 16 Python
Python分布式进程中你会遇到的问题解析
May 28 Python
python 设置xlabel,ylabel 坐标轴字体大小,字体类型
Jul 23 Python
浅谈python3中input输入的使用
Aug 02 Python
Flask框架学习笔记之表单基础介绍与表单提交方式
Aug 12 Python
keras模型可视化,层可视化及kernel可视化实例
Jan 24 Python
如何将tensorflow训练好的模型移植到Android (MNIST手写数字识别)
Apr 22 Python
Python selenium 加载并保存QQ群成员,去除其群主、管理员信息的示例代码
May 28 Python
python爬取youtube视频的示例代码
Mar 03 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漏洞全解(详细介绍)
2012/11/13 PHP
PHP数据对象PDO操作技巧小结
2016/09/27 PHP
比较全的JS checkbox全选、取消全选、删除功能代码
2008/12/19 Javascript
javascript 打开页面window.location和window.open的区别
2010/03/17 Javascript
javascript数组操作(创建、元素删除、数组的拷贝)
2014/04/07 Javascript
jQuery实现的原图对比窗帘效果
2014/06/15 Javascript
jQuery实现瀑布流的取巧做法分享
2015/01/12 Javascript
$.extend 的一个小问题
2015/06/18 Javascript
基于javascript实现动态显示当前系统时间
2016/01/28 Javascript
AngularJS入门教程中SQL实例详解
2016/07/27 Javascript
详解Angular-Cli中引用第三方库
2017/05/21 Javascript
vue.js实现标签页切换效果
2018/06/07 Javascript
vue-cli2.x项目优化之引入本地静态库文件的方法
2018/06/19 Javascript
jQuery实现表单动态添加数据并提交的方法
2018/07/19 jQuery
对angularJs中$sce服务安全显示html文本的实例
2018/09/30 Javascript
纯js实现无缝滚动功能代码实例
2020/02/21 Javascript
python使用点操作符访问字典(dict)数据的方法
2015/03/16 Python
整理Python中的赋值运算符
2015/05/13 Python
Python实现以时间换空间的缓存替换算法
2016/02/19 Python
Python socket实现简单聊天室
2018/04/01 Python
Python中使用pypdf2合并、分割、加密pdf文件的代码详解
2019/05/21 Python
Python Socket TCP双端聊天功能实现过程详解
2020/06/15 Python
python实现简单贪吃蛇游戏
2020/09/29 Python
美国著名的女性内衣零售商:Frederick’s of Hollywood
2018/02/24 全球购物
我的珠宝盒:Ma boîte à bijoux
2019/08/27 全球购物
创意爱尔兰礼物:Creative Irish Gifts
2020/01/29 全球购物
Chinti & Parker官网:奢华羊绒女装和创新针织设计
2021/01/01 全球购物
医院护士专业个人的求职信
2013/12/09 职场文书
《跨越海峡的生命桥》教学反思
2014/02/24 职场文书
安全生产月演讲稿
2014/05/09 职场文书
小学校长竞聘演讲稿
2014/05/16 职场文书
教师自荐信范文
2015/03/06 职场文书
2015年体育教学工作总结
2015/05/20 职场文书
详细介绍Next.js脚手架完整搭建封装
2022/04/26 Javascript
vue递归实现树形组件
2022/07/15 Vue.js
win10截图快捷键win+shift+s没有反应无法截图怎么解决?
2022/08/14 数码科技