matlab - 从分段指数分布生成数字

Question

我想从分段指数分布中生成一个随机数，它在不同的时间间隔具有不同的危险率。

我的理解是可以应用标准指数分布的无记忆特性。

任何人都知道以下代码是否适用于此目的？特别是，向量“S”包含与相应危险率相关的区间上限（包含在向量“lambda”中。

    function rand_PEXP=rand_PEXP(S,lambda)
        for j=1:numel(S)
             x=exprnd(lambda(j));
             if j==1 && x<S(j)
                 rand_PEXP=x;
                 break
             elseif x<S(j)
                 rand_PEXP=x+S(j-1);
                 break
             end
        end

我希望这个问题很清楚。