CN103281111A - 基于选择合并的解码转发协作系统的信号智能转发方法 - Google Patents

基于选择合并的解码转发协作系统的信号智能转发方法 Download PDF

Info

Publication number
CN103281111A
CN103281111A CN2013102164124A CN201310216412A CN103281111A CN 103281111 A CN103281111 A CN 103281111A CN 2013102164124 A CN2013102164124 A CN 2013102164124A CN 201310216412 A CN201310216412 A CN 201310216412A CN 103281111 A CN103281111 A CN 103281111A
Authority
CN
China
Prior art keywords
gamma
alpha
beta
integral
infin
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2013102164124A
Other languages
English (en)
Other versions
CN103281111B (zh
Inventor
赵洪林
赵大伟
张佳岩
马永奎
张中兆
贾敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201310216412.4A priority Critical patent/CN103281111B/zh
Publication of CN103281111A publication Critical patent/CN103281111A/zh
Application granted granted Critical
Publication of CN103281111B publication Critical patent/CN103281111B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

基于选择合并的解码转发协作系统的信号智能转发方法,涉及一种协作通信中的信号转发方法。为了解决对目的节点采用最大比合并器的信号转发方法极大地增加了目的节点接收机的复杂度的问题。所述方法所提出的智能中继适用于源节点S和中继节点R采用不同阶数的调制方式以提高频谱利用率的情况,能够充分利用瞬时的信道状态信息或平均的信道状态信息,采用功率分配、功率控制或基于信噪比的选择转发方法使最小化系统的误比特率,降低错误传播的影响,有效地提高解码转发协作传输的可靠性。它用于协作通信中的信号转发。

Description

基于选择合并的解码转发协作系统的信号智能转发方法
技术领域
本发明涉及一种协作通信中的信号转发方法,特别涉及一种基于选择合并的解码转发协作系统的信号智能转发方法。
背景技术
协作通信利用中继辅助源节点发送信息,不仅扩大了网络的覆盖范围,并且形成了虚拟的多天线系统,获得了空间分集增益。解码转发属于协作通信中的数字中继方案,中继在收到源节点发送的信号后,首先对信号进行解调、解码,之后再重新调制、转发。受源-中继信道非理想性的影响,中继的解码可能存在错误,固定的解码转发方案会导致错误传播而无法保证满分集增益。
针对这个问题,人们提出了基于信噪比门限的选择转发(即中继仅在源-中继的瞬时信噪比超过预设门限时才转发信号,否则不参与协作)以及智能中继(即中继能够根据信道状态信息连续地调整中继的发射功率),有效地改善了错误传播的影响。然而,目前的研究大多针对的是目的节点采用最大比合并器的情况,合并器需要知道各支路信道系数的幅度和相位信息,极大地增加了目的节点接收机的复杂度。
发明内容
本发明的目的是为了解决对目的节点采用最大比合并器的信号转发方法极大地增加了目的节点接收机的复杂度的问题,本发明提供一种基于选择合并的解码转发协作系统的信号智能转发方法。
基于选择合并的解码转发协作系统的信号智能转发方法,所述方法是基于单中继的解码转发协作系统和采用功率分配的智能中继方法实现的,所述系统包括源节点S、中继节点R和目的节点D,上述节点仅装配单根天线且工作于半双工模式,S→R支路、R→D支路和S→D支路均为符合独立非同分布的Rayleigh块衰落信道,
它包括如下步骤:
步骤A1:利用源节点S和中继节点R发送训练序列的方法,中继节点R获得信道状态信息T1(ΓSD,γSR,γRD)或T2(ΓSD,γSR,ΓRD);其中,ΓSD为S→D支路的平均信噪比,γSR为S→R支路的瞬时信噪比,γRD为R→D支路的瞬时信噪比,ΓRD为R→D支路的平均信噪比;
步骤A2:采用功率分配的智能中继方法,求得最优的功率修正系数αS和αR
步骤A3:源节点S和中继节点R根据功率修正系数αS和αR调整发射功率传输信号;
步骤A4:目的节点D收到各支路信号后,采用选择合并器输出合并后的信号:
Figure BDA00003295361900021
Figure BDA00003295361900022
Figure BDA00003295361900023
MS为源节点S的调制阶数,MR为中继的调制阶数,γSD为S→D支路的瞬时信噪比;调制方式采用基于格雷码的方形M-QAM调制;
所述采用功率分配的智能中继方法为:
当中继节点R获得信道状态信息为T1(ΓSD,γSR,γRD)时,此时的误比特率PT1S,αR)为:
P T 1 ( α S , α R ) = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD + ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD
其中,PSD为S→D支路的瞬时误比特率,PSD=PMSSγSD)
S→R→D支路的瞬时误比特率PSRD=PMSSγSR)+(1-2PMSSγSR))PMRRγRD),
P M ( γ ) = c M Q ( 2 d M 2 γ ) , Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M s = 2 M s - 1 M s log 2 M s ;
以中间变量
Figure BDA00003295361900028
表示 ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD ,
P T 1 1 = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD
= c MS Q ( 2 d MS 2 β α R γ RD ) exp ( - β α R γ RD α S Γ SD ) - c MS d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD Q ( 2 ( 1 + d MS 2 α S Γ SD ) β α R γ RD α S Γ SD ) ,
以中间变量
Figure BDA000032953619000212
表示 ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD ,
P T 1 2 = P SRD ∫ γ SD = 0 β α R γ RD / α S f ( γ SD ) d γ SD = P SRD ( 1 - exp ( - β α R γ RD α S Γ SD ) ) ;
使误比特率PT1S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS∈(0,2],αR∈(0,2],αSR=2;
当中继节点R获得信道状态信息为T2(ΓSD,γSR,ΓRD)时,此时的误比特率PT2S,αR)为:
P T 2 ( α S , α R ) = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD + ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
以中间变量
Figure BDA00003295361900032
表示 ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD ,
P T 2 1 = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD
= 1 2 c MS ( 1 - d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD )
- 1 2 c MS β α R Γ RD α S Γ SD + β α R Γ RD ( 1 - d MS 2 β α R α S Γ SD Γ RD α S Γ SD + β α R Γ RD + d MS 2 β α S α R Γ SD Γ RD )
以中间变量
Figure BDA00003295361900037
表示 ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD ,
P T 2 2 = ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
= c MS Q ( 2 d MS 2 α S γ SR ) ( β α R Γ RD α S Γ SD + β α R Γ RD ) + ( 1 - 2 c MS Q ( 2 d MS 2 α S γ SR ) ) ×
( 1 2 c MR ( 1 - d MR 2 α R Γ RD 1 + d MR 2 α R Γ RD ) - 1 2 c MS α S Γ SD α S Γ SD + β α R Γ RD ( 1 - d MR 2 α S α R Γ SD Γ RD α S Γ SD + β α R Γ RD + d MR 2 α S α R Γ SD Γ RD ) )
Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M R = 2 M R - 1 M R log 2 M R , c M s = 2 M s - 1 M s log 2 M s ;
使误比特率PT2S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS∈(0,2],αR∈(0,2],αSR=2。
本发明还提供一种基于选择合并的解码转发协作系统的信号智能转发方法,所述方法是基于单中继的解码转发协作系统和采用功率控制的智能中继方法实现的,所述系统包括源节点S、中继节点R和目的节点D,上述节点仅装配单根天线且工作于半双工模式,S→R支路、R→D支路和S→D支路均为符合独立非同分布的Rayleigh块衰落信道,
它包括如下步骤:
步骤B1:利用源节点S和中继节点R发送训练序列的方法,中继节点R获得信道状态信息T1(ΓSD,γSR,γRD)或T2(ΓSD,γSR,ΓRD);其中,ΓSD为S→D支路的平均信噪比,γSR为S→R支路的瞬时信噪比,γRD为R→D支路的瞬时信噪比,ΓRD为R→D支路的平均信噪比;
步骤B2:采用功率控制的智能中继方法,求得最优的功率修正系数αR
步骤B3:源节点S的发射功率不变,中继节点R根据功率修正系数αR调整发射功率传输信号;
步骤B4:目的节点D收到各支路信号后,采用选择合并器输出合并后的信号:
Figure BDA00003295361900041
Figure BDA00003295361900043
MS为源节点S的调制阶数,MR为中继的调制阶数,γSD为S→D支路的瞬时信噪比;调制方式采用基于格雷码的方形M-QAM调制;
所述采用功率控制的智能中继方法为:
当中继节点R获得信道状态信息为T1(ΓSD,γSR,γRD)时,此时的误比特率PT1S,αR)为:
P T 1 ( α S , α R ) = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD + ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD 其中,PSD为S→D支路的瞬时误比特率,PSD=PMSSγSD)
S→R→D支路的瞬时误比特率PSRD=PMSSγSR)+(1-2PMSSγSR))PMRRγRD),
P M ( γ ) = c M Q ( 2 d M 2 γ ) , Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M s = 2 M s - 1 M s log 2 M s ;
以中间变量
Figure BDA00003295361900048
表示 ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD ,
P T 1 1 = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD
= c MS Q ( 2 d MS 2 β α R γ RD ) exp ( - β α R γ RD α S Γ SD ) - c MS d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD Q ( 2 ( 1 + d MS 2 α S Γ SD ) β α R γ RD α S Γ SD ) ,
以中间变量
Figure BDA000032953619000412
表示 ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD ,
P T 1 2 = P SRD ∫ γ SD = 0 β α R γ RD / α S f ( γ SD ) d γ SD = P SRD ( 1 - exp ( - β α R γ RD α S Γ SD ) ) ;
使误比特率PT1S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS=1,αR∈[0,1];
当中继节点R获得信道状态信息为T2(ΓSD,γSR,ΓRD)时,此时的误比特率PT2S,αR)为:
P T 2 ( α S , α R ) = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD + ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
以中间变量
Figure BDA00003295361900052
表示 ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD ,
P T 2 1 = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD
= 1 2 c MS ( 1 - d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD ) ,
- 1 2 c MS β α R Γ RD α S Γ SD + β α R Γ RD ( 1 - d MS 2 β α R α S Γ SD Γ RD α S Γ SD + β α R Γ RD + d MS 2 β α S α R Γ SD Γ RD )
以中间变量
Figure BDA00003295361900057
表示 ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
P T 2 2 = ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
= c MS Q ( 2 d MS 2 α S γ SR ) ( β α R Γ RD α S Γ SD + β α R Γ RD ) + ( 1 - 2 c MS Q ( 2 d MS 2 α S γ SR ) ) ×
( 1 2 c MR ( 1 - d MR 2 α R Γ RD 1 + d MR 2 α R Γ RD ) - 1 2 c MS α S Γ SD α S Γ SD + β α R Γ RD ( 1 - d MR 2 α S α R Γ SD Γ RD α S Γ SD + β α R Γ RD + d MR 2 α S α R Γ SD Γ RD ) )
Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M R = 2 M R - 1 M R log 2 M R , c M s = 2 M s - 1 M s log 2 M s ;
使误比特率PT2S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS=1,αR∈[0,1]。
本发明还提供一种基于选择合并的解码转发协作系统的信号智能转发方法,所述方法是基于单中继的解码转发协作系统和采用于基于信噪比的选择转发的方法实现的,所述系统包括源节点S、中继节点R和目的节点D,上述节点仅装配单根天线且工作于半双工模式,S→R支路、R→D支路和S→D支路均为符合独立非同分布的Rayleigh块衰落信道,
它包括如下步骤:
步骤C1:利用源节点S和中继节点R发送训练序列的方法,中继节点R获得信道状态信息T1(ΓSD,γSR,γRD)或T2(ΓSD,γSR,ΓRD);其中,ΓSD为S→D支路的平均信噪比,γSR为S→R支路的瞬时信噪比,γRD为R→D支路的瞬时信噪比,ΓRD为R→D支路的平均信噪比;
步骤C2:采用于基于信噪比的选择转发的方法,求得最优的功率修正系数αR
步骤C3:源节点S的发射功率不变,中继节点R根据功率修正系数αR调整发射功率传输信号;
步骤C4:目的节点D收到各支路信号后,采用选择合并器输出合并后的信号:
Figure BDA00003295361900061
Figure BDA00003295361900062
Figure BDA00003295361900063
MS为源节点S的调制阶数,MR为中继的调制阶数,γSD为S→D支路的瞬时信噪比;调制方式采用基于格雷码的方形M-QAM调制;
所述采用于基于信噪比的选择转发的方法为:
当中继节点R获得信道状态信息为T1(ΓSD,γSR,γRD)时,此时的误比特率PT1S,αR)为:
P T 1 ( α S , α R ) = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD + ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD 其中,PSD为S→D支路的瞬时误比特率,PSD=PMSSγSD)
S→R→D支路的瞬时误比特率PSRD=PMSSγSR)+(1-2PMSSγSR))PMRRγRD),
P M ( γ ) = c M Q ( 2 d M 2 γ ) , Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M s = 2 M s - 1 M s log 2 M s ;
以中间变量
Figure BDA00003295361900068
表示 ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD ,
P T 1 1 = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD
= c MS Q ( 2 d MS 2 β α R γ RD ) exp ( - β α R γ RD α S Γ SD ) - c MS d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD Q ( 2 ( 1 + d MS 2 α S Γ SD ) β α R γ RD α S Γ SD ) ,
以中间变量
Figure BDA000032953619000612
表示 ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD ,
P T 1 2 = P SRD ∫ γ SD = 0 β α R γ RD / α S f ( γ SD ) d γ SD = P SRD ( 1 - exp ( - β α R γ RD α S Γ SD ) ) ;
当PT2(1,1)<PT2(1,0)时,αR=1,否则αR=0;
当中继节点R获得信道状态信息为T2(ΓSD,γSR,ΓRD)时,此时的误比特率PT2S,αR)为:
P T 2 ( α S , α R ) = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD + ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
以中间变量
Figure BDA00003295361900071
表示 ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD ,
P T 2 1 = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD
= 1 2 c MS ( 1 - d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD ) ,
- 1 2 c MS β α R Γ RD α S Γ SD + β α R Γ RD ( 1 - d MS 2 β α R α S Γ SD Γ RD α S Γ SD + β α R Γ RD + d MS 2 β α S α R Γ SD Γ RD )
以中间变量
Figure BDA00003295361900076
表示 ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD ,
P T 2 2 = ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
= c MS Q ( 2 d MS 2 α S γ SR ) ( β α R Γ RD α S Γ SD + β α R Γ RD ) + ( 1 - 2 c MS Q ( 2 d MS 2 α S γ SR ) ) ×
( 1 2 c MR ( 1 - d MR 2 α R Γ RD 1 + d MR 2 α R Γ RD ) - 1 2 c MS α S Γ SD α S Γ SD + β α R Γ RD ( 1 - d MR 2 α S α R Γ SD Γ RD α S Γ SD + β α R Γ RD + d MR 2 α S α R Γ SD Γ RD ) )
Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M R = 2 M R - 1 M R log 2 M R , c M s = 2 M s - 1 M s log 2 M s ;
当PT2(1,1)<PT2(1,0)时,αR=1,否则αR=0。
本发明的优点在于,采用选择合并的单中继解码转发协作通信系统,提供了三种方法,所述三种方法能够充分利用瞬时的信道状态信息或平均的信道状态信息,改善由中继可能出现的解码错误导致的错误传播问题,提高解码转发协作传输的可靠性。本发明选择合并器只需要知道各支路信道的瞬时信噪比信息,并且仅需要对信噪比最大的支路进行解调,接收机的结构相对简单。
附图说明
图1为本发明具体实施方式一所述的基于选择合并的解码转发协作系统的信号智能转发方法的流程示意图。
图2为本发明具体实施方式二所述的基于选择合并的解码转发协作系统的信号智能转发方法的流程示意图。
图3为本发明具体实施方式三所述的基于选择合并的解码转发协作系统的信号智能转发方法的流程示意图。
图4为具体实施方式一、具体实施方式二和具体实施方式三的功率修正系数αR在信道信息为T1时(ΓSD=10dB,γRD=10dB)随γSR变化的曲线示意图。A1表示具体实施方式一的曲线,B1表示具体实施方式二的曲线,C1表示具体实施方式三的曲线;MS=4,MR=16。
图5为具体实施方式一、具体实施方式二和具体实施方式三的瞬时误比特率在信道信息为T1时(ΓSD=10dB,γRD=10dB)随γSR变化的曲线示意图和。A2表示具体实施方式一的曲线,B2表示具体实施方式二的曲线,C2表示具体实施方式三的曲线;D表示源节点S和中继等分功率传输时,在信道信息为T1时(ΓSD=10dB,γRD=10dB)随γSR变化的曲线示意图;MS=4,MR=16。
图6为具体实施方式一、具体实施方式二和具体实施方式三的功率修正系数αR在信道信息为T2时(ΓSD=10dB,ΓRD=10dB)随γSR变化的曲线示意图。MS=4,MR=16,A3表示具体实施方式一的曲线,B3表示具体实施方式二的曲线,C3表示具体实施方式三的曲线;
图7为具体实施方式一、具体实施方式二和具体实施方式三的瞬时误比特率在信道信息为T2时(ΓSD=10dB,ΓRD=10dB)随γSR变化的曲线示意图和。MS=4,MR=16,A4表示具体实施方式一的曲线,B4表示具体实施方式二的曲线,C4表示具体实施方式三的曲线;E表示源节点S和中继等分功率传输时,在信道信息为T2时(ΓSD=10dB,ΓRD=10dB)随γSR变化的曲线示意图。
具体实施方式
具体实施方式一、结合图1说明本实施方式,本实施所述的基于选择合并的解码转发协作系统的信号智能转发方法,所述方法是基于单中继的解码转发协作系统和采用功率分配的智能中继方法实现的,所述系统包括源节点S、中继节点R和目的节点D,上述节点仅装配单根天线且工作于半双工模式,S→R支路、R→D支路和S→D支路均为符合独立非同分布的Rayleigh块衰落信道,
它包括如下步骤:
步骤A1:利用源节点S和中继节点R发送训练序列的方法,中继节点R获得信道状态信息T1(ΓSD,γSR,γRD)或T2(ΓSD,γSR,ΓRD);其中,ΓSD为S→D支路的平均信噪比,γSR为S→R支路的瞬时信噪比,γRD为R→D支路的瞬时信噪比,ΓRD为R→D支路的平均信噪比;
步骤A2:采用功率分配的智能中继方法,求得最优的功率修正系数αS和αR
步骤A3:源节点S和中继节点R根据功率修正系数αS和αR调整发射功率传输信号;
步骤A4:目的节点D收到各支路信号后,采用选择合并器输出合并后的信号:
Figure BDA00003295361900092
Figure BDA00003295361900093
MS为源节点S的调制阶数,MR为中继的调制阶数,γSD为S→D支路的瞬时信噪比;调制方式采用基于格雷码的方形M-QAM调制;
所述采用功率分配的智能中继方法为:
当中继节点R获得信道状态信息为T1(ΓSD,γSR,γRD)时,此时的误比特率PT1S,αR)为:
P T 1 ( α S , α R ) = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD + ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD
其中,PSD为S→D支路的瞬时误比特率,PSD=PMSSγSD)
S→R→D支路的瞬时误比特率PSRD=PMSSγSR)+(1-2PMSSγSR))PMRRγRD),
P M ( γ ) = c M Q ( 2 d M 2 γ ) , Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M s = 2 M s - 1 M s log 2 M s ;
以中间变量
Figure BDA00003295361900098
表示 ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD ,
P T 1 1 = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD
= c MS Q ( 2 d MS 2 β α R γ RD ) exp ( - β α R γ RD α S Γ SD ) - c MS d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD Q ( 2 ( 1 + d MS 2 α S Γ SD ) β α R γ RD α S Γ SD ) ,
以中间变量表示 ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD ,
P T 1 2 = P SRD ∫ γ SD = 0 β α R γ RD / α S f ( γ SD ) d γ SD = P SRD ( 1 - exp ( - β α R γ RD α S Γ SD ) ) ;
使误比特率PT1S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS∈(0,2],αR∈(0,2],αSR=2;
当中继节点R获得信道状态信息为T2(ΓSD,γSR,ΓRD)时,此时的误比特率PT2S,αR)为:
P T 2 ( α S , α R ) = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD + ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
以中间变量
Figure BDA00003295361900101
表示 ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD ,
P T 2 1 = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD
= 1 2 c MS ( 1 - d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD ) ,
- 1 2 c MS β α R Γ RD α S Γ SD + β α R Γ RD ( 1 - d MS 2 β α R α S Γ SD Γ RD α S Γ SD + β α R Γ RD + d MS 2 β α S α R Γ SD Γ RD )
以中间变量
Figure BDA00003295361900106
表示 ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD ,
P T 2 2 = ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
= c MS Q ( 2 d MS 2 α S γ SR ) ( β α R Γ RD α S Γ SD + β α R Γ RD ) + ( 1 - 2 c MS Q ( 2 d MS 2 α S γ SR ) ) ×
( 1 2 c MR ( 1 - d MR 2 α R Γ RD 1 + d MR 2 α R Γ RD ) - 1 2 c MS α S Γ SD α S Γ SD + β α R Γ RD ( 1 - d MR 2 α S α R Γ SD Γ RD α S Γ SD + β α R Γ RD + d MR 2 α S α R Γ SD Γ RD ) )
Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M R = 2 M R - 1 M R log 2 M R , c M s = 2 M s - 1 M s log 2 M s ;
使误比特率PT2s,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS∈(0,2],αR∈(0,2],αSR=2。
具体实施方式二、结合图2说明本实施方式,本实施方式所述的基于选择合并的解码转发协作系统的信号智能转发方法,所述方法是基于单中继的解码转发协作系统和采用功率控制的智能中继方法实现的,所述系统包括源节点S、中继节点R和目的节点D,上述节点仅装配单根天线且工作于半双工模式,S→R支路、R→D支路和S→D支路均为符合独立非同分布的Rayleigh块衰落信道,
它包括如下步骤:
步骤B1:利用源节点S和中继节点R发送训练序列的方法,中继节点R获得信道状态信息T1(ΓSD,γSR,γRD)或T2(ΓSD,γSR,ΓRD);其中,ΓSD为S→D支路的平均信噪比,γSR为S→R支路的瞬时信噪比,γRD为R→D支路的瞬时信噪比,ΓRD为R→D支路的平均信噪比;
步骤B2:采用功率控制的智能中继方法,求得最优的功率修正系数αR
步骤B3:源节点S的发射功率不变,中继节点R根据功率修正系数αR调整发射功率传输信号;
步骤B4:目的节点D收到各支路信号后,采用选择合并器输出合并后的信号:
Figure BDA00003295361900111
Figure BDA00003295361900112
Figure BDA00003295361900113
MS为源节点S的调制阶数,MR为中继的调制阶数,γSD为S→D支路的瞬时信噪比;调制方式采用基于格雷码的方形M-QAM调制;
所述采用功率控制的智能中继方法为:
当中继节点R获得信道状态信息为T1(ΓSD,γSR,γRD)时,此时的误比特率PT1S,αR)为:
P T 1 ( α S , α R ) = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD + ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD
其中,PSD为S→D支路的瞬时误比特率,PSD=PMSSγSD)
S→R→D支路的瞬时误比特率PSRD=PMSSγSR)+(1-2PMSSγSR))PMRRγRD),
P M ( γ ) = c M Q ( 2 d M 2 γ ) , Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M s = 2 M s - 1 M s log 2 M s ;
以中间变量
Figure BDA00003295361900118
表示 ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD ,
P T 1 1 = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD
= c MS Q ( 2 d MS 2 β α R γ RD ) exp ( - β α R γ RD α S Γ SD ) - c MS d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD Q ( 2 ( 1 + d MS 2 α S Γ SD ) β α R γ RD α S Γ SD ) ,
以中间变量
Figure BDA000032953619001112
表示 ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD ,
P T 1 2 = P SRD ∫ γ SD = 0 β α R γ RD / α S f ( γ SD ) d γ SD = P SRD ( 1 - exp ( - β α R γ RD α S Γ SD ) ) ;
使误比特率PT1S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS=1,αR∈[0,1];
当中继节点R获得信道状态信息为T2SD,γSR,ΓRD)时,此时的误比特率PT2S,αR)为:
P T 2 ( α S , α R ) = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD + ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
以中间变量
Figure BDA00003295361900121
表示 ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD ,
其中, P T 2 1 = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD = 1 2 c MS ( 1 - d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD ) - 1 2 c MS β α R Γ RD α S Γ SD + β α R Γ RD ( 1 - d MS 2 β α R α S Γ SD Γ RD α S Γ SD + β α R Γ RD + d MS 2 β α S α R Γ SD Γ RD ) ,
以中间变量表示 ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD ,
P T 2 2 = ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
= c MS Q ( 2 d MS 2 α S γ SR ) ( β α R Γ RD α S Γ SD + β α R Γ RD ) + ( 1 - 2 c MS Q ( 2 d MS 2 α S γ SR ) ) ×
( 1 2 c MR ( 1 - d MR 2 α R Γ RD 1 + d MR 2 α R Γ RD ) - 1 2 c MS α S Γ SD α S Γ SD + β α R Γ RD ( 1 - d MR 2 α S α R Γ SD Γ RD α S Γ SD + β α R Γ RD + d MR 2 α S α R Γ SD Γ RD ) )
Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M R = 2 M R - 1 M R log 2 M R , c M R = 2 M R - 1 M R log 2 M R ,
使误比特率PT2S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS=1,αR∈[0,1]。
具体实施方式三、结合图3说明本实施方式,本实施方式所述的基于选择合并的解码转发协作系统的信号智能转发方法,所述方法是基于单中继的解码转发协作系统和采用于基于信噪比的选择转发的方法实现的,所述系统包括源节点S、中继节点R和目的节点D,上述节点仅装配单根天线且工作于半双工模式,S→R支路、R→D支路和S→D支路均为符合独立非同分布的Rayleigh块衰落信道,
它包括如下步骤:
步骤C1:利用源节点S和中继节点R发送训练序列的方法,中继节点R获得信道状态信息T1(ΓSD,γSR,γRD)或T2(ΓSD,γSR,ΓRD);其中,ΓSD为S→D支路的平均信噪比,γSR为S→R支路的瞬时信噪比,γRD为R→D支路的瞬时信噪比,ΓRD为R→D支路的平均信噪比;
步骤C2:采用于基于信噪比的选择转发的方法,求得最优的功率修正系数αR
步骤C3:源节点S的发射功率不变,中继节点R根据功率修正系数αR调整发射功率传输信号;
步骤C4:目的节点D收到各支路信号后,采用选择合并器输出合并后的信号:
Figure BDA00003295361900132
Figure BDA00003295361900133
MS为源节点S的调制阶数,MR为中继的调制阶数,γSD为S→D支路的瞬时信噪比;调制方式采用基于格雷码的方形M-QAM调制;
所述采用于基于信噪比的选择转发的方法为:
当中继节点R获得信道状态信息为T1(ΓSD,γSR,γRD)时,此时的误比特率PT1S,αR)为:
P T 1 ( α S , α R ) = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD + ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD 其中,PSD为S→D支路的瞬时误比特率,PSD=PMSSγSD)
S→R→D支路的瞬时误比特率PSRD=PMSSγSR)+(1-2PMSSγSR))PMRRγRD),
P M ( γ ) = c M Q ( 2 d M 2 γ ) , Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M s = 2 M s - 1 M s log 2 M s ;
以中间变量
Figure BDA00003295361900138
表示 ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD ,
P T 1 1 = ∫ γ SD = β α R γ RD / α S ∞ P SD f ( γ SD ) d γ SD
= c MS Q ( 2 d MS 2 β α R γ RD ) exp ( - β α R γ RD α S Γ SD ) - c MS d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD Q ( 2 ( 1 + d MS 2 α S Γ SD ) β α R γ RD α S Γ SD ) ,
以中间变量表示 ∫ γ SD = 0 β α R γ RD / α S P SRD f ( γ SD ) d γ SD ,
P T 1 2 = P SRD ∫ γ SD = 0 β α R γ RD / α S f ( γ SD ) d γ SD = P SRD ( 1 - exp ( - β α R γ RD α S Γ SD ) ) ;
当PT2(1,1)<PT2(1,0)时,αR=1,否则αR=0;
当中继节点R获得信道状态信息为T2(ΓSD,γSR,ΓRD)时,此时的误比特率PT2S,αR)为:
P T 2 ( α S , α R ) = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD + ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
以中间变量
Figure BDA00003295361900141
表示 ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD ,
P T 2 1 = ∫ γ RD = 0 ∞ P T 1 1 f ( γ RD ) d γ RD
= 1 2 c MS ( 1 - d MS 2 α S Γ SD 1 + d MS 2 α S Γ SD ) ,
- 1 2 c MS β α R Γ RD α S Γ SD + β α R Γ RD ( 1 - d MS 2 β α R α S Γ SD Γ RD α S Γ SD + β α R Γ RD + d MS 2 β α S α R Γ SD Γ RD )
以中间变量
Figure BDA00003295361900146
表示 ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD ,
P T 2 2 = ∫ γ RD = 0 ∞ P T 1 2 f ( γ RD ) d γ RD
= c MS Q ( 2 d MS 2 α S γ SR ) ( β α R Γ RD α S Γ SD + β α R Γ RD ) + ( 1 - 2 c MS Q ( 2 d MS 2 α S γ SR ) ) ×
( 1 2 c MR ( 1 - d MR 2 α R Γ RD 1 + d MR 2 α R Γ RD ) - 1 2 c MS α S Γ SD α S Γ SD + β α R Γ RD ( 1 - d MR 2 α S α R Γ SD Γ RD α S Γ SD + β α R Γ RD + d MR 2 α S α R Γ SD Γ RD ) )
Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M R = 2 M R - 1 M R log 2 M R , c M s = 2 M s - 1 M s log 2 M s ;
当PT2(1,1)<PT2(1,0)时,αR=1,否则αR=0。
上述三种实施方式中,
当源节点S和中继等分功率发射信号时,用γJK和ΓJK,JK∈{SR,RD,SD},分别表示各支路的瞬时信噪比和平均信噪比。由于信道符合瑞利衰落,各瞬时信噪比符合如下指数分布:
f ( γ JK ) = 1 Γ JK exp ( - γ JK Γ JK ) , γ JK ≥ 0 , JK ∈ { SD , SR , RD } - - - 1
从而当采用上述三种实施方式产生源节点S和中继的功率修正系数αS和αR时,S=D、S→R和R→D支路实际的每比特瞬时信噪比分别记为γ′SD=αSγSD、γ′SRαSγSR、γ′RD=αRγRD。S→D、S→R和R→D支路实际的每比特平均信噪比分别记为Γ′SD=αSΓSD、Γ′SRαSΓSR、Γ′RD=αRΓRD
上述实施三种实施方式中的调制方式采用基于格雷码的方形M-QAM,源节点S和中继调制阶数分别用MS和MR表示,每符号包含的比特数分别为kS=log2MS,kR=log2MR。M-QAM信号星座点间的最小欧氏距离为
d M = 3 log 2 M 2 ( M - 1 ) . - - - ( 1 )
目的节点在收到各支路信号后,采用选择合并器对其进行合并,这里采用系数β=(dMR/dMS)2对各支路不同调制阶数信号的信噪比进行修正。从而目的节点将根据如下规则选择某一路信号作为合并器的输出:
上述三种实施方式中的误比特率分析:
对基于格雷码的方形M-QAM调制信号,当瞬时信噪比为γ时,其误比特率可近似为
P M ( γ ) = 1 2 c M erfc ( d M 2 γ ) = c M Q ( 2 d M 2 γ ) - - - ( 3 )
其中
Q ( x ) = 1 2 π ∫ x ∞ e - t 2 / 2 dt , c M = 2 M - 1 M log 2 M - - - ( 4 )
设PSD、PSRD分别为S→D、S→R→D支路的瞬时误比特率,则已知T0(γSD,γSR,γRD)瞬时误比特率时,选择合并器的瞬时误比特率可表示为
P T 0 ( &alpha; S , &alpha; R ) = P SD = P MS ( &alpha; S &gamma; SD ) , &alpha; S &gamma; SD &GreaterEqual; &beta; &alpha; R &gamma; RD P SRD = P MS ( &alpha; S &gamma; SR ) + ( 1 - 2 P MS ( &alpha; S &gamma; SR ) ) P MR ( &alpha; R &gamma; RD ) . &alpha; S &gamma; SD < &beta; &alpha; R &gamma; RD . - - - ( 5 )
根据信道反馈信息的不同,设中继已知两种信道状态信息:T1(ΓSD,γSR,γRD)和T2(ΓSD,γSR,ΓRD)。下面对这两种情况下的误比特率进行分析。
1、已知T1时智能中继的误比特率分析
对选择合并器的瞬时误比特率PT0在γSD上做平均可得已知T1时的误比特率为
P T 1 ( &alpha; S , &alpha; R ) = E &gamma; SD [ P T 0 ( &alpha; S , &alpha; R ) ]
= &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD + &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S P SRD f ( &gamma; SD ) d &gamma; SD - - - ( 6 )
计算第一部分积分可得
P T 1 1 = &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD
= c MS Q ( 2 d MS 2 &beta; &alpha; R &gamma; RD ) exp ( - &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD )
- c MS d MS 2 &alpha; S &Gamma; SD 1 + d MS 2 &alpha; S &Gamma; SD Q ( 2 ( 1 + d MS 2 &alpha; S &Gamma; SD ) &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) - - - ( 7 )
计算第二部分积分可得
P T 1 2 = P SRD &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S f ( &gamma; SD ) d &gamma; SD = P SRD ( 1 - exp ( - &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) ) - - - ( 8 )
将(8)和(9)代入(7)即可得已知T1时智能中继的误比特率。
2、已知T2时智能中继的误比特率分析
对PT0在γSD和γRD上做平均可得已知T2时的误比特率为
P T 2 ( &alpha; S , &alpha; R ) = E &gamma; SD , &gamma; RD [ P T 0 ( &alpha; S , &alpha; R ) ] = E &gamma; RD [ P T 1 ( &alpha; S , &alpha; R ) ]
= &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD + &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD - - - ( 9 )
计算第一部分积分可得
P T 2 1 = &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD
= 1 2 c MS ( 1 - d MS 2 &alpha; S &Gamma; SD 1 + d MS 2 &alpha; S &Gamma; SD ) ,
- 1 2 c MS &beta; &alpha; R &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ( 1 - d MS 2 &beta; &alpha; R &alpha; S &Gamma; SD &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD + d MS 2 &beta; &alpha; S &alpha; R &Gamma; SD &Gamma; RD ) - - - ( 10 )
计算第二部分积分可得
P T 2 2 = &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD
= c MS Q ( 2 d MS 2 &alpha; S &gamma; SR ) ( &beta; &alpha; R &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ) + ( 1 - 2 c MS Q ( 2 d MS 2 &alpha; S &gamma; SR ) ) &times;
( 1 2 c MR ( 1 - d MR 2 &alpha; R &Gamma; RD 1 + d MR 2 &alpha; R &Gamma; RD ) - 1 2 c MS &alpha; S &Gamma; SD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ( 1 - d MR 2 &alpha; S &alpha; R &Gamma; SD &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD + d MR 2 &alpha; S &alpha; R &Gamma; SD &Gamma; RD ) ) - - - ( 11 )
将(11)和(12)代入(10)即可得已知T2时智能中继的误比特率。
利用(7)和(10)给出的中继已知两种信道状态信息时的误比特率公式,可给出下面三种智能中继方案,以求解最优的αS和αR,最小化误比特率。表1所示为三种智能中继方案功率修正系数的取值范围。
表1上述三种实施方式中的的功率修正系数
αS αR 限制条件
具体实施方式一:功率分配的智能中继方法 (0,2] [0,2) αSR=2
具体实施方式二:功率控制的智能中继方法 1 [0,1]
具体实施方式三:基于信噪比的选择转发的方法 1 0或1
1、采用功率分配的智能中继方法:中继根据信道状态信息,寻找最优的功率分配系数,在源和中继之间分配总的发射功率,以最小化系统误比特率。该优化问题可表述为:
&alpha; R opt = arg min &alpha; R &Element; [ 0,2 ) P Ti ( 2 - &alpha; R , &alpha; R ) , i = 1,2 - - - ( 12 )
2、采用功率控制的智能中继方法:源节点发射功率不变,中继寻找最优的功率修正系数,连续地调整自身的发射功率。该优化问题可表述为:
&alpha; R opt = arg min &alpha; R &Element; [ 0,1 ] P Ti ( 1 , &alpha; R ) , i = 1,2 - - - ( 13 )
式(17)和(18)的最优解可通过黄金分割搜索法容易地得到。
3、基于信噪比的选择转发方法。中继仅在满足如下条件时才转发:
PTi(1,1)<PTi(1,0),i=1,2(14)
即中继仅在协作的误比特率小于直接传输时才参与协作转发信号。
从图4、图5、图6和图7中可以看出,与αS=αR=1的固定功率方案相比,当瞬时信噪比γSR较低,即S→R信道的传输不可靠时,上述三种实施方式都可以通过降低中继的发射功率(αR<1),降低系统的误比特率,消除错误传播的影响。而尽管在γSR较低时,具体实施方式一的方法与具体实施方式二的方法求得的系数αR非常接近,但由于在具体实施方式一的方法中若中继降低功率,源节点可以获得更大的发射功率,因此其性能要比具体实施方式二的方法好。而具体实施方式二的方法相比具体实施方式三的方法,性能差异很小,只能在基于信噪比的转发方案的信噪比门限附近获得一点可靠性的增益。在高信噪比情况下,对于具体实施方式一的方法,已知T2时,系数αR<1且性能与其它两种实施方式一致;而已知T1时,系数αR能接近2且误比特率大大下降了,即R→D支路的瞬时信噪比信息能给具体实施方式一的方法带来明显的增益。

