时间序列学习（2）：白噪声、随机游走

• • 1、白噪声
  • 2、对数收益率序列
  • 3、随机游走
  • 4、随机游走示例

1、白噪声

白噪声是非常简单的一种建模时间序列的模型。

对于时间序列 { w t } \{w_t\} {wt​}，若满足下面三个条件，该序列为一个离散的白噪声（white noise）：

• 每个时间点均值为0；
• 每个时间点方差为 σ 2 \sigma^2 σ2；
• 对于任意 k ≥ 1 k\ge1 k≥1，自相关系数 ρ k = 0 \rho_k=0 ρk​=0。

前两个条件并没有规定 w t w_t wt​一定是正态分布。如果是正态分布，那么这个时间序列又叫高斯白噪声（Gaussian white noise）。

回顾上一篇笔记中的平稳性要求，显然一个白噪声序列是符合弱平稳性的（其均值和自相关系数均为固定值，比弱平稳性更加严格一些）。

2、对数收益率序列

上一篇笔记中的日收益率序列，不严格地说，其实是一个差分序列。也就是每一个时间点的价格减去前一个时间点的价格（当然，还除以了前一时间点的价格）。

现在我们利用对数收益率处理成纯的差分序列，即每一个时间点价格的对数减去前一个时间点价格的对数。

对于价格的时间序列 { p t } \{p_t\} {pt​}，其对数收益率序列用下面的公式得到：

α t = ln ⁡ p t − ln ⁡ p t − 1 \alpha_t = \ln{p_t}-\ln{p_{t-1}} αt​=lnpt​−lnpt−1​

如果价格序列 { p t } \{p_t\} {pt​}体现出足够的自相关性（今天的价格和昨天强相关），那么对数收益率序列应该没有显著的自相关性了（自己减掉自己了）。也就是说对数收益率序列应该近似为一个白噪声。

我们将上一篇笔记的日收益率序列处理成对数收益率序列（即沪深300自2018年1月1日至2019年12月13日的对数收益率序列），如下：

import pandas as pd
import math


prices = get_price('000300.XSHG', start_date='2018-01-01', end_date='2019-12-13', frequency='daily', fields='close')

returns = {'return':{}}
for i in range(len(prices)):
    current_price = prices.iloc[i,0]
    if i == 0:
        current_return = 0
    else:
        last_price = prices.iloc[i-1,0]
        current_return = math.log(current_price)-math.log(last_price)
    date = list(prices.index)[i]
    returns['return'][date] = current_return
return_df = pd.DataFrame(returns, columns=['return'], index=list(prices.index))

fig = plt.figure(figsize=(10, 6))
ax = fig.add_axes([0.2, 0.2, 1.2, 1.2])
ax.plot(return_df, color="blue", linewidth=1.5, linestyle="-", label=r'hs300-daily-log-returns')
plt.legend(loc='upper right', frameon=False)

return_df.mean()
return_df.std()

它的均值和标准差分别是：-0.000062，0.01315。

然后看看其时间间隔在50以下的相关图（相关图是啥，见上一篇笔记）。

import matplotlib.pyplot as plt
from statsmodels.graphics.tsaplots import plot_acf

fig = plt.figure(figsize=(10, 6))
ax = fig.add_axes([0.2, 0.2, 1.2, 1.2])
plot_acf(return_df, ax, lags=50)
plt.show()

可以看到，上述收益率序列的均值-0.000062，非常接近于0；样本方差不能代表时间序列的方差，所以方差这一项依据我们暂时放在一边。看看自相关性，对于间隔为3的自相关系数超过了95%的置信区间，但其它都在区间内。

不管怎样，我们说上述对数收益率序列近似为一个白噪声。

3、随机游走

随机游走 random walk是对白噪声的简单延伸。

对于时间序列 { x t } \{x_t\} {xt​}，如果满足如下条件，该时间序列是一个随机游走：

x t = w t + w t − 1 + . . . + w 0 x_t=w_t+w_{t-1}+...+w_0 xt​=wt​+wt−1​+...+w0​

