%  Rayleigh distribution Figure 6.6

sig1=1;
sig2=2;
r=[0:.01:6];
r1=r.*exp(-r.^2/(2*sig1^2))/sig1^2;
r2=r.*exp(-r.^2/(2*sig2^2))/sig2^2;
figure;
plot(r,r1,'k',r,r2,'k:')
axis([0 6 0 1])
xlabel('Envelope r (mV)')
ylabel('p(r)')
legend('\sigma_1=.001','\sigma_2=.002')