Hurst指數(shù)與R/S(python代碼)

-- coding: utf-8 --

"""
Created on Mon Oct 31 21:39:30 2016

@author: Administrator
"""
import numpy as np
x=data
data2=[]
for i in range(1,len(x)):
data1= x[i] - x[i-1]
data2.append(data1)

mdata=[]
for i in range(len(data2)):
m=np.mean(data2[:i+1])
mdata.append(m)

r=[]
for i in range(len(mdata)):
print('.......')
c=[]
b=[]
for j in range(0,i+1):
p=data2[j]-mdata[i]
c.append(p)
print(c)
b.append(np.cumsum(c)[-1])
r.append(max(b)-min(b))

sdata=[]

for i in range(len(data2)):
std=np.std(data2[:i+1])
sdata.append(std)
v = list(map(lambda x: x[0]/x[1], zip(r, sdata)))
rs=np.log(v)

t=np.array(range(1,11))/2
g=np.polyfit(np.log(t)[1:],rs[1:],1)
rt=pow(2,(2*g[0]-1))-1

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容