CN103458521A - 一种基于鲁棒性设计的mimo传输功率分配优化方法 - Google Patents

一种基于鲁棒性设计的mimo传输功率分配优化方法 Download PDF

Info

Publication number
CN103458521A
CN103458521A CN2013104071730A CN201310407173A CN103458521A CN 103458521 A CN103458521 A CN 103458521A CN 2013104071730 A CN2013104071730 A CN 2013104071730A CN 201310407173 A CN201310407173 A CN 201310407173A CN 103458521 A CN103458521 A CN 103458521A
Authority
CN
China
Prior art keywords
tau
power
sigma
channel
gamma
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
CN2013104071730A
Other languages
English (en)
Other versions
CN103458521B (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201310407173.0A priority Critical patent/CN103458521B/zh
Publication of CN103458521A publication Critical patent/CN103458521A/zh
Application granted granted Critical
Publication of CN103458521B publication Critical patent/CN103458521B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Mobile Radio Communication Systems (AREA)

Abstract

本发明提出了一种基于鲁棒性设计的MIMO传输功率分配优化方法,首先利用矩阵函数凸凹性,将功率分配矩阵优化为对角阵,化简为标量函数,随后结合标量函数进行MIMO传输功率分配优化。在MIMO传输功率分配优化过程中,首先找到无限制条件的最优解,随后进行微调,保证算法的收敛性和低运算量。本发明提供的基于鲁棒性设计的MIMO传输功率分配优化方法,相对于普适的最优解获取算法,能够以相对较小的运算代价来获取最优解,并且在第一步就己将功率分配的基准值固定,不存在算法发散的问题。

Description

一种基于鲁棒性设计的MIMO传输功率分配优化方法
技术领域
本发明涉及一种在信道信息不完整时鲁棒性MIMO传输的具体功率分配方法,属于多天线干扰信道中从不完美的信道信息中寻找功率分配策略的问题。
背景技术
众所周知,MIMO(Multiple-Input Multiple-Out-put,多输入多输出)系统的性能很大程度上取决于信道状态信息(Channel State Information,简称CSI)。为实现MIMO信道的全部效益,可以获取在发送端的信道状态信息(Channel State Information at thetransmitter,简称CSIT)并采用恰当的预编码技术。如果能够获取完整的CSIT,在这种情况下已经有成熟的MIMO优化预编码技术。然而,实际上CSIT一般不完整或者有误差,这主要是由不精确地估计和量化信道信息,错误的或者迟滞的反馈信息,以及互信道延时和频偏。因此我们在研究预编码技术的时候要考虑到CSIT的不完整性。我们的目标是充分利用CSIT同时对不同的CSIT缺陷具有普适性或鲁棒性。
通常,CSIT缺陷一般用随机模型或者确定性模型来刻画。随机性模型指的是信道信息是一个随机量。它的即时信息未知但是符合统计规律,比如可以获取均值和协方差。在这种情况下,鲁棒性设计致力于优化长期的平均性能或者中断回复性能。确定性模型认为信道处在一个不确定域或者不确定集合中。它更适合刻画即时CSI的误差。这个集合的大小与信道的确定性密切相关,集合越大,信道的确定性越大。这样,我们可以把这样的预编码设计称为是鲁棒的:在不确定域内在最差信道获取最好性能,也就是在不确定域内对任何信道保证性能的不低于门限。实现这样的预编码设计可以抽象为优化最坏情况的性能,也就是解决一个最大最小(maximin)或者最小最大(minimax)问题。
一个经典的窄带MIMO信道的模型可以这样描述:
y=Hx+n
其中:x为N维发送向量,x∈CN,即有N个发射天线;y为M维接收向量,y∈CM,即有M个接收天线;H为M×N维信道矩阵,H∈CM×N;n为M维噪声向量,n∈CM,且服从高斯分布,即
Figure BDA0000378993870000011
Figure BDA0000378993870000012
为噪声功率;Q=E[xxH](对xxH求数学期望)为功率分配矩阵,Q可以分解为Q=FFH(F为传送符号的线性映射矩阵);定义传输符号的向量s,E[ssH]=I,那么,x=Fs;事实上,对于Q一般应满足Q∈Q(Q为Q所有可能取到的值的集合),其中:
Q={Q|Q≥0,λi(Q)≤Pi,tr(Q)≤Ps,i=1,2,…,N|}
其中:λi(Q)为矩阵Q的特征值且降序排列,即λ1(Q)≥λ2(Q)≥…≥λN(Q);tr(Q)为矩阵Q的迹;Pi为信道的功率上限值;Ps为总功率限制值;λi(Q)≤Pi是指第i个数据流上有功率限制Pi,tr(Q)≤Ps是指所有数据流有总功率限制Ps
在完整信道信息的情况下,信道矩阵H总是能够在接收端无误差地获取;然而实际情况是CSIT有缺陷,这就需要预编码设计能够利用CSIT并且能够对抗缺陷;可以假设H属于一个集合H,这个集合称为信道可行域;现在引入性能指标的目标函数Ψ(Q,H),即:
Figure BDA0000378993870000021
这是个最大最小(maxmin)问题。
在这里,我们采用接收信噪比作为性能的衡量指标。从上面的分析来看,接收信噪比的表达式为:
Figure BDA0000378993870000022
注意到
Figure BDA0000378993870000023
只是一个常数,那么接收信噪比正比于tr(HQHH)(接受信号功率,亦即归一化的接受信噪比),因此我们可以采用这样的性能函数:Ψ(Q,H)=tr(HQHH)。由于信道具有不确定性,定义
Figure BDA0000378993870000024
那么 Ψ ( Q , Δ ) = tr ( ( H ^ - Δ ) Q ( H ^ - Δ ) H ) , 所以目标函数可以化为:
max Q ∈ Q min Δ ∈ ϵ tr ( ( H ^ - Δ ) Q ( H ^ - Δ ) H )
其中:
Figure BDA0000378993870000027
为标称信道矩阵,Δ为实际信道矩阵H与标称信道矩阵
Figure BDA0000378993870000028
的差。为了将问题简化,保留其本质,我们不加证明地给出以下结论,当矩阵函数取最优解时:
Q = U q Λ q U q H = U h Λ q U h H
( H ^ - Δ ) H ( H ^ - Δ ) = U h ( Λ h - Λ δ ) H ( Λ h - Λ δ ) U h H
可以证明,这样的简化不会丢失最优解,不会影响功率分配策略,仅仅是去掉了关于信道的一些信息,再根据迹的酉不变性质有:
tr ( ( H ^ - Δ ) Q ( H ^ - Δ ) H ) = Σ i = 1 N ( γ i - δ i ) 2 p i ( γ i ≥ δ i ≥ 0 )
其中:Λq为Q的特征值pi所构成的对角阵,且pi按降序排列,即p1≥p2≥…≥pN;Λh为H的奇异值γi所构成的对角阵,且γi按降序排列,即γ1≥γ2≥…≥γN;Λδ为Δ的奇异值δi所构成的对角阵,且δi按降序排列,即δ1≥δ2≥…≥δN;Uq为功率分配矩阵特征值分解后的酉矩阵;
Figure BDA0000378993870000031
为功率分配矩阵特征值分解后的酉矩阵的共轭转置;Uh为信道矩阵矩阵特征值分解后的酉矩阵;
Figure BDA0000378993870000032
为信道矩阵矩阵特征值分解后的酉矩阵的共轭转置。
为了进一步求解,我们选择将集合不确定域ε具体化。这里我们选择平方误差(这也是常用的一种误差刻画方式),但是在这里,我们可以对平方误差做一个推广,其一般形式记为:
Figure BDA0000378993870000033
ε为允许误差,我们称之为加权平方误差。借助拉格朗日函数,求得此时的δi为:
Figure BDA0000378993870000034
优化问题转化为:
max { p i } , μ Σ i = 1 N μ τ i γ i 2 p i μ τ i + p i - μ ϵ 2
s . t . Σ i = 1 N p i = P s
0≤pi≤Pi(i=1,2,…,N)
μ>0
其中:τi为不确定域的加权系数,μ为拉格朗日乘子。由于前文提到的标称信道矩阵
Figure BDA0000378993870000037
不一定满秩,这里假设信道数目
Figure BDA0000378993870000038
当i>r时,γi=0,此时应该使pi=0以实现最优分配,所以优化问题修正为:
max { p i } , μ Σ i = 1 r μ τ i γ i 2 p i μ τ i + p i - μ ϵ 2
s . t . Σ i = 1 r p i = P s
0≤pi≤pi(i=1,2,…,r)
μ>0
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种在信道信息不完整时鲁棒性MIMO传输的具体功率分配方法,利用矩阵函数凸凹性,将功率矩阵优化为对角阵,化简为标量函数,最后结合标量函数给出算法,以获得功率分配的最优解。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种基于鲁棒性设计的MIMO传输功率分配优化方法,首先利用矩阵函数凹凸性,将功率矩阵优化为对角阵,化简为标量函数,最后结合标量函数进行MIMO传输功率的分配优化,在MIMO传输功率的分配优化过程中,首先找到无限制条件的最优解,随后进行微调,保证算法的收敛性和低运算量,具体包括如下步骤:
(1)将μ作为参数,对pi进行优化,即:
max { p i } , μ Σ i = 1 r μτ i γ i 2 p i μτ i + p i - μϵ 2 = max μ h ( μ )
经过计算:
p i = γ i τ i P s Σ i = 1 r γ i τ i + μτ i γ i Σ i = 1 r τ i - Σ i = 1 r τ i γ i Σ i = 1 r γ i τ i
h ( μ ) = Σ i = 1 r ( γ i 2 μ τ i ) - ( Σ i = 1 r γ i μ τ i ) 2 Σ i = 1 r ( μ τ i ) + P s - μϵ 2
其中:r为信道数目,
Figure BDA0000378993870000044
Figure BDA0000378993870000045
为标称信道矩阵;pi为功率分配矩阵Q的特征值,且pi按降序排列,即p1≥p2≥…≥pN;μ为拉格朗日乘子;τi为不确定域ε的加权系数;H为信道矩阵,γi为H的奇异值,且γi按降序排列,即γ1≥γ2≥…≥γN;Ps为所有数据流的总功率限制值。
(2)根据μ的限制条件,优化h(μ),即:
max μ h ( μ ) = h ( μ * )
其中μ*为使h(μ)取最大值时的拉格朗日乘子值;
定义 A = Σ i = 1 r τ i , B = Σ i = 1 r γ i τ i , C = Σ i = 1 r γ i 2 τ i , 则:
h ( μ ) = Cμ - μ 2 B 2 Aμ + P s - μϵ 2
对上式进行求导得到:
μ * = P s A ( B 2 B 2 - ( C - ϵ 2 ) A - 1 )
h ( μ * ) = P s A 2 ( B - B 2 - ( C - ϵ 2 ) A ) 2
(3)将μ*代入pi得到
Figure BDA0000378993870000056
验证
Figure BDA0000378993870000057
是否超过功率上下限:
若有k1个信道突破功率上限,称这类突破功率上限的信道为一型信道,并且这k1个信道所构成的集合称为K1;若有k2个信道突破功率下限,称这类突破功率下限的信道为二型信道,并且这k2个信道所构成的集合称为K2;剩下的信道既不突破功率上限也不突破功率下限,称这类信道为三型信道,这k3=r-k1-k2个信道所构成的集合称为K3;由假设容易想到,K1、K2、K3两两的交集为φ,三者的并集为全集(即所有的信道);对于所有的一型信道,
Figure BDA0000378993870000058
对于所有的二型信道,
Figure BDA0000378993870000059
对于所有的三型信道,
Figure BDA00003789938700000510
等于计算值,其中Pi为第i个数据流上的功率限制值;分两种情况进行讨论:
①如果经过计算所有的信道都是三型信道,那么输出结果
Figure BDA00003789938700000511
作为最终功率的功率分配方案,结束;
②如果经过计算存在一型信道或二型信道,那么一型信道分配的功率为功率上限,二型信道分配的功率为功率下限,三型信道分配的功率为计算值,将每一个信道的功率分配方案
Figure BDA00003789938700000512
保存,进入步骤(4);
(4)根据全局最优解进行剩余功率的再分配,记剩余功率
Figure BDA0000378993870000061
对Δp进行功率再分配:
Δ p i = γ i τ i Δp Σ i ∈ K 3 γ i τ i + μ τ i γ i Σ i ∈ K 3 τ i - Σ i ∈ K 3 τ i γ i Σ i ∈ K 3 γ i τ i
并且:
μ = Δp A ( B 2 B 2 - ( C - ϵ 2 ) A - 1 ) , A = Σ i ∈ K 3 τ i , B = Σ i ∈ K 3 γ i τ i
对于获得的Δpi进行验证:
①如果存在 p i * + &Delta; p i > P i 或者 p i * + &Delta; p i < 0 , 更新 p i * = p i * + &Delta; p i 或者 p i * = P i 或者 p i * = 0 , 再重新执行步骤(4);
②如果所有的
Figure BDA0000378993870000069
满足
Figure BDA00003789938700000610
那么确定最终的功率分配方案为: p i * * = p i * + &Delta; p i , 结束。
本算法的核心思想是先定调后微调,这样确保了算法的收敛,第一步分配的功率是Ps,第二部分配的功率是Δp……Δp,在数值上比Ps小很多,因此该算法一定是收敛的;而且可以人为控制收敛速度,比如当Δp<10-4时停止算法;该算法的复杂度为o(n),n为信道数量。
有益效果:本发明提供的基于鲁棒性设计的MIMO传输功率分配优化方法,相对于普适的最优解获取,能够以相对较小的运算代价来获取最优解,并且在第一步就已将功率分配的基准值固定不存在收敛问题。
附图说明
图1为本发明的流程框图;
图2为接收信噪比与总功率Ps的关系曲线;
图3为不同功率上限,接收信噪比与总功率Ps的关系曲线;
图4为接收信噪比与功率上限Pi的关系曲线;
图5为不同总功率,接收信噪比与功率上限Pi的关系曲线;
图6为接收信噪比与允许误差ε的关系曲线;
图7为接收信噪比与信道数目r的关系曲线。
具体实施方式
下面结合附图对本发明作更进一步的说明。
鲁棒性设计的目的是使信道的最小信噪比最大化。这是一个最大最小问题,其中信道的信息和功率分配的情况均以矩阵函数的形式体现。在优化过程中,我们利用矩阵函数的凸凹性,逐层优化。为不失一股性,我们可以证明,将功率分配矩阵按照信道矩阵的方式对角化时,可以取得最优解。本发明方法是基于这个结论而产生的。本方法的研究对象为对角化矩阵函数化简得的标量函数。针对这个标量函数,我们给出以下的信道功率分配方案,从而获得功率分配的最优解。由于普适的获取最优解需要很大的代价,本算法的优势在于能够以相对较小的运算代价来获取最优解。优化的目的是找到最佳的功率分配方案pi和μ值,使得目标函数h(p,μ)(p为分配的信道功率的向量)取最大值。
基于鲁棒性设计的MIMO传输功率分配优化方法首先利用矩阵函数凸凹性,将功率矩阵优化为对角阵,化简为标量函数,最后结合标量函数进行MIMO传输功率的分配优化,在MIMO传输功率的分配优化过程中,首先找到无限制条件的最优解,随后进行微调,保证算法的收敛性和低运算量,如图1所示,具体包括如下步骤:
(1)将μ作为参数,对pi进行优化,即:
max { p i } , &mu; &Sigma; i = 1 r &mu; &tau; i &gamma; i 2 p i &mu;&tau; i + p i - &mu;&epsiv; 2 = max &mu; h ( &mu; )
经过计算:
p i = &gamma; i &tau; i P s &Sigma; i = 1 r &gamma; i &tau; i + &mu;&tau; i &gamma; i &Sigma; i = 1 r &tau; i - &Sigma; i = 1 r &tau; i &gamma; i &Sigma; i = 1 r &gamma; i &tau; i
h ( &mu; ) = &Sigma; i = 1 r ( &gamma; i 2 &mu;&tau; i ) - ( &Sigma; i = 1 r &gamma; i &mu;&tau; i ) 2 &Sigma; i = 1 r ( &mu;&tau; i ) + P s - &mu;&epsiv; 2
其中:r为信道数目,
Figure BDA0000378993870000074
Figure BDA0000378993870000075
为标称信道矩阵;pi为功率分配矩阵Q的特征值,且pi按降序排列,即p1≥p2≥…≥PN;μ为拉格朗日乘子;τi为不确定域ε的加权系数;H为信道矩阵,γi为H的奇异值,且γi按降序排列,即γ1≥γ2≥…≥γN;Ps为所有数据流的总功率限制值;
(2)根据μ的限制条件,优化h(μ),即:
max &mu; h ( &mu; ) = h ( &mu; * )
其中μ*使h(μ)取最大值时的拉格朗日乘子值;
定义 A = &Sigma; i = 1 r &tau; i , B = &Sigma; i = 1 r &gamma; i &tau; i , C = &Sigma; i = 1 r &gamma; i 2 &tau; i , 则:
h ( &mu; ) = C&mu; - &mu; 2 B 2 A&mu; + P s - &mu;&epsiv; 2
对上式进行求导得到:
&mu; * = P s A ( B 2 B 2 - ( C - &epsiv; 2 ) A - 1 )
h ( &mu; * ) = P s A 2 ( B - B 2 - ( C - &epsiv; 2 ) A ) 2
(3)将μ*代入pi得到
Figure BDA0000378993870000087
验证是否超过功率上下限:
若有k1个信道突破功率上限,称这类突破功率上限的信道为一型信道,并且这k1个信道所构成的集合称为k1;若有k2个信道突破功率下限,称这类突破功率下限的信道为二型信道,并且这k2个信道所构成的集合称为k2;剩下的信道既不突破功率上限也不突破功率下限,称这类信道为三型信道,这k3=r-k1-k2个信道所构成的集合称为K3;由假设容易想到,K1、K2、K3两两的交集为φ,三者的并集为全集(即所有的信道);对于所有的一型信道,
Figure BDA0000378993870000089
对于所有的二型信道,
Figure BDA00003789938700000810
对于所有的三型信道,
Figure BDA00003789938700000811
等于计算值,其中Pi为第i个数据流上的功率限制值;分两种情况进行讨论:
①如果经过计算所有的信道都是三型信道,那么输出结果
Figure BDA00003789938700000812
作为最终功率的功率分配方案,结束;
②如果经过计算存在一型信道或二型信道,那么一型信道分配的功率为功率上限,二型信道分配的功率为功率下限,三型信道分配的功率为计算值,将每一个信道的功率分配方案
Figure BDA0000378993870000091
保存,进入步骤(4);
(4)根据全局最优解进行剩余功率的再分配,记剩余功率
Figure BDA00003789938700000915
对Δp进行功率再分配:
&Delta;p i = &gamma; i &tau; i &Delta;p &Sigma; i &Element; K 3 &gamma; i &tau; i + &mu;&tau; i &gamma; i &Sigma; i &Element; K 3 &tau; i - &Sigma; i &Element; K 3 &tau; i &gamma; i &Sigma; i &Element; K 3 &gamma; i &tau; i
并且:
&mu; = &Delta;p A ( B 2 B 2 - ( C - &epsiv; 2 ) A - 1 ) , A = &Sigma; i &Element; K 3 &tau; i , B = &Sigma; i &Element; K 3 &gamma; i &tau; i
对于获得的Δpi进行验证:
①如果存在 p i * + &Delta;p i > P i 或者 p i * + &Delta;p i < 0 , 更新 p i * = p i * + &Delta;p i 或者 p i * = P i 或者 p i * = 0 , 再重新执行步骤(4);
②如果所有的
Figure BDA00003789938700000910
满足
Figure BDA00003789938700000916
那么确定最终的功率分配方案为: p i * * = p i * + &Delta;p i 结束。
图2为在r=4,τ1234=1,γ1=4,γ2=3,γ3=2,γ4=1,ε=4,P1=P2=P3=P4=20dB的情况下,接收信噪比与总功率Ps的关系曲线;其中Pi,Ps均为对噪声
Figure BDA00003789938700000913
归一化的值,用dB;总功率越大,接收信噪比越高;此外,图中显示鲁棒性MIMO传输功率分配策略的接收信噪比优于非鲁棒性的设计(这里使用功率均匀分配的策略)。
图3为r=4,τ1234=1,γ1=4,γ2=3,γ3=2,γ4=1,ε=4,P1=P2=P3=P4=15dB,20dB,25dB的情况下,接收信噪比与总功率Ps的关系曲线;其中Pi,Ps均为对噪声
Figure BDA00003789938700000914
归一化的值,用dB;曲线的拐点为接收信噪比受到各信道功率限制的门限值;
图4为在为r=4,τ1=τ2=τ3=τ4=1,γ1=4,γ2=3,γ3=2,γ4=1,ε=4,Ps=20dB的情况下,接收信噪比与功率上限Pi的关系曲线;其中Pi,Ps均为对噪声
Figure BDA0000378993870000101
归一化的值,用dB;此外,图中显示鲁棒性MIMO传输功率分配策略的接收信噪比优于非鲁棒性的设计(这里使用功率均匀分配的策略)。
图5为在r=4,τ1=τ2=τ3=τ4=1,γ1=4,γ2=3,γ3=2,γ4=1,ε=4,Ps=15dB,20dB,25dB的情况下,接收信噪比与功率上限Pi的关系曲线;其中Pi,Ps均为对噪声
Figure BDA0000378993870000102
归一化的值,用dB;曲线的线性部分可以解释为接收信噪比主要由各信道功率上限决定,线性增加;曲线的水平部分可以解释为接收信噪比受到总功率的限制而不再提高。
图6为在r=4,τ1=τ2=τ3=τ4=1,γ1=4,γ2=3,γ3=2,γ4=1,P1=P2=P3=P4=15dB,ps=18dB的情况下,接收信噪比与允许误差ε的关系曲线,其中Pi,Ps均为对噪声
Figure BDA0000378993870000103
归一化的值,用dB;允许的误差越大,接收信噪比越低。
图7为在τi=1,γi=1,Pi=15dB,Ps=18dB的情况下,接收信噪比与信道数目r的关系曲线,其中Pi,Ps均为对噪声
Figure BDA0000378993870000104
归一化的值,用dB;相同的信道数目越多,接收信噪比越差,这是因为信道分配的功率较原来少的缘故。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (1)

1.一种基于鲁棒性设计的MIMO传输功率分配优化方法,其特征在于:包括如下步骤:
(1)将μ作为参数,对pi进行优化,即:
max { p i } , &mu; &Sigma; i = 1 r &mu; &tau; i &gamma; i 2 p i &mu;&tau; i + p i - &mu;&epsiv; 2 = max &mu; h ( &mu; )
经过计算:
p i = &gamma; i &tau; i P s &Sigma; i = 1 r &gamma; i &tau; i + &mu;&tau; i &gamma; i &Sigma; i = 1 r &tau; i - &Sigma; i = 1 r &tau; i &gamma; i &Sigma; i = 1 r &gamma; i &tau; i
h ( &mu; ) = &Sigma; i = 1 r ( &gamma; i 2 &mu;&tau; i ) - ( &Sigma; i = 1 r &gamma; i &mu;&tau; i ) 2 &Sigma; i = 1 r ( &mu;&tau; i ) + P s - &mu;&epsiv; 2
其中:r为信道数目,为标称信道矩阵;pi为功率分配矩阵Q的特征值,且pi按降序排列,即p1≥p2≥…≥pN;μ为拉格朗日乘子;τi为不确定域ε的加权系数;H为信道矩阵,γi为H的奇异值,且γi按降序排列,即γ1≥γ2≥…≥γN;Ps为所有数据流的总功率限制值;ε为允许误差;
(2)根据μ的限制条件,优化h(μ),即:
max &mu; h ( &mu; ) = h ( &mu; * )
其中μ*为使h(μ)取最大值时的拉格朗日乘子值;
定义 A = &Sigma; i = 1 r &tau; i , B = &Sigma; i = 1 r &gamma; i &tau; i , C = &Sigma; i = 1 r &gamma; i 2 &tau; i , 则:
h ( &mu; ) = C&mu; - &mu; 2 B 2 A&mu; + P s - &mu;&epsiv; 2
对上式进行求导得到:
&mu; * = P s A ( B 2 B 2 - ( C - &epsiv; 2 ) A - 1 )
h ( &mu; * ) = P s A 2 ( B - B 2 - ( C - &epsiv; 2 ) A ) 2
(3)将μ*代入pi得到
Figure FDA0000378993860000022
验证
Figure FDA0000378993860000023
是否超过功率上下限:
若有k1个信道突破功率上限,称这类突破功率上限的信道为一型信道,并且这k1个信道所构成的集合称为K1;若有k2个信道突破功率下限,称这类突破功率下限的信道为二型信道,并月.这k2个信道所构成的集合称为K2;剩下的信道既不突破功率上限也不突破功率下限,称这类信道为三型信道,这k3=r-k1-k2个信道所构成的集合称为K3;对于所有的一型信道,
Figure FDA0000378993860000024
对于所有的二型信道,
Figure FDA0000378993860000025
对于所有的三型信道,
Figure FDA0000378993860000026
等于计算值,其中Pi为第i个数据流上的功率限制值;分两种情况进行讨论:
①如果经过计算所有的信道都是三型信道,那么输出结果
Figure FDA0000378993860000027
作为最终功率的功率分配方案,结束;
②如果经过计算存在一型信道或二型信道,那么一型信道分配的功率为功率上限,二型信道分配的功率为功率下限,三型信道分配的功率为计算值,将每一个信道的功率分配方案
Figure FDA0000378993860000028
保存,进入步骤(4);
(4)根据全局最优解进行剩余功率的再分配,记剩余功率
Figure FDA0000378993860000029
对Δp进行功率再分配:
&Delta;p i = &gamma; i &tau; i &Delta;p &Sigma; i &Element; K 3 &gamma; i &tau; i + &mu;&tau; i &gamma; i &Sigma; i &Element; K 3 &tau; i - &Sigma; i &Element; K 3 &tau; i &gamma; i &Sigma; i &Element; K 3 &gamma; i &tau; i
并且:
&mu; = &Delta;p A ( B 2 B 2 - ( C - &epsiv; 2 ) A - 1 ) , A = &Sigma; i &Element; K 3 &tau; i , B = &Sigma; i &Element; K 3 &gamma; i &tau; i
对于获得的Δpi进行验证:
①如果存在 p i * + &Delta;p i > P i 或者 p i * + &Delta; p i < 0 , 更新 p i * = p i * + &Delta;p i 或者 p i * = P i 或者 p i * = 0 , 再重新执行步骤(4);
②如果所有的满足
Figure FDA0000378993860000032
那么确定最终的功率分配方案为: p i * * = p i * + &Delta;p i , 结束。
CN201310407173.0A 2013-09-09 2013-09-09 一种基于鲁棒性设计的mimo传输功率分配优化方法 Expired - Fee Related CN103458521B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310407173.0A CN103458521B (zh) 2013-09-09 2013-09-09 一种基于鲁棒性设计的mimo传输功率分配优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310407173.0A CN103458521B (zh) 2013-09-09 2013-09-09 一种基于鲁棒性设计的mimo传输功率分配优化方法

Publications (2)

Publication Number Publication Date
CN103458521A true CN103458521A (zh) 2013-12-18
CN103458521B CN103458521B (zh) 2016-03-23

Family

ID=49740364

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310407173.0A Expired - Fee Related CN103458521B (zh) 2013-09-09 2013-09-09 一种基于鲁棒性设计的mimo传输功率分配优化方法

Country Status (1)

Country Link
CN (1) CN103458521B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103997775A (zh) * 2014-06-03 2014-08-20 东南大学 频分复用多用户mimo能效优化方法
CN104393956A (zh) * 2014-11-26 2015-03-04 北京邮电大学 一种用于无线携能通信系统最大化和速率预编码方法
CN106941375A (zh) * 2017-02-09 2017-07-11 东南大学 多色可见光通信中滤光片参数的鲁棒性联合优化方法
CN106982085A (zh) * 2016-01-15 2017-07-25 上海交通大学 多小区mimo系统中块对角化辅助的鲁棒收发机设计方法
CN107306244A (zh) * 2016-04-21 2017-10-31 中兴通讯股份有限公司 一种对角加载方法和装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102098244A (zh) * 2011-01-05 2011-06-15 东南大学 一种显式延时反馈下多用户预编码的鲁棒构造方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102098244A (zh) * 2011-01-05 2011-06-15 东南大学 一种显式延时反馈下多用户预编码的鲁棒构造方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
TADILO ENDESHAW BOGALE等: "Robust Transceiver Optimization for Downlink Coordinated Base Station Systems:Distributed Algorithm", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 *
许威等: "MIMO-OFDM系统中的比特和功率联合分配算法", 《电子与信息学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103997775A (zh) * 2014-06-03 2014-08-20 东南大学 频分复用多用户mimo能效优化方法
CN103997775B (zh) * 2014-06-03 2017-03-15 东南大学 频分复用多用户mimo能效优化方法
CN104393956A (zh) * 2014-11-26 2015-03-04 北京邮电大学 一种用于无线携能通信系统最大化和速率预编码方法
CN104393956B (zh) * 2014-11-26 2017-12-05 北京邮电大学 一种用于无线携能通信系统最大化和速率预编码方法
CN106982085A (zh) * 2016-01-15 2017-07-25 上海交通大学 多小区mimo系统中块对角化辅助的鲁棒收发机设计方法
CN106982085B (zh) * 2016-01-15 2020-08-04 上海交通大学 多小区mimo系统中块对角化辅助的鲁棒收发机设计方法
CN107306244A (zh) * 2016-04-21 2017-10-31 中兴通讯股份有限公司 一种对角加载方法和装置
CN107306244B (zh) * 2016-04-21 2021-08-31 中兴通讯股份有限公司 一种对角加载方法和装置
CN106941375A (zh) * 2017-02-09 2017-07-11 东南大学 多色可见光通信中滤光片参数的鲁棒性联合优化方法
CN106941375B (zh) * 2017-02-09 2019-04-30 东南大学 多色可见光通信中滤光片参数的鲁棒性联合优化方法

Also Published As

Publication number Publication date
CN103458521B (zh) 2016-03-23

Similar Documents

Publication Publication Date Title
US20200044899A1 (en) Method for automatically identifying modulation mode for digital communication signal
CN103458521B (zh) 一种基于鲁棒性设计的mimo传输功率分配优化方法
CN101399798B (zh) 一种ofdma无线通信系统的稳健信号传输方法及装置
CN102111197B (zh) 预编码矩阵索引信息上报方法和设备
CN102340784B (zh) 选择用户终端以增强上下行互逆误差校准的方法和装置
CN104702378A (zh) 混合高斯分布的参数估计方法和装置
Lin et al. On secrecy capacity of fast fading MIMOME wiretap channels with statistical CSIT
CN101345592B (zh) 应用于mimo的自适应信号的检测器及检测方法
US20100234053A1 (en) Systems and method for coordinated multipoint downlink transmissions
CN105577329A (zh) 基于空间调制的物理层安全传输方法
CN101951678B (zh) 多用户mimo系统中的功率分配方法和装置
CN102347820A (zh) 一种多小区协作无线通信系统联合编解码方法
CN103036601A (zh) 一种确定秩指示和预编码矩阵索引的方法和装置
CN105871434A (zh) Mimo干扰信道的拓扑干扰对齐方法
Coon et al. Combined bulk and per-tone transmit antenna selection in OFDM systems
CN114362855A (zh) 一种基于lstm的信道状态预测方法及系统
CN106027184A (zh) 一种基于lte-a的链路级到系统级仿真的接口方法
CN102487309B (zh) 一种mimo系统下的信号检测方法和装置
CN104283634A (zh) 一种数据的发送方法、接收方法、系统及装置
CN103607234B (zh) 一种中继合并接收方法和系统
CN104836605A (zh) 一种基于空间复用的新型收发天线联合选择方法
CN104539309B (zh) 基于极化失配的功放非线性影响下的全双工自干扰消除方法
CN105049107B (zh) 基于证据理论与物理层网络编码结合的满分集多天线双向中继接收方法
CN105634579A (zh) 基于解码转发网络编码的多天线双向中继证据理论接收方法
CN105099530B (zh) 认知无线电mimo-ofdm系统中基于认知用户漏功率的干扰抑制预编码方法

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160323

Termination date: 20190909