CN109861793A - 一种基于二阶统计量的空频分组码信号盲识别方法 - Google Patents
一种基于二阶统计量的空频分组码信号盲识别方法 Download PDFInfo
- Publication number
- CN109861793A CN109861793A CN201910090771.7A CN201910090771A CN109861793A CN 109861793 A CN109861793 A CN 109861793A CN 201910090771 A CN201910090771 A CN 201910090771A CN 109861793 A CN109861793 A CN 109861793A
- Authority
- CN
- China
- Prior art keywords
- ofdm
- signal
- order statistic
- peak value
- determined
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 238000010276 construction Methods 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 abstract description 3
- 238000000926 separation method Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 201000004569 Blindness Diseases 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
Landscapes
- Radio Transmission System (AREA)
Abstract
本发明公开了一种基于二阶统计量的空频分组码信号盲识别方法。该方法具体包括以下步骤:构造不同接收天线接收信号的二阶统计量,然后进行傅里叶变换(DFT)得到频谱图;通过检测频谱图是否存在峰值区分SM‑OFDM和AL‑OFDM信号。其中峰值检测方法具体是根据二阶统计量的周期性,可以推测有16个尖峰,因此可以通过检测峰值的数量,若检测出16个峰值,则判定为AL‑OFDM信号,否则判定为SM‑OFDM信号;或者通过检测峰值间距离,若距离为二阶统计量序列的1/16,则判定为AL‑OFDM信号,否则判定为SM‑OFDM信号。
Description
技术领域
本发明属于信号处理领域中非协作通信信号处理技术,具体是指一种基于二阶统计量的空频分组码信号盲识别方法。
背景技术
信号盲识别技术一直是近几年研究的热点和难点,特别是在频谱资源日益紧张的信息时代。随着频谱资源的紧缺,多载波技术由于其频带利用率高和抗多径干扰能力强的优点而应用的越来越广泛,因此多载波多天线识别技术也成为研究的热点。近些年多载波OFDM系统与空间编码(Space Block Coding,SBC)结合,已经成为IEEE 806.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(0)和s(1)分别为:
s(0)=[x(0,0),x(1,0),x(2,0),…]
s(1)=[x(0,1),x(1,1),x(2,1),…]
因此,发射信号经过信道后第i根天线接收信号可以表示为:
其中,i=0,1,…nr-1,hfi(l)是第f个发射天线和第i个接收天线的信道脉冲响应,Lh代表未知路径数,n(i)(m)是高斯白噪声序列,s(f)(m)是s(f)的第m个序列。
所述的构造不同接收天线的接收信号二阶统计量,然后对二阶统计量进行DFT得到频谱图具体为:构造第k时刻和(k+N/2)时刻二阶统计量,其中N为OFDM子载波个数,k为接收信号时刻,然后对二阶统计量进行DFT变换。
为了进一步说明问题,本发明所述的构造不同接收天线的接收信号二阶统计量,具体为:将两个不同接收天线的接收信号ri(k)和ri′(k),构造二阶统计量形式为y(k)=E[r(i)(k)r(i')(k+N/2)],i≠i′,ri(k)表示第i根接收天线第k个时刻接收信号,ri′(k+N/2)表示第i′根接收天线第(k+N/2)个时刻接收信号。
由于随机变量可以表示为其均值和另一个代表偏差的零均值的变量的和,因此SM和AL信号的y(k)可以表示为:
ySM(k)=E{ySM(k)}+ψSM(k)
yAL(k)=E{yAL(k)}+ψAL(k)
其中,E{yξ(k)}是yξ(k)均值,ψξ(k)代表均值的偏差,ξ∈{AL,SM}。
对于yξ(k),ψξ(k)认为是隐藏在E{yξ(k)}中的噪声项,当不存在噪声时由于发射信号是独立同分布的,很容易得到:
ySM=[0,0,0,…]
对于AL信号,不考虑噪声,E{yAL(k)}表达式为:
根据发射端信号相关性,E[x(b,0)(n)x(b,1)(n+N/2)]表达式为:
结合发射端的相关性,E{yAL(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)}中的噪声项,当不存在噪声时:
yAL=[C1,0,0…C2,0,0…C3,0,0…C4,0,0…]
其中,C1=-1、C2=-j、C3=1和C4=j,中间是15个连零,因此序列yAL是一个周期函数,其周期为16。
所述的对二阶统计量进行离散傅里叶变换(DFT)得到频谱图,具体为:假设Y=[Y(0),Y(1),…Y(K-1)]是向量y的DFT变换:
因此,YSM和YAL分别为:
YSM(n)=ΨSM,n=0,1,…K-1
所述的通过检测频谱图是否存在峰值,如果存在峰值则判定为AL-OFDM,如果没有峰值则判定是SM-OFDM。具体为:如果有16个尖峰,则判定为AL-OFDM信号,如果没有峰值则判定是SM-OFDM,具体为:检测频谱图主尖峰的数量,如果主尖峰为16个,则判定为AL-OFDM信号,否则判定为SM-OFDM信号。
算法流程如下:
输入:接收信号r(i)(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(i)(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所示。从图中曲线可以发现,调制方式对算法基本无影响。
Claims (5)
1.一种基于二阶统计量的空频分组码信号盲识别方法,其特征在于:构造不同接收天线的接收信号的二阶统计量,然后对二阶统计量进行离散傅里叶变换得到频谱图,通过检测频谱图是否存在峰值,如果存在峰值则判定为AL-OFDM信号,如果没有峰值则判定是SM-OFDM。
2.根据权利要求1所述的盲识别方法,其特征在于,所述的构造不同接收天线的接收信号二阶统计量,然后对二阶统计量进行离散傅里叶变换得到频谱图,具体为:构造第k个时刻和第(k+N/2)个时刻二阶统计量,其中N为OFDM子载波个数,k为接收信号时刻,然后对二阶统计量进行离散傅里叶变换。
3.根据权利要求2所述的盲识别方法,其特征在于,所述的构造不同接收天线的接收信号二阶统计量,具体为:将两个不同接收天线的接收信号ri(k)和ri′(k),构造二阶统计量形式为y(k)=E[r(i)(k)r(i')(k+N/2)],i≠i′,ri(k)表示第i根接收天线第k个时刻接收信号,ri′(k+N/2)表示第i′根接收天线第(k+N/2)个时刻接收信号。
4.根据权利要求1至3中任一所述的盲识别方法,其特征在于,所述的通过检测频谱图是否存在峰值,如果存在峰值则判定为AL-OFDM信号,如果没有峰值则判定是SM-OFDM的具体方法为:检测频谱图主尖峰的数量,如果主尖峰数量为16个,则判定为AL-OFDM信号,否则判定为SM-OFDM信号。
5.根据权利要求1至3中任一所述的盲识别方法,其特征在于,所述的通过检测频谱图是否存在峰值,如果存在峰值则判定为AL-OFDM信号,如果没有峰值则判定是SM-OFDM的具体方法为:检测频谱图尖峰间距离,如果相邻尖峰间距离为则判定为AL-OFDM信号,否则判定为SM-OFDM信号,K为y(k)序列数量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910090771.7A CN109861793A (zh) | 2019-01-29 | 2019-01-29 | 一种基于二阶统计量的空频分组码信号盲识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910090771.7A CN109861793A (zh) | 2019-01-29 | 2019-01-29 | 一种基于二阶统计量的空频分组码信号盲识别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109861793A true CN109861793A (zh) | 2019-06-07 |
Family
ID=66896905
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910090771.7A Pending CN109861793A (zh) | 2019-01-29 | 2019-01-29 | 一种基于二阶统计量的空频分组码信号盲识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109861793A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108512793A (zh) * | 2017-02-28 | 2018-09-07 | 畅想科技有限公司 | Ofdm信号存在检测 |
CN114244659A (zh) * | 2021-11-17 | 2022-03-25 | 西安电子科技大学 | 一种非高斯噪声干扰下的mimo系统空时模式识别方法及系统 |
CN114244660A (zh) * | 2021-11-17 | 2022-03-25 | 西安电子科技大学 | 一种脉冲噪声干扰下mimo-ofdm系统sfbc类型识别方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103763222A (zh) * | 2014-01-16 | 2014-04-30 | 西安电子科技大学 | 一种mimo信号盲检测中的信道模糊度去除方法 |
CN108416674A (zh) * | 2018-02-12 | 2018-08-17 | 上海翌固数据技术有限公司 | 时频谱的应用方法和设备 |
CN108923829A (zh) * | 2018-07-11 | 2018-11-30 | 中国人民解放军海军航空大学 | 一种基于相关函数的空频分组码盲识别方法 |
-
2019
- 2019-01-29 CN CN201910090771.7A patent/CN109861793A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103763222A (zh) * | 2014-01-16 | 2014-04-30 | 西安电子科技大学 | 一种mimo信号盲检测中的信道模糊度去除方法 |
CN108416674A (zh) * | 2018-02-12 | 2018-08-17 | 上海翌固数据技术有限公司 | 时频谱的应用方法和设备 |
CN108923829A (zh) * | 2018-07-11 | 2018-11-30 | 中国人民解放军海军航空大学 | 一种基于相关函数的空频分组码盲识别方法 |
Non-Patent Citations (2)
Title |
---|
YAHIA A. ELDEMERDASH 等: "Blind Identification of SM and Alamouti STBC-OFDM Signals", 《IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS》 * |
闫文君 等: "空时分组码的二阶循环平稳特性分析", 《海军航空工程学院学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108512793A (zh) * | 2017-02-28 | 2018-09-07 | 畅想科技有限公司 | Ofdm信号存在检测 |
CN108512793B (zh) * | 2017-02-28 | 2022-08-26 | 畅想科技有限公司 | 用于ofdm信号存在检测的方法、检测器及介质 |
CN114244659A (zh) * | 2021-11-17 | 2022-03-25 | 西安电子科技大学 | 一种非高斯噪声干扰下的mimo系统空时模式识别方法及系统 |
CN114244660A (zh) * | 2021-11-17 | 2022-03-25 | 西安电子科技大学 | 一种脉冲噪声干扰下mimo-ofdm系统sfbc类型识别方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102279390B (zh) | 一种低信噪比雷达辐射源信号脉内调制识别方法 | |
CN109861793A (zh) | 一种基于二阶统计量的空频分组码信号盲识别方法 | |
CN101984612B (zh) | 基于压缩感知的非连续正交频分复用信道估计方法 | |
CN104410590A (zh) | 一种基于压缩感知的短波ofdm抑制干扰联合信道估计方法 | |
CN106713190B (zh) | 基于随机矩阵理论和特征阈值估计的mimo发射天线数目盲估计算法 | |
CN106059639B (zh) | 基于矩阵盖尔圆的发射天线数目盲估计方法 | |
CN104393963A (zh) | 基于循环平稳特性的空时分组码mc-cdma信号盲识别方法 | |
CN106357574A (zh) | 基于顺序统计量的bpsk/qpsk信号调制盲识别方法 | |
CN112347871A (zh) | 一种用于通信载波监视系统的干扰信号调制识别方法 | |
CN102546061A (zh) | 基于小波变换的自适应时频空穴检测方法 | |
CN105187352A (zh) | 一种基于ofdm前导的整数频偏估计方法 | |
CN106506427A (zh) | 一种基于folp的stbc‑ofdm信号盲识别方法 | |
CN108923829A (zh) | 一种基于相关函数的空频分组码盲识别方法 | |
CN104486284B (zh) | 基于增强型六维64psk星座的正交频分复用方法 | |
CN103368700B (zh) | 基于特征量预估计的时滞相关空时码模式盲识别方法 | |
CN101848183B (zh) | 多入多出正交频分复用系统中信道估计方法及装置 | |
CN103475603B (zh) | 基于序参量非正交变换的通信系统盲信道估计方法 | |
CN103138856B (zh) | 一种检测干扰的方法及装置 | |
US8374229B2 (en) | Method for the detection and generation of a useful signal and associated devices and communications system | |
CN103117821A (zh) | 一种基于瑞利商的加权协作频谱感知方法 | |
Wang et al. | Indoor fingerprinting localization based on fine-grained CSI using principal component analysis | |
CN110351211A (zh) | 一种短波通信信道估计中导频图案位置的智能搜索方法 | |
CN103152790B (zh) | 基于相关性分簇的两级融合调制识别方法 | |
CN113726710B (zh) | Ofdm信号的发送方法、接收方法、装置、存储介质及电子设备 | |
CN110011744A (zh) | 端口检测方法、系统和终端 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190607 |
|
RJ01 | Rejection of invention patent application after publication |