数字通信中的Rice和Nakagami特性研究

%Matlab仿真：

     %1.%%Nakagami源代码

x=0:0.01:3;

omega=1;

m1=3;

m2=2;

m3=1.5;

m4=1;

m5=0.75;

m6=0.5;

p1=exp(-(m1*x.*x)/omega).*x.^(2*m1-1)*(m1/omega)^m1*2/gamma(m1)+eps;

p2=exp(-(m2*x.*x)/omega).*x.^(2*m2-1)*(m2/omega)^m2*2/gamma(m2)+eps;

p3=exp(-(m3*x.*x)/omega).*x.^(2*m3-1)*(m3/omega)^m3*2/gamma(m3)+eps;

p4=exp(-(m4*x.*x)/omega).*x.^(2*m4-1)*(m4/omega)^m4*2/gamma(m4)+eps;

p5=exp(-(m5*x.*x)/omega).*x.^(2*m5-1)*(m5/omega)^m5*2/gamma(m5)+eps;

p6=exp(-(m6*x.*x)/omega).*x.^(2*m6-1)*(m6/omega)^m6*2/gamma(m6)+eps;

figure(1);

plot(x,p1,'r',x,p2,'g',x,p3,'b',x,p4,'k',x,p5,'y',x,p6,'m');

xlabel('x');

ylabel('PDF');

legend('m1=3','m2=2','m3=1.5','m4=1','m5=0.75','m6=0.5');

title('Nakagami m PDF');

hold on;

P1=cumtrapz(x,p1)+eps;

P2=cumtrapz(x,p2)+eps;

P3=cumtrapz(x,p3)+eps;

P4=cumtrapz(x,p4)+eps;

P5=cumtrapz(x,p5)+eps;

P6=cumtrapz(x,p6)+eps;

figure(2);

plot(x,P1,'r',x,P2,'g',x,P3,'b',x,P4,'k',x,P5,'y',x,P6,'m');

xlabel('x');

ylabel('CDF');

legend('m1=3','m2=2','m3=1.5','m4=1','m5=0.75','m6=0.5');

title('Nakagami m CDF' );

hold on;

%2.%%Rice源代码

x=0:0.01:3;

A=1;

k1=10;

k2=1;

k3=0.1;

p1=exp(-(k1+1)/A*(x.^2+A*k1/(k1+1))).*besseli(0,2*x.*sqrt(k1*(k1+1)/A)).*x*2*(k1+1)/A+eps;

p2=exp(-(k2+1)/A*(x.^2+A*k2/(k2+1))).*besseli(0,2*x.*sqrt(k2*(k2+1)/A)).*x*2*(k2+1)/A+eps;

p3=exp(-(k3+1)/A*(x.^2+A*k3/(k3+1))).*besseli(0,2*x.*sqrt(k3*(k3+1)/A)).*x*2*(k3+1)/A+eps;

figure(1);

plot(x,p1,'r',x,p2,'g',x,p3,'b');

legend('k1=10','k2=1','k3=0.1');

xlabel('x');

ylabel('PDF');

title('Rice PDF');

hold on;

P1=cumtrapz(x,p1)+eps;

P2=cumtrapz(x,p2)+eps;

P3=cumtrapz(x,p3)+eps;

figure(2);

plot(x,P1,'r',x,P2,'g',x,P3,'b');

legend('k1=10','k2=1','k3=0.1');

xlabel('x');

ylabel('CDF');

title('Rice CDF');

hold on;