PYTHON绘制强度为0.1的泊松过程500条样本轨道，并标注均值函数

随机过程的作业

# -*- coding: utf-8 -*-
"""
Spyder Danning's Work
-2022-9-21
"""

import matplotlib.pyplot as plt
import numpy as np

MX = [0 for i in range(1001)]
i=0
while i < 500:
     Y = [0] 
     X = [0]
     MX[0]=0
     j=0
     while j < 1000:
         X.append(X[j]+np.random.poisson(lam = 0.1))
         j=j+1
         Y.append(j)
         MX[j]=MX[j]+X[j]
         print(MX[j])
     i=i+1
     plt.plot(Y,X)

j=0
while j < 1001:
    MX[j] = MX[j] / 500
    j=j+1
    #print (MX[i])
plt.plot(Y, MX, 'ro',label='mean', linewidth=20)
plt.show()

最终效果如图所示：