Claims (3)

1.基于选择合并的解码转发协作系统的信号智能转发方法,其特征在于,所述方法是基于单中继的解码转发协作系统和采用功率分配的智能中继方法实现的,所述系统包括源节点S、中继节点R和目的节点D,上述节点仅装配单根天线且工作于半双工模式,S→R支路、R→D支路和S→D支路均为符合独立非同分布的Raψleigh块衰落信道,
它包括如下步骤:
步骤A1:利用源节点S和中继节点R发送训练序列的方法,中继节点R获得信道状态信息T1(ΓSD,γSR,γRD)或T2(ΓSD,γSR,ΓRD);其中,ΓSD为S→D支路的平均信噪比,γSR为S→R支路的瞬时信噪比,γRD为R→D支路的瞬时信噪比,ΓRD为R→D支路的平均信噪比;
步骤A2:采用功率分配的智能中继方法,求得最优的功率修正系数αS和αR
步骤A3:源节点S和中继节点R根据功率修正系数αS和αR调整发射功率传输信号;
步骤A4:目的节点D收到各支路信号后,采用选择合并器输出合并后的信号:
Figure FDA00003295361800011
Figure FDA00003295361800012
Figure FDA00003295361800013
MS为源节点S的调制阶数,MR为中继的调制阶数,γSD为S→D支路的瞬时信噪比;调制方式采用基于格雷码的方形M-QAM调制;
所述采用功率分配的智能中继方法为:
当中继节点R获得信道状态信息为T1(ΓSD,γSR,γRD)时,此时的误比特率PT1S,αR)为:
P T 1 ( &alpha; S , &alpha; R ) = &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD + &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S P SRD f ( &gamma; SD ) d &gamma; SD
其中,PSD为S→D支路的瞬时误比特率,PSD=PMSSγSD)
S→R→D支路的瞬时误比特率PSRD=PMSSγSR)+(1-2PMSSγSR))PMRRγRD),
P M ( &gamma; ) = c M Q ( 2 d M 2 &gamma; ) , Q ( x ) = 1 2 &pi; &Integral; x &infin; e - t 2 / 2 dt , c M s = 2 M s - 1 M s log 2 M s ;
以中间变量
Figure FDA00003295361800018
表示 &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD ,
P T 1 1 = &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD
= c MS Q ( 2 d MS 2 &beta; &alpha; R &gamma; RD ) exp ( - &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) - c MS d MS 2 &alpha; S &Gamma; SD 1 + d MS 2 &alpha; S &Gamma; SD Q ( 2 ( 1 + d MS 2 &alpha; S &Gamma; SD ) &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) ,
以中间变量表示 &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S P SRD f ( &gamma; SD ) d &gamma; SD ,
P T 1 2 = P SRD &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S f ( &gamma; SD ) d &gamma; SD = P SRD ( 1 - exp ( - &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) ) ;
使误比特率PT1S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS∈(0,2],αR∈(0,2],αSR=2;
当中继节点R获得信道状态信息为T2(ΓSD,γSR,ΓRD)时,此时的误比特率PT2S,αR)为:
P T 2 ( &alpha; S , &alpha; R ) = &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD + &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD
以中间变量
Figure FDA00003295361800027
表示 &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD ,
P T 2 1 = &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD
= 1 2 c MS ( 1 - d MS 2 &alpha; S &Gamma; SD 1 + d MS 2 &alpha; S &Gamma; SD ) ,
- 1 2 c MS &beta; &alpha; R &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ( 1 - d MS 2 &beta; &alpha; R &alpha; S &Gamma; SD &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD + d MS 2 &beta; &alpha; S &alpha; R &Gamma; SD &Gamma; RD )
以中间变量
Figure FDA000032953618000212
表示 &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD ,
P T 2 2 = &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD
= c MS Q ( 2 d MS 2 &alpha; S &gamma; SR ) ( &beta; &alpha; R &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ) + ( 1 - 2 c MS Q ( 2 d MS 2 &alpha; S &gamma; SR ) ) &times;
( 1 2 c MR ( 1 - d MR 2 &alpha; R &Gamma; RD 1 + d MR 2 &alpha; R &Gamma; RD ) - 1 2 c MS &alpha; S &Gamma; SD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ( 1 - d MR 2 &alpha; S &alpha; R &Gamma; SD &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD + d MR 2 &alpha; S &alpha; R &Gamma; SD &Gamma; RD ) )
Q ( x ) = 1 2 &pi; &Integral; x &infin; e - t 2 / 2 dt , c M s = 2 M s - 1 M s log 2 M s ; c M s = 2 M s - 1 M s log 2 M s ;
使误比特率PT2S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS∈(0,2],αR∈(0,2],αSR=2。
2.基于选择合并的解码转发协作系统的信号智能转发方法,其特征在于,所述方法是基于单中继的解码转发协作系统和采用功率控制的智能中继方法实现的,所述系统包括源节点S、中继节点R和目的节点D,上述节点仅装配单根天线且工作于半双工模式,S→R支路、R→D支路和S→D支路均为符合独立非同分布的Rayleigh块衰落信道,
它包括如下步骤:
步骤B1:利用源节点S和中继节点R发送训练序列的方法,中继节点R获得信道状态信息T1(ΓSD,γSR,γRD)或T2(ΓSD,γSR,ΓRD);其中,ΓSD为S→D支路的平均信噪比,γSR为S→R支路的瞬时信噪比,γRD为R→D支路的瞬时信噪比,ΓRD为R→D支路的平均信噪比;
步骤B2:采用功率控制的智能中继方法,求得最优的功率修正系数αR
步骤B3:源节点S的发射功率不变,中继节点R根据功率修正系数αR调整发射功率传输信号;
步骤B4:目的节点D收到各支路信号后,采用选择合并器输出合并后的信号:
Figure FDA00003295361800031
Figure FDA00003295361800032
Figure FDA00003295361800033
MS为源节点S的调制阶数,MR为中继的调制阶数,γSD为S→D支路的瞬时信噪比;调制方式采用基于格雷码的方形M-QAM调制;
所述采用功率控制的智能中继方法为:
当中继节点R获得信道状态信息为T1(ΓSD,γSR,γRD)时,此时的误比特率PT1S,αR)为:
P T 1 ( &alpha; S , &alpha; R ) = &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD + &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S P SRD f ( &gamma; SD ) d &gamma; SD
中,PSD为S→D支路的瞬时误比特率,PSD=PMSSγSD)
S→R→D支路的瞬时误比特率PSRD=PMSSγSR)+(1-2PMSSγSR))PMRRγRD),
P M ( &gamma; ) = c M Q ( 2 d M 2 &gamma; ) , Q ( x ) = 1 2 &pi; &Integral; x &infin; e - t 2 / 2 dt , c M s = 2 M s - 1 M s log 2 M s ;
以中间变量
Figure FDA00003295361800038
表示 &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD ,
P T 1 1 = &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD
= c MS Q ( 2 d MS 2 &beta; &alpha; R &gamma; RD ) exp ( - &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) - c MS d MS 2 &alpha; S &Gamma; SD 1 + d MS 2 &alpha; S &Gamma; SD Q ( 2 ( 1 + d MS 2 &alpha; S &Gamma; SD ) &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) ,
以中间变量
Figure FDA00003295361800043
表示 &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S P SRD f ( &gamma; SD ) d &gamma; SD ,
P T 1 2 = P SRD &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S f ( &gamma; SD ) d &gamma; SD = P SRD ( 1 - exp ( - &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) ) ;
使误比特率PT1S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS=1,αR∈[0,1];
当中继节点R获得信道状态信息为T2(ΓSD,γSR,ΓRD)时,此时的误比特率PT2S,αR)为:
P T 2 ( &alpha; S , &alpha; R ) = &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD + &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD
以中间变量
Figure FDA00003295361800047
表示 &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD ,
P T 2 1 = &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD
= 1 2 c MS ( 1 - d MS 2 &alpha; S &Gamma; SD 1 + d MS 2 &alpha; S &Gamma; SD ) ,
- 1 2 c MS &beta; &alpha; R &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ( 1 - d MS 2 &beta; &alpha; R &alpha; S &Gamma; SD &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD + d MS 2 &beta; &alpha; S &alpha; R &Gamma; SD &Gamma; RD )
以中间变量
Figure FDA000032953618000412
表示 &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD ,
P T 2 2 = &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD
= c MS Q ( 2 d MS 2 &alpha; S &gamma; SR ) ( &beta; &alpha; R &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ) + ( 1 - 2 c MS Q ( 2 d MS 2 &alpha; S &gamma; SR ) ) &times;
( 1 2 c MR ( 1 - d MR 2 &alpha; R &Gamma; RD 1 + d MR 2 &alpha; R &Gamma; RD ) - 1 2 c MS &alpha; S &Gamma; SD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ( 1 - d MR 2 &alpha; S &alpha; R &Gamma; SD &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD + d MR 2 &alpha; S &alpha; R &Gamma; SD &Gamma; RD ) )
Q ( x ) = 1 2 &pi; &Integral; x &infin; e - t 2 / 2 dt , c M R = 2 M R - 1 M R log 2 M R , c M R = 2 M R - 1 M R log 2 M R ,
使误比特率PT2S,αR)最小,采用黄金分割搜索法求得最优的功率修正系数αS和αR,其中αS=1,αR∈[0,1]。
3.基于选择合并的解码转发协作系统的信号智能转发方法,其特征在于,所述方法是基于单中继的解码转发协作系统和采用于基于信噪比的选择转发的方法实现的,所述系统包括源节点S、中继节点R和目的节点D,上述节点仅装配单根天线且工作于半双工模式,S→R支路、R→D支路和S→D支路均为符合独立非同分布的Rayleigh块衰落信道,
它包括如下步骤:
步骤C1:利用源节点S和中继节点R发送训练序列的方法,中继节点R获得信道状态信息T1(ΓSD,γSR,γRD)或T2(ΓSD,γSR,ΓRD);其中,ΓSD为S→D支路的平均信噪比,γSR为S→R支路的瞬时信噪比,γRD为R→D支路的瞬时信噪比,ΓRD为R→D支路的平均信噪比;
步骤C2:采用于基于信噪比的选择转发的方法,求得最优的功率修正系数αR
步骤C3:源节点S的发射功率不变,中继节点R根据功率修正系数αR调整发射功率传输信号;
步骤C4:目的节点D收到各支路信号后,采用选择合并器输出合并后的信号:
Figure FDA00003295361800051
Figure FDA00003295361800052
Figure FDA00003295361800053
MS为源节点S的调制阶数,MR为中继的调制阶数,γSD为S→D支路的瞬时信噪比;调制方式采用基于格雷码的方形M-QAM调制;
所述采用于基于信噪比的选择转发的方法为:
当中继节点R获得信道状态信息为T1(ΓSD,γSR,γRD)时,此时的误比特率PT1S,αR)为:
P T 1 ( &alpha; S , &alpha; R ) = &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD + &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S P SRD f ( &gamma; SD ) d &gamma; SD
其中,PSD为S→D支路的瞬时误比特率,PSD=PMSSγSD)
S→R→D支路的瞬时误比特率PSRD=PMSSγSR)+(1-2PMSSγSR))PMRRγRD),
P M ( &gamma; ) = c M Q ( 2 d M 2 &gamma; ) , Q ( x ) = 1 2 &pi; &Integral; x &infin; e - t 2 / 2 dt , c M s = 2 M s - 1 M s log 2 M s ;
以中间变量表示 &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD ,
P T 1 1 = &Integral; &gamma; SD = &beta; &alpha; R &gamma; RD / &alpha; S &infin; P SD f ( &gamma; SD ) d &gamma; SD
= c MS Q ( 2 d MS 2 &beta; &alpha; R &gamma; RD ) exp ( - &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) - c MS d MS 2 &alpha; S &Gamma; SD 1 + d MS 2 &alpha; S &Gamma; SD Q ( 2 ( 1 + d MS 2 &alpha; S &Gamma; SD ) &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) ,
以中间变量
Figure FDA00003295361800063
表示 &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S P SRD f ( &gamma; SD ) d &gamma; SD ,
P T 1 2 = P SRD &Integral; &gamma; SD = 0 &beta; &alpha; R &gamma; RD / &alpha; S f ( &gamma; SD ) d &gamma; SD = P SRD ( 1 - exp ( - &beta; &alpha; R &gamma; RD &alpha; S &Gamma; SD ) ) ;
当PT1(1,1)<PT1(1,0)时,αR=1,否则αR=0;
当中继节点R获得信道状态信息为T2(ΓSD,γSR,ΓRD)时,此时的误比特率PT2S,αR)为:
P T 2 ( &alpha; S , &alpha; R ) = &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD + &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD
以中间变量
Figure FDA00003295361800067
表示 &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD ,
P T 2 1 = &Integral; &gamma; RD = 0 &infin; P T 1 1 f ( &gamma; RD ) d &gamma; RD
= 1 2 c MS ( 1 - d MS 2 &alpha; S &Gamma; SD 1 + d MS 2 &alpha; S &Gamma; SD ) ,
- 1 2 c MS &beta; &alpha; R &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ( 1 - d MS 2 &beta; &alpha; R &alpha; S &Gamma; SD &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD + d MS 2 &beta; &alpha; S &alpha; R &Gamma; SD &Gamma; RD )
以中间变量表示 &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD ,
P T 2 2 = &Integral; &gamma; RD = 0 &infin; P T 1 2 f ( &gamma; RD ) d &gamma; RD
= c MS Q ( 2 d MS 2 &alpha; S &gamma; SR ) ( &beta; &alpha; R &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ) + ( 1 - 2 c MS Q ( 2 d MS 2 &alpha; S &gamma; SR ) ) &times;
( 1 2 c MR ( 1 - d MR 2 &alpha; R &Gamma; RD 1 + d MR 2 &alpha; R &Gamma; RD ) - 1 2 c MS &alpha; S &Gamma; SD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD ( 1 - d MR 2 &alpha; S &alpha; R &Gamma; SD &Gamma; RD &alpha; S &Gamma; SD + &beta; &alpha; R &Gamma; RD + d MR 2 &alpha; S &alpha; R &Gamma; SD &Gamma; RD ) )
Q ( x ) = 1 2 &pi; &Integral; x &infin; e - t 2 / 2 dt , c M R = 2 M R - 1 M R log 2 M R , c M s = 2 M s - 1 M s log 2 M s ;
当PT2(1,1)<PT2(1,0)时,αR=1,否则αR=0。
CN201310216412.4A 2013-06-04 2013-06-04 基于选择合并的解码转发协作系统的信号智能转发方法 Active CN103281111B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310216412.4A CN103281111B (zh) 2013-06-04 2013-06-04 基于选择合并的解码转发协作系统的信号智能转发方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310216412.4A CN103281111B (zh) 2013-06-04 2013-06-04 基于选择合并的解码转发协作系统的信号智能转发方法

