python-3.x - 根据 Python 中的给定方程模拟股票价格

Question

如何使用以下等式生成价格时间序列：

p(t) = p0(1+A * sin(ωt +0.5η(t)))

其中t在1000 个时间步长范围内从0到1 ， p0 = 100，A = 0.1和ω = 100。η(t)是具有零均值和单位方差的iid 高斯随机变量序列。

我已经使用如下代码生成价格，但似乎不是必需的。所以我需要社区的帮助。提前致谢。

from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
mu = 0
sigma = 1
np.random.seed(2020)

dist = norm(loc = mu,scale=sigma)
sample = dist.rvs(size = 1000)
stock_price = np.exp(sample.cumsum())
print(stock_price)
plt.plot(stock_price)
plt.xlabel("Day")
plt.ylabel("Price")
plt.title("Simulated Stock price")
plt.show()

Answer 1

假设我没有错过任何东西，这应该可以解决问题

import numpy as np
import matplotlib.pyplot as plt

n_t = np.random.normal(0, 1, 1000)
t = np.arange(0, 1, 1/1000)

p_0, A, w = 100, 0.1, 100
ts = p_0 * (1 + A * np.sin(w * t + 0.5 * n_t))

plt.plot(t, ts)
plt.xlabel("Day")
plt.ylabel("Price")
plt.show()

这给出了情节

Answer 2

我的试用，不知道是否正确，欢迎给我一些意见。

import numpy as np
import math
np.random.seed(2020)

mu = 0
sigma = 1
dt = 0.01
p0 = 100
A = 0.1
w = 100
N = 1000
for t in np.linspace(0, 1, 1000):
    X = np.random.normal(mu * dt, sigma* np.sqrt(dt), N)
    X = np.cumsum(X)
    pt = p0 * (1+ A*np.sin(w*t + 0.5*X))
#     print(pt)
plt.plot(pt)
plt.xlabel("Day")
plt.ylabel("Price")
plt.title("Simulated Stock price")
plt.show()

出去：