raypdf.m

rice信道在matlab下的仿真

close all
clear all
clc

N=1000000;
for Kdb=[0,1,2,4,10,16]
r=rician(Kdb,N);
rms=sqrt(sum(r.^2)/length(r));
step=0.05;
ran=0.001:step:3;
h=hist(r,ran);
fr_appr=h/(step*sum(h));
fr=(ran/0.5).*exp(-ran.^2);
ran=20*log10(ran/rms);
figure(1)
plot(ran,fr_appr)%,ran,fr)
hold on
end

xlabel('amplitude/RMS')
ylabel('PDF')
title('The PDF of Rician Channel with sigma=0.5')
legend('K=0','K=1','K=2','K=4','K=10','K=16');
axis([-30 10 0 2.5])
grid on