poisson过程——随机模拟(Python和R实现)

Python实现

exponential()使用,自动poisson过程实现。

python 复制代码
import numpy as np
import matplotlib.pyplot as plt

# Parameters
lambda_rate = 5  # rate parameter (events per time unit)
T = 10  # total time

# Generate Poisson process
times = np.random.exponential(1/lambda_rate, int(lambda_rate * T * 1))
arrival_times = np.cumsum(times)
arrival_times = arrival_times[arrival_times < T]

# Plot Poisson process
plt.step(np.concatenate([[0], arrival_times]), np.arange(len(arrival_times) + 1), where='post')
plt.xlabel('Time')
plt.xticks(np.arange(0, T + 1, 1))
plt.ylabel('Number of events')
plt.title('Poisson Process')
plt.show()

def poisson函数使用。

python 复制代码
#正常显示中文标签
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False
#%%
import random
import numpy as np
import matplotlib.pyplot as plt

def poisson(T,lambdas):
    t = 0
    i = 0
    s = [t]
    while True:
        u = random.expovariate(lambdas)
        t = t + u
        if t <= T:
            i = i + 1#计数
            s.append(np.round(t,2))#记录发生事件的时间,保留两位小数
        else:
            break
    print(f'{T}时间内事件发生的次数:',i,'\n 各个事件发生的时刻:',s)

    #绘制样本轨迹
    n = i + 1
    x = np.arange(n)
    plt.step(s,x, where='post')#绘制阶梯图
    plt.xticks(s)#设置x轴刻度
    plt.xlabel('时间t')
    plt.ylabel('N(t)')
    plt.title('Poisson Process')
    plt.show()
    return s
#生成图像
poisson(4,1)#参数分别为T和lambda,T为时间长度,lambda为单位时间内事件发生的次数

R语言实现

R 复制代码
#R语言实现poisson过程
# Parameters
lambda_rate <- 5  # rate parameter (events per time unit)
T <- 10  # total time

# Generate Poisson process
set.seed(123)  # for reproducibility
times <- rexp(n = ceiling(lambda_rate * T * 1.5), rate = lambda_rate)
arrival_times <- cumsum(times)
arrival_times <- arrival_times[arrival_times < T]

# Plot Poisson process
plot(stepfun(arrival_times, 0:length(arrival_times)), do.points = FALSE, xlab = "Time", ylab = "Number of events", main = "Poisson Process")
相关推荐
起床气2335 分钟前
C++海战棋开发日记(序)
开发语言·c++
APItesterCris17 分钟前
TypeScript 与淘宝 API:构建类型安全的商品数据查询前端 / Node.js 服务
开发语言·php
ftpeak21 分钟前
《Cargo 参考手册》第二十一章:Cargo 包命令
开发语言·rust
可触的未来,发芽的智生22 分钟前
触摸未来2025-10-18:生成文字的小宇宙矩阵溯源
人工智能·python·神经网络·程序人生·自然语言处理
_码力全开_23 分钟前
P1005 [NOIP 2007 提高组] 矩阵取数游戏
java·c语言·c++·python·算法·矩阵·go
陈一Tender27 分钟前
JavaWeb后端实战(登录认证 & 令牌技术 & 拦截器 & 过滤器)
java·开发语言·spring boot·mysql
Camel卡蒙27 分钟前
红黑树详细介绍(五大规则、保持平衡操作、Java实现)
java·开发语言·算法
jerryinwuhan31 分钟前
机器人模拟器(python)
开发语言·python·机器人
AhriProGramming39 分钟前
Flask-SQLAlchemy精读-双语精选文章
python·算法·flask
列兵阿甘42 分钟前
知微传感Dkam系列3D相机SDK例程篇:Python获取内外参
python·数码相机·3d