以上步骤统一制作
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
S0=100 #初始价格
mean=0.15 #收益率均值
sigma=0.3 #收益率标准差
T=1.0 #时间1年
S1=S0*np.exp(mean*dt+\
sigma*np.sqrt(dt)*np.random.standard_normal(10000))
#模拟股票价格走势的公式
plt.subplot(211)
plt.hist(S1,bins=50)
plt.rcParams['font.sans-serif']=['SimHei'] #显示中文
plt.xlabel('价格')
plt.ylabel('频数')
plt.title('模拟价格的区间价格的个数')
M=250 #1年250个交易日
dt=T/M #1的250个等分时间间隔
I=3 #制作模拟价格走势
S=np.zeros((M+1,I))
S[0]=S0 #定义S[0]=S0
for t in range(1,M+1):
S[t]=S[t-1]*np.exp(mean*dt+\
sigma*np.sqrt(dt)*np.random.standard_normal(I))
#模拟股票价格走势的公式
plt.subplot(212)
plt.plot(S[:,:],lw=1.5)
plt.rcParams['font.sans-serif']=['SimHei'] #显示中文
plt.xlabel('时间')
plt.ylabel('价格')
plt.title('模拟价格走势')