pdf.asv

基于调和平均数方法分析瑞利信道多点两跳中继合作分集性能

% tingting yang
% 2008-11-12
% multi-node Dual-Hop
% MPSK; SER
% nonregenerative relays
% 最大比合并
% PDF expressions
clc
clear
m=4;  % the number of relay terminals
M=4;    % the orders of modulation
g_MPSK=(sin(pi/M))^2;
r=0:0.2:4;
rh=[0.5 1 2 4];
rg=[0.5 1 2 4];
pdf=linspace(0,0,21);
%no realistic sense here, the aim is just to make the length of pdf the same as
%that of r.
%for j=1:length(r);

    for i=1:m
        va1=-(r)*(rh(i)^(-1)+rg(i)^(-1));
        f=exp(va1);
        va2=2.*r.*f;
        va3=rh(i)*rg(i);
        va4=rh(i)+rg(i);
        va5=va4*besselk(1,2*r/sqrt(va3));
        va6=2*K0(r,va3);
        va7=va2/va3;
        sum2=va7.*(va5+va6);
        
        pdf=pdf+sum2;
    end 
 m=4;  % the number of relay terminals
M=4;    % the orders of modulation
g_MPSK=(sin(pi/M))^2;
r=0:0.2:4;
rh2=[0.1 0.3 1 2.5];
rg2=[2 2.5 2.7 4];
pdf2=linspace(0,0,21);
%no realistic sense here, the aim is just to make the length of pdf the same as
%that of r.
%for j=1:length(r);

    for i=1:m
        va1=-(r)*(rh2(i)^(-1)+rg2(i)^(-1));
        f=exp(va1);
        va2=2.*r.*f;
        va3=rh2(i)*rg2(i);
        va4=rh2(i)+rg2(i);
        va5=va4*besselk(1,2*r/sqrt(va3));
        va6=2*K0(r,va3);
        va7=va2/va3;
        sum2=va7.*(va5+va6);
        
        pdf2=pdf2+sum2;
    end    
   m=10;  % the number of relay terminals
M=8;    % the orders of modulation
g_MPSK=(sin(pi/M))^2;
r=0:0.2:4;
rh2=[0.1 0.3 0.5 0.7 0.9 1.2 1.4 1.6 1.8 2.5];
rg2=[2.1 2.3 2.5 2.7 2.9 3.2 3.4 3.6  2.5 2.7];
pdf3=linspace(0,0,21);
%no realistic sense here, the aim is just to make the length of pdf the same as
%that of r.
%for j=1:length(r);

    for i=1:m
        va1=-(r)*(rh2(i)^(-1)+rg2(i)^(-1));
        f=exp(va1);
        va2=2.*r.*f;
        va3=rh2(i)*rg2(i);
        va4=rh2(i)+rg2(i);
        va5=va4*besselk(1,2*r/sqrt(va3));
        va6=2*K0(r,va3);
        va7=va2/va3;
        sum2=va7.*(va5+va6);
        
        pdf3=pdf3+sum2;
    end    
    
       m=10;  % the number of relay terminals
M=8;    % the orders of modulation
g_MPSK=(sin(pi/M))^2;
r=0:0.2:4;
rg3=[0.1 0.3 0.5 0.7 0.9 1.2 1.4 1.6 1.8 2.5];
rh3=[2.1 2.3 2.5 2.7 2.9 3.2 3.4 3.6  2.5 2.7];
pdf4=linspace(0,0,21);
%no realistic sense here, the aim is just to make the length of pdf the same as
%that of r.
%for j=1:length(r);

    for i=1:m
        va1=-(r)*(rh3(i)^(-1)+rg3(i)^(-1));
        f=exp(va1);
        va2=2.*r.*f;
        va3=rh3(i)*rg3(i);
        va4=rh3(i)+rg3(i);
        va5=va4*besselk(1,2*r/sqrt(va3));
        va6=2*K0(r,va3);
        va7=va2/va3;
        sum2=va7.*(va5+va6);
        
        pdf4=pdf4+sum2;
    end 
    
hold on
plot(r,pdf)
%mesh(r,pdf)
%grid on
plot(r,pdf2,':')
plot(r,pdf3,'o')
plot(r,pdf4,'x')
xlabel('')