タイプ のシステムの解決に関与していますAx = b。ここで、A は正方スパース行列、x は未知数のベクトル (計算する必要があります)、b は最後の要素を除くすべてゼロのベクトルです。は 1 です。行列 A の最後の行は正規化に使用されるため、1 で満たされます。
このシステムの解は確率であるため、条件0<x(i)<1を尊重する必要があります。
システムを解決するために、Matlab コマンドx = A \ b;が使用されます。
この方法はうまく機能しているように見えますが、ベクトル x にも負の値が含まれる特殊なケースがあります。マトリックス A の任意の要素に非常に小さな値 (10^-6) を追加すると、解像度が条件を満たすように戻ります。
私は数学者ではないので、コードの問題なのか、解がすべて 0 と 1 の間にあることを保証するために行列 A がいくつかのプロパティを尊重する必要があるのかどうかはわかりません。