一种基于二阶统计量的空频分组码信号盲识别方法与流程

本发明属于信号处理领域中非协作通信信号处理技术，具体是指一种基于二阶统计量的空频分组码信号盲识别方法。

背景技术：

信号盲识别技术一直是近几年研究的热点和难点，特别是在频谱资源日益紧张的信息时代。随着频谱资源的紧缺，多载波技术由于其频带利用率高和抗多径干扰能力强的优点而应用的越来越广泛，因此多载波多天线识别技术也成为研究的热点。近些年多载波ofdm系统与空间编码(spaceblockcoding,sbc)结合，已经成为ieee806.12n和ieee806.12e核心技术，广泛应用在第四代移动通信中。总的来说，sbc与ofdm结合有两种方式，一种是stbc-ofdm，它是把数据符号先ofdm编码然后再sbc编码，而另一种是sfbc-ofdm，是先进行sbc编码，然后再ofdm编码。目前大多数识别算法的研究对象是stbc-ofdm信号，主要针对stbc-ofdm信号盲识别发射端编码方式，主要的方法有统计特性、循环谱、k-s检测和folp等，从识别效果来看，识别效果较理想，算法的鲁棒性也较强。但是sfbc-ofdm盲识别问题相关的文献还处于空白。本发明的算法是针对sfbc-ofdm信号中空间复用信号空频分组码(sm-ofdm)和al空频分组码(al-ofdm)信号进行识别。

技术实现要素：

本发明的内容是在全盲条件下，对sfbc-ofdm信号中空间复用信号空频分组码(sm-ofdm)和al空频分组码(al-ofdm)信号进行识别。

为解决上述技术问题，本发明是通过以下技术方案实现的：构造不同接收天线的接收信号的二阶统计量，然后对二阶统计量进行离散傅里叶变换(dft)得到频谱图，通过检测频谱图是否存在峰值，如果存在峰值则判定为al-ofdm信号，如果没有峰值则判定是sm-ofdm。

sm和al信号编码方式详细的方案如下：

对sm通信系统，将整个数据集d以2n为单位分块，以块为单位进行信号发射，第b块数据集通过空间复用生成两个新的数据向量，并通过天线0和天线1发射：

q^(b,0)＝[d^(b)(0),d^(b)(2),…d^(b)(2n-2)]

q^(b,1)＝[d^(b)(1),d^(b)(3),…d^(b)(2n-1)]

其中，q^(b,i)分别为第i(i＝0,1)根天线第b个长度为n的发射信号。

对于al通信系统，将整个数据集d以n为单位分块，以块为单位进行信号发射，第b块数据向量通过空间编码为新的数据向量u为：

u^(b,0)＝[d^(b)(0),-d^(b)*(1),d^(b)(2),…d^(b)(n-2),-d^(b)*(n-1)]

u^(b,1)＝[d^(b)(1),d^(b)*(0),d^(b)(3),d^(b)*(2)…d^(b)(n-1),d^(b)*(n-2)]

其中，u^(b,i)分别为第i(i＝0,1)根天线第b个长度为n的发射信号。

在sbc编码后，进行ofdm变换，具体方法是：对得到的新数据集u和q做ifft变换，并且加上前缀，可得到ofdm块，当信号为al编码时x^(b,f)为：

当信号为sm时x^(b,f)为：

其中，n＝-v,…n-1，f＝0,1代表发射天线的编号，u^(b,f)(k)和q^(b,f)(k)分别表示u^(b,f)和q^(b,f)第k个数据符号。

s由不同块发射信号x^(b)组合而成，因此s⁽⁰⁾和s⁽¹⁾分别为：

s⁽⁰⁾＝[x^(0,0),x^(1,0),x^(2,0),…]

s⁽¹⁾＝[x^(0,1),x^(1,1),x^(2,1),…]

因此，发射信号经过信道后第i根天线接收信号可以表示为：

其中，i＝0,1,…nr-1，hfi(l)是第f个发射天线和第i个接收天线的信道脉冲响应，lh代表未知路径数，n⁽ⁱ⁾(m)是高斯白噪声序列，s^(f)(m)是s^(f)的第m个序列。

所述的构造不同接收天线的接收信号二阶统计量，然后对二阶统计量进行dft得到频谱图具体为：构造第k时刻和(k+n/2)时刻二阶统计量，其中n为ofdm子载波个数，k为接收信号时刻，然后对二阶统计量进行dft变换。

