【移动通信】瑞利衰落通道上的差分MPSK（Matlab代码实现）

👨‍🎓个人主页：研学社的博客 

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

时变瑞利衰落通道上的差分MPSK。通信是点对点的。源通过无线瑞利信道与目的地通信。通道是时变的。移动通信就是这种情况。自相关值决定了通道速度的大小。

📚2 运行结果

 

部分代码：

% differential decoder for time-varying channels
function vh= diff_decod_2Antenna(y1,y2,M)
N=length(y1);
if M==2
    V=[-1,1];   
elseif M==4
    V=[1,1i,-1,-1i];
end

for k=3:2:N
    for row=1:length(V)
        v_km1=V(row);
        for col=1:length(V)
            v_k=V(col);
            costF1=(y1(k)-v_k*y1(k-1))*(y2(k-1)-v_km1*y2(k-2));
            costF2=(y1(k-1)-v_km1*y1(k-2))*(y2(k)-v_k*y2(k-1));
            costF(row,col)=costF1-costF2;
        end
    end
    [nLl_rm,r_min]=min(abs(costF));
    [~,col_min]=min(nLl_rm);
    vh(k-2)=V(r_min(col_min));
    vh(k-1)=V(col_min);
end

end

% total symbol error probability for D-AF with R-rleys

function Ps=Ps_multi_node(P,M,fsr,frd,R)
% inputs
% P : source power
% M : M-PSK constelastion
% fsr & frd : normalized doppler for source-relay and relay-destination
% R : number of relays
% output
% Ps : total symbol error probability

% amplification factor
global A2;

% autocorelations of individual and cascaded channel
alfa1=besselj(0,2*pi*fsr*(R+1));
alfa2=besselj(0,2*pi*frd*(R+1));
alfa=alfa1*alfa2;

% minimum distance in M-PSK
dmin=2*sin(pi/M);

% compute theoritical SER
% c=alfa^2*A2.*P;
% b=4;
% a=2*(1-alfa^2)*A2.*P+4*A2;
% beta=2*b./(c*dmin^2+2*a);
% d=c*dmin^2./(c*dmin^2+2*a);
%Ex_F_gama=d.*ex_fgama(beta,1);

temp=4*(1-alfa^2)+alfa^2*dmin^2;
K=P.*A2*temp+8*A2;
H=8./K;
Ex_F_gama=alfa^2*dmin^2*H.*P.*A2.*(1-exp(H).*H.*expint(H))/8;
Ex=1-Ex_F_gama;
sumEx=0;
for l=0:R-1
   sumEx=sumEx+nchoosek(2*l,l)*(Ex/4).^l;
end
Ps=.5*(1-sqrt(Ex_F_gama).*sumEx);

% nearest neighbour approximation when M>2
if M>2
    Ps=2*Ps;
end

end
%% subfunction integral x^R/(x+beta) exp(-x) dx, x=0 to inf
function out=ex_fgama(beta,R)

temp=0;
for k=1:R
    temp=temp+factorial(R-1)*(-beta).^(R-k);
end
out=(temp+(-1)^R*beta.^R.*exp(beta).*expint(beta))/factorial(R-1);

end

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]Differential MPSK over Rayleigh fading channels

🌈4 Matlab代码实现