使用python實(shí)現(xiàn)畫AR模型時(shí)序圖
更新時(shí)間：2019年11月20日 09:29:56 作者：潛水的飛魚baby
今天小編就為大家分享一篇使用python實(shí)現(xiàn)畫AR模型時(shí)序圖，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過來看看吧
背景：
用python畫AR模型的時(shí)序圖。
結(jié)果：
代碼：
import numpy as np
import matplotlib.pyplot as plt
"""
AR(1)的時(shí)序圖：x[t]=a*x[t-1]+e
"""
num = 2000
e = np.random.rand(num)
x = np.empty(num)
"""
平穩(wěn)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")
"""
非平穩(wěn)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")
"""
平穩(wěn)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")
"""
非平穩(wěn)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")
"""
非平穩(wěn)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實(shí)現(xiàn)畫AR模型時(shí)序圖就是小編分享給大家的全部?jī)?nèi)容了，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。
您可能感興趣的文章:
相關(guān)文章
Python項(xiàng)目文件中安裝?setup.py的步驟
大家都知道setup.py是python模塊分發(fā)與安裝的指導(dǎo)文件，這篇文章主要介紹了Python項(xiàng)目文件中如何安裝?setup.py,需要的朋友可以參考下2023-05-05
想學(xué)python，這5本書籍你必看！本文為大家推薦了學(xué)習(xí)python的5本書籍，5本經(jīng)典書籍，感興趣的小伙伴們可以參考一下2018-12-12
Python編程使用有限狀態(tài)機(jī)識(shí)別地址有效性
這篇文章主要介紹了Python編程中如何使用有限狀態(tài)機(jī)識(shí)別地址有效性，如何識(shí)別一個(gè)地址是否有效，確切的講，如何編程識(shí)別一個(gè)中國地址是否有效2021-09-09![Pandas庫中iloc[]函數(shù)的使用方法](http://img.jbzj.com/images/xgimg/bcimg3.png)
在數(shù)據(jù)分析過程中,很多時(shí)候需要從數(shù)據(jù)表中提取出相應(yīng)的數(shù)據(jù),而這么做的前提是需要先“索引”出這一部分?jǐn)?shù)據(jù),下面這篇文章主要給大家介紹了關(guān)于Pandas庫中iloc[]函數(shù)的使用方法,需要的朋友可以參考下2023-01-01
Python中的類型提示(Type Hints)總結(jié)
Python3.5 版本引入了類型提示（Type Hints），它允許開發(fā)者在代碼中顯式地聲明變量、函數(shù)、方法等的類型信息，下面小編就來帶大家一起看看Python類型提示的初步使用吧2023-05-05
python保留格式匯總各部門excel內(nèi)容的實(shí)現(xiàn)思路
這篇文章主要介紹了python保留格式匯總各部門excel內(nèi)容,本文通過實(shí)例代碼給大家介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值，需要的朋友可以參考下2020-06-06
python-opencv 中值濾波{cv2.medianBlur(src, ksize)}的用法
這篇文章主要介紹了python-opencv 中值濾波{cv2.medianBlur(src, ksize)}的用法，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。2021-06-06
python實(shí)現(xiàn)超市管理系統(tǒng)(后臺(tái)管理)
這篇文章主要為大家詳細(xì)介紹了python實(shí)現(xiàn)超市管理系統(tǒng)，增加后臺(tái)管理，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2019-10-10
Pytorch計(jì)算網(wǎng)絡(luò)參數(shù)的兩種方法
PyTorch是一個(gè)流行的深度學(xué)習(xí)框架,它允許研究人員和開發(fā)者快速構(gòu)建和訓(xùn)練神經(jīng)網(wǎng)絡(luò),計(jì)算一個(gè)PyTorch網(wǎng)絡(luò)的參數(shù)量通常涉及兩個(gè)步驟,本文給大家介紹了在PyTorch中計(jì)算網(wǎng)絡(luò)參數(shù)量的一般方法,需要的朋友可以參考下2024-05-05
python CMD命令行傳參實(shí)現(xiàn)方法(argparse、click、fire)
這篇文章主要介紹了python CMD命令行傳參實(shí)現(xiàn)方法(argparse、click、fire),本文通過實(shí)例代碼給大家介紹的非常詳細(xì),感興趣的朋友跟隨小編一起看看吧2024-07-07