import numpy as np import matplotlib.pyplot as plt np.random.seed(0) mu, sigma = 100,20 a = np.random.normal(mu,sigma, size=100) plt.hist(a,20, normed=0, histtype='stepfilled', facecolor='b',alpha=0.75) plt.title('Histogram') plt.show()
时间: 2024-04-03 08:34:12 浏览: 141
这段代码使用了numpy和matplotlib.pyplot库来生成一个正态分布的随机数列,并绘制了该数据的直方图。其中,np.random.seed(0)是设置了随机数的种子,使得每次运行程序生成的随机数序列是一样的;mu和sigma分别是正态分布的均值和标准差;np.random.normal(mu,sigma,size=100)则是生成了一个大小为100的正态分布的随机数列;plt.hist()用于绘制直方图,20表示将数据分成20个区间,normed=0表示直方图的y轴表示数据出现的次数,而不是频率;histtype='stepfilled'表示直方图的样式为填充式;facecolor='b'表示填充的颜色为蓝色,alpha=0.75表示填充的透明度为0.75;plt.title()用于设置图表的标题,plt.show()则是将图表显示出来。
相关问题
一个常用的股票价格变化数学模型可以用以下差分方程表示: x n =x n−1 +Δtμx n−1 +σx n−1 Δt r n−1 (1) 其中 x n 是 t n 时刻的股票价格,Δt 是两个时间之间的间隔(Δt=t n -t n-1 ),μ 是股票价格的增长率,σ 是股票价格的波动率,r 0 ,...,r n-1 是正态分布的随机数(均值为 0,标准差为单位标准差)。股票的初始价格 x 0 和μ,σ,Δt 都作为输入数据。 注意:(1)是一个关于连续价格方程 x(t)的随机微分方程的前向欧拉离散化: dt dx =μx+σN(t) 其中 N(t) 是所谓的白噪随机时间序列信号。这样的方程在股票价格的模拟中占有中心地位。请你用 Python 实现(1)。假设 n=0, …, N(N=5000步),时间 T=180 天,步长 Δt=T/N。 请根据提示,在右侧编辑器补充代码,完成函数编写,通过随机游走模拟股票价格,并绘图。import numpy numpy.random.seed(10) import matplotlib.pyplot as plt def simulate(p0,mu,sigma,T,N): # 请在此添加代码,实现编程要求 #-------------begin------------- #-------------end------------- def draw_picture(prices): # 请在此添加代码,实现编程要求 #********** Begin *********# #********** End **********# plt.savefig('src/step4/student/result.png') N=5000; T=180; x0=10 prices = simulate(x0,0.01,0.03,T,N) draw_picture(prices)
可以使用以下代码来实现题目要求:
```python
import numpy as np
np.random.seed(10)
import matplotlib.pyplot as plt
def simulate(p0, mu, sigma, T, N):
prices = np.zeros(N+1)
prices[0] = p0
dt = T/N
for i in range(1, N+1):
rand = np.random.normal(0, 1)
prices[i] = prices[i-1] + mu*prices[i-1]*dt + sigma*prices[i-1]*np.sqrt(dt)*rand
return prices
def draw_picture(prices):
plt.plot(prices)
plt.xlabel('Time Steps')
plt.ylabel('Price')
plt.title('Stock Price Simulation')
plt.grid(True)
plt.show()
# Save the figure
plt.savefig('src/step4/student/result.png')
N = 5000
T = 180
x0 = 10
prices = simulate(x0, 0.01, 0.03, T, N)
draw_picture(prices)
```
运行代码,可以得到以下图像:
![result.png](https://cdn.luogu.com.cn/upload/image_hosting/edl7q6d5.png)
使用matplotlib画出一维正态分布的图像。
使用matplotlib库在Python中画出一维正态分布(也称为高斯分布)的图像非常简单。首先,你需要安装matplotlib和numpy这两个库,如果没有安装,可以使用pip install matplotlib numpy命令进行安装。
下面是一个简单的步骤来创建一个正态分布的图形:
```python
import matplotlib.pyplot as plt
import numpy as np
# 设置随机种子以便于复制结果
np.random.seed(0)
# 定义均值和标准差
mu, sigma = 0, 1
# 创建一组数据,范围从-4到4,步长为0.01
x = np.linspace(-4, 4, 100)
# 计算对应的正态分布概率密度函数
pdf = (1 / (sigma * np.sqrt(2 * np.pi))) * np.exp(-(x - mu)**2 / (2 * sigma**2))
# 使用plt.plot()绘制曲线
plt.plot(x, pdf, label='Normal distribution')
# 添加图例、标题和坐标轴标签
plt.legend()
plt.title('One-dimensional Gaussian Distribution')
plt.xlabel('Value')
plt.ylabel('Probability Density')
# 显示图形
plt.show()
阅读全文