从分段指数分布生成数字

Question

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

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

有人知道以下代码是否正确用于此目的？特别是，矢量＆＃39; S＆＃39;包含间隔＆＃39;与相应危险率相关的上限（载于向量＆＃39; lambda＆＃39;。

    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

我希望问题很明确。