摘要：隨著數(shù)字經(jīng)濟(jì)的快速發(fā)展，比特幣成為當(dāng)代重要的金融討論話(huà)題。比特幣是一種基于區(qū)塊鏈技術(shù)，保證整個(gè)支付體系和虛擬貨幣產(chǎn)生和使用的點(diǎn)對(duì)點(diǎn)協(xié)議，比特幣具有去中心化、匿名性和內(nèi)置激勵(lì)性等特征，缺點(diǎn)是缺少信用背書(shū)、幣值不穩(wěn)定和發(fā)行者過(guò)多占據(jù) “鑄幣稅”。由于比特幣主要扮演的是資產(chǎn)角色而非充當(dāng)貨幣角色，其市場(chǎng)投機(jī)性強(qiáng)和幣值波動(dòng)性大的特點(diǎn)，因而和其他投資市場(chǎng)相比，因而比特幣存在市場(chǎng)效率低下的現(xiàn)象，故研究比特幣價(jià)格的預(yù)測(cè)與分析具有重要意義。本文主要選用2019年3月到2021年7月的比特幣收盤(pán)價(jià)，采用了ARIMA模型，通過(guò)使用R語(yǔ)言對(duì)數(shù)據(jù)進(jìn)行平穩(wěn)化處理并進(jìn)行單位根檢驗(yàn)，對(duì)比各參數(shù)，建立合理的ARIMA模型，同時(shí)，與自動(dòng)生成的ARIMA模型進(jìn)行比較，選擇較為優(yōu)良的模型進(jìn)行比特幣價(jià)格的短期預(yù)測(cè)。

關(guān)鍵詞：比特幣;ARIMA模型;時(shí)間序列

引言

2008年10月，一位名為中本聰密碼學(xué)者發(fā)布了《比特幣： 一種點(diǎn)對(duì)點(diǎn)的電子現(xiàn)金系統(tǒng)》的文章，文內(nèi)描述了一種被他稱(chēng)為“比特幣”的電子貨幣及其算法。隨后，中本聰挖出第一個(gè)區(qū)塊，即所謂的創(chuàng)世區(qū)塊，由此比特幣正式誕生。比特幣是基于區(qū)塊鏈技術(shù)，具有分布式賬本、可追溯、不可篡改、去中心化等特點(diǎn)，為金融創(chuàng)新帶來(lái)了新的可能。在比特幣出現(xiàn)后，還有以太坊、瑞波幣等基于區(qū)塊鏈技術(shù)的數(shù)字貨幣大量涌現(xiàn)。

比特幣的底層技術(shù)是區(qū)塊鏈。區(qū)塊鏈，顧名思義是由“區(qū)塊”和 “鏈”組成的。區(qū)塊簡(jiǎn)單而言就是一個(gè)信息塊，通過(guò)加密算法記錄交易信息，如果區(qū)塊的交易信息被認(rèn)證、接受，那就可以寫(xiě)入鏈條中，并發(fā)布給所有的節(jié)點(diǎn)，形成區(qū)塊鏈。所以說(shuō)，區(qū)塊鏈可以理解為一個(gè)嚴(yán)格按照時(shí)間順序排列的交易信息記錄簿，而且分布在所有節(jié)點(diǎn)上。

2009年10月，1美元相當(dāng)于1309. 03比特幣;2010年11月6日，第一個(gè)比特幣交易所成立，1比特幣相當(dāng)于0.5美元;此后，比特幣的價(jià)格出現(xiàn)飛躍式上漲，在2013年11月漲至超過(guò)1000美元的歷史高點(diǎn);2014年2月，當(dāng)時(shí)世界最大的比特幣交易所 Mt. Gox 被盜，引發(fā)虛擬貨幣歷史上的第二次大熊市;到2020年，比特幣的價(jià)格暴漲，又回到超過(guò) 1 萬(wàn)美元的高位。可見(jiàn)比特幣的價(jià)格是極度波動(dòng)的，而隨機(jī)的、大幅的波動(dòng)性，使得以比特幣為載體的商業(yè)交易面臨著極大的不確定性，同時(shí)比特幣主要用于資產(chǎn)而非貨幣，在金融市場(chǎng)和投資組合管理中發(fā)揮重要作用，因此，對(duì)比特幣波動(dòng)性的研究具有重要的現(xiàn)實(shí)意義。

1.時(shí)間序列模型——ARIMA模型

