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

質問が明確であることを願っています。