Publications (2)

Publication Number Publication Date
CN103281111A true CN103281111A (zh) 2013-09-04
CN103281111B CN103281111B (zh) 2016-08-10

Family

ID=49063576

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310216412.4A Active CN103281111B (zh) 2013-06-04 2013-06-04 基于选择合并的解码转发协作系统的信号智能转发方法

Country Status (1)

Country Link
CN (1) CN103281111B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103561447A (zh) * 2013-11-08 2014-02-05 山东大学 基于机会中继的增量混合式译码放大转发协作方法
WO2018210095A1 (zh) * 2017-05-19 2018-11-22 华为技术有限公司 一种功率控制方法及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101282199A (zh) * 2008-04-14 2008-10-08 北京邮电大学 用于多中继协作通信的中继策略的自适应选择方法
WO2012148256A1 (en) * 2011-04-25 2012-11-01 Mimos Berhad Improved macro diversity handover in wireless multi-hop relay- networks
US20120300666A1 (en) * 2010-02-08 2012-11-29 Ji Woong Jang Coordinated communication method using multiple terminals
CN103095360A (zh) * 2013-02-27 2013-05-08 哈尔滨工业大学 一种自适应解码转发协作系统的信号选择合并方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101282199A (zh) * 2008-04-14 2008-10-08 北京邮电大学 用于多中继协作通信的中继策略的自适应选择方法
US20120300666A1 (en) * 2010-02-08 2012-11-29 Ji Woong Jang Coordinated communication method using multiple terminals
WO2012148256A1 (en) * 2011-04-25 2012-11-01 Mimos Berhad Improved macro diversity handover in wireless multi-hop relay- networks
CN103095360A (zh) * 2013-02-27 2013-05-08 哈尔滨工业大学 一种自适应解码转发协作系统的信号选择合并方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103561447A (zh) * 2013-11-08 2014-02-05 山东大学 基于机会中继的增量混合式译码放大转发协作方法
CN103561447B (zh) * 2013-11-08 2016-08-17 山东大学 基于机会中继的增量混合式译码放大转发协作方法
WO2018210095A1 (zh) * 2017-05-19 2018-11-22 华为技术有限公司 一种功率控制方法及设备

