一种提高频谱利用率的无线通信系统的设计方法
【专利摘要】一种提高频谱利用率的无线通信系统的设计方法，由源信号、信道、DSP分离系统和输出信号组成。N路未知源信号通过无线信道后，由M个天线接收混合信号，混合信号由DSP分离系统进行分离，恢复出源信号，实现无线通信。DSP分离系统是此系统的核心部分，它通过时频分析估计出源信号的数目和信道矩阵，然后判断源信号数目N与接收天线数目M之间的关系，根据不同的关系选择不同的分离算法或者估计算法。本发明系统能进一步提高频谱利用率，放宽对信号时隙、频率、码字等的限制，实现系统对用户或目标的灵活接入。
【专利说明】
-种提高频谱利用率的无线通信系统的设计方法
技术领域
[0001] 本发明属于无线通信技术领域，具体设及线性均匀天线接收条件下，源信号数目 未知的无线通信系统的设计，是一种适用面更广的高效频谱利用率的无线通信方法。
【背景技术】
[0002] 在无线通信系统中，一般假设信号源的数目是已知的、固定的，然而在实际应用 中，无线通信环境是复杂的，源信号数目是未知的，甚至是随时间变化的，例如移动通讯系 统中，在某一个特定的区域内，用户接入或者断开的数目在不断的变化，因此在某一时间段 内的用户数目是无法预知的；深海无源探测网络中，某一有效面积内的目标数目也是未知 的。在运种实际情况下，传统的假设源信号数目已知的无线通信系统模型可能不再适用，本 发明提出一种新的无线通信系统的设计，该系统适用于在线性均匀天线接收条件下，源信 号未知的情况。
[0003] 常用的无线通信系统模型一般假设信号源数目已知，且对信号的时间、频率、码字 等加 W-些限制，在接收端接收运些信号的混合信号，通过信号在时间、频率、码字等的特 征下，提取或者恢复源信号，实现无线信道的复用，运样的系统采用的复用技术包括频分复 用技术(FDM)、时分复用技术(TDM)和码分复用技术(CDM)等，它们能较好的解决频谱匿乏的 问题，但在时隙上，或是在频率上，或是在码字上对信号有所限制，采用运些复用技术的无 线通信系统对用户或目标的接入数量的灵活性有较严重的限制。2009年提出的统计复用无 线通信系统(WSDM)能较好的解决传统复用技术在时间、频率、码字上的限制，提高频谱利用 率，但是WSDM适用于源信号数目已知且与接收天线数目相同的情况，极大的限制了WSDM的 应用环境。
【发明内容】
[0004] 本发明要解决的问题是:现有无线通信系统如WSDM都是假设源信号数目已知的情 况，而实际无线通信系统是复杂的，源信号数目往往是未知的，新的无线通信系统能较好的 解决运一问题，且能进一步的提高无线频谱利用率，保证对源信号数目估计的准确性和对 源信号恢复的精确性，实现无线通信信道复用。
[0005] 本发明的技术方案为:一种提高频谱利用率的无线通信系统的设计方法，源信号 数目N未知，接收天线线性均匀排列，设定A)源信号空间时频分布矩阵存在单源点;B)信道 矩阵的任意两个列矢量线性独立;源信号通过无线信道后，由Μ个天线接收混合信号，混合 信号由DSP分离系统进行分离，恢复出源信号，其中DSP分离系统进行时频分析，首先求取混 合信号的空间时频分布STFD矩阵，筛选出特殊的时频点，所述特殊的时频点处的STFD矩阵 仅有一个对角元素非零，其余元素均为零，运些时频点为单源点，单源点上的能量由一路源 信号形成，其他的源信号的贡献为零，所有单源点对应的波达角呈现直线聚类特性，通过峰 值检测和聚类算法估计出聚类中屯、的数量和方向，即为源信号的数目和波达角的值，由此 得到源信号数目和信道矩阵，根据源信号的数目和接收天线的数目的关系选择相应的恢复 算法，恢复出源信号，实现无线通信。
[0006] S(t) = (Sl(t)，S2(t)，，SN(t))T 为 N 路源信号组成的矢量，X(t) = (Xl(t)，X2(t)，，XM (t))T为Μ路混合信号组成的矢量其中t=l，2，，Ts，Ts为通信信号的持续时间，源信号矩阵S =[3(1)，3(2)，，3巧5)]，混合信号矩阵乂=^(1)，义(2)，，义巧5)]，接收天线为1个阵元组成的 直线型线性均匀天线阵，设Θι，θ2,，Θν分别为N路源信号的波达角，则天线阵接收的混合信号 表示如下：
[0007] X=AS (1)
[000引其中A e C"x、为信道矩阵，具有范德蒙德结构，有如下形式：
[0009]
(2),
[0010] 采用线性时频变换和二次时频分布相结合的方法来进行时频分析，并根据通信环 境的实际情况选择二次时频分布的类型，W保证源信号在时频域具有稀疏性，
[0011] DSP分离系统的处理具体步骤为：
[0012] 1)根据二次时频分布类型，求混合信号的空间时频分布矩阵，即ST抑矩阵，对噪声 点进行去除：
[OOU] 对源信号矢量34) = (3如），32(〇，，3如）八其51'抑矩阵胖3山門为：
[0020] 降噪方法如下：
[0021] 设时频点(ta，时)为噪声点，则有：
[0022]
(3)
[002引其中I I · II表示F-范数，ε功接近于0的正阔值，maxf I I Wx(ta,f川表示对时间ta选 择混合矢量的STFD矩阵的F-范数在所有频点f中的最大值，如果时频点满足公式(3)噪声点 的条件，则将其去除；
[0024] 2)单源点筛选：
[0025] 通过分析混合信号的STFD矩阵来间接的分析源信号的STFD矩阵，依次通过两个筛 选定理来对单源点进行筛选：
[00%]定理1:如果时频点(t，f)属于单源点域，则满足：
[0027]
W
[0028] 其中I · I表示实数的绝对值或者复数的模值，62为接近于零的正阔值，若此点（t, f)为第η个源信号的单源点，则此时第η个信号的波达角为：
[0029]
(5)
[0030] 其中 m，ke{l，2，，M};
[0031] 定理2:如果时频点(t，f)属于单源点域，则满足：
[00创
辟
[0033] 其中E3为接近于零的正阔值，知:1'((，/)=片戸^化/).知,化/)为一种降低交叉项干扰 的时频分布，><"「(/，/)为第i路混合信号Xi(t)的谱图，谱图指Xi(t)的短时傅里叶变换的模 值的平方，即
[0034]
[0035] 其中*表示取共辆，丫（t)为窗函数；
[0036] 对时频点先经定理1进行初选，初选结果再经定理2进行再选，得到单源点；
[0037] 3)峰值检测和聚类算法：
[0038] 筛选出单源点后，计算出波达角，所有单源点对应的波达角呈现直线聚类特性，通 过峰值检测和聚类算法估计出聚类中屯、的数量和方向，即为源信号的数目和波达角的值， 再通过波达角计算出信道矩阵；
[0039] 4)通过上述步骤估计出源信号的数目N和信道矩阵A，根据源信号的数目N和接收 天线的数目Μ的关系选择相应的恢复算法恢复源信号。
[0040] 所述恢复算法分Ξ种情况：
[0041] a)、当Ν<Μ时，即为超定的情况，设定了信道矩阵的任意两个列矢量线性独立，则信 道矩阵A列满秩，根据关系式Χ=AS，源信号的估计为：
[0042]
(7)
[0043] b)、当N=M时，即为适定的情况，信道矩阵A列满秩，源信号的估计为：
[0044]
(8)
[0045] C)、当N〉M时，即为欠定的情况，对信道矩阵A已知的欠定情况采用欠定盲源分离算 法，包括基于最小均方误差的波束赋形法和基于活跃信源数估计的ISSR算法。
[0046] 本发明提出一种新的无线通信系统的设计方案，先对源信号的数目进行检测，同 时对信道矩阵进行估计，然后利用源信号数目和信道矩阵信息对源信号进行恢复，实现信 道复用。本发明的设计方案对源信号只要求其空间时频分布矩阵存在单源点，由于只要源 信号的频率之间有少量的不重合，即可保证单源点的存在，因此本发明方案允许源信号之 间有大量的频率交叉，即该方案能同时传输多路频率交叠的源信号，相比于传统的复用方 式如FDM、TDM、CDM，能进一步提高频谱利用率，放宽对信号时隙、频率、码字等的限制，实现 系统对用户或目标的灵活接入。
【附图说明】
[0047] 图1为本发明提出的无线通信系统设计方案的原理图。
[0048] 图2为本发明接收天线的线性均匀阵列模型。
[0049] 图3为本发明方法中源信号数目估计和信道矩阵估计的流程图。
[0050] 图4为本发明实施例中对波达角的示意图，（a)为源信号的波达角；（b)为波达角的 估计。
[0051] 图5为本发明实施例中源信号和输出信号的时频分布图，（a)、（b)、（c)为3路源信 号实部幅度值，（d)、（e)、（f)为源信号的wigner-ville分布，（g)为混合信号的wigner- ville分布，化)为降低交叉项干扰后的混合信号的时频分布，（i)为筛选出的单源点的时频 分布，（j)、化)、（1)为输出信号的时频分布，(m)、(n)、（o)为分别为输出信号的实部幅度值。
[0052] 图6为本发明设计的无线通信系统在不同信噪比环境下的信道矩阵估计的MSE (地）。
【具体实施方式】
[0053] 图1是本发明提出的无线通信系统设计方案示意图，无线通信系统由源信号、信 道、DSP分离系统和输出信号组成。源信号未知，设有N路，N是事先不可知的，源信号通过无 线信道后，由Μ个天线接收混合信号，混合信号由DSP分离系统进行分离，恢复出源信号，实 现无线通信。DSP分离系统是此系统的核屯、部分，它首先通过时频分析估计出源信号的数目 和信道矩阵，也就是混合矩阵，然后判断源信号数目Ν与接收天线数目Μ之间的关系，根据不 同的关系选择不同的分离算法或者估计算法，来估计恢复出源信号。
[0054] 设S(t) = (Sl(t)，S2(t)，，SN(t))T，X(t) = (Xl(t)，X2(t)，，XM(t))T，y(t) = (yi(t)， y2(t)，，yN(t))T分别为N路源信号组成的矢量、M路混合信号组成的矢量、N路输出信号组成 的矢量，输出信号即源信号的估计，其中t = l，2，，Ts，Ts为信号的持续时间；源信号矩阵S = [3(1)，3(2)，，3(了5)]，混合信号矩阵乂=^(1)，义(2)，，义(了5)]，输出信号矩阵￥=[7(1)，7 (2)，，y(Ts)]。本发明的接收天线考虑线性均匀阵列（uniform linear array，ULA)，如图2 所示，Μ个阵元组成直线型天线阵，相邻的两个阵元之间距离为d，源信号Sn(t)位于天线阵 列的远场区域，则可认为源信号W平面波的形式到达天线阵列，此时不考虑接收信号间的 幅度差异，仅考虑时延差，设θ?，02,，θΝ(θηΕ [-31/2,31/2])分别为N路源信号的波达角 (Direction of Arrival，D0A)，即源信号的入射方向，则天线阵列接收的混合信号可W表 示如下：
[005引 X=AS (1)
[0化6]其中AsC心》为信道矩阵（或混合矩阵），具有范德蒙德(Vandermonde)结构，有如 下形式：
[0057]
(2)
[0058] 由信道矩阵A的表示式可知，对源信号数目的估计和对信道矩阵的估计等价于估 计波达角θι，02，，Θν的数目及其数值，由于本发明假设信道矩阵的任意两个列矢量线性独 立，则有气:;鳴(V/V 7'），即达波角互不相等。
[0059] 下面详细介绍本发明设计的DSP分离系统，DSP分离系统采用时频分析方法对信道 矩阵及源信号数目进行估计，常用的时频分析有线性时频变换和二次时频分布，线性时频 变换包括短时傅里叶变换（STFT)、小波变换等，二次时频分布包括Cohen类时频分布、 wigner-ville分布等。二次时频分布较线性时频变换在时频域上具有更高的能量聚集性， 使得一些非平稳信号在时频域上更加稀疏，但是二次时频分布存在交叉项的干扰会影响估 计性能，而线性时频变换不存在交叉项的问题，本发明DSP分离系统采用线性时频变换和二 次时频分布相结合的方法来进行时频分析，结合两者的优点，既保持了信号在时频域上能 量聚集性，又消除了交叉项的干扰。
[0060] 实际应用中应根据实际情况和场景合理的选择二次时频分布的类型，W保证源信 号在时频域具有一定的稀疏性。本发明Wwigner-Ville分布为例进行说明。对第η个源信号 sn(t)，自wi即er-ville分布为：
[0061]
(3)
[00创由上式易看出而'托/) = A。,、,化/)'，'即Sn(t)的自wi即er-villeA,,,、;,化/)分布 一定是实值。对两个不同的源信号sn(t)、sw(t)，n声W，它们的互wigner-ville分布为：
[0063]
(4)
[0064] 可W看出不同源信号的互wigner-ville分布满足A 化./ ) = 化./ )。结合式 (3)和(4)，对源信号矢量34) = (3加），32(〇，，3^*))了，可^得到矩阵：
[00化]
(5)
[0066] Ws (t，f)称为源信号的空间时频分布矩阵，空间时频分布矩阵简称STFD矩阵，它的 对角线元素为源信号的自wigner-ville分布，非对角线元素为互wigner-ville分布。通过 胖3(*，門的结构，可^把￥3(*，門对应的时频点分为^类：
[0067] 1、自源点：
[006引定义1:如果(te，fe)称为自源点，则需满足Ws(te，fe)对角元素不全为零且非对角元 素为零，即Ws(te，fe)为对角阵，此时Ws(te，fe)反映出源信号自身的能量聚集性。此处为现有 技术定义，不再详述。
[0069] 2、互源点：
[0070] 定义2:如果(心，片)称为互源点，则需满足Ws(tn，f。)对角元素为零且非对角元素不 全为零，此时Ws(tn，f。)反映出源信号之间的能量聚集性。此处为现有技术定义，不再详述。
[0071] 3、一般时频点：
[0072] 定义3:若在(tY，fv)处，Ws(tY，fY)不满足上述两个性质，则称为一般时频点，即在 (tY，fy)处源信号的自身能量和信号间的能量是重叠的。
[0073] 进一步的，对于混合矢量X(t) = (Xl(t)，X2(t)，，XM(t))T，可W用同样的方式定义 其ST抑矩阵，得到混合信号的ST抑矩阵Wx(t，f):
[0074]
巧)
[0075] 则通过混合矢量与源信号矢量的关系和wigner-ville分布的定义可得：
[0076]
(7)
[0077] 由于在通信系统的接收端只有混合信号是已知的，因此唯一可W获取的是混合信 号的STFD矩阵Wx(t，f)，由公式(6)结合公式(3)(4)得到。然而不是所有的时频点处的Wx(t, f)都可W直接使用，因为一些时频点处的能量太小对于估计信道矩阵没有明显的作用，甚 至会对估计产生干扰，因此首先应该去除运些能量太小的时频点。为了简便，称运些能量太 小的点为噪声点，一种常见的降噪方法如下：
[007引如果时频点（ta，时)为噪声点，贝IJ
[0079]
(8)
[0080] 其中II · II表示F-范数(Frobenius范数），ει为接近于0的正阔值，仿真试验中通 常设为0.1，maxf II Wx(ta，f川表示对时间点ta，选择混合矢量的STFD矩阵的F-范数在所有 频点f中的最大值。如果时频点满足噪声点的条件，则将其去除。
[0081 ]如前所述，Wx(t，f)是唯一可获取的STFD矩阵，需要通过分析混合矢量的STFD矩阵 来间接的分析源信号的STFD矩阵，在一些特殊的时频点处能找出Wx(t，f)与Ws(t，f)的联系。
[0082] 定义4:若自源点（ta，fa)满足在此点上的能量主要由一路源信号形成，其他的源信 号的贡献为零，则称此点为单源点，否则称为多源点。本发明的单源点与在"盲源分离"领域 中的"单源点"的含义一致，而与在其他领域例如"最短路径"中的"单源点"的含义不一致， 特此用定义4标明。
[0083] 容易得出，在单源点（ta，fa)处，Ws(ta，fa)的主对角线只有一个元素非零，非对角元 素都为零，则可W得到如下等式：
[0084]
(9)
[0085] 其中源信号si(t)为能量的主要贡献者，ill表示信道矩阵A的第i列矢量，Ωι为所有 满足上式的时频点组成的集合，为了简便，称之为第i个信号的单源点域。实际应用中很难 有信号严格满足单源点的条件，因此只要有一个对角元素明显大于其他元素，则认为此时 频点为单源点，式(9)可W近似成立。
[0086] 下面为对单源点的筛选方法进行具体说明。
[0087] 定理1:如果时频点(t，f)属于单源点域，则必须满足：
[00则
（…）
[0089] 其中I · I表示实数的绝对值或者复数的模值，62为接近于零的正阔值，例如0.1。 公式（10)适用于所有二次时频分布。若此点（t，f)为第η个源信号的单源点，则此时第η个信 号的波达角为：
[0090]
〇\)
[0091] 其中 m，ke{l，2，，M}。
[0092] 证明：假设(t，f)为第C个信号的单源点，设ac = -3isin(0c)，其中Θ。为第C个信号的 波达角，则根据式(9)，有
[0096] 上述等式即为时频点（t，f)成为单源点的一个必要条件，然而实际应用中无法严 格满足上述等式，因此引入一个接近于零的正阔值62,上述必要条件修改为式（10)。通过式 (12)还可得：
[0097]
(14)
[009引其中111，4￡{1，2，，1}。由于假设曰。=-313山(目。），则可得达波角如式（11)所示。
[0099] 定理1证毕。
[0100] 定理1可W作为筛选单源点的条件，由于只考虑了STFD矩阵对角元素，因此计算量 较小，下面验证定理1的条件对非单源点的严格性。
[0101] 设非单源点，则由式(7)可得：
[0102]
[0103] 若点处也满足定理1中单源点的必要条件，贝U
[0104]
[0105] 如果上式对任意的/\,、(^/|)吨，9￡{1，2，，《，口声9都满足，则可得(^ =日。，呵^， 由此可得ep = 0q，Vp，g，与前述的假设波达角互不相等相矛盾。一些特殊的非单源点也可能 满足定理1中单源点的必要条件，例如一些非单源点若满足
[010引其中Vp，qe{l，2，，N}，p声q，贝lJ(t'，f'）即可W满足定理l中单源点的必要条件，但 是上述条件的满足等价于满足
[0109]
[0110] 式（17)对非单源点的要求是十分严格的，因此非单源点很难满足定理1的条件，定 理1中条件可W很好的用作单源点的筛选。但由于此条件仅为必要条件，且选择精度受阔值 的影响，所W可将定理1的条件作为单源点的一个粗筛选。
[0111] 由于线性时频变换没有交叉项的干扰，因此本发明考虑线性时频变换与二次时频 分布的结合，W短时傅里叶变换(STFT)为例，混合信号Xi(t)的短时傅里叶变换为：
[0112]
(㈱
[0113] 其中*表示取共辆，丫（t)为窗函数。短时傅里叶变换模值的平方为谱图，记作
则线性时频变换与二次时频分布的结合得到的一种 降低交叉项干扰的时频分布3
i于时频分布，此 公式是现有的"减小交叉项干扰的分布(RID)"中的一个，易知扣:了 (?，/')在一定程度上保持 了 wigner-ville分布的能量聚集性，降低了交叉项的干扰。
[0114] 下面构造基于短时傅里叶变换的空间时频分布矩阵，记为Wfn"，/)，短时傅里叶 变换不存在交叉项，因此可令Wf-(/，/)非对角线元素为0,对角线元素为信号的短时傅里叶 变换，即
[0115]
(19}
[0116] 其中(wffT化/成表示矩阵wrT(/，/)的第i行第j列元素。降低交叉项干扰的时频分 布h /)的空间时频分布矩阵记为(^，/)，表示为
[0117]
(20)
[011引其中0为哈德蒙德化adamard)积。则Wfw典乃矩阵的对角线元素为降低交叉项干 扰的时频分布，非对角线元素为0。
[0119] 定理2:如果时频点(t，f)属于单源点域，则必须满足：
[0120]
(21)
[0121] 其中E3为接近于零的正阔值。
[0122] 证明：同定理1的证明，假设(t，f)为第C个信号的单源点，设ac = -3isin(0c)，则
[0123]
[0124] 由上述矩阵结构可得
[0125]
隣）
[0126] 式(23)即为单源点应满足的必要条件。在实际应用中，式(23)很难严格满足，因此 引入一个接近于0的正阔值63,将式(23)的条件修改为式(21)。
[0127] 定理2证毕。
[01%]定理2同定理1，也只需考虑对角元素，因此计算量较小。下面验证定理2对非单源 点的严格性。
[0129] 设非单源点，则由式(20)可得：
[0130]
[0131] 若点处也满足定理2中单源点的必要条件，则等价于满足式（16)，由于前 面已经验证了式（16)对非单源点的严格性，所W定理2的条件对非单源点也具有同样的严 格性。
[0132] 单源点的筛选结合定理1和定理2,先经定理1进行初选，初选结果再经定理2进行 再选，确定单源点后通过式（11)，可W计算出波达角，所有单源点对应的波达角呈现直线聚 类特性，通过峰值检测和聚类算法，如K-means可W估计出聚类中屯、的数量和方向，即为源 信号的数目和波达角的值，通过波达角可W计算出信道矩阵。图3所示为源信号数目估计和 信道矩阵估计的流程图。
[0133] 通过上述流程，可W有效的估计出源信号的数目N和信道矩阵A，为了恢复出源信 号，需要根据源信号的数目N和接收天线的数目Μ的关系选择相应的恢复算法，如图1，分Ξ 种情况。
[0134] 1、当Ν<Μ时，即为超定的情况。由于假设信道矩阵的任意两个列矢量线性独立，则A 列满秩，根据关系式X=AS，可根据下式给出源信号的估计：
[0135]
(24)
[0136] 2、当N=M时，即为适定的情况。同样，A列满秩，可由下式给出源信号的估计：
[0137]
(25)
[013引 3、当N〉M时，即为欠定的情况。对A已知的欠定问题通常为欠定盲源分离问题的第 二步，已存在很多有效的算法，如基于最小均方误差的波束赋形法(MMSE Beamforming)、基 于活跃信源数估计的ISSR算法等，本发明后续仿真分析WISSR算法为例。
[0139] W上是本发明提出的无线通信系统的详细实施方法和步骤，下面通过仿真实验来 验证本系统的有效性和可靠性。
[0140] 考虑更复杂的欠定的情况，设定3个源信号，2个接收天线的情景，即N=3，M = 2。源 信号为线性调制化FM)信号，其中两路为单分量LFM信号、一路为多分量LFM信号，Ξ路源信 号分别W波达角-15°、-45°、60°入射到天线阵列。在单源点的筛选中，设参数61 = 0.1、62 = 0.5、ε3 = 0.015。
[0141] 将源信号波达角θη表示成复数:e^·，则在复平面上形成的散点与原点的连线如图4 (a)所示，运Ξ条直线与实部正轴的夹角即为波达角，图4(b)表示波达角的估计矣构成的复 数6心《在复平面上的散点图，为了便于观察，每个点均乘W-个(0,3)上均匀分布的随机数。 由图4可知看出，波达角的估计组成的散点图自然形成Ξ组，每组对应一条直线，与左图原 始波达角的方向基本一致，后续采用峰值检测及K-means聚类算法能准确的估计出波达角 的数量W及达波角的值，并W此来确定源信号数目N和信道矩阵A。此时，可W确定源信号数 目N=3,则判断系统属于欠定的情况，采用ISSR算法对信号进行恢复。如图5所示，（a)、（b)、 (C)分别为Ξ路源信号实部幅度值，（d)、（e)、（f)为相应的wigner-ville分布，（g)为混合信 号的wigner-ville分布（包括自源点、互源点和一般时频点处的时频分布），化）为采用式 (20)进行降低交叉项干扰后的混合信号的时频分布，（i)为通过粗筛选和再筛选后单源点 的时频分布，（j)、化）、（1)分别为Ξ路输出信号的时频分布，（m)、（η)、（0)分别为Ξ路输出 信号的实部幅度值。
[0142] 比较图5中源信号和输出信号的时域和时频域可知，输出信号实现了对源信号的 恢复，尽管幅值不一致，即存在幅度不确定性，但是不影响系统正常通信。观察(g)、化）、（i) 可W看到，化)实现了对交叉项干扰的抑制，保证了系统的正确性和稳定性，（i)表示单源点 筛选的两个相关定理能很好的提取出单源点，如图（i)所示，图中的点（除了少部分干扰点） 均取自某一源信号的自源点处，且运些点处有且仅有一路源信号的wigner-ville分布非 零。
[0143] 在存在噪声的环境中，验证本发明系统的可行性和实用性。仿真试验中，选用高斯 白噪声进行测试，仅考虑估计信道矩阵的准确性，采用原始信道矩阵与估计信道矩阵的均 方误差作为评价指标，即：
[0144]
(26)：
[0145] 其中A为A的估计矩阵，此处假设i中的列矢量排列已调整为与A-致(一种有效的 调整方法可W分别计算A的每一列矢量与A的第i列矢量的均方误差，均方误差最小对应的 列矢量即为A的第i列矢量的估计）。设参数ei = 0.1、E2 = 0.5,分别取63 = 0.01,0.015,0.02， 进行100次蒙特卡洛实验，则在不同信噪比环境下信道矩阵估计的MSE(地)如图6所示。
[0146] 由图6可W看出当63 = 0.01,0.015,0.02时信道矩阵的估计准确度都很高，在信噪 比比较低的时候也实现了信道矩阵的精确估计，因此验证了本系统在存在高斯白噪声环境 下的可行性和实用性。从图6中还可W得出有关阔值选择的结论，由于E3是系统进行再筛选 时使用的阔值，对系统的性能影响较大，因此固定ει、ε2,考虑不同63对信道矩阵估计性能的 影响，图6反映出Ε3的取值可W在[0.01,0.02]区间选取，在低信噪比情况下，阔值越小越 好；当信噪比较高时(大于36地），63 = 0.015是最佳选择，运也是前面仿真中取63 = 0.015的 依据。
[0147]综上所述，本发明提出的无线通信系统设计方案是可行的、可靠的，而且是高效 的，适用于在线性均匀天线接收条件下，源信号未知的情况。在系统实现过程中，放宽了对 源信号的时隙、频率、码字等的限制，实现系统对用户或目标的灵活接入，进一步提高了频 谱利用率。
【主权项】
1. 一种提高频谱利用率的无线通信系统的设计方法，其特征是源信号数目N未知，接收 天线线性均匀排列，设定A)源信号空间时频分布矩阵存在单源点;B)信道矩阵的任意两个 列矢量线性独立;源信号通过无线信道后，由Μ个天线接收混合信号，混合信号由DSP分离系 统进行分离，恢复出源信号，其中DSP分离系统进行时频分析，首先求取混合信号的空间时 频分布STro矩阵，筛选出特殊的时频点，所述特殊的时频点处的STro矩阵仅有一个对角元 素非零，其余元素均为零，这些时频点为单源点，单源点上的能量由一路源信号形成，其他 的源信号的贡献为零，所有单源点对应的波达角呈现直线聚类特性，通过峰值检测和聚类 算法估计出聚类中心的数量和方向，即为源信号的数目和波达角的值，由此得到源信号数 目和信道矩阵，根据源信号的数目和接收天线的数目的关系选择相应的恢复算法，恢复出 源信号，实现无线通信。2. 根据权利要求1所述的一种提高频谱利用率的无线通信系统的设计方法，其特征是s (t) = (Sl(t)，S2(t)，···，SN(t))T为N路源信号组成的矢量，X(t) = (Xl(t)，X2(t)，···，XM(t))T 为Μ路混合信号组成的矢量，其中t = 1，2，…，Ts，Ts为通信信号的持续时间，源信号矩阵S = [8(1)，8(2)，~，8(1'5)]，混合信号矩阵乂=[以1)，奴2)，~，以1' 5)]，接收天线为1个阵元组成 的直线型线性均匀天线阵，设θ^θ%···，ΘΝ分别为N路源信号的波达角，则天线阵接收的混合 信号表不如下： X=AS (1) 其中4SCT*为信道矩阵，具有范德蒙德结构，有如下形式：采用线性时频变换和二次时频分布相结合的方法来进行时频分析，并根据通信环境的 实际情况选择二次时频分布的类型，以保证源信号在时频域具有稀疏性， DSP分离系统的处理具体步骤为： 1)根据二次时频分布类型，求混合信号的空间时频分布矩阵，即STro矩阵，对噪声点进 行去除： 对源信号矢量8(〇 =(幻(〇,82(认~，洲(〇)7，其51?0矩阵1(扒〇为：混合信号的STFD矩阵：降噪方法如下： 设时频点（ta，fe)为噪声点，则有：其中| | . | |表示F-范数，￡1为接近于0的正阈值，maxf I |Wx(U，f) I I表示对时间ta选择混 合矢量的SIFD矩阵的F-范数在所有频点f中的最大值，如果时频点满足公式(3)噪声点的条 件，则将其去除； 2) 单源点筛选： 通过分析混合信号的STro矩阵来间接的分析源信号的STro矩阵，依次通过两个筛选定 理来对单源点进行筛选： 定理1:如果时频点(t，f)属于单源点域，则满足：其中I · I表示实数的绝对值或者复数的模值，ε2为接近于零的正阈值，若此点（t，f)为 第η个源信号的单源点，则此时第η个信号的波达角为：其中 m，ke{l，2，'"，M}; 定理2:如果时频点(t，f)属于单源点域，则满足：其中ε3为接近于零的正阈值，Ρ:1/) = <￡%/)·Α^,/)为一种降低交叉项干扰 的时频分布，<、，/)为第i路混合信号Xl(t)的谱图，谱图指^⑴的短时傅里叶变换的模 值的平方，即其中*表示取共辄，γ (t)为窗函数； 对时频点先经定理1进行初选，初选结果再经定理2进行再选，得到单源点； 3) 峰值检测和聚类算法： 筛选出单源点后，计算出波达角，所有单源点对应的波达角呈现直线聚类特性，通过峰 值检测和聚类算法估计出聚类中心的数量和方向，即为源信号的数目和波达角的值，再通 过波达角计算出信道矩阵； 4)通过上述步骤估计出源信号的数目N和信道矩阵A，根据源信号的数目N和接收天线 的数目Μ的关系选择相应的恢复算法恢复源信号。3.根据权利要求1或2所述的一种提高频谱利用率的无线通信系统的设计方法，其特征 是所述恢复算法分三种情况： a) 、当Ν〈Μ时，即为超定的情况，设定了信道矩阵的任意两个列矢量线性独立，则信道矩 阵Α列满秩，根据关系式X=AS，源信号的估计为： 8: = (^?)^1 AaX (7) b) 、当N=M时，即为适定的情况，信道矩阵A列满秩，源信号的估计为： §: = Α-4Χ (8) c) 、当Ν>Μ时，即为欠定的情况，对信道矩阵Α已知的欠定情况采用欠定盲源分离算法， 包括基于最小均方误差的波束赋形法和基于活跃信源数估计的ISSR算法。
【文档编号】H04W72/04GK106060834SQ201610593962
【公开日】2016年10月26日
【申请日】2016年7月26日 公开号201610593962.1, CN 106060834 A, CN 106060834A, CN 201610593962, CN-A-106060834, CN106060834 A, CN106060834A, CN201610593962, CN201610593962.1
【发明人】苏巧, 魏以民, 沈越泓
【申请人】中国人民解放军理工大学