为了进一步说明问题，本发明所述的构造不同接收天线的接收信号二阶统计量，具体为：将两个不同接收天线的接收信号rⁱ(k)和r^i′(k)，构造二阶统计量形式为y(k)＝e[r⁽ⁱ⁾(k)r⁽ⁱ'⁾(k+n/2)],i≠i′，rⁱ(k)表示第i根接收天线第k个时刻接收信号，r^i′(k+n/2)表示第i′根接收天线第(k+n/2)个时刻接收信号。

由于随机变量可以表示为其均值和另一个代表偏差的零均值的变量的和，因此sm和al信号的y(k)可以表示为：

y^sm(k)＝e{y^sm(k)}+ψ^sm(k)

y^al(k)＝e{y^al(k)}+ψ^al(k)

其中，e{y^ξ(k)}是y^ξ(k)均值，ψ^ξ(k)代表均值的偏差，ξ∈{al,sm}。

对于y^ξ(k)，ψ^ξ(k)认为是隐藏在e{y^ξ(k)}中的噪声项，当不存在噪声时由于发射信号是独立同分布的，很容易得到：

y^sm＝[0,0,0,…]

对于al信号，不考虑噪声，e{y^al(k)}表达式为：

根据发射端信号相关性，e[x^(b,0)(n)x^(b,1)(n+n/2)]表达式为：

结合发射端的相关性，e{y^al(k)}表达式为：

其中，π(k,k+n/2-l-l′)代表指示函数，其含义是s^(b,f)(k-l)和s^(b,f′)(k+n/2-l′)在同一个ofdm块中。

对于y^ξ(k)，ψ^ξ(k)认为是隐藏在e{y^ξ(k)}中的噪声项，当不存在噪声时：

y^al＝[c1,0,0…c2,0,0…c3,0,0…c4,0,0…]

其中，c1＝-1、c2＝-j、c3＝1和c4＝j，中间是15个连零，因此序列y^al是一个周期函数，其周期为16。

所述的对二阶统计量进行离散傅里叶变换(dft)得到频谱图，具体为：假设y＝[y(0),y(1),…y(k-1)]是向量y的dft变换：

因此，y^sm和y^al分别为：

y^sm(n)＝ψ^sm,n＝0,1,…k-1

所述的通过检测频谱图是否存在峰值，如果存在峰值则判定为al-ofdm，如果没有峰值则判定是sm-ofdm。具体为：如果有16个尖峰，则判定为al-ofdm信号，如果没有峰值则判定是sm-ofdm，具体为：检测频谱图主尖峰的数量，如果主尖峰为16个，则判定为al-ofdm信号，否则判定为sm-ofdm信号。

算法流程如下：

输入：接收信号r⁽ⁱ⁾(k)和r^(i′)(k)；

1)计算接收信号的二阶统计量估计值

2)对二阶统计量进行dft变换，得到频域图；

3)如果主尖峰为16个，则判定为al-ofdm信号(h1为真)，否则判定为sm-ofdm信号(h0为真)。

所述的通过检测频谱图是否存在峰值，如果存在峰值则判定为al-ofdm，如果没有峰值则判定是sm-ofdm。具体为：如果尖峰之间距离为则判定为al-ofdm信号，如果没有峰值则判定是sm-ofdm，k为y(k)序列数量。

输入：接收信号r⁽ⁱ⁾(k)和r^(i′)(k)；

1)计算接收信号的二阶统计量估计值

2)对二阶统计量进行dft变换，得到频域图；

3)如果尖峰之间距离为则判定为al-ofdm信号(h1为真)，如果没有峰值则判定是sm-ofdm(h0)。

本发明方法能识别在全盲条件下al-ofdm信号和sm-ofdm信号，且在低信噪比下识别效果较好；且不受调制方式等影响，方法的鲁棒性较好。

附图说明

图1是本发明所述方法的总体流程图。

图2是本发明所述的sfbc-ofdm信号传输框图。

图3是sm-ofdm和al-ofdm二阶统计量的频谱图

图4是不同信噪比下sm-ofdm和al-ofdm识别概率。

图5实施案例中不同调制方式对识别性能对比。

具体实施方式

本实施方式的总体流程图如图1所示，发射sfbc-ofdm信号传输框图如图2所示。

本实施例所述方法实现过程如下：

1)采样，初始化数据；

2)分别计算sm-ofdm和al-ofdm信号的二阶统计量；

3)对二阶统计量进行dft变换，得到频域图；

4)采用识别算法中方法一或者方法二峰值检测算法检测峰值；

5)重复步骤1～步骤4进行1000次蒙特卡洛仿真，计算正确识别概率。

实例中无特殊说明，仿真参数如下设置：子载波n＝64，循环前缀v＝5，采用qpsk调制方式，对于每一条传输路径，信道假定为频率选择信道，其中包括3个统计独立的抽头，每一个抽头建模为均值为0的复高斯随机变量，其功率服从指数pdf，其中l＝0,1,…,lh-1。信噪比定义为其中为每根发射天线的功率，为高斯噪声功率。正确识别概率pc(λ|λ)衡量算法性能。每个信噪比做1000次蒙特卡洛实验。

sm-ofdm和al-ofdm二阶统计量的频谱图如图3所示。对于al信号，|y(n)|确实每隔有一个尖峰，主峰值的数量为16个，在每组主峰值旁边还有两个次峰值，数量也是16个，且相邻主峰之间的间隔也是

图4为p(λ|λ)与snr关系，图中列出了不同的ofdm块数目下正确识别概率，从图4分析，增大ofdm块数目，al识别性能得到很大提高，主要是由于块数目增大，的估计值更准确，其峰值更明显，另外，sm信号不随着ofdm块数目变化，主要是由于即使块数目增大，其互相关函数也没有展现峰值。

图5为sm-ofdm和al-ofdm在不同调制方式下下识别概率图。其中调制方式为qpsk、8-psk、16-psk和64-qam。其余仿真条件采用默认条件，实验仿真结果如图5所示。从图中曲线可以发现，调制方式对算法基本无影响。