Also Published As

Publication number Publication date
CN103281111B (zh) 2016-08-10

Similar Documents

Publication Publication Date Title
CN103763015B (zh) 一种有直连链路的多天线中继网络中发送天线选择方法
CN107040310B (zh) 基于ihdaf协议的fso通信系统的工作方法
CN105871506B (zh) 无线通信自适应传输方法
CN106413099B (zh) 全双工协作中继网络的最优功率分配方法
Ikki et al. Performance of decode-and-forward cooperative diversity networks over Nakagami-m fading channels
CN103906198A (zh) 放大转发协作网络的中继选择与功率分配方法
CN106533514B (zh) 基于ihdaf的协作空间调制系统的工作方法
Annamalai et al. Tight bounds on the ergodic capacity of cooperative analog relaying with adaptive source transmission techniques
CN105517096A (zh) 一种多中继放大转发协作网络的中继选择方法
CN102769486B (zh) 双向多跳中继系统中的中继端信号处理方法
CN103281111A (zh) 基于选择合并的解码转发协作系统的信号智能转发方法
Jia et al. Two-way denoise-and-forward network coding opportunistic relaying with jointing adaptive modulation relay selection criterions
CN106535278A (zh) 基于完全信道状态信息的全双工中继选择方法
CN103441787A (zh) 双向中继系统中基于天线和用户联合选择的信息传输方法
CN103095360B (zh) 一种自适应解码转发协作系统的信号选择合并方法
CN105553536A (zh) 一种重叠码分复用混合转发协作通信方法
CN105007113A (zh) 一种信息和能量互换的双向中继通信方法
CN102780545B (zh) 放大转发分布式协作系统发射天线选择方法
Michalopoulos et al. Distributed transmit antenna selection (DTAS) under performance or energy consumption constraints
CN110278019B (zh) 基于自干扰最小化准则的全双工能量采集中继传输方法
CN107197500A (zh) 一种用于空间调制解码转发中继系统的中继节点选择方法
CN102111357B (zh) 一种克服信号畸变的中继解调转发系统及方法
CN104994043B (zh) 基于节点选择的卫星移动通信自适应协作传输方法
CN102665260A (zh) 结合功率分配的多天线终端协作式mimo中继选择方法
CN103368692A (zh) 一种双向中继系统中自适应变时隙模拟网络编码策略

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant