MA過程を描いてみる。 import matplotlib.pyplot as plt import numpy as np def ma(n, mu, theta, sigma): """ 1次の移動平均過程を返す """ epsilon = (np.random.rand(n) * 2 - 1) * sigma y = [] for i in range(n): if i == 0: yi = mu + epsilon[i] el…