其中， w t w_t wt​是一个均值为0，方差为 σ 2 \sigma^2 σ2的白噪声。

所以随机游走是白噪声的一个累积。比如，非严格地说，如果把股价看成随机游走，它其实就是每日收益率（如果看成白噪声）的不断累积。

根据均值和方差的可加性可以求得随机游走的的均值和方差：

μ x t = 0 ; V a r ( x t ) = t σ 2 \mu_{x_t}=0; Var(x_t)=t\sigma^2 μxt​​=0;Var(xt​)=tσ2

进一步，随机游走的自协方差也可以计算出来：

C o v ( x t , x t + k ) = C o v ( x t , x t + w t + 1 + . . . + w k ) Cov(x_t,x_{t+k})=Cov(x_t,x_t+w_{t+1}+...+w_k) Cov(xt​,xt+k​)=Cov(xt​,xt​+wt+1​+...+wk​)
= C o v ( x t , x t ) + ∑ i = t + 1 k C o v ( x t , w i ) =Cov(x_t,x_t)+\sum^k_{i=t+1}Cov(x_t,w_i) =Cov(xt​,xt​)+i=t+1∑k​Cov(xt​,wi​)
= C o v ( x t , x t ) + 0 = t σ 2 =Cov(x_t,x_t)+0=t\sigma^2 =Cov(xt​,xt​)+0=tσ2

有了自协方差，自相关系数也可以算出来：

ρ k ( t ) = C o v ( x t , x t + k ) V a r ( x t ) V a r ( x t + k ) \rho_k(t)=\frac{Cov(x_t,x_{t+k})}{\sqrt{Var(x_t)}\sqrt{Var(x_{t+k})}} ρk​(t)=Var(xt​) ​Var(xt+k​) ​Cov(xt​,xt+k​)​
= t σ 2 t σ 2 ( t + k ) σ 2 =\frac{t\sigma^2}{\sqrt{t\sigma^2}\sqrt{(t+k)\sigma^2}} =tσ2 ​(t+k)σ2 ​tσ2​
= 1 1 + k t =\frac{1}{\sqrt{1+\frac{k}{t}}} =1+tk​ ​1​

显然，上面的自相关系数和时间间隔有关，所以不具有平稳性。

另外，当k很小时，自相关系数接近于1。如果把股价走势看成随机游走，那么它的当日价格和昨日价格的相关性理论上非常高。

很多人把上述结论看成是一个重要的工具，于是以今日价格作为次日价格的预测值（因为相关性很高嘛）。看似很有道理，但是策略的期望收益大概率为0。

这里要注意的是，当k很小时，自相关系数接近于1，但不是1！

4、随机游走示例

下面我们用随机游走来生成一些序列看看。就用第二节中的均值（-0.000062）和标准差（0.01315）构建一个高斯白噪声。

设定初始值为3000（昨天，上证第44次突破3000点），然后，我们生成8条随机游走曲线如下：

import numpy as np
import pandas as pd
import math

mu = 0  
sigma = 0.01315  
num = 500
offsets = [i for i in range(500)]
fig = plt.figure(figsize=(10, 6))
ax = fig.add_axes([0.2, 0.2, 1.2, 1.2])
colors = ['grey', 'brown', 'red', 'sienna', 'olive', 'green', 'blue', 'purple']

for i in range(0, 8):
    rand_returns = np.random.normal(mu, sigma, num)
    prices = []
    for j in range(500):
        if j == 0:
            current_price = 3000
        else:
            last_price = prices[j-1]
            current_return = rand_returns[j]
            current_price = math.pow(math.e, math.log(last_price)+current_return)
        prices.append(current_price)
    ax.plot(offsets, prices, color=colors[i], linewidth=1.5, linestyle="-", label=r'random walk %d' % i)

plt.legend(loc='upper left', frameon=False)

可以看到，这些数值走势毫无相关性和规律性，在时间推进过程中，波动由于单点白噪声的累积在远端越来越大，但是没有明显的方向性。