Python实现Paralyzable系统模型的代码问题求助

首先得指出你原代码的核心问题：嵌套循环的逻辑没抓住Paralyzable模型的核心——死时间是累加延长的，而不是固定标记一段死时间就完事。原代码的嵌套循环会导致死时间跟踪混乱，甚至错误地提前终止死时间标记。

正确的实现思路（高效单循环） #

我们只需要跟踪一个死时间截止时刻（deadline），遍历每个时间步即可：

• 初始时探测器处于活跃状态，deadline设为-1（表示当前没有死时间限制）
• 对于每个时间点：
  1. 如果当前时间在deadline之后，探测器是活跃的：
     • 若此时有粒子入射（值为1），保留这个探测信号，同时更新deadline为当前时间+死时间d
  2. 如果当前时间在deadline之内，探测器处于死状态：
     • 不管有没有粒子，这个时间点的信号都设为0
     • 若此时有粒子入射，需要延长死时间，把deadline更新为当前时间+死时间d（这是Paralyzable模型的关键）

修正后的代码 #

import numpy as np
np.random.seed(2)
random_set = np.random.poisson(lam=1, size=500)
# 模拟饱和效应：同一时刻最多探测1个粒子
random_set[random_set > 1] = 1
d = 2  # 死时间长度

deadline = -1  # 初始无死时间限制
for i, val in enumerate(random_set):
    if i >= deadline:
        # 探测器活跃状态
        if val == 1:
            # 记录信号，同时设置新的死时间截止时刻
            deadline = i + d
    else:
        # 探测器处于死时间
        random_set[i] = 0
        if val == 1:
            # 新粒子入射，延长死时间
            deadline = i + d

代码逻辑解释 #

• 用deadline变量精准跟踪探测器恢复活跃的时间点，避免了嵌套循环的冗余操作
• 单循环遍历一次数组，时间复杂度是O(n)，比原代码的O(n*d)高效得多，尤其是当d很大时
• 严格遵循Paralyzable模型规则：任何时刻的粒子入射都会重置死时间，不管探测器当前是否处于死状态

你可以运行这段代码，对比原代码的输出，就能看到正确的死时间累加效果了——死时间内的粒子不仅会被屏蔽，还会把后续的死时间往后推。

内容的提问来源于stack exchange，提问作者Nikki