技术领域
[0001] 本
发明属于无线通信技术领域，尤其涉及一种基于空间复用的新型收发天 线联合选择方法。
背景技术
[0002] 随着无线通信的快速发展，微型电脑、智能手机的广泛应用，人们更希望 任何时间任何地点都能够接入互联网来获取信息。这要求通信网络具有更高的 传输速率和更好的服务
质量。传统的单输入单输出系统(SISO)远远不能满足实 际需求。大规模的多输入多输出系统(又称多发射天线多接收天线系统，MIMO) 是无线通信系统重大的技术突破，它在发射端和接收端配置大量的天线。与传 统的SISO系统相比，在不增加系统带宽和发射功率的情况下，MIMO系统的 容量和可靠性有很大的改善。由于MIMO通信系统有着较大的性能优势，已经 引起了该领域专家、学者和运营商的广泛兴趣。更重要的是第三代移动通信系 统(3GPP)已经在标准中已部分引入了MIMO技术。目前，MIMO技术也已广泛 应用于超3G和4G标准中。可以预见，在未来移动通信系统势必大规模地采用 MIMO技术。
[0003] MIMO系统可以通过空间复用或者空间分集来提高系统的容量和降低误比 特率，改善通信系统的性能。采用MIMO技术的通信系统设备需要与天线数相 同的射频链路(主要包括
放大器、上下频转换器、模/数/模转换器)，导致系统的
硬件成本和复杂度较大。因此，如何保持MIMO系统性能优越性同时降低系统 的复杂度是一个重要的研究课题。天线选择是一种可以保持系统性能优势，同 时降低系统复杂度和硬件成本的有效方法。天线选择的思想是用有限的射频链 路来匹配较多的天线，它可以有效地降低通信设备的射频链路数目，同时保持
[0004] 基于综上所述原因，天线选择技术已经引起了广泛的研究兴趣。天线选择 可以有效的降低MIMO系统的复杂度和硬件成本，同时有效地改善MIMO通 信系统的性能。为了选择适当的天线子集，近年来MIMO系统中的天线选择算 法在大量的文献中被深入研究。很大一部分的天线选择技术都是基于发射端或 接收端，收发联合天线选择
算法很少被研究。给定发射和接收端天线数目，为 了找到最优的发射和接收天线子集，最优的天线子集选择算法(本发明简称， 最优算法)计算所有可能的天线子集组合，从中选择出能够使系统性能最优的 天线子集。虽然最优算法能够取得最好的系统性能，但是，由于最优算法计算 复杂度过大，算法的复杂度随着天线数目的增加成指数增长，特别是在快变信 道中，天线选择需要快速估计出信道，而最优算法的复杂度大，很难达到实时 性要求。
发明内容
[0005] 本发明的目的在于提供一种基于空间复用的新型收发天线联合选择方法， 旨在解决现有的收发天线联合选择方法存在的操作复杂，工作量较大，实用性 较低的问题。
[0006] 本发明是这样实现的，一种基于空间复用的新型收发天线联合选择方法， 该基于空间复用的新型收发天线联合选择方法在起始时对应用于天线选择的发 射接收天线信道集合为满集，在每一次选择过程中选择出对信道容量贡献最小 的一对收发天线。
[0007] 进一步，该基于空间复用的新型收发天线联合选择方法包括以下步骤：
[0008] 第一步：天线选择初始化，设置多天线通信系统的参数ρ、NR、NT、Lt、 Lr大小,按图2所示描述系统生成起始满集信道矩阵H，清零选择步长n＝0；
[0009] 第二步：计算
[0010] 第三步：从发射天线集 中选择某一个列向量 使得 最小,将其从发射天线集中删除
[0011] for i＝1:NT-n
[0012]
[0013] end
[0014] Jn＝arg maxJn,i
[0015]
[0016] 第四步：对发射天线选择后，信道矩阵更新为
[0017] 第五步：计算
[0018] 第六步：并从接收天线集 中选择某一个行向量 使得 最小，并将其从接收天线集中删除
[0019] for j＝1:NR-n
[0020]
[0021] end
[0022] λn＝argmaxλn,j
[0023]
[0024] 第七步：对接收天线选择后，信道矩阵更新为
[0025] 第八步：如果n<＝min{NR-Lr,NT-Lt}则返回至第2步；如果n>min{NR-Lr, NT-Lt},那么发射端和接收端中具有最小值的一端将停止天线选择，而另外一端 将继续执行天线选择，直到n＝max{NR-Lr,NT-Lt}时天线选择结束。
[0026] 进一步，在发射端和接收端分别配置NT根和NR根天线；发射端和 接收端分别配置Lt个和Lr个射频链路；假设所有的发射天线和接收天线都参与 天线选择；接收端从NR天线中选择Lr接收天线，发射端从NT天线中选择Lt发 射天线；在发射端天线和接收端天线之间的所有无线信道是独立的、且服从平 坦的
瑞利衰落，用随机信道矩阵H表示；这样，经过匹配
滤波器后，接收端接 收到的复包络
信号为：
[0027] y＝HX+z；
[0028] 其中发送信号X是选自于复
星座数字调制(如多进制
相移键控 (MPSK)或多进制
正交幅度调制(MQAM))的NT维符号，每发送符号的平 均功率满足E[XHX]＝ES；z是一个具有NR维的复高斯噪声矢量，均值为0，方 差为N0/2，N0是单边高斯噪声的
功率谱密度。信道矩阵H定义为：
[0029]
[0030] 其中元素hij(i＝1,2,…NR；j＝1,2,…NT)代表第i根接收天线和第j根发射天线 之间衰落信道的增益大小，它是服从均值为0，
实部和
虚部每维的方差为0.5 的循环对称的复高斯随机变量。假设接收端已知信道状态信息(CSI)，且发射 端所有天线采用等功率分配。那么，NR×NT MIMO系统的信道容量大小为：
[0031] 当NR≥NT
[0032] 或者 当NR≤NT
[0033] 其中，ρ＝ES/N0表示接收信号的平均
信噪比。天线选择时，接收端 已知CSI，而发射端根据接收边反馈通道取到CSI，这样，发射端的从NT天线 中选择Lt发射天线与Lt射频链路相连接，接收端的从NR天线中选择Lr发射天 线与Lr射频链路相连接。定义Θ作为所有可能的联合收发天线选择集，大小为 设Hsel为接收端选择Lr根天线发射端选择Lt根天线后的矩阵；那 么，MIMO系统经过天线选择后的信道容量大小为：
[0034]
[0035] Hsel的维数是由发射端和接收端的射频链路数决定；由于各信道不相关， 那么，联合收发天线选择的信道容量可表示为：
[0036]
[0037] 其中，Hsel定义为被选信道矩阵子集，Lt和Lr分别是发射端和接收端的被 选天线数目。
[0038] 进一步，该于空间复用的新型收发天线联合选择方法的容量分析包括：
[0039] 发射天线满集为 接收天线满集为 执行到第 n步天线选择时，若第x根发射天线对信道容量贡献最小的，hx是对应到信道 矩阵中的第x列向量；因此，hx从发射天线集合 中删除，被选发射天线子集 表示为ΓNT-(n+1)＝ΓNT-n-hx；
在接收端也采用类似天线选择策略；若第y根接收天 线是对信道容量贡献最小的，gxy是对应到信道矩阵中的第y行向量；因此，gxy从接收天线集合 中删除，被选接收天线子集表示为 选择后 的信道矩阵Hn表示为一个(NR-n)×(NT-n)的矩阵，其中，n＝
min{NR-Lr,NT-Lt}， 表示已删除了n根发射天线和n根接收天线；这样，选择后的信道矩阵Hn保留 信道矩阵H中的NR-n行向量和NT-n列向量；在每一步天线选择执行中，对信 道容量C贡献最小的一对发射接收天线被挑选出来，hx是对信道容量C贡献最 小的一个(NR-n)×1的列向量，gxy是对信道容量C贡献最小的一个1×(NT-n-1)行 向量；那么，在第(n+1)步执行中，被选择后的信道矩阵Hn+1表示为：
[0040] H′n(H′n)H+hx(hx)H＝Hn(Hn)H；
[0041] (Hn+1)HHn+1+(gxy)Hgxy＝(H′n)HH′n；
[0042] 其中，H＇n是一个(NR-n)×(NT-n-1)矩阵，Hn+1是一个(NR-n-1)×(NT-n-1)矩阵；
[0043] 那么，第(n+1)步信道容量C(Hn+1)表示为：
[0044]
[0045] 应用谢尔莫行列式公式，上式可写成：
[0046]
[0047]
[0048] 又C(H′n)写成:
[0049]
[0050]
[0051] 第(n+1)步信道容量C(Hn+1)重新表示为：
[0052]
[0053] 其中
[0054]
[0055]
[0056]
[0057] 由于C(Hn)是由之前的第n步决定的，因此，信道容量C(Hn+1)的减少量是 由等式的后两项决定，而这两项都是递减函数；那么天线选择算法重新写为：
[0058]
[0059] 在第n+1步执行中，通过联合更新对信道容量贡献最小的发射和接收天线 对gxy和hx，并删除之。
[0060] 本发明提供的基于空间复用的新型收发天线联合选择方法，该选择算法起 始于满集，即包含所有用于天线选择的发射端多天线与接收端多天线间的信道 集合；相比于基于范数的天线选择(Norm-Based Selection，NBS)算法仅考虑向 量的范数，在每一次选择过程中挑选对信道容量贡献最小的一对收发天线进行 删除；建议的新型算法在被选天线数目Lt
附图说明
[0061] 图1是本发明
实施例提供的基于空间复用的新型收发天线联合选择方法
流程图；图中，flag＝0代表最小值p在接收端，故剩下的q-p次天线选择发生 在发射端；flag＝1代表最小值p在发射端，故剩下的q-p次天线选择发生在接 收端；
[0062] 图2是本发明实施例提供的MIMO天线选择系统
框图；
[0063] 图3是本发明实施例提供的了8发8收的MIMO系统，比较新型收发天线 联合选择算法(本发明简称：建议算法)、CJAS算法、及最优算法的计算复杂 度随被选择天线数目的变化曲线；
[0064] 图4是本发明实施例提供的了12发12收的MIMO系统，比较建议算法、 CJAS算法、及最优算法的计算复杂度随被选择天线数目的变化曲线；
[0065] 图5是本发明实施例提供的是在发射端和接收端天线数目都是5，而发射 端和接收端被选择的天线数目都是2的情况下，比较了建议算法、CJAS算法、 NBS算法、随机算法及最优算法的信道容量随信噪比变化曲线；
[0066] 图6是本发明实施例提供的是在发射端和接收端天线数目都是5，而发射 端和接收端被选择的天线数目都是3的情况下，比较了建议算法、CJAS算法、 NBS算法、随机算法及最优算法的信道容量随信噪比变化曲线；
[0067] 图7是本发明实施例提供的是在发射端和接收端天线数目都是7，而发射 端和接收端被选择的天线数目都是2的情况下，比较了建议算法、CJAS算法、 NBS算法、随机算法及最优算法的信道容量随信噪比变化曲线；
[0068] 图8是本发明实施例提供的是在发射端和接收端天线数目都是7，而发射 端和接收端被选择的天线数目都是3的情况下，比较了建议算法、CJAS算法、 NBS算法、随机算法及最优算法的信道容量随信噪比变化曲线；
[0069] 图9是本发明实施例提供的是在发射端和接收端天线数目都是5，而发射 端和接收端的被选择的天线数目都是3，信噪比为8dB的情况下，比较了建议 算法、CJAS算法、NBS算法、随机算法及最优算法的
累积分布函数随信道容 量变化曲线；
[0070] 图10是本发明实施例提供的是在发射端和接收端天线数目都是5，而发射 端和接收端的被选择的天线数目都是3，信噪比为20dB的情况下，比较了建议 算法、CJAS算法、NBS算法、随机算法及最优算法的累积分布函数随信道容 量变化曲线；
[0071] 图11是本发明实施例提供的是在发射端和接收端天线数目都是7，而发射 端和接收端的被选择的天线数目都是3，信噪比为8dB的情况下，比较了建议 算法、CJAS算法、NBS算法、随机算法及最优算法的累积分布函数随信道容 量变化曲线；
[0072] 图12是本发明实施例提供的是在发射端和接收端天线数目都是7，而发射 端和接收端的被选择的天线数目都是3，信噪比为20dB的情况下，比较了建议 算法、CJAS算法、NBS算法、随机算法及最优算法的累积分布函数随信道容 量变化曲线。
具体实施方式
[0073] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例， 对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以 解释本发明，并不用于限定本发明。
[0074] 下面结合附图及具体实施例对本发明的应用原理作进一步描述。
[0075] 本发明建议的新型收发联合天线选择算法很好地在MIMO系统的信道容 量性能与系统复杂度之间寻找一个平衡点；本发明建议算法是基于信道容量最 大化准则的收发联合天线选择；建议算法是一种递减算法，在起始时，所有用 于天线选择的发射端天线与接收端天线间的信道矩阵是满集；每一步选择执行 删除对系统信道容量贡献最小的某一行或者某一列，这种删除对系统的总信道 容量损失非常微小；建议算法计算复杂度可简化为等同于计算一个矩阵二项展 开式；这种基于已知的共轭对称矩阵的矩阵二项展开式，使天线选择算法的计 算复杂度大为下降；本发明的仿真实验结果表明：在相同的信道容量要求下， 建议的新型算法所需的中断概率远小于NBS算法，也比CJAS算法来得小；本 发明的仿真实验结果也表明：建议的新型收发联合天线选择算法信道容量性能 优于NBS算法，与CJAS算法相当，且接近于最优算法；总之，本发明的建议 的算法更充分地利用了多天线的空间复用增益，基本保留了MIMO通信系统容 量；
[0076] 建议的新型收发联合天线选择算法的一些符号变量解释如下：
[0077] ρ＝ES/N0表示接收符号的平均信噪比，ES＝E[xHx]表示每发送符号的平均功 率，N0是单边高斯噪声的功率谱密度；NR为接收端天线数目；NT为发射端天线 数目；Lr为接收端射频连路数(即接收端被选天线数目)；Lt为发射端射频链路数 (即发射端被选天线数目)；H为满集时系统随机信道矩阵；
[0078] 代表发射端天线满集，它意味着发射天线还未进行天线选 择；
[0079] 代表hx所对应的发射天线对信道容量贡献最小；
[0080] 代表hx所对应的发射天线对信道容量贡献最小,将从发射天 线集中删除；
[0081] 代表接收端天线满集，它意味着发射天线还未进行天线选 择；
[0082] 代表gxy所对应的接收天线对信道容量贡献最小；
[0083] 代表gxy所对应的接收天线对信道容量贡献最小,将从接收 天线集中删除；
[0084] H＇n是一个(NR-n)×(NT-n-1)矩阵，Hn+1是一个(NR-n-1)×(NT-n-1)矩阵；
[0085] hx是一个(NR-n)×1的列向量，gxy是一个1×(NT-n-1)行向量；
[0086] 是一个(NT-n-1)×(NT-n-1)矩阵；
[0087] 是一个(NR-n)×(NR-n)矩阵。
[0088] 如图1所示，本发明实施例的基于空间复用的新型收发天线联合选择方法 包括以下步骤：
[0089] 第一步：天线选择初始化，设置多天线通信系统的参数ρ、NR、NT、Lt、 Lr大小,按图2所示描述系统生成起始满集信道矩阵H，清零选择步长n＝0；
[0090] 第二步：计算
[0091] 第三步：从发射天线集 中选择某一个列向量 使得 最小,将其从发射天线集中删除
[0092] for i＝1:NT-n
[0093]
[0094] end
[0095] Jn＝arg maxJn,i
[0096]
[0097] 第四步：对发射天线选择后，信道矩阵更新为
[0098] 第五步：计算
[0099] 第六步：并从接收天线集 中选择某一个行向量 使得 最小，并将其从接收天线集中删除
[0100] for j＝1:NR-n
[0101]
[0102] end
[0103] λn＝arg maxλn,j
[0104]
[0105] 第七步：对接收天线选择后，信道矩阵更新为
[0106] 第八步：如果n<＝min{NR-Lr,NT-Lt}则返回至第2步；如果n>min{NR-Lr, NT-Lt},那么发射端和接收端中具有最小值的一端将停止天线选择，而另外一端 将继续执行天线选择，直到n＝max{NR-Lr,NT-Lt}时天线选择结束。
[0107] 如图2所示，本发明实施例的基于空间复用的新型收发天线联合选择方法 包括以下假设：
[0108] 第一：假设图2系统的所有发射天线和接收天线都能用于通信，接收端收 到的
射频信号已经过解调至基带、时间已同步、
相位和
频率已校准，
输出信号 已经过
匹配滤波器处理；
[0109] 第二：假设在发射端天线和接收端天线间所经历的信道是独立不相关的， 它们均服从相同的平坦瑞利衰落信道，信道矩阵H是一个随机矩阵，每一组符 号处理期间保持不变；
[0110] 第三：假设在每一接收符号时刻，接收端的每根天线接收信号串入加性高 斯白噪声(AWGN)，这些噪声匀服从0均值、实部和虚部的方差为0.5的循环 对称复高斯随机变量；
[0111] 第四：假设接收端已知CSI，从接收端到发射机的信道反馈通道模
块不产 生差错，因此，发射端的CSI是完美的，系统对选择的发射天线进行等功率分 配。
[0112] 如图3、图4所示，本发明实施例的基于空间复用的新型收发天线联合选 择方法包括以下说明：
[0113] 本发明实施例提供的了8发8收和12发12收的两个MIMO系统，比较建 议的新型算法(本发明所有曲线图中匀用“建议算法”表示)、CJAS算法、及 最优算法的计算复杂度随被选择天线数目的变化曲线。当建议算法在收发端被 选天线数目Lt
[0114] 如图5、图6、图7、图8所示，本发明实施例的基于空间复用的新型收发 天线联合选择方法包括以下解释：
[0115] 第一：图5和图6说明了建议算法的性能比CJAS好。同时，建议算法的 性能与最优算法相接近。这是因为建议算法每次删除对系统信道容量贡献最小 的某一行或者某一列，对信道容量的损失量可忽略不计；
[0116] 第二：图5和图6也说明了系统的信道容量随着发射端和接收端被选天线 数目的增多而增加，在发射端和接收端配置多天线可以改善无线通信系统的性 能；
[0117] 第三：图5和图7说明了当射频链路数不变的情况下，在发射端和接收端 配置更多的天线，通过天线选择后系统的性能更好；
[0118] 第四：图6和图8说明了当射频链路数不变的情况下，在发射端和接收端 配置更多的天线，通过天线选择后系统的性能更好；
[0119] 第五：图7和图8说明了随着被选择的发射和接收天线数目增多，MIMO 无线通信系统的性能也得到改善，这是因为MIMO的空间复用能
力在多径传播 环境下具有改善通信系统容量的优势。
[0120] 如图9、图10、图11、图12所示，本发明实施例的基于空间复用的新型 收发天线联合选择方法包括以下解释：
[0121] 第一：图9说明了NBS算法虽然具有较低的复杂度，但是其性能较差，在 Lr＝Lt＝3,NT＝NR＝5，SNR＝8dB时NBS算法信道容量低于7bit/s/Hz的中断概率大 达0.85，而建议算法信道容量低于7bit/s/Hz的中断概率为0.5。建议算法的中 断概率性能比CJAS算法和NBS算法好；
[0122] 第二：图9和图10说明了当信噪比变大的时，MIMO天线选择系统的信 道容量性能更好；
[0123] 第三：图9和图11说明了当射频链路数不变的情况下，在发射端和接收端 配置更多的天线，通过天线选择后系统的性能更好；
[0124] 第四：图11和图12说明了建议算法的性能优于CJAS算法和NBS算法， 可以获得更多的信道容量。
[0125] 本发明的具体实现方法：
[0126] 本发明的矩阵和向量用
黑体，大写字母和小写字母分别代表矩阵和向量， IN和det(·)分别代表N×N的单位矩阵和N×N行列式。上标[.]H表示共轭转置。 上标[.]+表示转置。
[0127] 如图2所示，MIMO天线选择系统框图；在发射端和接收端分别 配置NT根和NR根天线；发射端和接收端分别配置Lt个和Lr个射频链路；假设 所有的发射天线和接收天线都参与天线选择；接收端从NR天线中选择Lr接收 天线，发射端从NT天线中选择Lt发射天线；在发射端天线和接收端天线之间 的所有无线信道是独立的、且服从平坦的瑞利衰落，用随机信道矩阵H表示； 这样，经过匹配滤波器后，接收端接收到的复包络信号为：
[0128] y＝HX+z.(1)
[0129] 其中发送信号X是选自于复星座数字调制(如多进制相移键控 (MPSK)或多进制正交幅度调制(MQAM))的NT维符号，每发送符号的平 均功率满足E[XHX]＝ES；z是一个具有NR维的复高斯噪声矢量，均值为0，方 差为N0/2，N0是单边高斯噪声的功率谱密度。信道矩阵H定义为：
[0130]
[0131] 其中元素hij(i＝1,2,…NR；j＝1,2,…NT)代表第i根接收天线和第j根发射天线 之间衰落信道的增益大小，它是服从均值为0，实部和虚部每维的方差为0.5 的循环对称的复高斯随机变量。假设接收端已知信道状态信息(CSI)，且发射 端所有天线采用等功率分配。那么，NR×NT的MIMO系统的信道容量大小为：
[0132] 当NR≥NT
[0133] 或者 当NR≤NT (3)
[0134] 其中，ρ＝ES/N0表示接收信号的平均信噪比。天线选择时，接收端 已知CSI，而发射端根据接收边反馈通道取到CSI，这样，发射端的从NT天线 中选择Lt发射天线与Lt射频链路相连接，接收端的从NR天线中选择Lr发射天 线与Lr射频链路相连接。定义Θ作为所有可能的联合收发天线选择集，大小为 设Hsel为接收端选择Lr根天线发射端选择Lt根天线后的矩阵；那 么，MIMO系统经过天线选择后的信道容量大小为：
[0135]
[0136] Hsel的维数是由发射端和接收端的射频链路数决定；由于各信道不相关， 那么，联合收发天线选择的信道容量可表示为：
[0137]
[0138] 其中，Hsel定义为被选信道矩阵子集，Lt和Lr分别是发射端和接收端的被 选天线数目。
[0139] 本发明的新型收发联合天线选择方法：
[0140] 一种低复杂度收发联合天线选择算法，基于容量最大化的最优收发联合天 线选择算法的基本思想是，分别从NT根发射天线和NR根接收天线中选择Lt根发射天线和Lr根接收天线，但最优算法的计算复杂度很大，很难满足实时性 要求高的MIMO通信系统；为了减小MIMO系统实现的复杂度，如NBS、CJAS 等一些次最优的快速天线选择算法相继被提出；尽管NBS算法简单，不过NBS 具有较大的信道容量损失；虽然CJAS天线选择算法的容量性能接近于最优算 法，但是当Lt≥NT/2，Lr≥NR/2时，该算法仍存在较大的计算复杂度；为此， 本发明提出了一种低复杂度的新型天线选择算法，该算法复杂度在Lt≥NT/2， Lr≥NR/2时比CJAS要小很多；该选择算法起始于满集，即包含所有用于天线 选择的发射端多天线与接收端多天线间的信道集合，在每一次选择过程中选择 出对信道容量贡献最小的一对收发天线进行删除；这种建议的新型天线选择算 法的性能与最优算法非常接近；
[0141] 收发联合天线选择算法的容量分析：
[0142] 建议的新型天线选择算法的思想是在每一个步中计算各行、列的信道容量， 并从删除对信道容量贡献最小的某一行和某一列，该算法替代了对所有可能组 合的天线子集的信道容量计算；发射天线满集为 接收天线满集 为执行到第n步天线选择时，若第x根发射天线对信道容量贡 献最小的，hx是对应到信道矩阵中的第x列向量；因此，hx从发射天线集合ΓNT中删除，被选发射天线子集表示为ΓNT-(n+1)＝ΓNT-n-hx；在接收端也采用类似天线 选择策略；若第y根接收天线是对信道容量贡献最小的，gxy是对应到信道矩阵 中的第y行向量；因此，gxy从接收天线集合 中删除，被选接收天线子集表 示为 选择后的信道矩阵Hn表示为一个(NR-n)×(NT-n)的矩
阵， 其中，n＝min{NR-Lr,NT-Lt}，表示已删除了n根发射天线和n根接收天线； 这样，选择后的信道矩阵Hn保留信道矩阵H中的NR-n行向量和NT-n列向量； 在每一步天线选择执行中，对信道容量C贡献最小的一对发射接收天线被挑选 出来，hx是对信道容量C贡献最小的一个(NR-n)×1的列向量，gxy是对信道容量 C贡献最小的一个1×(NT-n-1)行向量；那么，在第(n+1)步执行中，被选择后的 信道矩阵Hn+1表示为：
[0143] H′n(H′n)H+hx(hx)H＝Hn(Hn)H； (6)
[0144] (Hn+1)HHn+1+(gxy)Hgxy＝(H′n)HH′n； (7)
[0145] 其中，H＇n是一个(NR-n)×(NT-n-1)矩阵，Hn+1是一个(NR-n-1)×(NT-n-1)矩阵；
[0146] 那么，第(n+1)步信道容量C(Hn+1)表示为：
[0147]
[0148] 应用谢尔莫行列式公式，公式(8)可写成：
[0149]
[0150] 又C(H′n)可以写成:
[0151]
[0152] 第(n+1)步信道容量C(Hn+1)可以重新表示为：
[0153]
[0154] 其中
[0155]
[0156]
[0157]
[0158] 由于C(Hn)是由之前的第n步决定的，因此，信道容量C(Hn+1)的减少量是 由等式(11)的后两项决定，而这两项都是递减函数；那么天线选择算法重新 写为：
[0159]
[0160] 在第n+1步执行中，通过联合更新对信道容量贡献最小的发射和接收天线 对gxy和hx，并删除之。
[0161] 因此，建议新型收发联合天线选择算法可以总结归纳为：
[0162] 第一步：天线选择初始化，设置多天线通信系统的参数ρ、NR、NT、Lt、 Lr大小,按图2所示描述系统生成起始满集信道矩阵H，清零选择步长n＝0；
[0163] 第二步：计算
[0164] 第三步：从发射天线集 中选择某一个列向量 使得 最小,将其从发射天线集中删除
[0165] for i＝1:NT-n
[0166]
[0167] end
[0168] Jn＝arg maxJn,i
[0169]
[0170] 第四步：对发射天线选择后，信道矩阵更新为
[0171] 第五步：计算
[0172] 第六步：并从接收天线集 中选择某一个行向量 使得 最小，并将其从接收天线集中删除
[0173] for j＝1:NR-n
[0174]
[0175] end
[0176] λn＝arg maxλn,j
[0177]
[0178] 第七步：对接收天线选择后，信道矩阵更新为
[0179] 第八步：如果n<＝min{NR-Lr,NT-Lt}则返回至第2步；如果n>min{NR-Lr, NT-Lt},那么发射端和接收端中具有最小值的一端将停止天线选择，而另外一端 将继续执行天线选择，直到n＝max{NR-Lr,NT-Lt}时天线选择结束。
[0180] 对算法进行进一步的解释；为什么建议的天线选择算法可以比NBS算法有 着较好的性能优势；建议算法的容量是由等式(11)的后两项决定，而这两项为 负数，对信道容量呈递减作用；为了满足第(n+1)步天线选择的信道容量C(Hn+1) 最大，就要求这两项之和最小；建议的新型天线选择算法可以写为：
[0181]
[0182] 等式(13)意味着要计算某一行或者某一列对信道容量贡献最小就等同于 gxyBn-1gxyH和hxHDn-1hx同时最大，对信道容量贡献最小的某一行或者某一列将 会被删除，仅剩下的信道矩阵向量就是对信道容量贡献较大的向量；然而，NBS 仅考虑向量的范数，没有考虑某一行或者某一列向量对信道容量的贡献；因此 建议的新型天线选择算法的性能会比NBS优越，主要是因为建议的天线选择算 法是基于信道容量最大化策略，而NBS仅考虑信道矩阵向量的范数，不考虑系 统的信道容量；
[0183] 本发明的天线选择算法的复杂度主要是矩阵gxyBn-1gxyH和hxHDn-1hx的复乘 和复加运算及矩阵Bn-1和Dn-1的求逆运算，对于任意的zAz+，这里上标“+” 表示矩阵或矢量的转置。由二项式公式展开得：
[0184]
[0185] 如果 那么 这意味着当i>j时，只需计算 而不需要计算 因此，zAz+的计算复杂度是n2/2复乘，其中 已知；
[0186] 又由于Bn-1和Dn-1是共轭对称矩阵，根据二项式展开公式(14)更新gxyBn-1gxyH和H -1 -1 -1 -1 H H -1hxDn hx时总的计算复杂度减半；因此，当Bn 和Dn 已知时，计算gxyBn gxy 和hx Dn hx的计算复杂度可以表示为
[0187]
[0188] 在每一步执行天线选择时Bn-1和Dn-1都已知，因此，只需要计算NT-n次 hxHhx和NR-n次gxygxyH；根据二项式公式，计算hxHhx和gxygxyH的复杂度可以表 示为：
[0189]
[0190] 除此之外，在每一步天线选择执行中都要对矩阵Dn和Bn进行求逆运算， 而矩阵求逆运算的复杂度为O(n3)，其中n是不可逆矩阵的秩；因此，计算Bn-1和Dn-1的复杂度可以表示为：
[0191]
[0192] 因此，建议的新型天线选择算法总的计算复杂度可以表示为：
[0193] vC＝Vhg+VBD (18)
[0194] 经过进一步简化，等式(18)可以近似为：
[0195]
[0196] 图3和4所示为CJAS算法，建议算法和最优算法的三种天线选择算法复 杂度比较；当Lt＝1,Lr＝1，建议的新型天线选择算法的计算复杂度比最优算法高， 这是因为建议算法需要对等式(12)做(NT-1)×(NR-1)次的比较，而最优算法只需 要选取信道矩阵中最大的数，没有涉及到行列式运算；随着被选天线数目Lt和 Lr的逐渐增大，可以发现建议算法的计算复杂度比最优算法小；当Lt≥NT/2，Lr≥ NR/2时，建议算法的计算复杂度也比CJAS小；这是因为当Lt≥NT/2，Lr≥NR/2 时，建议算法只需经过执行max{NR-Lr,NT-Lt}步就可以选出符合条件的天线 数，而CJAS需经过执行max{Lr,Lt}步才能选出符合条件的天线数；如，当NT＝ NR＝N＝8，Lt＝Lr＝L＝6时，CJAS需经过6步，而建议算法只需要2步就能选 出所需的天线数目。在这种情况下，建议算法只需要进行2次的复乘和复加运 算，而CJAS需要经过6次的复乘和复加运算；因此，建议算法在Lt≥NT/2，Lr≥ NR/2时的计算复杂度会比CJAS小，该算法特别适用于被选天线数目Lt≥NT/2， Lr≥NR/2的情况。
[0197] 通过以下的仿真试验对本发明的应用效果做进一步的说明：
[0198] 通过仿真实验验证建议算法的性能，仿真中同时比较最优算法，NBS，CJAS 及本发明建议的新型收发联合天线选择算法的系统性能曲线；仿真中使用的信 道是准静态平坦瑞利衰落信道，通过仿真给出信道容量随信噪比变化曲线及给 定信噪比的情况下累积分布函数(PDF)随信道容量的变化曲线图；为了能够 更好说明各个算法的性能，本发明采用蒙特卡洛仿真实验，每张曲线图都是由 1000个独立信道矩阵H的随机
种子数仿真生成的结果；
[0199] 图5和图6分别为Lr＝Lt＝2，NT＝NR＝5和Lr＝Lt＝3，NT＝NR＝5条件下，4种算法的 信道容量随信噪比变化曲线，信噪比SNR从0dB到20dB；图5和图6说明了建议 算法的性能比CJAS好；同时，建议算法的性能与最优算法相接近，这是因为建 议算法每次删除对系统信道容量贡献最小的某一行或者某一列，对信道容量的 损失量可忽略不计；图5和图6也说明了系统的信道容量随着发射和接收天线数 目的增多而增加，在发射端和接收端配置多根天线可以改善无线通信系统的性 能；图5和图7及图6和图8都说明了随着被选发射和接收天线数目的增多，无线 通信系统的信道容量性能能得到改善，这是因为MIMO无线通信系统在多径传 播环境下具有改善通信系统容量的潜质；
[0200] 图9是在Lr＝Lt＝3,NT＝NR＝5，SNR＝8dB条件下，各天线选择算法的累积分布 函数随信道容量的变化曲线图；图10是在Lr＝Lt＝3,NT＝NR＝5，SNR＝20dB条件下， 各天线选择算法的累积分布函数随信道容量的变化曲线图；图9和图10说明了 NBS算法虽然具有较低的复杂度，但是其性能较差，在SNR＝8dB时，NBS算法 信道容量低于7bit/s/Hz的中断概率达0.85，CJAS算法信道容量低于7bit/s/Hz的 中断概率为0.5，而建议算法信道容量低于7bit/s/Hz的中断概率仅为0.4；显然， 建议算法的性能优于CJAS算法和NBS算法；
[0201] 在Lr＝Lt＝3,NT＝NR＝7条件下，图11和图12分别描绘了SNR＝8dB和SNR＝20dB 时各算法的累积分布函数随信道容量的变化曲线图，可以清楚的看出，建议算 法的性能同样优于CJAS算法和NBS算法；因此，建议的新型天线选择算法可以 获得更多的天线空间复用增益，提高通信系统容量；
[0202] 天线选择技术是一种降低MIMO系统复杂度的有效方法；为了更好的在 MIMO系统的性能与系统复杂度之间寻找一个平衡点，本发明提出了一种新型 的收发联合天线选择算法；本发明建议的新型收发联合天线选择算法是基于信 道容量最大化准则；仿真结果说明了建议的天线选择算法的性能与最优天线选 择算法的性能非常接近，同时，建议算法的计算复杂度远小于最优算法；该算 法尤其适合于被选天线数目Lt≥NT/2，Lr≥NR/2的场合应用，能很好的获取多 天线的空间复用增益；
[0203] 在本发明的实施例中，天线选择中Bn-1和Dn-1都已知，因此，只需要计算 NT-n次H H H Hhxhx和NR-n次gxygxy；根据二项式展开公式，计算hxhx和gxygxy的复杂度可以表示为：
[0204]
[0205] 将公式(A3～A5)代入到(A1)，可以得到：
[0206]
[0207] 公式(A3～A5)分别为：
[0208]
[0209]
[0210]
[0211] 在每一步天线选择中都要对矩阵Dn和Bn进行求逆运算，而矩阵求逆运算 的复杂度为O(n3)，其中n是可逆矩阵的秩；因此，计算Bn-1和Dn-1的复杂度 可以表示为：
[0212]
[0213] 同理，将公式(A3～A5)代入到(A6)，可以得到：
[0214]
[0215] 因此，建议天线选择算法总的计算复杂度可以表示为：
[0216]
[0217] 以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发 明的精神和原则之内所作的任何
修改、等同替换和改进等，均应包含在本发明 的保护范围之内。
