【采用差分二相移键控DBPSK的D-AF中继网络中】选择合并技术在时变信道上的差分放大转发中继性能】研究（Matlab代码实现）

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

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

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

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

目录

 ⛳️赠与读者

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、文章下载

---

 ⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑。哲学是科学之母，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

     或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

摘要：对于时变信道上的差分放大转发中的目的地进行选择合并（SC）是具有吸引力的，因为与半最大比合并（semi-MRC）相比，它不需要通道状态信息，同时提供接近的性能。最近有关SC方案的性能分析仅针对缓慢衰落信道的情况进行了报告。本文提供了SC方案在一般情况下的时变瑞利衰落信道上的精确平均比特误码率（BER），并且当DBPSK调制与目的地的非相干检测一起使用时。所呈现的分析经过了各种衰落场景的模拟结果的彻底验证。结果表明，系统的性能与信道的自相关值有关。还表明，SC方法的性能非常接近于半MRC方法，并且在高信噪比区域存在误码率底线在这两种方法中是不可避免的。对于SC方法得到的BER分析也可用于近似MRC方法的BER性能，后者在时变信道中的确切分析评估似乎很困难。

合作通信现在已成为一个成熟的研究课题。目前，一种特殊类型的协同通信（借助一个中继器）已经在3GPP LTE技术中标准化，以解决蜂窝网络的覆盖问题，并且预期LTE-高级版本将包含协同中继特性，以克服容量和干扰等其他限制[1]。合作中继网络在无线局域网、车对车通信和无线传感器网络中的应用也在[2]-[5]及其中引用的文献中有讨论。

在合作通信中，网络中的一个用户充当中继，接收来自源节点的信号，处理后并重新广播给目的节点。通过这种方式，除了源节点到目的节点的直接链路之外，还可以通过中继构建额外的链路，从而系统整体的空间多样性得到增加。根据中继所采用的信号处理策略，中继网络通常被分为解码转发（DF）和放大转发（AF）两大类[6]。

在采用差分二相移键控（DBPSK）的D-AF中继网络中，研究了时变瑞利衰落信道中目的端接收到的信号的选择合并。得益于差分编码和选择合并器，目的端进行信息检测时无需信道状态信息。系统的确切误比特率被推导出来。在不同衰落速率和信道方差条件下的仿真结果验证了分析，并显示选择合并器的性能非常接近于更复杂的半最大比合并方法（该方法需要目的端所有信道的二阶统计信息）。

📚2 运行结果

部分代码：

M=2;% MPSK modulation
Ns=1E5;% number of symbols
N0=1;% noise variance
Ptot_dB=0:5:30;% SNR scan
Ptot=10.^(Ptot_dB/10)*N0;

% ------------------------------------% 
%choose channel scenario 1,2,3
ch_scn=1;
%choose fading scenario, 1,2,3 
fade_scn=1;

% ------------- channel variances
vsig_sd=[1,1,1];
vsig_sr=[1,10,1];
vsig_rd=[1,1,10];
sig_sd=vsig_sd(ch_scn);
sig_sr=vsig_sr(ch_scn);
sig_rd=vsig_rd(ch_scn);

% distance between two channel uses
ch_dis=1; % 1 for block-by-block and R+1 for symbol-by-symbol

% normalized Dopplers frequencies
vfsd=[.001,.02,.05];
vfsr=[.001,.02,.01];
vfrd=[.001,.001,.05];

% scenario number
fsd=vfsd(fade_scn);
fsr=vfsr(fade_scn);
frd=vfrd(fade_scn);

% auto-correlation values
alfa0=besselj(0,2*pi*fsd*ch_dis);
alfa1=besselj(0,2*pi*fsr*ch_dis);
alfa2=besselj(0,2*pi*frd*ch_dis);
alfa=alfa1*alfa2;

% power allocation
vpaf=[.66,.54,.8];
paf=vpaf(ch_scn);

% source power
P0=paf*Ptot;
% relay power
P1=(1-paf)*Ptot;
% amplification factor
Ai2= P1./(P0*sig_sr+N0);

% this loop scans the SNR range
for ind=1:length(Ptot)

nbits=0;%total number of info sent
err_mrc=0;% error counter MRC method
err_sc=0;% error counter Selection Combining

clc
Ptot_dB(ind)

% this loop keeps going to get a certain amount of errors
ERR_TH=100;
while err_mrc<ERR_TH || err_sc<ERR_TH 

% info bits
xb=bits(log2(M)*Ns);

%binary  to MPSK
v=bin2mpsk(xb,M);

% DPSK modulation
s=diff_encoder(v);
Nd=length(s);

% S-D channel
hsd=sqrt(sig_sd)*flat_cos(Nd,fsd,ch_dis);

% S-R and R-D channels
hsr=sqrt(sig_sr)*flat_cos(Nd,fsr,ch_dis);
hrd=sqrt(sig_rd)*flat_cos(Nd,frd,ch_dis);

% AWGN noise CN(0,N0)
z_sd=cxn(Nd,N0);
z_sr=cxn(Nd,N0);
z_rd=cxn(Nd,N0);

% received signals at destination
y_sd=sqrt(P0(ind))*hsd.*s+z_sd;
y_sr=sqrt(P0(ind))*hsr.*s+z_sr;
y_rd=sqrt(Ai2(ind))*hrd.*y_sr+z_rd;

% decision variables
zeta_sd=real(diff_detector(y_sd));
zeta_rd=real(diff_detector(y_rd));

% MRC combining
a0_mrc=1/N0;
a1_mrc=1/((1+Ai2(ind)*sig_rd)*N0);
vh_mrc=a0_mrc*zeta_sd+a1_mrc*zeta_rd;

% selection combining
vh_sc=zeros(1,length(zeta_sd));
 for k=1:length(zeta_sd)
     if abs(zeta_sd(k))> abs(zeta_rd(k))
         vh_sc(k)=zeta_sd(k);
     else
         vh_sc(k)=zeta_rd(k);
     end
 end

% binary detection
bh_mrc=mpsk2bin(vh_mrc,M);
bh_sc=mpsk2bin(vh_sc,M);

% error count
err_mrc=err_mrc+sum(abs(xb-bh_mrc));
err_sc=err_sc+sum(abs(xb-bh_sc));
nbits=log2(M)*Ns+nbits;

end

% compute practical BER 
BER_mrc(ind)=err_mrc/nbits;
BER_sc(ind)=err_sc/nbits;

end

% theoretical BER, numerical average
Pb_sc=Pb_sc_dpsk(P0,N0,Ai2,M,fsd,fsr,frd,ch_dis,sig_sr,sig_rd,sig_sd);

% BER, analytic average
%Pb_sc2=Pb_sc_dpsk2(P0,N0,Ai2,M,fsd,fsr,frd,ch_dis,sig_sr,sig_rd,sig_sd);

% error floor
Pbf=Pe_floor(paf,fsd,fsr,frd,M,ch_dis,sig_sd,10,sig_rd)*ones(1,length(P0));
%% plot 
lcm=['b-+';'r-X';'k:>';'g-.';'y-*'];
mrk_sz=8;

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

🌈4 Matlab代码、文章下载

资料获取，更多粉丝福利，MATLAB|Simulink|Python资源获取