sourcecodes/bnt-master/KPMtools/softeye.m


1
2
3
4
5
6
7
8
9
10
11

function M = softeye(K, p)
% SOFTEYE Make a stochastic matrix with p on the diagonal, and the remaining mass distributed uniformly
% M = softeye(K, p)
%
% M is a K x K matrix.

M = p*eye(K);
q = 1-p;
for i=1:K
  M(i, [1:i-1  i+1:K]) = q/(K-1);
end