CN104506223B - 任意时间偏差下鲁棒的分布式波束形成器设计方法 - Google Patents
任意时间偏差下鲁棒的分布式波束形成器设计方法 Download PDFInfo
- Publication number
- CN104506223B CN104506223B CN201410842278.3A CN201410842278A CN104506223B CN 104506223 B CN104506223 B CN 104506223B CN 201410842278 A CN201410842278 A CN 201410842278A CN 104506223 B CN104506223 B CN 104506223B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mtd
- eta
- 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.)
- Expired - Fee Related
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0426—Power distribution
- H04B7/043—Power distribution using best eigenmode, e.g. beam forming or beam steering
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/022—Site diversity; Macro-diversity
- H04B7/024—Co-operative use of antennas of several sites, e.g. in co-ordinated multipoint or co-operative multiple-input multiple-output [MIMO] systems
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Radio Relay Systems (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种任意时间偏差下鲁棒的分布式波束形成器设计方法,将整数倍的时间偏差假设为频率选择性信道的影响,将从发射节点到中继节点以及从中继节点到接收节点的信道都假设为频率选择性信道;N个协作中继节点以AF方式将接收到的信号同步转发给接收节点,利用复系数{ωn}进行波束形成加权设计;假设脉冲成型函数为升余弦函数,用线段来拟合升余弦函数曲线,将得到的近似升余弦函数代入接收信号中,化简得到接收信干噪比,利用柯西施瓦茨不等式以及三角不等式求出信干噪比的下限值;采用凸优化的方法,最大化接收信干噪比的下限值,分别在总的中继发射功率和单个中继发射功率的约束下;求出对任意偏差下时间异步鲁棒的分布式波束形成方法的权值。
Description
【技术领域】
本发明涉及协作通信系统中的分布式波束形成技术领域,特别涉及一种任意偏差时间异步即同时存在整数和分数符号间隔时间偏差时鲁棒的分布式波束形成器设计方法。
【背景技术】
无线通信技术在近年内取得了突飞猛进的发展,早已成为人们生活和工作中必不可少的部分。而随着无线通信技术的发展,人们对于移动通信质量尤其是数据速率的要求越来越高。
MIMO系统通过在收发两端配置多天线,可以有效地利用空间信道的独立衰落特性,获得分集和复用增益,已经成为移动通信系统提高系统性能和频谱效率的重要方式。然而,实际中的无线终端往往由于尺寸、功耗、成本或者硬件复杂度的限制,很难配置多天线。此时,协作通信的概念被提出。
协作通信系统通过多个单天线中继节点的协作来组成虚拟的天线阵列,可以取得相当于MIMO系统的性能。在协作通信系统中,各个中继用户以一定的方式共享天线,就构成了一个虚拟的MIMO系统。目前常见的应用协作通信技术的网络系统有:协作中继网络(Relay)、协作分布式天线系统(Cooperative Distributed Antenna Systems,Cooperative DAS)、协作多点传输网络(Coordinated Multiple Point transmission andreception,CoMP)等。
协作波束形成技术允许各个中继节点对接收信号进行复加权,调整各个信号的幅度和相位使得各个中继节点通过相互协作形成期望的波束,从而通过用户间协作获得空间分集增益和空间阵列增益。协作通信技术给我们带来性能提升的同时,在实际应用中也面临着许多挑战,一个突出的难题就是时间异步问题。由于协作通信中的各个协作节点是分布式的,处于不同位置的协作节点的处理延时以及它们到接收节点的传播延时的不同,使得它们发射的中继信号到达接收节点处时存在不同的时偏,即时间异步问题。由于接收端的信号是多个不同时偏信号的叠加,将会导致系统性能的严重恶化。特别地,当各个协作中继节点相距很远时,它们到达接收端的相对时间偏差很大,可能会超过一个符号周期,此时就会存在整数和分数倍符号间隔的时间偏差,系统性能下降更加明显。
目前大多数分布式波束形成的研究都假设理想时间同步,而针对时间异步下的情况研究很少。对时间异步的分布式波束形成鲁棒的设计至今仍很不完善。文献“The effectof time synchronization errors on the performance of cooperative MISOsystems”分析了时间异步对协作MIMO系统的性能影响,但是没有提出解决方案。已公开的中国发明专利“CN103402213A”提出了一种协作通信系统中对小尺度偏差时间异步鲁棒的分布式波束形成方法,但是假设时间偏差不超过一个符号周期,而对于偏差较大即同时存在整数和分数倍符号间隔时间偏差的情况,该方法也将失效。本发明将对会对同时存在整数和分数倍符号间隔时间偏差时鲁棒的分布式波束形成技术进行研究。
【发明内容】
本发明的目的是提供一种任意时间偏差下鲁棒的分布式波束形成器设计方法,能对任意偏差时间异步即同时存在整数和分数倍符号间隔时间偏差的协作通信系统的性能进行显著提高。
为了达到以上目的,本发明是采用如下的技术方案实现的:
一种任意时间偏差下鲁棒的分布式波束形成器设计方法,其特征在于,该方法应用于采用一个发射节点,一个接收节点,N个中继节点的协作通信系统模型中;其中发射节点、N个中继节点和接收节点均配置单天线;中继节点采用协作地方式共享彼此的天线以形成一个虚拟的多天线系统,实现空域分集;采用放大转发的中继策略进行通信;
所述任意时间偏差下鲁棒的分布式波束形成器设计方法,包括下述步骤:
(1)获取发射端所发射的信号序列sl,经调制后所发射的基带信号获取发射成型脉冲g(t),符号周期Ts,发射信号功率Ps;
(2)信道模型为从发射节点到N个中继节点的1×N的SIMO信道和从N个中继节点到接收节点的N×1的MISO信道都是频率选择性信道,获取信道阶数M1和M2,获取对应的复衰落系数和将从发射端到中继的信道和从中继到接收端信道建模为:
其中,是发射节点与第n个协作中继节点之间信道的第m1个抽头的复衰落系数,为相应的延时;是第n个协作中继节点与接收节点之间信道的第m2个抽头的复衰落系数,为相应的延时;
(3)获取第n个中继节点处的高斯白噪声υn的方差值将第n个中继节点处的接收信号建模为:
(4)中继节点对它接收到的信号rn乘以后转发给接收端,其中是加权系数,bn是伸缩系数,计算伸缩系数
(5)获取接收机处的高斯白噪声nk的功率对接收信号进一步简化并将其写成矩阵的形式:
其中,d=(k-l)-(m1+m2),η=d+m1+m2;
(6)表示出接收信干噪比,利用柯西施瓦茨不等式和三角不等式给出接收信干噪比的最小值;
(7)采用凸优化的方法求出对时间异步鲁棒的分布式波束形成的加权值;
(8)在各个中继处应用所求出的加权值进行波束形成。
优选的,所述步骤(1)中,采用升余弦函数来替代脉冲成型函数;升余弦函数的主瓣两边各有L0个旁瓣,将升余弦函数用线段进行拟合,这些线段的选取与偏差的最大值有关,若偏差是在(-ξ,ξ)之间均匀分布的,那么这些线段就是以整数倍的符号周期lTs为起点,lTs分别加减ξ为终点拟合得到的;当采样的时候往左偏的话,升余弦函数用lTs点左边的线段进行拟合,往右偏的话就用lTs点右边的线段进行拟合;令ξn=δ-εn,将g(-m1-m2+ξn)和g(d+ξn)进行泰勒展开:
g′(mTs)用线段的斜率来近似:
其中,g′n(m)表示经第n个中继信道延时后的等效脉冲信号对应曲线的导数值g′(mTs),和分别为mTs点左右两边线段的斜率。
优选的,所述步骤(5)中,为了将接收机处的信号表示成矩阵形式,定义了以下矩阵:对于第n个中继节点,
——发射端到该中继的信道参数;
——该中继到接收端的信道参数;
——该中继处的噪声向量;
——该中继处等效成型脉冲函数各点的导数值;
——经该中继的收发端等效信道参数;
为了进一步化简接收信号,记定义以下三个辅助矩阵:
并记
其中并记
第k个采样时刻的接收信号为:
优选的,所述步骤(6)中,接收信干噪比为:
其中,Ps为发射功率,为平均信噪比;
显然:
假设hη为范数有界的定值,且满足||hη||≤μη,利用柯西施瓦茨不等式和三角不等式:
接收信干噪比的最小值:
其中,
优选的,所述步骤(6)中,采用凸优化的方法求出权值,其中凸优化的目标函数为最大化最小的接收信干噪比,约束总的中继发射功率:
重新定义矩阵并借助一些辅助变量,将该优化问题等效转化为下面的SOCP凸优化问题:
这是一个标准的凸优化问题,通过求解该优化问题即可得到对时间异步鲁棒的分布式波束形成加权系数。
优选的,所述协作通信系统为Relay系统;从源节点到中继节点的SIMO信道为非频率选择性信道,从中继节点到目的节点的MISO信道为频率选择性信道。
优选的,所述协作通信系统为DAS系统;DAS系统的分布式天线RAU与基带处理单元BPU之间通过高速光纤等有线方式连接,在分布式天线RAU端能够准确地获得并解码所传输的信号。
与现有的波束形成方法相比,本发明的有益效果是:
1、利用协作通信技术可以利用仅配置单天线的终端构成虚拟的MIMO系统,有效地利用空域信道的独立衰落特性,获得分集和复用增益,可以有效地扩大系统容量,提高传输质量。
2、可以克服由于整数和分数时间间隔的时间异步所造成的系统性能的恶化,使得协作通信技术更加实用。
3、本发明适用于多种协作通信系统,适用场景广泛。
【附图说明】
图1是本发明方法所涉及的通用信号传输模型。
图2是本发明对升余弦函数进行拟合的示意图。
图3是现有方法(最大比合并和非鲁棒波束形成)和本发明方法仿真结果对比图。其中:图3a为偏差的最大值为0.4,中继功率约束为总功率约束时接收信干噪比随发射总功率变化曲线图;图3b为偏差的最大值为0.4,中继功率约束为单个功率约束时接收信干噪比随发射总功率变化曲线图;图3c为偏差的最大值为0.2,中继功率约束为总功率约束时接收信干噪比随发射总功率变化曲线图;图3d为偏差的最大值为0.2,中继功率约束为单个功率约束时接收信干噪比随发射总功率变化曲线图。
图4是本发明方法所涉及的实施例1中协作Relay系统下的信号传输模型。
图5是本发明方法所涉及的实施例2中协作DAS系统下的信号传输模型。
【具体实施方式】
下面结合附图与具体实施例对本发明做进一步的详细说明。
本发明一种任意时间偏差下鲁棒的分布式波束形成器设计方法,应用于图1所示协作通信系统模型中,该模型有一个发射节点和一个接收节点,以及N=8个协作中继节点。假设从发射节点到协作中继节点的SIMO信道和从协作中继节点到接收节点的MISO信道都是频率选择性信道,信道阶数分别为M1和M2。我们假设整数倍的时间偏差是由于频率选择性信道的影响,因此频率选择性信道的阶数就相当于偏差的整数倍数。从发射端到中继的信道和从中继到接收端信道可以表示为:
其中,是发射节点与第n个协作中继节点之间信道的第m1个抽头的复衰落系数,是相应的延时;是第n个协作中继节点与接收节点之间信道的第m2个抽头的复衰落系数,是相应的延时。
假设发射端的所发射的信号序列为sl,经过调制后所发射的基带信号为sl是发送序列,g(t)是发射成型脉冲,采用升余弦函数,Ts是符号周期。
第n个协作中继节点处的接收信号为:
其中,υn是第n个中继节点处的高斯白噪声,且假设
协作中继节点以AF(amplify and forward)方式同步工作,并利用复系数对接收到的信号进行转发。为波束形成加权系数,为伸缩系数,将中继处的发射功率归一化为||ωn||2,即
接收端所接收到的信号是N个中继节点所发射信号在接收端的叠加:
其中,n(t)是接收端处的噪声,且假设其为高斯白噪声,满足且其与中继节点处噪声具有相同的方差值,记为
对延时τn和采样偏差Δ进行归一化,其中,m1,m2表示对应偏差的整数部分;表示小数部分,且令δ表示符号归一化后的采样偏差。本发明把由于整数倍的偏差导致的时间异步假设为频率选择性信道的影响,即当偏差为(m1+m2)Ts+εnTs时,认为频率选择性信道有(m1+m2)阶,此时判决时刻k的离散接收信号可以表示为:
其中,nk是接收端处噪声的采样值,满足
记d=(k-l)-(m1+m2),ξn=δ-εn,若假设升余弦函数的主瓣两旁各有L0个旁瓣,则此时有:
应用此公式简化后的接收信号为:
将g(-(m1+m2)+ξn)和g(d+ξn)进行泰勒展开,当ξn很小时:
所采用的升余弦函数如图2所示,当旁瓣数大于3时,旁瓣的幅度变得很小,因此令L0=3。另外整个升余弦函数曲线可以用12条线段来拟合,这12条线段表示为sg6,sg5,…,sg1,sg1',…,sg5',sg6'。其中sg1',…,sg6'分别与sg1,…,sg6相对称。因此g′n(m)可以用这些线段的斜率来近似:
其中,g′n(m)表示经第n个中继信道延时后的等效脉冲信号对应曲线的导数值g′(mTs),和分别为对应曲线中mTs点左右两边线段的斜率;
下面将接收信号表示为矩阵形式,对于第n个中继节点,
——发射端到该中继的信道参数;
——该中继到接收端的信道参数;
——该中继处的噪声向量;
——该中继处等效成型脉冲函数各点的导数值;
——经该中继的收发端等效信道参数;
为了进一步化简接收信号,记定义以下三个辅助矩阵:
并记
其中并记
则第k个采样时刻接收端的接收信号可以表示为如下三个部分:
信号部分:
干扰部分:
噪声部分:
接收信干噪比为:
其中,Ps为发射功率,ρs=Ps/σ2是平均符号信噪比。
由此可以看到接收信干噪比的值与hη的取值密切相关。hη可以被看做一个欧几里得范数有界的定值。即||hη||≤μη,η∈[-L0,M1+M2+L0]。
采用凸优化的方法来求得权值,凸优化问题可以表述为:
其中P是中继处的总发射功率。
显然:
利用柯西施瓦茨不等式和三角不等式进行化简:
代换之后得:
其中,
此优化问题可以表示为:
首先我们令将上面的优化问题等效转化为下面的形式:
引入正实数辅助变量t和{kη},优化问题可以转化为:
假设并假设|xHp0|>μ0t,则上述优化问题等价于下面的形式:
接下来考虑:
不难看出,此优化问题在第一个约束条件取等号时才是最优的,因此上面两个优化问题是等价的。
通过对x进行相位旋转使xHp0的虚部为零,这并不影响优化结果,可以将上面的优化问题等效变换为:
为了进一步化简,可以定义矩阵:
就变成了最终的优化问题形式:
这是一个标准的SOCP优化问题,用此凸优化问题的最优解的前N个元素得到x,然后令ω=x/y作为原问题的解ω,即得到对时间异步鲁棒的分布式波束形成的加权值;利用获得的复加权值对各中继处的信号进行加权发射,即可在整数和分数倍符号偏差的时间异步环境下实现鲁棒的波束形成器,稳定地提高系统性能,使得协作通信技术更加实用。
为了验证在相同发射总功率下,本发明方法可以获取比参考方法(最大比合并方法和非鲁棒的波束形成方法)更加好的性能,此处的性能用接收信干噪比来表征。图3给出了偏差的最大值分别为0.4和0.2以及单个中继发射功率约束和各中继总发射功率约束的情况下,相同条件下本发明方法与参考方法接收信干噪比随发射功率变化曲线的对比图。从图中可以看出,使用本发明方法能得到比参考方法更好的性能。
至此,从技术方案和仿真结果都可以验证本发明方法在克服时间异步影响方面的积极效果。
下面将介绍两个具体实施例:
实施例1:
实施例1考虑一种如图4所示的协作Relay系统,有一个发射节点和一个接收节点,以及N=8个协作中继节点。当协作中继节点相对于发射节点分布较为均匀时,从发射节点到各个协作中继节点的相对时间偏差较小,假设它们之间的SIMO信道为非频率选择性信道;同时仍假设从协作中继节点到接收节点的相对时偏较大,假设它们之间的MISO信道为频率选择性信道,信道阶数为M。
我们仍假设整数倍的时间偏差是由于频率选择性信道的影响,因此频率选择性信道的阶数就相当于偏差的整数倍数。信道可以表示为:
其中,αn是发射节点与第n个协作中继节点之间信道的复衰落系数,为信道延时;βn,m是第n个协作中继节点与接收节点之间信道的第m个抽头的复衰落系数,为相应的延时。
假设发射端的所发射的信号序列为sl,经调制后所发射的基带信号为sl是发送序列,g(t)是发射成型脉冲,Ts是符号周期。
第n个协作中继节点处的接收信号为:
其中,υn是第n个中继节点处的高斯白噪声,且假设
协作中继节点以AF(amplify and forward)方式同步工作,并利用复系数对接收到的信号进行转发。将中继处的发射功率归一化为|ωn|2,为波束形成加权系数。接收端所接收到的信号是N个中继节点所发射信号在接收端的叠加:
其中,n(t)是接收机处的高斯白噪声,且假设并假设
对延时τn和采样偏差Δ进行归一化,其中,ηn和m分别表示和符号归一化后的整数部分;和表示对应的小数部分,令δ表示符号归一化后的采样偏差。对于从发射基站到中继的非频率选择性信道,当协作中继节点之间的间隔不太远时,相对时间延迟的整数部分是相同的,只有小数部分是不同的。不失一般性,我们假设本发明把由于整数倍的偏差导致的时间异步假设为频率选择性信道的影响,即当偏差为mTs+εnTs时,认为从中继到接收节点的频率选择性信道有m阶,采样时刻k时刻的离散接收信号可以表示为:
假设g(t)为升余弦函数,主瓣两旁各有L0个旁瓣,令d=(k-l)-m,ξn=δ-εn,则:
应用此公式简化后的接收信号可以表示为:
将g(-m+ξn)和g(d+ξn)进行泰勒展开,当ξn很小时:
仍然假设所用升余弦函数的旁瓣数目L0=3,按照上文的方法用12条线段来拟合,因此g'(mTs)可以用这些线段的斜率来近似:
其中,g′n(m)表示经第n个中继信道延时后的等效脉冲信号对应曲线的导数值g′(mTs),和分别为对应曲线中mTs点左右两边线段的斜率;
下面将接收信号表示为矩阵形式,对于第n个中继节点,
——该中继到接收端的信道参数;
——该中继处的噪声向量;
——该中继处等效成型脉冲函数各点的导数值;
——经该中继的收发端等效信道参数;
为了进一步化简接收信号,记定义以下三个辅助矩阵:
并记
其中并记
则第k个采样时刻接收端的接收信号可以表示为如下三个部分:
信号部分:
干扰部分:
噪声部分:
接收信干噪比为:
其中,Ps是发射功率,ρs=Ps/σ2是平均符号信噪比。
由此可以看到接收信干噪比的值与hη的取值密切相关。hη可以被看做一个欧几里得范数有界的定值。即||hη||≤μη,η∈[-L0,M+L0]。
采用凸优化的方法来求得权值,凸优化问题可以表述为:
其中P是中继处的总发射功率。
利用柯西施瓦茨不等式和三角不等式可以求出:
其中,
凸优化问题可以表示为:
此优化问题与上一节中的优化问题形式完全相同,利用上一节中的变换方法即可将其等效变换为标准的SOCP优化问题,通过解此凸优化问题求出对时间异步鲁棒的分布式波束形成的权值。
实施例2:
实施例2考虑一种协作DAS系统。在协作DAS系统中,通常分为远端天线单元(Remote Antenna Units,RAUs)和基带处理单元(Baseband Processing Unit,BPU)。RAUs只进行信号发射,信号的基带处理在BPU中完成,通常远端天线单元RAUs是通过射频光纤与基带处理单元BPU连接的,因此RAUs处可以准确地获取所发射的信号,类似于采用DF中继策略的协作中继系统。此时只需要考虑从天线单元组到接收节点之间的时间异步问题。
如图5所示,假设DAS系统有N=8个天线单元RAUs,有一个接收节点。假设从天线单元组到接收节点的相对时间偏差较大,包含整数和分数倍符号间隔,因此假设为它们之间的MISO信道为频率选择性信道,阶数假设为M。假设整数倍的时间偏差是频率选择性信道的影响,偏差的整数倍数就是频率选择性信道的阶数。
从天线单元组到接收节点的MISO信道可以表示为:
其中,M为此频率选择性信道的阶数,αn,m是第n个天线单元与接收节点之间信道的第m个抽头的复衰落系数,τn,m为相应的延时。
假设发射端的所发射的信号序列为sl,经调制后所发射的基带信号为sl是发送序列,g(t)是发射成型脉冲,Ts是符号周期。
多个天线站协同工作,利用复加权系数{ωn}进行分布式波束形成传输,接收端所接收到的信号是N个天线站所发射信号在接收端的叠加:
其中,n(t)是接收机处的高斯白噪声,且假设
在接收端对接收到的信号在t=kTs+Δ时刻进行采样,其中Δ为采样偏差。并对相对时延和采样偏差Δ进行归一化,τn=mTs+εnTs,Δ=δTs。其中,m是偏差的整数倍数,也就是频率选择性信道的阶数,εn和δ为小数部分。
此时可以得到k时刻离散的接收信号表达式:
假设g(t)为升余弦函数,且主瓣两旁各有L0个旁瓣,令d=(k-l)-m,ξn=δ-εn,则:
简化后的公式为:
将g(-m+ξn)和g(d+ξn)进行泰勒展开,当ξn很小时:
假设旁瓣数目为L0=3。整个升余弦函数曲线用12条线段来拟合,g'(mTs)可以用这些线段的斜率来近似:
其中,g′n(m)表示经第n个中继信道延时后的等效脉冲信号对应曲线的导数值g′(mTs),和分别为对应曲线中mTs点左右两边线段的斜率;
记η=d+m,定义以下矩阵:
——各天线到接收端的信道参数;
——各天线处等效成型脉冲函数导数值;
并记
其中并记
则此时的接收信号可以表示为三个部分:
信号部分:
干扰部分:
噪声部分:yk3=nk,
接收信干噪比为:
其中,Ps为发射功率,ρs=Ps/σn 2是平均符号信噪比。
由此可以看到接收信干噪比的值与hη的取值密切相关。hη可以被看做一个欧几里得范数有界的定值。即||hη||≤μη,η∈[-L0,M1+M2+L0]。
采用凸优化的方法来求得权值,凸优化问题可以表述为:
其中P是中继处的总发射功率。
利用柯西施瓦茨不等式和三角不等式可以得到:
其中,
上面的优化问题可以转化为:
仿照上一节的变换方式,即可将此优化问题等效变换为标准的SOCP凸优化问题,通过求解SOCP凸优化问题求出对时间异步鲁棒的分布式波束形成的权值。
上述两种具体实施例最后得到的优化问题只是在本发明方法所涉及的通用模型的一种特例。
Claims (5)
1.一种任意时间偏差下鲁棒的分布式波束形成器设计方法,其特征在于,该方法应用于采用一个发射节点,一个接收节点,N个中继节点的协作通信系统模型中;其中发射节点、N个中继节点和接收节点均配置单天线;中继节点采用协作地方式共享彼此的天线以形成一个虚拟的多天线系统,实现空域分集;采用放大转发的中继策略进行通信;
所述任意时间偏差下鲁棒的分布式波束形成器设计方法,包括下述步骤:
(1)获取发射端所发射的信号序列sl,经调制后所发射的基带信号l为序列号码;获取发射成型脉冲g(t),符号周期Ts,发射信号功率Ps;t为采样时刻;
(2)信道模型为从发射节点到N个中继节点的1×N的SIMO信道和从N个中继节点到接收节点的N×1的MISO信道都是频率选择性信道,获取信道阶数M1和M2,获取对应的复衰落系数和将从发射端到中继的信道和从中继到接收端信道建模为:
<mrow>
<msubsup>
<mi>c</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
</munderover>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</mrow>
</msub>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msubsup>
<mi>&tau;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msubsup>
<mi>c</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</munderover>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msubsup>
<mi>&tau;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中,是发射节点与第n个协作中继节点之间信道的第m1个抽头的复衰落系数,为相应的延时;是第n个协作中继节点与接收节点之间信道的第m2个抽头的复衰落系数,为相应的延时;
(3)获取第n个中继节点处的高斯白噪声υn的方差值将第n个中继节点处的接收信号rn(t)建模为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CircleTimes;</mo>
<msubsup>
<mi>c</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&upsi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</munder>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</mrow>
</msub>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msubsup>
<mi>&tau;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&upsi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
其中,υn(t)表示第n个中继节点处的噪声;
(4)中继节点对它接收到的信号rn乘以后转发给接收端,其中是加权系数,bn是伸缩系数,计算伸缩系数
(5)获取接收机处的高斯白噪声nk的功率对接收信号进一步简化并将其写成矩阵的形式,则第k个采样时刻的接收信号采样值yk为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<mi>n</mi>
</munder>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
<msubsup>
<mi>&omega;</mi>
<mi>n</mi>
<mo>*</mo>
</msubsup>
<munder>
<mi>&Sigma;</mi>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<munder>
<mi>&Sigma;</mi>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</munder>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</mrow>
</msub>
<msub>
<mi>g</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<munder>
<mi>&Sigma;</mi>
<mi>n</mi>
</munder>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
<msubsup>
<mi>&omega;</mi>
<mi>n</mi>
<mo>*</mo>
</msubsup>
<munder>
<mi>&Sigma;</mi>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<munder>
<mi>&Sigma;</mi>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</munder>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>d</mi>
<mo>&NotEqual;</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<msub>
<mi>g</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>d</mi>
<mo>+</mo>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>d</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<munder>
<mi>&Sigma;</mi>
<mi>n</mi>
</munder>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
<msubsup>
<mi>&omega;</mi>
<mi>n</mi>
<mo>*</mo>
</msubsup>
<munder>
<mi>&Sigma;</mi>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</munder>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<msub>
<mi>&upsi;</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</munderover>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>&eta;</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
<msub>
<mi>s</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>&eta;</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
</munderover>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
<msub>
<mi>s</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>&eta;</mi>
</mrow>
</msub>
<mo>+</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mi>u</mi>
<mo>+</mo>
<msub>
<mi>n</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
1
其中,d=(k-l)-(m1+m2),η=d+m1+m2,gn表示第n个中继处的脉冲成型函数,ξn表示第n个中继处传播时延的小数部分,k表示采样时刻,sk表示k采样时刻的信号采样值,ωH表示中继加权向量的Hermitian变换,p0,h0,pη,hη,u是为了变换定义的向量,定义如下:
η=0,…,M1+M2,并记
η=-L0,…,L0+M1+M2,其中并记
(6)表示出接收信干噪比,利用柯西施瓦茨不等式和三角不等式给出接收信干噪比的最小值;
(7)采用凸优化的方法求出对时间异步鲁棒的分布式波束形成的加权值;
(8)在各个中继处应用所求出的加权值进行波束形成;
所述步骤(6)中,接收信干噪比γ为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&gamma;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>|</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mi>u</mi>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>|</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mi>Z</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,Ps为发射功率,为平均信噪比;ω=x/y;ωH为ω转置共轭;
显然:
<mrow>
<msub>
<mi>&gamma;</mi>
<mi>k</mi>
</msub>
<mo>&GreaterEqual;</mo>
<mfrac>
<mrow>
<munder>
<mi>min</mi>
<msub>
<mi>h</mi>
<mn>0</mn>
</msub>
</munder>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munder>
<mi>max</mi>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</munder>
<mrow>
<mo>(</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mi>Z</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
假设hη为范数有界的定值,且满足||hη||≤μη,利用柯西施瓦茨不等式和三角不等式:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>|</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&GreaterEqual;</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>|</mo>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&GreaterEqual;</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>|</mo>
<mo>-</mo>
<munder>
<mi>max</mi>
<msub>
<mi>h</mi>
<mn>0</mn>
</msub>
</munder>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>|</mo>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>&omega;</mi>
<mo>|</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&le;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>p</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&le;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>p</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>+</mo>
<munder>
<mi>max</mi>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</munder>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>p</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mi>&eta;</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>&omega;</mi>
<mo>|</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
接收信干噪比的最小值:
<mrow>
<msup>
<mi>&gamma;</mi>
<mi>l</mi>
</msup>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<msub>
<mi>h</mi>
<mi>&eta;</mi>
</msub>
</munder>
<msub>
<mi>&gamma;</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>|</mo>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>&omega;</mi>
<mo>|</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>p</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mi>&eta;</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>&omega;</mi>
<mo>|</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>&mu;</mi>
<mn>2</mn>
</msup>
<mo>|</mo>
<mo>|</mo>
<mi>&omega;</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mi>Z</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
2
其中,
所述步骤(7)中,采用凸优化的方法求出权值,其中凸优化的目标函数为最大化最小的接收信干噪比,约束总的中继发射功率:
<mrow>
<munder>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mi>W</mi>
</munder>
<msup>
<mi>&gamma;</mi>
<mi>l</mi>
</msup>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
<mo>|</mo>
<mo>|</mo>
<mi>&omega;</mi>
<mo>|</mo>
<mo>|</mo>
<mo>&le;</mo>
<msqrt>
<mi>P</mi>
</msqrt>
<mo>,</mo>
</mrow>
即
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<munder>
<mi>max</mi>
<mi>&omega;</mi>
</munder>
</mtd>
<mtd>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>|</mo>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>&omega;</mi>
<mo>|</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>p</mi>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mi>&eta;</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>&omega;</mi>
<mo>|</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>&mu;</mi>
<mn>2</mn>
</msup>
<mo>|</mo>
<mo>|</mo>
<mi>&omega;</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<mi>Z</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
s.t.||ω||2≤P,
P是中继处的总发射功率;重新定义矩阵并借助一些辅助变量,将该优化问题等效转化为下面的SOCP凸优化问题:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<munder>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
</munder>
</mtd>
<mtd>
<mrow>
<msup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>p</mi>
<mo>~</mo>
</mover>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<mi>d</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>H</mi>
</msup>
<mi>D</mi>
<mo>|</mo>
<mo>|</mo>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>A</mi>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mo>|</mo>
<mo>|</mo>
<mo>&le;</mo>
<msup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>H</mi>
</msup>
<mi>e</mi>
<msqrt>
<mi>P</mi>
</msqrt>
<mo>,</mo>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>A</mi>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mo>|</mo>
<mo>|</mo>
<mo>&le;</mo>
<msup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>H</mi>
</msup>
<mi>d</mi>
<mo>,</mo>
</mrow>
<mrow>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>H</mi>
</msup>
<msub>
<mover>
<mi>p</mi>
<mo>~</mo>
</mover>
<mi>&eta;</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mo>&le;</mo>
<msup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>H</mi>
</msup>
<msub>
<mi>c</mi>
<mi>&eta;</mi>
</msub>
<mo>,</mo>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>H</mi>
</msup>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>H</mi>
</msup>
<msub>
<mover>
<mi>p</mi>
<mo>~</mo>
</mover>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>H</mi>
</msup>
<msub>
<mi>c</mi>
<mi>&eta;</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>H</mi>
</msup>
<mi>e</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
d=(k-l)-(m1+m2);
其中,
<mrow>
<msub>
<mi>c</mi>
<mi>&eta;</mi>
</msub>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mn>0</mn>
<mrow>
<mn>1</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>+</mo>
<mi>&eta;</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msub>
<mn>0</mn>
<mrow>
<mn>1</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>2</mn>
<mo>-</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>H</mi>
</msup>
<mo>,</mo>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<mi>d</mi>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mn>0</mn>
<mrow>
<mn>1</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>H</mi>
</msup>
<mo>,</mo>
</mrow>
<mrow>
<mi>e</mi>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mn>0</mn>
<mrow>
<mn>1</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>H</mi>
</msup>
<mo>,</mo>
</mrow>
<mrow>
<mi>D</mi>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<mo>&lsqb;</mo>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<mi>d</mi>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mi>d</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&mu;</mi>
<mi>d</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mover>
<mi>Z</mi>
<mo>~</mo>
</mover>
<mo>/</mo>
<msqrt>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
</msqrt>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>e</mi>
<mo>/</mo>
<msqrt>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>p</mi>
<mo>~</mo>
</mover>
<mi>&eta;</mi>
</msub>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>p</mi>
<mi>&eta;</mi>
<mi>H</mi>
</msubsup>
</mtd>
<mtd>
<msub>
<mn>0</mn>
<mrow>
<mn>1</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>H</mi>
</msup>
<mo>,</mo>
<mi>&eta;</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>x</mi>
<mi>H</mi>
</msup>
</mtd>
<mtd>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>k</mi>
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mi>t</mi>
</mtd>
<mtd>
<mi>y</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>H</mi>
</msup>
<mo>,</mo>
</mrow>
3
<mrow>
<msub>
<mover>
<mi>p</mi>
<mo>~</mo>
</mover>
<mi>&eta;</mi>
</msub>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>p</mi>
<mi>&eta;</mi>
<mi>H</mi>
</msubsup>
</mtd>
<mtd>
<msub>
<mn>0</mn>
<mrow>
<mn>1</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>H</mi>
</msup>
<mo>,</mo>
<mi>&eta;</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<munder>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<msub>
<mi>h</mi>
<mn>0</mn>
</msub>
</munder>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>&omega;</mi>
<mi>H</mi>
</msup>
<msub>
<mi>h</mi>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>&omega;</mi>
<mo>|</mo>
<mo>|</mo>
<mo>,</mo>
</mrow>
这是一个标准的凸优化问题,通过求解该优化问题得到对时间异步鲁棒的分布式波束形成加权系数。
2.如权利要求1所述的一种任意时间偏差下鲁棒的分布式波束形成器设计方法,其特征在于,所述步骤(1)中,采用升余弦函数来替代脉冲成型函数;升余弦函数的主瓣两边各有L0个旁瓣,将升余弦函数用线段进行拟合,这些线段的选取与偏差的最大值有关,若偏差是在(-ξ,ξ)之间均匀分布的,那么这些线段就是以整数倍的符号周期lTs为起点,lTs分别加减ξ为终点拟合得到的;当采样的时候往左偏的话,升余弦函数用lTs点左边的线段进行拟合,往右偏的话就用lTs点右边的线段进行拟合;对第n个中继处的传播延时τn和采样偏差Δ进行归一化,Δ=δTs;其中,m1,m2表示对应偏差的整数部分;表示小数部分,且令表示经过第n个中继的两级传播的总延时的归一化的小数部分;δ表示符号归一化后的采样偏差;令ξn=δ-εn,将g(-m1-m2+ξn)和g(d+ξn)进行泰勒展开:
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>&le;</mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>&NotEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>></mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>+</mo>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>d</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mi>d</mi>
<mo>&NotEqual;</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>d</mi>
<mo>&NotEqual;</mo>
<mn>0</mn>
<mo>,</mo>
<mi>d</mi>
<mo>&NotEqual;</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
g′(mTs)用线段的斜率来近似:
<mrow>
<msubsup>
<mi>g</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>a</mi>
<mi>m</mi>
<mo>+</mo>
</msubsup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>a</mi>
<mi>m</mi>
<mo>-</mo>
</msubsup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中,g′n(m)表示经第n个中继信道延时后的等效脉冲信号对应曲线的导数值g′(mTs),和分别为mTs点左右两边线段的斜率。
3.如权利要求1所述的一种任意时间偏差下鲁棒的分布式波束形成器设计方法,其特征在于,所述步骤(5)中,为了将接收机处的信号表示成矩阵形式,定义了以下矩阵:
对于第n个中继节点,
——发射端到该中继的信道参数;
——该中继到接收端的信道参数;
——该中继处的噪声向量;
——该中继处等效成型脉冲函数各点的导数值;
——经该中继的收发端等效信道参数;
为了进一步化简接收信号,记定义以下三个辅助矩阵:
η=0,…,M1+M2,并记
η=-L0,…,L0+M1+M2,其中并记
<mrow>
<mi>u</mi>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<msup>
<mrow>
<mo>[</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&beta;</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>v</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&beta;</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>v</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>.</mo>
<mo>.</mo>
<mo>.</mo>
<mo>,</mo>
<msub>
<mi>b</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&beta;</mi>
<mi>N</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>v</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>]</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>;</mo>
</mrow>
第k个采样时刻的接收信号为:
4.如权利要求1所述的一种任意时间偏差下鲁棒的分布式波束形成器设计方法,其特征在于,所述协作通信系统模型为Relay系统;从源节点到中继节点的SIMO信道为非频率选择性信道,从中继节点到目的节点的MISO信道为频率选择性信道。
5.如权利要求1所述的一种任意时间偏差下鲁棒的分布式波束形成器设计方法,其特征在于,所述协作通信系统模型为DAS系统;DAS系统的分布式天线RAU与基带处理单元BPU之间通过高速光纤有线方式连接,在分布式天线RAU端能够准确地获得并解码所传输的信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410842278.3A CN104506223B (zh) | 2014-12-30 | 2014-12-30 | 任意时间偏差下鲁棒的分布式波束形成器设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410842278.3A CN104506223B (zh) | 2014-12-30 | 2014-12-30 | 任意时间偏差下鲁棒的分布式波束形成器设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104506223A CN104506223A (zh) | 2015-04-08 |
CN104506223B true CN104506223B (zh) | 2017-10-20 |
Family
ID=52947946
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410842278.3A Expired - Fee Related CN104506223B (zh) | 2014-12-30 | 2014-12-30 | 任意时间偏差下鲁棒的分布式波束形成器设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104506223B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103402213B (zh) * | 2013-07-24 | 2016-04-27 | 西安交通大学 | 对时间异步鲁棒的分布式波束形成方法 |
CN104993859B (zh) * | 2015-08-05 | 2018-07-06 | 中国电子科技集团公司第五十四研究所 | 一种适用于时间异步环境下的分布式波束赋形方法 |
EP3403339B1 (en) * | 2016-01-13 | 2020-11-04 | Telefonaktiebolaget LM Ericsson (publ) | Practical hybrid precoding scheme for multi-user massive mimo systems |
CN106647253B (zh) * | 2016-09-28 | 2019-10-11 | 南京航空航天大学 | 航空发动机分布式控制系统多性能鲁棒跟踪控制方法 |
CN107888254B (zh) * | 2016-09-30 | 2021-02-12 | 电信科学技术研究院 | 一种信息传输方法及通信节点 |
CN108983201B (zh) * | 2018-07-23 | 2022-06-17 | 南京航空航天大学 | 一种基于概率约束的稳健宽带波束形成器设计方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102106162B (zh) * | 2008-06-23 | 2013-06-05 | 上海贝尔股份有限公司 | 用于上行链路的低复杂度多基站mimo方法和设备 |
CN103402213B (zh) * | 2013-07-24 | 2016-04-27 | 西安交通大学 | 对时间异步鲁棒的分布式波束形成方法 |
CN104052580B (zh) * | 2014-06-25 | 2017-08-01 | 西安交通大学 | 无线传感器网络中的多节点协同信号发射和接收方法 |
-
2014
- 2014-12-30 CN CN201410842278.3A patent/CN104506223B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104506223A (zh) | 2015-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104506223B (zh) | 任意时间偏差下鲁棒的分布式波束形成器设计方法 | |
CN104601209B (zh) | 一种适用于3d‑mimo系统的协作多点传输方法 | |
CN109302250A (zh) | 能量传输全双工中继进行中继选择与发送功率分配方法 | |
CN107113052A (zh) | 用于中继器小小区的频分双工(fdd)大规模mimo回程 | |
CN102710309B (zh) | 适用于大规模天线阵列的同步信号发射方法 | |
CN106301457B (zh) | 双向全双工mimo中继通信系统中的自干扰消除方法 | |
WO2009107314A1 (ja) | マルチホップ無線通信システム | |
CN102694628B (zh) | 基于多用户mimo协同中继系统的干扰抑制方法 | |
CN105897319A (zh) | 一种mimo全双工中继系统信源中继联合预编码方法 | |
KR20160146692A (ko) | Mimo 송신기의 전처리 필터를 이용한 전송 신호 생성 방법 | |
CN105429688A (zh) | 大规模分布天线系统中抑制导频污染的多小区预编码方法 | |
CN107852614A (zh) | 具有基于fir的信道均衡器的无线中继器 | |
CN106231665A (zh) | 数能一体化网络中基于rrh动态模式切换的资源分配方法 | |
CN103516410A (zh) | 一种利用有源天线系统实现下行数据发射方法及装置 | |
CN103269241B (zh) | 一种多跳多中继协作网络的波束成形方法 | |
CN103297108B (zh) | 一种多源多中继协作网络的上行波束成形方法 | |
CN106452536A (zh) | 用于多输入多输出通信的长期反馈的方法和装置 | |
CN108880633A (zh) | 一种波束成形天线选择分组算法联合设计优化方法 | |
CN104821840B (zh) | 一种大规模多输入多输出下行系统的抗干扰方法 | |
CN114745032B (zh) | 一种无蜂窝大规模mimo智能分布式波束选择方法 | |
CN104113367B (zh) | 混合协作多点传输系统中的反馈方法及相关的用户设备 | |
CN103402213A (zh) | 对时间异步鲁棒的分布式波束形成方法 | |
CN101783699B (zh) | 多用户mimo系统中下行链路的信号发射方法和装置 | |
CN107104715B (zh) | 一种基于天线选择的干扰对齐方法 | |
CN105227224A (zh) | 3d-mimo系统中基于最小slnr最大准则的分布式干扰协调方法 |
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: 20171020 Termination date: 20201230 |
|
CF01 | Termination of patent right due to non-payment of annual fee |