CN105282814A - 一种预测csi的中继选择方法和系统 - Google Patents

一种预测csi的中继选择方法和系统 Download PDF

Info

Publication number
CN105282814A
CN105282814A CN201510581857.1A CN201510581857A CN105282814A CN 105282814 A CN105282814 A CN 105282814A CN 201510581857 A CN201510581857 A CN 201510581857A CN 105282814 A CN105282814 A CN 105282814A
Authority
CN
China
Prior art keywords
gamma
rho
prime
overbar
lsqb
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
CN201510581857.1A
Other languages
English (en)
Other versions
CN105282814B (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.)
Guangxi Normal University
Original Assignee
Guangxi Normal University
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 Guangxi Normal University filed Critical Guangxi Normal University
Priority to CN201510581857.1A priority Critical patent/CN105282814B/zh
Publication of CN105282814A publication Critical patent/CN105282814A/zh
Application granted granted Critical
Publication of CN105282814B publication Critical patent/CN105282814B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W40/00Communication routing or communication path finding
    • H04W40/02Communication route or path selection, e.g. power-based or shortest path routing
    • H04W40/22Communication route or path selection, e.g. power-based or shortest path routing using selective relaying for reaching a BTS [Base Transceiver Station] or an access point
    • 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

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Picture Signal Circuits (AREA)
  • Data Exchanges In Wide-Area Networks (AREA)

Abstract

本发明涉及一种预测CSI的中继选择方法和系统,其中方法包括步骤S1,计算系统模型中所有链路的预测信噪比;步骤S2,根据步骤S1的计算结果,计算被选择的中继节点;在执行步骤S1之前,还可以包括步骤S0:计算系统模型中源节点和所有中继节点的发射功率。本发明采用预测CSI进行中继选择,有效克服了过时CSI对系统性能的影响,具有科学意义和应用价值;本发明在测量信噪比基础上加上特定分布的随机数生成预测信噪比,预测精度比直接预测信噪比高,操作简单、计算量少、精度较高;本发明提出的考虑功率资源公平分配的功率分配方法,并用于中继选择方法中,为解决协作通信中的公平问题提供了有益的参考。

Description