假設(shè)時(shí)間序列{X_t}具有如下形式：

X_t=φ_0+φ_1 X_（t-1）+φ_2 X_（t-2）+...+φ_p X_（t-p）+ε_(tái)t-θ_1 ε_(tái)（t-1）-θ_2 ε_(tái)（t-2）-...-θ_q ε_(tái)（t-q） ?（1）則稱(chēng)式（1）為自回歸移動(dòng)平均模型，記為 ARMA（ p，q） 模型，稱(chēng){X_t}為ARMA（ p，q）過(guò)程，其中，E（ε_(tái)t）=0，var（ε_(tái)t）=σ_ε^2，cov（ε_(tái)t，ε_(tái)t）=0（s≠t），cov（x_s，ε_(tái)t）=0（?s<t）

如果一個(gè)時(shí)間序列{X_t}的d次差分Wt = ?^d X_t時(shí)ARMA（p，q）過(guò)程，

即 W_t=φ_1 W_（t-1）+φ_2 W_（t-2）+...+φ_p W_（t-p）+ε_(tái)t-θ_1 ε_(tái)（t-1）-θ_2 ε_(tái)（t-2）-...θ_q ε_(tái)（t-1） ? （2）

則稱(chēng)上式（2）為自回歸滑動(dòng)平均求和模型，記為ARIMA（p，d，q）模型，稱(chēng){X_t}為ARIMA（p，d，q）過(guò)程。當(dāng) d = 0，ARIMA（p，0，q） 模型實(shí)際上就是 ARMA（p，q）;當(dāng) p = 0，ARMA（0，d，q） 模型可以簡(jiǎn)記為IMA（d，q） 模型;當(dāng) q = 0 時(shí)，ARMA（p，d，0）模型可以簡(jiǎn)記為ARI（p，d）模型。

ARIMA模型是差分整合移動(dòng)平均自回歸模型，又稱(chēng)整合移動(dòng)平均自回歸模型，是時(shí)間序列預(yù)測(cè)分析方法之一。通常對(duì)線(xiàn)性趨勢(shì)可以用一階差分可以使之平穩(wěn)化，對(duì)二階曲線(xiàn)使用二階差分。在實(shí)際應(yīng)用過(guò)程中，可能由于過(guò)差分使數(shù)據(jù)失真，因此產(chǎn)生了分?jǐn)?shù)階差分來(lái)進(jìn)行優(yōu)化，本文使用一階差分，因此對(duì)分?jǐn)?shù)階差分不再贅述。

2.實(shí)證分析

2.1建立時(shí)間序列

本文所使用的數(shù)據(jù)主要是2019年10月至2021年1月比特幣交易平臺(tái)指數(shù)的每日收盤(pán)價(jià)原始數(shù)據(jù)。這些數(shù)據(jù)在http：//www.coindesk上公開(kāi)，選取數(shù)據(jù)包中的收盤(pán)價(jià)作時(shí)間序列。

2.2平穩(wěn)性檢驗(yàn)

從該圖1可知，該序列是非平穩(wěn)的，由于建立ARIMA模型需要平穩(wěn)性序列，因此要對(duì)該序列進(jìn)行差分，經(jīng)過(guò)R語(yǔ)言返回的結(jié)果是一階差分，因此對(duì)其進(jìn)行一階差分處理，再繪制圖形。

從圖2可大體看出，一階差分后的序列是平穩(wěn)的，同時(shí)對(duì)差分后的序列做ADF單位根檢驗(yàn)，檢驗(yàn)結(jié)果顯示p值是0.01，所以我們拒絕原假設(shè)，即認(rèn)為一階差分處理后的序列是平穩(wěn)的。

2.3模型的定階及擬合

接下來(lái)對(duì)序列作一階差分后的自相關(guān)圖和偏自相關(guān)圖。

從圖3自相關(guān)圖可知，在滯后階數(shù)逐漸增加時(shí)，自相關(guān)系數(shù)逐漸減小到0.

從偏自相關(guān)圖可發(fā)現(xiàn)較難確定模型中的q值，因此用forecast包中的auto.arima（） 函數(shù)來(lái)建立最優(yōu)模型。

2.4 自動(dòng)建模

利用R語(yǔ)言自帶的forecast包中的auto.arima（） 函數(shù)進(jìn)行自動(dòng)選擇模型。通過(guò) R 程序自動(dòng)選擇了 ARIMA（4，1，4） 模型，且AIC值為17149.56，標(biāo)準(zhǔn)差方差估計(jì)值為39749012，比上述建模后的AIC值小，說(shuō)明模型ARIMA（4，1，4）比模型ARIMA（2，1，2）效果要好。

2.5模型診斷

進(jìn)行模型的診斷需要看兩方面，一方面是殘差是否服從正態(tài)分布，另一方面是看殘差之間是否相關(guān)。首先我們繪制Q-Q圖觀(guān)察殘差是否服從正態(tài)分布。

從圖5可看出，它處于45度分位線(xiàn)上，所以認(rèn)為是服從正態(tài)分布的。接下來(lái)用Ljung-Bo函數(shù)來(lái)檢驗(yàn)殘差之間是否相關(guān)，結(jié)果如下：

從結(jié)果可知p值為0.9253大于0.05，因此接受原假設(shè)，即認(rèn)為殘差之間不相關(guān)，也就是殘差是平穩(wěn)的。

2.6用ARIMA模型進(jìn)行預(yù)測(cè)

由前文可知ARIMA（4，1，4） 模型的AIC值比模型ARIMA（2，1，2）的AIC值小，故選擇ARIMA（4，1，4） 模型使用forecast函數(shù)來(lái)向前預(yù)測(cè)五期的收盤(pán)價(jià)，結(jié)果如下：

該函數(shù)預(yù)測(cè)了接下來(lái)三個(gè)月的比特幣收盤(pán)價(jià)，第一列表示的是點(diǎn)估計(jì)值，第二列和第三列表示的是百分之80的置信下限和百分之80的置信上限，第四列和第五列表示的是百分之九十五的置信下限和百分之95的置信上限，最后畫(huà)下模型預(yù)測(cè)圖。

由圖6可知，圖中三個(gè)點(diǎn)表示點(diǎn)預(yù)測(cè)值，藍(lán)色的區(qū)間表示80%的置信區(qū)間，灰色的表示95%的置信區(qū)間。

3.結(jié)論

本文利用時(shí)間序列分析的建模思想，對(duì)比特幣的收盤(pán)價(jià)進(jìn)行預(yù)測(cè)，通過(guò)模型識(shí)別、建立和檢驗(yàn)，確定ARIMA（4，1，4）模型為比特幣收盤(pán)價(jià)的預(yù)測(cè)模型。通過(guò)模型的預(yù)測(cè)發(fā)現(xiàn)，ARIMA 模型在實(shí)際中有著廣泛的適用性，可以推廣到其他領(lǐng)域的預(yù)測(cè)，不過(guò)需要根據(jù)所要解決的問(wèn)題和問(wèn)題的特點(diǎn)等因素來(lái)綜合考慮，選擇出相對(duì)最優(yōu)的預(yù)測(cè)模型。

綜上所述，ARMA模型充分利用有限的數(shù)據(jù)集對(duì)未來(lái)的發(fā)展趨勢(shì)進(jìn)行了預(yù)測(cè)分析。通過(guò) ARMA模型進(jìn)行建模和實(shí)證分析得到的短期比特幣的預(yù)測(cè)值較為理想，為學(xué)者分析價(jià)格波動(dòng)提供了一定的參考意義。

參考文獻(xiàn)：

[1]劉剛，劉娟，唐婉容.比特幣價(jià)格波動(dòng)與虛擬貨幣風(fēng)險(xiǎn)防范——基于中美政策信息的事件研究法[J].廣東財(cái)經(jīng)大學(xué)學(xué)報(bào)，2015，30（03）：30-40.

[2]段格格.基于A(yíng)RIMA模型的比特幣價(jià)格預(yù)測(cè)研究[J].現(xiàn)代營(yíng)銷(xiāo)（下旬刊），2021（01）：27-29.

[3]蔡承智，黃軍結(jié)，梁穎.基于A(yíng)RIMA模型的世界大豆總產(chǎn)預(yù)測(cè)分析[J].大豆科學(xué)，2019，38（02）：298-303.

[4]Saralees Nadarajah，Jeffrey Chu. On the inefficiency of Bitcoin[J]. Economics Letters，2017，150：

[5]Yonghong Jiang，He Nie，Weihua Ruan. Time-varying long-term memory in Bitcoin market[J]. Finance Research Letters，2018，25.

作者簡(jiǎn)介：朱曉婕（1997--） ，女，漢族，籍貫：江蘇南通，碩士，單位：福建師范大學(xué)，研究方向：金融