% generate a rician distribution for 4 K values

sig2=1;
r=0:.005:12;
figure;
% Create Rician distribution
K=[-10 10 15];
for j=1:3
    A=sqrt(2*sig2*10^(K(j)/10))
    ex=exp(-(r.^2+A^2)/(2*sig2));
    I=besseli(0,A*r/sig2);
    p=r.*ex.*I/sig2;
    plot(r,p,'k')
    hold on
end
xlabel('r')
ylabel('p(r)')
legend('K=-10dB','K=10dB','K=15dB')