一种预测CSI的中继选择方法和系统
技术领域
本发明涉及无线协作通信领域,尤其涉及一种预测CSI的中继选择方法和系统。
背景技术
协作分集是无线通信中对抗多径效应的有效方法,通过引入中继节点协助发送端与接收端之间的信息通信,使得传输路径得到共享,可以显著提高网络吞吐量和信息传输的可靠性。当可选的中继节点增多时,如何从众多的中继节点中选择合适的中继节点转发信息是非常重要的研究问题。中继选择不仅要考虑误码率、中断概率等性能指标,也要考虑资源的优化利用如功率分配等。
已有的中继选择和功率分配方法一般都建立在过时的信道状态信息(CSI,ChannelStatusInformation)基础上,与中继转发时的CSI可能存在显著的差异。因此在中继选择和功率分配中引入预测信噪比能有效保障系统的性能。但是目前基于预测信噪比的中继选择和功率分配方法尚不多见,而且存在预测误差大、预测算法复杂等缺陷。
发明内容
本发明所要解决的技术问题,是提供一种预测精度高的预测CSI的中继选择方法和系统。
本发明解决上述技术问题的技术方案如下:一种预测CSI的中继选择方法,包括以下步骤:
步骤S1,计算系统模型中所有链路的预测信噪比;
所述系统模型包括一个源节点S、一个目的节点D和Nrelay个中继节点,其中Nrelay≥1,记第i个中继节点为中继节点i,其中i=1,2,...,Nrelay,Srelay表示中继节点集,且Srelay={1,2,...,Nrelay},所述链路包括源节点S到所有中继节点之间的链路和所有中继节点到目的节点D之间的链路;源节点S到中继节点i之间的链路和中继节点i到目的节点D之间的链路构成中继链路i;a∈{1,2},当a=1时表示源节点S到中继节点i之间的链路,当a=2时表示中继节点i到目的节点D之间的链路;
步骤S2,根据步骤S1的计算结果,计算被选择的中继节点。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步地,在执行步骤S1之前,还包括步骤S0:计算系统模型中源节点和所有中继节点的发射功率。
进一步地,步骤S1中,记源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的预测信噪比为γai,则
其中为源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的测量信噪比,Δai为服从概率密度函数为的分布的随机数,
f Δ a i | γ ^ a i ( Δ a i | γ ^ a i ) = 1 γ ‾ a i ( 1 - ( ρ a i ) 2 ) exp ( - [ Δ a i + ( 1 + ( ρ a i ) 2 ) γ ^ a i ] γ ‾ a i ( 1 - ( ρ a i ) 2 ) ) × I 0 ( 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ‾ a i ( 1 - ( ρ a i ) 2 ) ) ;
其中表示源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的平均信噪比,ρai表示源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的预测信道系数和测量信道系数的相关系数, I 0 ( 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ‾ a i ( 1 - ( ρ a i ) 2 ) ) 表示自变量为 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ‾ a i ( 1 - ( ρ a i ) 2 ) 的第一类零阶变形贝塞尔函数。
进一步地,步骤S2中,用i*表示被选择的中继节点,则
i * = arg m a x i ∈ S r e l a y { γ i } ;
其中i*∈Srelay,γi表示中继链路i的预测信噪比,γi=min{γ1i2i},γ1i、γ2i分别为a=1及a=2时的γai
进一步地,步骤S0中,用P0表示源节点S的发射功率,P1,P2,...,PNrelay分别表示Nrelay个中继节点的发射功率,记发射功率集合为P,则
P = { P 0 , P 1 , P 2 , ... , P N r e l a y } ;
在特定的能量限制下最小化系统中断概率以分配源节点和所有中继节点的发射功率,即得出发射功率集合所述特定能量限制为:
s . t . P 0 + Σ i = 1 N r e l a y P i p s e l e c t i ( P ) ≤ E T ,
E min i ≤ P i p s e l e c t i ( P ) ≤ E m a x i ,
E min 0 ≤ P 0 ≤ E m a x 0 ;
其中, J ( P ) = Σ i = 1 N r e l a y F i ( γ T ) p s e l e c t i ( P ) ;
p s e l e c t i ( P ) = Π j = 1 j ≠ i N r e l a y p j i ;
p j i = p 1 i 2 i p 1 j 2 j p 1 j 1 i + p 2 i 1 i p 1 j 2 j p 1 j 2 i + p 1 i 2 i p 2 j 1 j p 2 j 1 i + p 2 i 1 i p 2 j 1 j p 2 j 2 i ;
p 1 i 2 i = exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 1 = 0 ∞ Σ k 2 = 0 ∞ Σ k 3 = 0 k 2 ( k 1 + k 3 ) ! ( k 1 ! ) 2 k 2 ! k 3 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 1 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 2 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 1 + k 3 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 1 - k 2 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 1 + k 3 ;
p 1 j 2 j = exp ( - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ ) Σ k 4 = 0 ∞ Σ k 5 = 0 ∞ Σ k 6 = 0 k 5 ( k 4 + k 6 ) ! ( k 4 ! ) 2 k 5 ! k 6 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 4 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 5 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 4 + k 6 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 4 - k 5 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 4 + k 6 ;
p 1 j 1 i = exp ( - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 7 = 0 ∞ Σ k 8 = 0 ∞ Σ k 9 = 0 k 8 ( k 7 + k 9 ) ! ( k 7 ! ) 2 k 8 ! k 9 !
× [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 7 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 8 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 7 + k 9 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 7 - k 8 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 7 + k 9 ;
p 2 i 1 i = exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 10 = 0 ∞ Σ k 11 = 0 ∞ Σ k 12 = 0 k 11 ( k 10 + k 12 ) ! ( k 10 ! ) 2 k 11 ! k 12 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 11 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 10 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 10 - k 11 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 10 + k 12 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 10 + k 12 ;
p 1 j 2 i = exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ ) Σ k 13 = 0 ∞ Σ k 14 = 0 ∞ Σ k 15 = 0 k 14 ( k 13 + k 15 ) ! ( k 13 ! ) 2 k 14 ! k 15 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 13 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 14 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 13 + k 15 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 13 - k 14 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 13 + k 15 ;
p 2 j 1 j = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ ) Σ k 16 = 0 ∞ Σ k 17 = 0 ∞ Σ k 18 = 0 k 17 ( k 16 + k 18 ) ! ( k 16 ! ) 2 k 17 ! k 18 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 17 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 16 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] 1 + k 16 - k 17 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 16 + k 18 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 16 + k 18 ;
p 2 j 1 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 19 = 0 ∞ Σ k 20 = 0 ∞ Σ k 21 = 0 k 20 ( k 19 + k 21 ) ! ( k 19 ! ) 2 k 20 ! k 21 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 20 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 19 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 19 - k 20 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 19 + k 21 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 19 + k 21 ;
p 2 j 2 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 22 = 0 ∞ Σ k 23 = 0 ∞ Σ k 24 = 0 k 23 ( k 22 + k 24 ) ! ( k 22 ! ) 2 k 23 ! k 24 ! × [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 23 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 22 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 22 - k 23 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 22 + k 24 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 22 + k 24 ;
F i ( γ T ) = exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 25 = 0 ∞ 1 ( k 25 ! ) 2 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 25 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 25 γ ( k 25 + 1 , γ T ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) + exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 26 = 0 ∞ 1 ( k 26 ! ) 2 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 26 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 26 γ ( k 26 + 1 , γ T ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ )
- exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 27 = 0 ∞ Σ k 28 = 0 ∞ 1 ( k 27 ! ) 2 ( k 28 ! ) 2 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 28 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 27 × [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 28 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 27 γ ( k 28 + 1 , γ T ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) γ ( k 27 + 1 , γ T ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) ;
其中Pi表示中继节点i的发射功率;为中继节点i被选择的概率;ET为归一化的可用的总能量;为归一化的可分配给中继节点i的最小能量,为归一化的可分配给中继节点i的最大能量,为归一化的可分配给源节点S的最小能量,为归一化的可分配给源节点S的最大能量;J(P)是系统中断概率;FiT)为中继链路i发生中断的概率,γT为中继链路发生中断的信噪比阈值;j∈Srelay,即j=1,2,...,Nrelay,且j≠i;为中继链路i的预测信噪比大于等于中继链路j的预测信噪比的概率;为源节点S到中继节点i之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于中继节点j到目的节点D之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率;为中继节点i到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率;为中继节点j到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点j之间链路的预测信噪比的概率,为中继节点j到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率;为中继节点j到目的节点D之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率;为分配功率前源节点S到中继节点i之间链路的平均信噪比,为源节点S到中继节点i之间链路的信道系数的方差,为中继节点i处源节点S的噪声方差,为分配功率前中继节点i到目的节点D之间链路的平均信噪比,为中继节点i到目的节点D之间链路的信道系数的方差,为目的节点D处中继节点i的噪声方差;为分配功率前源节点S到中继节点i之间链路的测量信噪比,为分配功率前中继节点i到目的节点D之间链路的测量信噪比;为分配功率前源节点S到中继节点j之间链路的平均信噪比,为源节点S到中继节点j之间链路的信道系数的方差,为中继节点j处源节点S的噪声方差,为分配功率前中继节点j到目的节点D之间链路的平均信噪比,Pj表示中继节点j的发射功率,为中继节点j到目的节点D之间链路的信道系数的方差,为目的节点D处中继节点j的噪声方差;为分配功率前源节点S到中继节点j之间链路的测量信噪比,为分配功率前中继节点j到目的节点D之间链路的测量信噪比;k1,k2,k3,……,k27,k28为求和下标指示,均为大于或等于0的整数;ρ1i为源节点S到中继节点i之间链路的预测信道系数和测量信道系数的相关系数,ρ1j为源节点S到中继节点j之间链路的预测信道系数和测量信道系数的相关系数,ρ2i为中继节点i到目的节点D之间链路的预测信道系数和测量信道系数的相关系数,ρ2j为中继节点j到目的节点D之间链路的预测信道系数和测量信道系数的相关系数;γ()表示不完全伽马函数。
本发明解决上述技术问题的另一种技术方案如下:一种预测CSI的中继选择系统,包括预测信噪比计算模块和被选择的中继节点计算模块;
所述预测信噪比计算模块用于计算系统模型中所有链路的预测信噪比;
所述系统模型包括一个源节点S、一个目的节点D和Nrelay个中继节点,其中Nrelay≥1,记第i个中继节点为中继节点i,其中i=1,2,...,Nrelay,Srelay表示中继节点集,且Srelay={1,2,...,Nrelay},所述链路包括源节点S到所有中继节点之间的链路和所有中继节点到目的节点D之间的链路;源节点S到中继节点i之间的链路和中继节点i到目的节点D之间的链路构成中继链路i;a∈{1,2},当a=1时表示源节点S到中继节点i之间的链路,当a=2时表示中继节点i到目的节点D之间的链路;
所述被选择的中继节点计算模块用于根据所述预测信噪比计算模块的计算结果,计算被选择的中继节点。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步地,所述预测CSI的中继选择系统,还包括发射功率计算模块,所述发射功率计算模块用于计算系统模型中源节点和所有中继节点的发射功率,并将计算结果发送给相应的节点,相应节点按照所计算的发射功率进行发射。
进一步地,所述预测信噪比计算模块的计算过程具体为:
记源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的预测信噪比为γai,则
其中为源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的测量信噪比,Δai为服从概率密度函数为的分布的随机数,
f Δ a i | γ ^ a i ( Δ a i | γ ^ a i ) = 1 γ ‾ a i ( 1 - ( ρ a i ) 2 ) exp ( - [ Δ a i + ( 1 + ( ρ a i ) 2 ) γ ^ a i ] γ ‾ a i ( 1 - ( ρ a i ) 2 ) ) × I 0 ( 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ‾ a i ( 1 - ( ρ a i ) 2 ) ) ;
其中表示源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的平均信噪比,ρai表示源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的预测信道系数和测量信道系数的相关系数, I 0 ( 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ‾ a i ( 1 - ( ρ a i ) 2 ) ) 表示自变量为 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ‾ a i ( 1 - ( ρ a i ) 2 ) 的第一类零阶变形贝塞尔函数。
进一步地,所述被选择的中继节点计算模块的计算过程具体为:
用i*表示被选择的中继节点,则
i * = arg m a x i ∈ S r e l a y { γ i } ;
其中i*∈Srelay,γi表示中继链路i的预测信噪比,γi=min{γ1i2i},γ1i、γ2i分别为a=1及a=2时的γai
进一步地,所述发射功率计算模块的计算过程具体为:
用P0表示源节点S的发射功率,P1,P2,...,PNrelay分别表示Nrelay个中继节点的发射功率,记发射功率集合为P,则
在特定的能量限制下最小化系统中断概率以分配源节点和所有中继节点的发射功率,即得出发射功率集合所述特定能量限制为:
s . t . P 0 + Σ i = 1 N r e l a y P i p s e l e c t i ( P ) ≤ E T ,
E min i ≤ P i p s e l e c t i ( P ) ≤ E m a x i ,
E min 0 ≤ P 0 ≤ E m a x 0 ;
其中, J ( P ) = Σ i = 1 N r e l a y F i ( γ T ) p s e l e c t i ( P ) ;
p s e l e c t i ( P ) = Π j = 1 j ≠ i N r e l a y p j i ;
p j i = p 1 i 2 i p 1 j 2 j p 1 j 1 i + p 2 i 1 i p 1 j 2 j p 1 j 2 i + p 1 i 2 i p 2 j 1 j p 2 j 1 i + p 2 i 1 i p 2 j 1 j p 2 j 2 i ;
p 1 i 2 i = exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 1 = 0 ∞ Σ k 2 = 0 ∞ Σ k 3 = 0 k 2 ( k 1 + k 3 ) ! ( k 1 ! ) 2 k 2 ! k 3 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 1 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 2 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 1 + k 3 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 1 - k 2 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 1 + k 3 ;
p 1 j 2 j = exp ( - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ ) Σ k 4 = 0 ∞ Σ k 5 = 0 ∞ Σ k 6 = 0 k 5 ( k 4 + k 6 ) ! ( k 4 ! ) 2 k 5 ! k 6 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 4 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 5 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 4 + k 6 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 4 - k 5 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 4 + k 6 ;
p 1 j 1 i = exp ( - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 7 = 0 ∞ Σ k 8 = 0 ∞ Σ k 9 = 0 k 8 ( k 7 + k 9 ) ! ( k 7 ! ) 2 k 8 ! k 9 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 7 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 8 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 7 + k 9 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 7 - k 8 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 7 + k 9 ;
p 2 i 1 i = exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 10 = 0 ∞ Σ k 11 = 0 ∞ Σ k 12 = 0 k 11 ( k 10 + k 12 ) ! ( k 10 ! ) 2 k 11 ! k 12 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 11 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 10 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 10 - k 11 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 10 + k 12 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 10 + k 12 ;
p 1 j 2 i = exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ ) Σ k 13 = 0 ∞ Σ k 14 = 0 ∞ Σ k 15 = 0 k 14 ( k 13 + k 15 ) ! ( k 13 ! ) 2 k 14 ! k 15 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 13 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 14 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 13 + k 15 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 13 - k 14 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 13 + k 15 ;
p 2 j 1 j = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ ) Σ k 16 = 0 ∞ Σ k 17 = 0 ∞ Σ k 18 = 0 k 17 ( k 16 + k 18 ) ! ( k 16 ! ) 2 k 17 ! k 18 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 17 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 16 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] 1 + k 16 - k 17 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 16 + k 18 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 16 + k 18 ;
p 2 j 1 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 19 = 0 ∞ Σ k 20 = 0 ∞ Σ k 21 = 0 k 20 ( k 19 + k 21 ) ! ( k 19 ! ) 2 k 20 ! k 21 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 20 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 19 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 19 - k 20 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 19 + k 21 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 19 + k 21 ;
p 2 j 2 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 22 = 0 ∞ Σ k 23 = 0 ∞ Σ k 24 = 0 k 23 ( k 22 + k 24 ) ! ( k 22 ! ) 2 k 23 ! k 24 ! × [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 23 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 22 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 22 - k 23 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 22 + k 24 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 22 + k 24 ;
F i ( γ T ) = exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 25 = 0 ∞ 1 ( k 25 ! ) 2 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 25 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 25 γ ( k 25 + 1 , γ T ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) + exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 26 = 0 ∞ 1 ( k 26 ! ) 2 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 26 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 26 γ ( k 26 + 1 , γ T ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) - exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 27 = 0 ∞ Σ k 28 = 0 ∞ 1 ( k 27 ! ) 2 ( k 28 ! ) 2 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 28 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 27 × [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 28 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 27 γ ( k 28 + 1 , γ T ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) γ ( k 27 + 1 , γ T ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) ;
其中Pi表示中继节点i的发射功率;为中继节点i被选择的概率;ET为归一化的可用的总能量;为归一化的可分配给中继节点i的最小能量,为归一化的可分配给中继节点i的最大能量,为归一化的可分配给源节点S的最小能量,为归一化的可分配给源节点S的最大能量;J(P)是系统中断概率;FiT)为中继链路i发生中断的概率,γT为中继链路发生中断的信噪比阈值;j∈Srelay,即j=1,2,...,Nrelay,且j≠i;为中继链路i的预测信噪比大于等于中继链路j的预测信噪比的概率;为源节点S到中继节点i之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于中继节点j到目的节点D之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率;为中继节点i到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率;为中继节点j到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点j之间链路的预测信噪比的概率,为中继节点j到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率;为中继节点j到目的节点D之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率;为分配功率前源节点S到中继节点i之间链路的平均信噪比,为源节点S到中继节点i之间链路的信道系数的方差,为中继节点i处源节点S的噪声方差,为分配功率前中继节点i到目的节点D之间链路的平均信噪比,为中继节点i到目的节点D之间链路的信道系数的方差,为目的节点D处中继节点i的噪声方差;为分配功率前源节点S到中继节点i之间链路的测量信噪比,为分配功率前中继节点i到目的节点D之间链路的测量信噪比;为分配功率前源节点S到中继节点j之间链路的平均信噪比,为源节点S到中继节点j之间链路的信道系数的方差,为中继节点j处源节点S的噪声方差,为分配功率前中继节点j到目的节点D之间链路的平均信噪比,Pj表示中继节点j的发射功率,为中继节点j到目的节点D之间链路的信道系数的方差,为目的节点D处中继节点j的噪声方差;为分配功率前源节点S到中继节点j之间链路的测量信噪比,为分配功率前中继节点j到目的节点D之间链路的测量信噪比;k1,k2,k3,……,k27,k28为求和下标指示,均为大于或等于0的整数;ρ1i为源节点S到中继节点i之间链路的预测信道系数和测量信道系数的相关系数,ρ1j为源节点S到中继节点j之间链路的预测信道系数和测量信道系数的相关系数,ρ2i为中继节点i到目的节点D之间链路的预测信道系数和测量信道系数的相关系数,ρ2j为中继节点j到目的节点D之间链路的预测信道系数和测量信道系数的相关系数;γ()表示不完全伽马函数。
本发明的有益效果是:本发明采用预测CSI进行中继选择,有效克服了过时CSI对系统性能的影响,具有科学意义和应用价值;本发明在测量信噪比基础上加上特定分布的随机数生成预测信噪比,预测精度比直接预测信噪比高,操作简单、计算量少、精度较高;本发明提出的考虑功率资源公平分配的功率分配方法,并用于中继选择方法中,为解决协作通信中的公平问题提供了有益的参考。
附图说明
图1为本发明实施例一所述的预测CSI的中继选择方法的流程图;
图2为本发明实施例二所述的预测CSI的中继选择方法的流程图;
图3为本发明实施例三所述的预测CSI的中继选择系统的结构图;
图4为本发明实施例四所述的预测CSI的中继选择系统的结构图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
图1为本发明所述实施例一的预测CSI的中继选择方法的流程图;
如图1所示,一种预测CSI的中继选择方法,包括以下步骤:
步骤S1,计算系统模型中所有链路的预测信噪比;所述系统模型包括一个源节点S、一个目的节点D和Nrelay个中继节点,其中Nrelay≥1,所述中继节点类型为DF(Decode-and-Forward,译码转发)类型,用Srelay表示中继节点集,则有Srelay={1,2,...,Nrelay},由于相距很远或被阻挡,所述目的节点接收不到来自所述源节点的信号,所述系统模型的信道为瑞利信道。记第i个中继节点为中继节点i,其中i∈Srelay,即i=1,2,...,Nrelay,所述链路包括源节点S到所有中继节点之间的链路和所有中继节点到目的节点D之间的链路,源节点S到中继节点i之间的链路和中继节点i到目的节点D之间的链路构成中继链路i。a∈{1,2},当a=1时表示源节点S到中继节点i之间的链路,当a=2时表示中继节点i到目的节点D之间的链路。
记源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的预测信噪比为γai,则
其中为源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的测量信噪比,Δai为服从概率密度函数(pdf)为的分布的随机数,
f Δ a i | γ ^ a i ( Δ a i | γ ^ a i ) = 1 γ ‾ a i ( 1 - ( ρ a i ) 2 ) exp ( - [ Δ a i + ( 1 + ( ρ a i ) 2 ) γ ^ a i ] γ ‾ a i ( 1 - ( ρ a i ) 2 ) ) × I 0 ( 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ‾ a i ( 1 - ( ρ a i ) 2 ) ) ;
其中表示源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的平均信噪比,ρai表示源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的预测信道系数和测量信道系数的相关系数, I 0 ( 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ‾ a i ( 1 - ( ρ a i ) 2 ) ) 表示自变量为 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ‾ a i ( 1 - ( ρ a i ) 2 ) 的第一类零阶变形贝塞尔函数;根据特定的分布产生随机数的方法已有现成算法,本发明不再赘述。
步骤S2,根据步骤S1的计算结果,计算被选择的中继节点;
用i*表示被选择的中继节点,则
i * = arg m a x i ∈ S r e l a y { γ i } ;
其中i*∈Srelay,γi表示中继链路i的预测信噪比,γi=min{γ1i2i},γ1i、γ2i分别为a=1及a=2时的γai,分别表示源节点S到中继节点i之间链路和中继节点i到目的节点D之间链路的预测信噪比。
图2为本发明所述实施例二的预测CSI的中继选择方法的流程图。
如图2所示,在实施例一所述预测CSI的中继选择方法的步骤S1之前,还可以包括步骤S0:计算系统模型中源节点和所有中继节点的发射功率;
用P0表示源节点S的发射功率,P1,P2,...,PNrelay分别表示Nrelay个中继节点的发射功率,记发射功率集合为P,则在特定的能量限制下最小化系统中断概率以分配源节点和所有中继节点的发射功率,即得出发射功率集合所述特定能量限制为:
s . t . P 0 + Σ i = 1 N r e l a y P i p s e l e c t i ( P ) ≤ E T ,
E min i ≤ P i p s e l e c t i ( P ) ≤ E m a x i ,
E min 0 ≤ P 0 ≤ E m a x 0 ;
其中, J ( P ) = Σ i = 1 N r e l a y F i ( γ T ) p s e l e c t i ( P ) ;
p s e l e c t i ( P ) = Π j = 1 j ≠ i N r e l a y p j i ;
p j i = p 1 i 2 i p 1 j 2 j p 1 j 1 i + p 2 i 1 i p 1 j 2 j p 1 j 2 i + p 1 i 2 i p 2 j 1 j p 2 j 1 i + p 2 i 1 i p 2 j 1 j p 2 j 2 i ;
p 1 i 2 i = exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 1 = 0 ∞ Σ k 2 = 0 ∞ Σ k 3 = 0 k 2 ( k 1 + k 3 ) ! ( k 1 ! ) 2 k 2 ! k 3 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 1 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 2 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 1 + k 3 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 1 - k 2 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 1 + k 3 ;
p 1 j 2 j = exp ( - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ ) Σ k 4 = 0 ∞ Σ k 5 = 0 ∞ Σ k 6 = 0 k 5 ( k 4 + k 6 ) ! ( k 4 ! ) 2 k 5 ! k 6 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 4 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 5 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 4 + k 6 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 4 - k 5 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 4 + k 6 ;
p 1 j 1 i = exp ( - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 7 = 0 ∞ Σ k 8 = 0 ∞ Σ k 9 = 0 k 8 ( k 7 + k 9 ) ! ( k 7 ! ) 2 k 8 ! k 9 !
× [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 7 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 8 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 7 + k 9 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 7 - k 8 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 7 + k 9 ;
p 2 i 1 i = exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 10 = 0 ∞ Σ k 11 = 0 ∞ Σ k 12 = 0 k 11 ( k 10 + k 12 ) ! ( k 10 ! ) 2 k 11 ! k 12 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 11 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 10 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 10 - k 11 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 10 + k 12 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 10 + k 12 ;
p 1 j 2 i = exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ ) Σ k 13 = 0 ∞ Σ k 14 = 0 ∞ Σ k 15 = 0 k 14 ( k 13 + k 15 ) ! ( k 13 ! ) 2 k 14 ! k 15 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 13 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 14 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 13 + k 15 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 13 - k 14 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 13 + k 15 ;
p 2 j 1 j = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ ) Σ k 16 = 0 ∞ Σ k 17 = 0 ∞ Σ k 18 = 0 k 17 ( k 16 + k 18 ) ! ( k 16 ! ) 2 k 17 ! k 18 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 17 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 16 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] 1 + k 16 - k 17 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 16 + k 18 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 16 + k 18 ;
p 2 j 1 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 19 = 0 ∞ Σ k 20 = 0 ∞ Σ k 21 = 0 k 20 ( k 19 + k 21 ) ! ( k 19 ! ) 2 k 20 ! k 21 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 20 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 19 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 19 - k 20 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 19 + k 21 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 19 + k 21 ;
p 2 j 2 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 22 = 0 ∞ Σ k 23 = 0 ∞ Σ k 24 = 0 k 23 ( k 22 + k 24 ) ! ( k 22 ! ) 2 k 23 ! k 24 ! × [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 23 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 22 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 22 - k 23 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 22 + k 24 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 22 + k 24 ;
F i ( γ T ) = exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 25 = 0 ∞ 1 ( k 25 ! ) 2 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 25 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 25 γ ( k 25 + 1 , γ T ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) + exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 26 = 0 ∞ 1 ( k 26 ! ) 2 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 26 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 26 γ ( k 26 + 1 , γ T ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ )
- exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 27 = 0 ∞ Σ k 28 = 0 ∞ 1 ( k 27 ! ) 2 ( k 28 ! ) 2 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 28 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 27 × [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 28 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 27 γ ( k 28 + 1 , γ T ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) γ ( k 27 + 1 , γ T ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) ;
其中Pi表示中继节点i的发射功率;为中继节点i被选择的概率;ET为归一化的可用的总能量;为归一化的可分配给中继节点i的最小能量,为归一化的可分配给中继节点i的最大能量,为归一化的可分配给源节点S的最小能量,为归一化的可分配给源节点S的最大能量;J(P)是系统中断概率;FiT)为中继链路i发生中断的概率,γT为中继链路发生中断的信噪比阈值;j∈Srelay,即j=1,2,...,Nrelay,且j≠i;为中继链路i的预测信噪比大于等于中继链路j的预测信噪比的概率;为源节点S到中继节点i之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于中继节点j到目的节点D之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率;为中继节点i到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率;为中继节点j到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点j之间链路的预测信噪比的概率,为中继节点j到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率;为中继节点j到目的节点D之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率;为分配功率前源节点S到中继节点i之间链路的平均信噪比,为源节点S到中继节点i之间链路的信道系数的方差,为中继节点i处源节点S的噪声方差,为分配功率前中继节点i到目的节点D之间链路的平均信噪比,为中继节点i到目的节点D之间链路的信道系数的方差,为目的节点D处中继节点i的噪声方差;为分配功率前源节点S到中继节点i之间链路的测量信噪比,为分配功率前中继节点i到目的节点D之间链路的测量信噪比;为分配功率前源节点S到中继节点j之间链路的平均信噪比,为源节点S到中继节点j之间链路的信道系数的方差,为中继节点j处源节点S的噪声方差,为分配功率前中继节点j到目的节点D之间链路的平均信噪比,Pj表示中继节点j的发射功率,为中继节点j到目的节点D之间链路的信道系数的方差,为目的节点D处中继节点j的噪声方差;为分配功率前源节点S到中继节点j之间链路的测量信噪比,为分配功率前中继节点j到目的节点D之间链路的测量信噪比;k1,k2,k3,……,k27,k28为求和下标指示,均为大于或等于0的整数;ρ1i为源节点S到中继节点i之间链路的预测信道系数和测量信道系数的相关系数,ρ1j为源节点S到中继节点j之间链路的预测信道系数和测量信道系数的相关系数,ρ2i为中继节点i到目的节点D之间链路的预测信道系数和测量信道系数的相关系数,ρ2j为中继节点j到目的节点D之间链路的预测信道系数和测量信道系数的相关系数;γ()表示不完全伽马函数。
在特定的能量限制下最小化系统中断概率可算出功率集合P,具体的解法可参照几何规划问题的解法,已有现成算法,此处不赘述。
在执行完步骤S0之后,再执行步骤S1时,在步骤S1中使用到的即为中继节点被分配发射功率后的源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的测量信噪比。
图3为本发明实施例三所述的预测CSI的中继选择系统的结构图。
将上述实施例一所述预测CSI的中继选择方法的技术特征一一对应为系统,即为如图3所示的一种预测CSI的中继选择系统,包括预测信噪比计算模块和被选择的中继节点计算模块;
所述预测信噪比计算模块用于计算系统模型中所有链路的预测信噪比;
所述被选择的中继节点计算模块用于根据所述预测信噪比计算模块的计算结果,计算被选择的中继节点。
预测信噪比计算模块和被选择的中继节点计算模块的具体计算过程同于上述实施例一所述预测CSI的中继选择方法中的具体计算过程,因此此处不再赘述。
图4为本发明实施例四所述的预测CSI的中继选择系统的结构图。
将上述实施例二所述预测CSI的中继选择方法的技术特征一一对应为系统,即为如图4所示的一种预测CSI的中继选择系统,包括发射功率计算模块、预测信噪比计算模块和被选择的中继节点计算模块;所述发射功率计算模块用于计算系统模型中源节点和所有中继节点的发射功率,并将计算结果发送给相应的节点,相应节点按照所计算的发射功率进行发射。发射功率计算模块、预测信噪比计算模块和被选择的中继节点计算模块的具体计算过程同于上述实施例二所述预测CSI的中继选择方法中的具体计算过程,因此此处不再赘述。
下述为一个具体实施例:
系统模型中有3个DF中继节点,即Nrelay=3。
一、采用实施例一所述的预测CSI的中继选择方法
步骤S1,计算系统模型中所有链路的预测信噪比γai
通过测量可得到如下:
γ ^ 11 = 13.5412209650588 , γ ^ 12 = 28.8749209459298 , γ ^ 13 = 14.0367748592063 ,
γ ^ 21 = 80.7307313375716 , γ ^ 22 = 7.09938020630797 , γ ^ 23 = 121578191357054 ;
通过测量可得到如下:
γ ‾ 11 = 22.0996361918168 , γ ‾ 12 = 24.7767266893373 , γ ‾ 13 = 23.6148620334904 ,
γ ‾ 21 = 75.7391602146676 , γ ‾ 22 = 11.0294895193293 , γ ‾ 23 = 6.75463244952510 ;
通过测量得到ρai如下:
ρ11=0.5,ρ12=0.55,ρ13=0.45,ρ21=0.65,ρ22=0.42,ρ23=0.65;
随机生成Δai如下,其中Δai为服从概率密度函数(pdf)为的分布的随机数;
Δ11=-1.55030879086126,Δ12=-2.3906676234756,Δ13=3.2764104792587,
Δ21=-5.2251844450911,Δ22=-1.20109726505910,Δ23=-1.82972045248100;
则根据可得γai如下:
γ11=11.9909121741975,γ12=26.4842533224542,γ13=17.3131853384650,
γ21=75.5055468924805,γ22=5.89828294124887,γ23=10.3280986832244;
步骤S2,根据步骤S1的计算结果,计算被选择的中继节点。
由步骤S1可得:
γ1=11.9909121741975,γ2=5.89828294124887,γ3=10.3280986832244,
中继节点1即为被选择的中继节点,用于转发信息。
二、采用实施例二所述的预测CSI的中继选择方法
步骤S0,,在特定的能量限制下最小化系统中断概率进行分配各节点的发射功率,得出发射功率集合P,
E T = 25 ; E min 0 = 1 , E min 1 = 1 , E min 2 = 1 , E m i n 3 = 1 ; E m a x 0 = 12.5 , E m a x 1 = 3.75 , E m a x 2 = 5 , E m a x 3 = 3.75 ;
通过测量可得到如下:
γ ^ 11 ′ = 20.7371499797968 , γ ^ 12 ′ = 85.8960116815784 , γ ^ 13 ′ = 31.0683391436942 ,
γ ^ 21 ′ = 39.4526226936800 , γ ^ 22 ′ = 16.4505763345343 , γ ^ 23 ′ = 14.7443336767864 ;
通过测量可得到i=1,2,3,如下:
通过测量可得到i=1,2,3,如下:
通过测量可得到i=1,2,3,如下:
σ 11 2 = 4.87885660949638 , σ 12 2 = 2.38683916951300 , σ 13 2 = 5.84848429575290 ;
通过测量可得到i=1,2,3,如下:
σ 21 2 = 3.25176228126154 , σ 22 2 = 4.61601898883682 , σ 23 2 = 5.42090755197626 ;
通过测量得到:
ρ11=0.5,ρ12=0.6,ρ13=0.5,ρ21=0.6,ρ22=0.5,ρ23=0.6;
通过以上数据解上述功率分配问题,可以得出发射功率集合P={P0,P1,P2,P3}={12.5,3.75,5,3.75}。
步骤S1,测量分配功率作用下的计算各条链路的预测信噪比γai
在分配功率作用下通过测量可得到如下:
γ ^ 11 = 14.1241485611585 , γ ^ 12 = 75.3943566377932 , γ ^ 13 = 32.8448797267051 ,
γ ^ 21 = 31.8661027325756 , γ ^ 22 = 10.4942806149574 , γ ^ 23 = 13.6673985657023 ;
通过测量可得到如下:
γ ‾ 11 = 33.3821381032940 , γ ‾ 12 = 107.057735784745 , γ ‾ 13 = 28.2036200065261 ,
γ ‾ 21 = 23.9159692106506 , γ ‾ 22 = 20.7566799928244 , γ ‾ 23 = 11.7918835294520 ;
通过测量得到ρai,如下:
ρ11=0.5,ρ12=0.55,ρ13=0.45,ρ21=0.65,ρ22=0.42,ρ23=0.65;
随机生成Δai如下,其中Δai为服从概率密度函数(pdf)为的分布的随机数;
Δ11=-1.77158333716048,Δ12=3.1130826764402,Δ13=2.3770957295203,
Δ21=-1.2337705644191,Δ22=-1.8567672610689,Δ23=2.09170378664851;
则根据可得γai如下:
γ11=12.3525652239981,γ12=78.5074393142334,γ13=35.2219754562254,
γ21=30.6323321681565,γ22=8.63751335388847,γ23=15.7591023523508;
步骤S2,根据步骤S1的计算结果,计算被选择的中继节点。
γ1=12.3525652239981,γ2=8.63751335388847,γ3=15.7591023523508;
中继节点为被选择的中继节点,用于转发信息。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种预测CSI的中继选择方法,其特征在于,包括以下步骤:
步骤S1,计算系统模型中所有链路的预测信噪比;
所述系统模型包括一个源节点S、一个目的节点D和Nrelay个中继节点,其中Nrelay≥1,记第i个中继节点为中继节点i,其中i=1,2,...,Nrelay,Srelay表示中继节点集,且Srelay={1,2,...,Nrelay},所述链路包括源节点S到所有中继节点之间的链路和所有中继节点到目的节点D之间的链路;源节点S到中继节点i之间的链路和中继节点i到目的节点D之间的链路构成中继链路i;a∈{1,2},当a=1时表示源节点S到中继节点i之间的链路,当a=2时表示中继节点i到目的节点D之间的链路;
步骤S2,根据步骤S1的计算结果,计算被选择的中继节点。
2.根据权利要求1所述的预测CSI的中继选择方法,其特征在于,在执行步骤S1之前,还包括步骤S0:计算系统模型中源节点和所有中继节点的发射功率。
3.根据权利要求1或2所述的预测CSI的中继选择方法,其特征在于,步骤S1中,记源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的预测信噪比为γai,则
其中为源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的测量信噪比,Δai为服从概率密度函数为的分布的随机数,
f Δ a i | γ ^ a i ( Δ a i | γ ^ a i ) = 1 γ ^ a i ( 1 - ( ρ a i ) 2 ) exp [ - [ Δ a i + ( 1 - ( ρ a i ) 2 ) γ ^ a i ] γ ^ a i ( 1 - ( ρ a i ) 2 ) ] × I 0 ( 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ^ a i ( 1 - ( ρ a i ) 2 ) ) ;
其中表示源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的平均信噪比,ρai表示源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的预测信道系数和测量信道系数的相关系数,表示自变量为的第一类零阶变形贝塞尔函数。
4.根据权利要求3所述的预测CSI的中继选择方法,其特征在于,步骤S2中,用i*表示被选择的中继节点,则
i * = arg m a x i ∈ S r e l a y { γ i } ;
其中i*∈Srelay,γi表示中继链路i的预测信噪比,γi=min{γ1i2i},γ1i、γ2i分别为a=1及a=2时的γai
5.根据权利要求2所述的预测CSI的中继选择方法,其特征在于,步骤S0中,用P0表示源节点S的发射功率,分别表示Nrelay个中继节点的发射功率,记发射功率集合为P,则
在特定的能量限制下最小化系统中断概率以分配源节点和所有中继节点的发射功率,即得出发射功率集合所述特定能量限制为:
s . t . P 0 + Σ i = 1 N r e l a y P i p s e l e c t i ( P ) ≤ E T ,
E m i n i ≤ P i p s e l e c t i ( P ) ≤ E m a x i ,
E m i n 0 ≤ P 0 ≤ E m a x 0 ;
其中, J ( P ) = Σ i = 1 N r e l a y F i ( γ T ) p s e l e c t i ( P ) ;
p s e l e c t i ( P ) = Π j = 1 j ≠ 1 N r e l a y p j i ;
p j i = p 1 i 2 i p 1 j 2 j p 1 j 1 i + p 2 i 1 i p 1 j 2 j p 1 j 2 i + p 1 i 2 i p 2 j 1 j p 2 j 1 i + p 2 i 1 i p 2 j 1 j p 2 j 2 i ;
p 1 i 2 i = exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2i ′ ) Σ k 1 = 0 ∞ Σ k 2 = 0 ∞ Σ k 3 = 0 k 2 ( k 1 + k 3 ) ! ( k 1 ! ) 2 k 2 ! k 3 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 1 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 2 [ γ ‾ 1i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 1 + k 3 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 1 - k 2 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 1 + k 3 ;
p 1 j 2 j = exp ( - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ ) Σ k 4 = 0 ∞ Σ k 5 = 0 ∞ Σ k 6 = 0 k 5 ( k 4 + k 6 ) ! ( k 4 ! ) 2 k 5 ! k 6 !
× [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 4 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 5 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 4 + k 6 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 4 - k 5 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 4 + k 6 ;
p 1 j 1 i = exp ( - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1i ′ ) Σ k 7 = 0 ∞ Σ k 8 = 0 ∞ Σ k 9 = 0 k 8 ( k 7 + k 9 ) ! ( k 7 ! ) 2 k 8 ! k 9 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 7 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 8 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 7 + k 9 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 7 - k 8 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 7 + k 9 ;
p 2 i 1 i = exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1i ′ ) Σ k 10 = 0 ∞ Σ k 11 = 0 ∞ Σ k 12 = 0 k 11 ( k 10 + k 12 ) ! ( k 10 ! ) 2 k 11 ! k 12 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 11 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 10 [ γ ‾ 1i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 10 - k 11 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 10 + k 12 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 10 + k 12 ;
p 1 j 2 i = exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1j ′ ) Σ k 13 = 0 ∞ Σ k 14 = 0 ∞ Σ k 15 = 0 k 14 ( k 13 + k 15 ) ! ( k 13 ! ) 2 k 14 ! k 15 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 13 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 14 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 13 + k 15 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 13 - k 14 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 13 + k 15 ;
p 2 j 1 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1j ′ ) Σ k 16 = 0 ∞ Σ k 17 = 0 ∞ Σ k 18 = 0 k 17 ( k 16 + k 18 ) ! ( k 16 ! ) 2 k 17 ! k 18 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 17 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 16 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] 1 + k 16 - k 17 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 16 + k 18 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 16 + k 18 ;
p 2 j 1 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1i ′ ) Σ k 19 = 0 ∞ Σ k 20 = 0 ∞ Σ k 21 = 0 k 20 ( k 19 + k 21 ) ! ( k 19 ! ) 2 k 20 ! k 21 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 20 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 19 [ γ ‾ 1i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 19 - k 20 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 19 + k 21 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 19 + k 21 ;
p 2 j 2 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 22 = 0 ∞ Σ k 23 = 0 ∞ Σ k 24 = 0 k 23 ( k 22 + k 24 ) ! ( k 22 ! ) 2 k 23 ! k 24 !
× [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 23 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 22 [ γ ‾ 2i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 22 - k 23 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 22 + k 24 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 22 + k 24 ;
F i ( γ T ) = exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 25 = 0 ∞ 1 ( k 25 ! ) 2 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 25 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 25 γ ( k 25 + 1 , γ T ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) + exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 26 = 0 ∞ 1 ( k 26 ! ) 2 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 26 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 26 γ ( k 26 + 1 , γ T ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) - exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 2 7 = 0 ∞ Σ k 2 8 = 0 ∞ 1 ( k 2 7 ! ) 2 ( k 2 8 ! ) 2 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 28 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 2 7 × [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 28 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 2 7 γ ( k 2 8 + 1 , γ T ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) γ ( k 2 7 + 1 , γ T ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) ;
其中Pi表示中继节点i的发射功率;为中继节点i被选择的概率;ET为归一化的可用的总能量;为归一化的可分配给中继节点i的最小能量,为归一化的可分配给中继节点i的最大能量,为归一化的可分配给源节点S的最小能量,为归一化的可分配给源节点S的最大能量;J(P)是系统中断概率;FiT)为中继链路i发生中断的概率,γT为中继链路发生中断的信噪比阈值;j∈Srelay,即j=1,2,...,Nrelay,且j≠i;为中继链路i的预测信噪比大于等于中继链路j的预测信噪比的概率;为源节点S到中继节点i之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于中继节点j到目的节点D之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率;为中继节点i到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率;为中继节点j到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点j之间链路的预测信噪比的概率,为中继节点j到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率;为中继节点j到目的节点D之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率;为分配功率前源节点S到中继节点i之间链路的平均信噪比, 为源节点S到中继节点i之间链路的信道系数的方差,为中继节点i处源节点S的噪声方差,为分配功率前中继节点i到目的节点D之间链路的平均信噪比, 为中继节点i到目的节点D之间链路的信道系数的方差,为目的节点D处中继节点i的噪声方差;为分配功率前源节点S到中继节点i之间链路的测量信噪比,为分配功率前中继节点i到目的节点D之间链路的测量信噪比;为分配功率前源节点S到中继节点j之间链路的平均信噪比, 为源节点S到中继节点j之间链路的信道系数的方差,为中继节点j处源节点S的噪声方差,为分配功率前中继节点j到目的节点D之间链路的平均信噪比,Pj表示中继节点j的发射功率,为中继节点j到目的节点D之间链路的信道系数的方差,为目的节点D处中继节点j的噪声方差;为分配功率前源节点S到中继节点j之间链路的测量信噪比,为分配功率前中继节点j到目的节点D之间链路的测量信噪比;k1,k2,k3,……,k27,k28为求和下标指示,均为大于或等于0的整数;ρ1i为源节点S到中继节点i之间链路的预测信道系数和测量信道系数的相关系数,ρ1j为源节点S到中继节点j之间链路的预测信道系数和测量信道系数的相关系数,ρ2i为中继节点i到目的节点D之间链路的预测信道系数和测量信道系数的相关系数,ρ2j为中继节点j到目的节点D之间链路的预测信道系数和测量信道系数的相关系数;γ()表示不完全伽马函数。
6.一种预测CSI的中继选择系统,其特征在于,包括预测信噪比计算模块和被选择的中继节点计算模块;
所述预测信噪比计算模块用于计算系统模型中所有链路的预测信噪比;
所述系统模型包括一个源节点S、一个目的节点D和Nrelay个中继节点,其中Nrelay≥1,记第i个中继节点为中继节点i,其中i=1,2,...,Nrelay,Srelay表示中继节点集,且Srelay={1,2,...,Nrelay},所述链路包括源节点S到所有中继节点之间的链路和所有中继节点到目的节点D之间的链路;源节点S到中继节点i之间的链路和中继节点i到目的节点D之间的链路构成中继链路i;a∈{1,2},当a=1时表示源节点S到中继节点i之间的链路,当a=2时表示中继节点i到目的节点D之间的链路;
所述被选择的中继节点计算模块用于根据所述预测信噪比计算模块的计算结果,计算被选择的中继节点。
7.根据权利要求6所述的预测CSI的中继选择系统,其特征在于,还包括发射功率计算模块,所述发射功率计算模块用于计算系统模型中源节点和所有中继节点的发射功率,并将计算结果发送给相应的节点,相应节点按照所计算的发射功率进行发射。
8.根据权利要求6或7所述的预测CSI的中继选择系统,其特征在于,所述预测信噪比计算模块的计算过程具体为:
记源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的预测信噪比为γai,则
其中为源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的测量信噪比,Δai为服从概率密度函数为的分布的随机数,
f Δ a i | γ ^ a i ( Δ a i | γ ^ a i ) = 1 γ ^ a i ( 1 - ( ρ a i ) 2 ) exp [ - [ Δ a i + ( 1 - ( ρ a i ) 2 ) γ ^ a i ] γ ^ a i ( 1 - ( ρ a i ) 2 ) ] × I 0 ( 2 ( ρ a i ) 2 γ ^ a i ( Δ a i + γ ^ a i ) γ ^ a i ( 1 - ( ρ a i ) 2 ) ) ;
其中表示源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的平均信噪比,ρai表示源节点S到中继节点i之间链路或中继节点i到目的节点D之间链路的预测信道系数和测量信道系数的相关系数,表示自变量为的第一类零阶变形贝塞尔函数。
9.根据权利要求8所述的预测CSI的中继选择系统,其特征在于,所述被选择的中继节点计算模块的计算过程具体为:
用i*表示被选择的中继节点,则
i * = arg m a x i ∈ S r e l a y { γ i } ;
其中i*∈Srelay,γi表示中继链路i的预测信噪比,γi=min{γ1i2i},γ1i、γ2i分别为a=1及a=2时的γai
10.根据权利要求7所述的预测CSI的中继选择系统,其特征在于,所述发射功率计算模块的计算过程具体为:
用P0表示源节点S的发射功率,分别表示Nrelay个中继节点的发射功率,记发射功率集合为P,则
在特定的能量限制下最小化系统中断概率以分配源节点和所有中继节点的发射功率,即得出发射功率集合所述特定能量限制为:
s . t . P 0 + Σ i = 1 N r e l a y P i p s e l e c t i ( P ) ≤ E T ,
E m i n i ≤ P i p s e l e c t i ( P ) ≤ E m a x i ,
E m i n 0 ≤ P 0 ≤ E m a x 0 ;
其中, J ( P ) = Σ i = 1 N r e l a y F i ( γ T ) p s e l e c t i ( P ) ;
p s e l e c t i ( P ) = Π j = 1 j ≠ i N r e l a y p j i ;
p j i = p 1 i 2 i p 1 j 2 j p 1 j 1 i + p 2 i 1 i p 1 j 2 j p 1 j 2 i + p 1 i 2 i p 2 j 1 j p 2 j 1 i + p 2 i 1 i p 2 j 1 j p 2 j 2 i ;
p 1 i 2 i = exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2i ′ ) Σ k 1 = 0 ∞ Σ k 2 = 0 ∞ Σ k 3 = 0 k 2 ( k 1 + k 3 ) ! ( k 1 ! ) 2 k 2 ! k 3 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 1 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 2 [ γ ‾ 1i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 1 + k 3 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 1 - k 2 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 1 + k 3 ;
p 1 j 2 j = exp ( - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ ) Σ k 4 = 0 ∞ Σ k 5 = 0 ∞ Σ k 6 = 0 k 5 ( k 4 + k 6 ) ! ( k 4 ! ) 2 k 5 ! k 6 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 4 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 5 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 4 + k 6 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 4 - k 5 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 4 + k 6 ;
p 1 j 1 i = exp ( - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1 j ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1i ′ ) Σ k 7 = 0 ∞ Σ k 8 = 0 ∞ Σ k 9 = 0 k 8 ( k 7 + k 9 ) ! ( k 7 ! ) 2 k 8 ! k 9 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 7 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 8 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 7 + k 9 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 7 - k 8 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 7 + k 9 ;
p 2 i 1 i = exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 10 = 0 ∞ Σ k 11 = 0 ∞ Σ k 12 = 0 k 11 ( k 10 + k 12 ) ! ( k 10 ! ) 2 k 11 ! k 12 !
× [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 11 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 10 [ γ ‾ 1i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 10 - k 11 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 10 + k 12 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 10 + k 12 ;
p 1 j 2 i = exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1j ′ ) Σ k 13 = 0 ∞ Σ k 14 = 0 ∞ Σ k 15 = 0 k 14 ( k 13 + k 15 ) ! ( k 13 ! ) 2 k 14 ! k 15 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 13 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 14 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] - k 13 + k 15 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 13 - k 14 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 13 + k 15 ;
p 2 j 1 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 1 j ) 2 γ ^ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) γ ‾ 1j ′ ) Σ k 16 = 0 ∞ Σ k 17 = 0 ∞ Σ k 18 = 0 k 17 ( k 16 + k 18 ) ! ( k 16 ! ) 2 k 17 ! k 18 ! × [ γ ^ 1 j ′ ( ρ 1 j ) 2 ] k 17 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 16 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) ] 1 + k 16 - k 17 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 16 + k 18 [ γ ‾ 1 j ′ ( 1 - ( ρ 1 j ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 16 + k 18 ;
p 2 j 1 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1i ′ ) Σ k 19 = 0 ∞ Σ k 20 = 0 ∞ Σ k 21 = 0 k 20 ( k 19 + k 21 ) ! ( k 19 ! ) 2 k 20 ! k 21 ! × [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 20 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 19 [ γ ‾ 1i ′ ( 1 - ( ρ 1 i ) 2 ) ] 1 + k 19 - k 20 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 19 + k 21 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 19 + k 21 ;
p 2 j 2 i = exp ( - ( ρ 2 j ) 2 γ ^ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) γ ‾ 2 j ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 22 = 0 ∞ Σ k 23 = 0 ∞ Σ k 24 = 0 k 23 ( k 22 + k 24 ) ! ( k 22 ! ) 2 k 23 ! k 24 ! × [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 23 [ γ ^ 2 j ′ ( ρ 2 j ) 2 ] k 22 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] 1 + k 22 - k 23 [ γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] - k 22 + k 24 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) + γ ‾ 2 j ′ ( 1 - ( ρ 2 j ) 2 ) ] 1 + k 22 + k 24 ;
F i ( γ T ) = exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) Σ k 25 = 0 ∞ 1 ( k 25 ! ) 2 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 25 [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 25 γ ( k 25 + 1 , γ T ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) + exp ( - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 26 = 0 ∞ 1 ( k 26 ! ) 2 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 26 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 26 γ ( k 26 + 1 , γ T ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) - exp ( - ( ρ 1 i ) 2 γ ^ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ - ( ρ 2 i ) 2 γ ^ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) Σ k 2 7 = 0 ∞ Σ k 2 8 = 0 ∞ 1 ( k 2 7 ! ) 2 ( k 2 8 ! ) 2 [ γ ^ 1 i ′ ( ρ 1 i ) 2 ] k 28 [ γ ^ 2 i ′ ( ρ 2 i ) 2 ] k 2 7 × [ γ ‾ 1 i ′ ( 1 - ( ρ 1 i ) 2 ) ] - k 28 [ γ ‾ 2 i ′ ( 1 - ( ρ 2 i ) 2 ) ] - k 2 7 γ ( k 2 8 + 1 , γ T ( 1 - ( ρ 1 i ) 2 ) γ ‾ 1 i ′ ) γ ( k 2 7 + 1 , γ T ( 1 - ( ρ 2 i ) 2 ) γ ‾ 2 i ′ ) ;
其中Pi表示中继节点i的发射功率;为中继节点i被选择的概率;ET为归一化的可用的总能量;为归一化的可分配给中继节点i的最小能量,为归一化的可分配给中继节点i的最大能量,为归一化的可分配给源节点S的最小能量,为归一化的可分配给源节点S的最大能量;J(P)是系统中断概率;FiT)为中继链路i发生中断的概率,γT为中继链路发生中断的信噪比阈值;j∈Srelay,即j=1,2,...,Nrelay,且j≠i;为中继链路i的预测信噪比大于等于中继链路j的预测信噪比的概率;为源节点S到中继节点i之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于中继节点j到目的节点D之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率;为中继节点i到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率,为源节点S到中继节点j之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率;为中继节点j到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点j之间链路的预测信噪比的概率,为中继节点j到目的节点D之间链路的预测信噪比小于等于源节点S到中继节点i之间链路的预测信噪比的概率;为中继节点j到目的节点D之间链路的预测信噪比小于等于中继节点i到目的节点D之间链路的预测信噪比的概率;为分配功率前源节点S到中继节点i之间链路的平均信噪比, 为源节点S到中继节点i之间链路的信道系数的方差,为中继节点i处源节点S的噪声方差,为分配功率前中继节点i到目的节点D之间链路的平均信噪比, 为中继节点i到目的节点D之间链路的信道系数的方差,为目的节点D处中继节点i的噪声方差;为分配功率前源节点S到中继节点i之间链路的测量信噪比,为分配功率前中继节点i到目的节点D之间链路的测量信噪比;为分配功率前源节点S到中继节点j之间链路的平均信噪比, 为源节点S到中继节点j之间链路的信道系数的方差,为中继节点j处源节点S的噪声方差,为分配功率前中继节点j到目的节点D之间链路的平均信噪比,Pj表示中继节点j的发射功率,为中继节点j到目的节点D之间链路的信道系数的方差,为目的节点D处中继节点j的噪声方差;为分配功率前源节点S到中继节点j之间链路的测量信噪比,为分配功率前中继节点j到目的节点D之间链路的测量信噪比;k1,k2,k3,……,k27,k28为求和下标指示,均为大于或等于0的整数;ρ1i为源节点S到中继节点i之间链路的预测信道系数和测量信道系数的相关系数,ρ1j为源节点S到中继节点j之间链路的预测信道系数和测量信道系数的相关系数,ρ2i为中继节点i到目的节点D之间链路的预测信道系数和测量信道系数的相关系数,ρ2j为中继节点j到目的节点D之间链路的预测信道系数和测量信道系数的相关系数;γ()表示不完全伽马函数。
CN201510581857.1A 2015-09-14 2015-09-14 一种预测csi的中继选择方法和系统 Expired - Fee Related CN105282814B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510581857.1A CN105282814B (zh) 2015-09-14 2015-09-14 一种预测csi的中继选择方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510581857.1A CN105282814B (zh) 2015-09-14 2015-09-14 一种预测csi的中继选择方法和系统

Publications (2)

Publication Number Publication Date
CN105282814A true CN105282814A (zh) 2016-01-27
CN105282814B CN105282814B (zh) 2019-02-26

Family

ID=55150950

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510581857.1A Expired - Fee Related CN105282814B (zh) 2015-09-14 2015-09-14 一种预测csi的中继选择方法和系统

Country Status (1)

Country Link
CN (1) CN105282814B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106060892A (zh) * 2016-07-07 2016-10-26 广西师范大学 一种基于过时信道状态信息的中继选择及功率分配方法
CN108271206A (zh) * 2017-12-08 2018-07-10 中山大学 一种基于阈值检测的中继节点选择方法
CN108600991A (zh) * 2018-05-03 2018-09-28 南通大学 车联网协作通信系统、路边单元的功率分配方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100098045A1 (en) * 2008-10-21 2010-04-22 Fujitsu Limited Apparatus, System, And Method For Communication
CN101932064A (zh) * 2010-07-27 2010-12-29 北京大学 一种双向中继网络中基于联合中继选择的通信方法
CN103281718A (zh) * 2013-06-07 2013-09-04 广西师范大学 移动通信系统中基于信道预测和公平性原则的中继节点选取方法
CN104105157A (zh) * 2014-08-17 2014-10-15 广西师范大学 协同通信系统中基于构造函数的中继选择方法
CN104754606A (zh) * 2015-03-26 2015-07-01 广西师范大学 一种基于信道预测的中继选择方法和系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100098045A1 (en) * 2008-10-21 2010-04-22 Fujitsu Limited Apparatus, System, And Method For Communication
CN101932064A (zh) * 2010-07-27 2010-12-29 北京大学 一种双向中继网络中基于联合中继选择的通信方法
CN103281718A (zh) * 2013-06-07 2013-09-04 广西师范大学 移动通信系统中基于信道预测和公平性原则的中继节点选取方法
CN104105157A (zh) * 2014-08-17 2014-10-15 广西师范大学 协同通信系统中基于构造函数的中继选择方法
CN104754606A (zh) * 2015-03-26 2015-07-01 广西师范大学 一种基于信道预测的中继选择方法和系统

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106060892A (zh) * 2016-07-07 2016-10-26 广西师范大学 一种基于过时信道状态信息的中继选择及功率分配方法
CN106060892B (zh) * 2016-07-07 2019-08-23 广西师范大学 一种基于过时信道状态信息的中继选择及功率分配方法
CN108271206A (zh) * 2017-12-08 2018-07-10 中山大学 一种基于阈值检测的中继节点选择方法
CN108271206B (zh) * 2017-12-08 2020-12-01 中山大学 一种基于阈值检测的中继节点选择方法
CN108600991A (zh) * 2018-05-03 2018-09-28 南通大学 车联网协作通信系统、路边单元的功率分配方法及系统
CN108600991B (zh) * 2018-05-03 2019-08-02 南通大学 车联网协作通信系统、路边单元的功率分配方法及系统

Also Published As

Publication number Publication date
CN105282814B (zh) 2019-02-26

Similar Documents

Publication Publication Date Title
Ding et al. Power allocation strategies in energy harvesting wireless cooperative networks
Prior et al. Network coding protocols for smart grid communications
Ko et al. Outage probability and channel capacity for the Nth best relay selection AF relaying over INID Rayleigh fading channels
Diamant et al. Robust spatial reuse scheduling in underwater acoustic communication networks
Ding et al. Multi-user SWIPT cooperative networks: Is the max–min criterion still diversity-optimal?
Stamatiou et al. Delay characterization of multihop transmission in a Poisson field of interference
Efazati et al. Effective capacity maximization in multirelay networks with a novel cross-layer transmission framework and power-allocation scheme
Ahmed et al. Relay subset selection and fair power allocation for best and partial relay selection in generic noise and interference
Beres et al. Outage probability of selection cooperation in the low to medium SNR regime
Hwang et al. Outage probability of cooperative diversity systems with opportunistic relaying based on decode-and-forward
CN105282814A (zh) 一种预测csi的中继选择方法和系统
CN104601473A (zh) 一种带约束的多目标路径的路由生成方法及系统
CN106332259A (zh) 一种双向中继网络增加安全能量效率的最优功率分配方法
Zhu et al. On latency reductions in vehicle-to-vehicle networks by random linear network coding
Mao et al. Joint link selection and relay power allocation for energy harvesting relaying systems
CN105050155A (zh) 基于csi预测的多用户中继选择方法
Dao et al. Uplink non-orthogonal multiple access protocol in two-way relaying networks: realistic operation and performance analysis
Islam et al. High-throughput link-channel selection and power allocation in wireless mesh networks
Liang et al. A game‐theoretic approach for relay assignment over distributed wireless networks
Zhang et al. Design and analysis of multi-relay selection for cooperative spatial multiplexing
Fedorenko et al. Criteria for the deployment of a heterogeneous linear WSN: Operability vs energy efficiency
CN101764678A (zh) 基于空间映射的中继方法和装置
Lee et al. Performance analysis of cooperative multi-hop networks with hop selection in Nakagami-m fading channels
Osman Empowering internet-of-everything (IoE) networks through synergizing Lagrange optimization and deep learning for enhanced performance
Bithas et al. Transmit antenna selection in vehicle-to-vehicle time-varying fading channels

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190226

Termination date: 20210914

CF01 Termination of patent right due to non-payment of annual fee