CN107046512B - 一种simo-fir信道阶数估计方法 - Google Patents
一种simo-fir信道阶数估计方法 Download PDFInfo
- Publication number
- CN107046512B CN107046512B CN201710126715.5A CN201710126715A CN107046512B CN 107046512 B CN107046512 B CN 107046512B CN 201710126715 A CN201710126715 A CN 201710126715A CN 107046512 B CN107046512 B CN 107046512B
- Authority
- CN
- China
- Prior art keywords
- matrix
- detection statistic
- simo
- channel
- statistic
- 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.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/0238—Channel estimation using blind estimation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/0212—Channel estimation of impulse response
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
本发明公开了一种SIMO‑FIR信道阶数估计方法,包括:(1)构造接收信号矩阵和信道矩阵;(2)计算接收信号自协方差矩阵;(3)对协方差矩阵做特征分解;(4)构造检测统计量;(5)求解判决门限值;(6)确定检测统计量属性;(7)确定信道阶数。本发明的有益之处在于:在本发明的方法中,采用随机理论中的方法构造检测统计量,相比之前基于概率论的传统方法,该检测统计量具有精确的分布特性,具有快速收敛特性,因此可以利用有限的采样点得到可靠的估计性能;在本发明的方法中,判决方式是基于假设检验的,当虚警概率接近于0时,其正确识别概率趋近于1。
Description
技术领域
本发明涉及一种信道阶数的估计方法,具体涉及一种频率选择性衰落的SIMO-FIR(Single Input Multiple Output-Finite Impulse Response,单输入多输出-有限冲激响应)信道的多径数目的估计方法,可用于通信信号盲识别和认知无线电,属于无线通信技术领域。
背景技术
下一代无线通信系统提供更高的传输速率、更好的服务质量和更高的频谱效率,多天线技术可以为实现这一目标提供技术支撑。在通信信号盲识别场景中,需要对非合作通信信号参数进行盲估计和识别,信道阶数估计是其中所涉及到的盲信道估计和信号解码等关键技术的前提。
针对SIMO-FIR信道的阶数估计问题,最早的方法为基于信息论准则的模式选择方法。M.Wax和T.Kailath在“Detection of signals by information theoretic criteria”中提出了基于赤池信息论准则(Akaike Information Theoretic Criterion,AIC)和基于最小描述长度(Minimum Description Length,MDL)准则的算法。此类方法假设采集的信号必须是零均值独立同分布的高斯随机变量,且加性噪声是与信道输出信号不相关的高斯白噪声。但是由于实际的通信系统很难满足这个条件,往往会导致过估或者欠估的问题。Liavas A.P在“Blind channel approximation:effective channel orderdetermination”中提出了一种基于数值分析的信道阶数检测准则(Liavas Criteria,LC)。该准则可以对接收信号的协方差矩阵进行稳定的信号和噪声子空间的分解,但是其不足在于,在低信噪比的条件下识别性能差,不适用于有用信号功率较低的通信场景。王玉红等在“利用特征极值比的盲信道阶数估计方法”中提出了一种特征极值比定理(Maximumeigenvalue Minimum eigenvalue Ratio Ratio,MMRR),利用子空间方法中噪声子空间矢量构成特殊矩阵的奇异性与信道阶数之间的关系实现信道阶数的估计。但是该方法需要计算特殊矩阵,计算量大,不适用于实际无线通信系统。
发明内容
针对现有技术的不足,本发明旨在提供一种SIMO-FIR信道阶数估计方法,采用随机理论中的方法构造检测统计量,相比之前基于概率论的传统方法,该检测统计量具有精确的分布特性,具有快速收敛特性,因此可以利用有限的采样点得到可靠的估计性能;另外,判决方式是基于假设检验的,当虚警概率接近于0时,其正确识别概率趋近于1。
为了实现上述目的,本发明采用如下技术方案:
一种SIMO-FIR信道阶数估计方法,所述方法使用的系统模型为1×Nr的SIMO系统,包括一个发射机、一个接收机和一个侦听机,每个子信道表示为一个FIR滤波器,其中,发射机的天线数目为1,侦听机的天线数目为Nr,信道的实际阶数为L,发送信号采用4QAM调制方式,接收窗长为N,N≥L;所述方法包括如下步骤:
S1构造接收信号矩阵和信道矩阵:
侦听机第i根接收天线在时刻n接收到的观测输出为yi(n),将接收窗长内N个连续的观测输出向量堆叠在一起,构成NrN×1维的接收矩阵YN(n);
S2计算接收信号自协方差矩阵R:
R=E{YNYN Η};
S3对R进行特征值分解:R=UΛUΗ,然后将R的特征值从大到小排列;其中,U为NNr×NNr的酉矩阵,Λ是R的特征值构成的对角矩阵;
S4构造检验统计量:
lk表示步骤S3中经过排列的R的第k个特征值;
S5求解判决门限
其中,表示第一类或第二类Tracy-Widom分布的逆累积分布函数,Pfa表示预先设定的虚警概率;参数和是Nr和N的函数;
S6确定检测统计量属性:
如果检测统计量Uk大于判决门限γk,则将检测统计量Uk判定为信号统计量;
如果检测统计量Uk小于或等于判决门限γk,则将对应检测统计量Uk判定为噪声统计量;
统计得到信号统计量的总数K;
S7确定SIMO系统信道阶数估计值:
需要说明的是,步骤S5中,参数和通过下式确定:
当β=1时代表实数噪声:
当β=2时代表复数噪声:
本发明的有益效果在于:
1、在本发明的方法中,采用随机理论中的方法构造检测统计量,相比之前基于概率论的传统方法,该检测统计量具有精确的分布特性,具有快速收敛特性,因此可以利用有限的采样点得到可靠的估计性能。
2、在本发明的方法中,判决方式是基于假设检验的,当虚警概率接近于0时,其正确识别概率趋近于1。
附图说明
图1是本发明的方法适用的系统模型示意图;
图2是本发明的方法的流程图;
图3是本发明的方法在不同采样点下正确识别概率仿真效果图;
图4是本发明的方法在不同信噪比下正确识别概率仿真效果图。
具体实施方式
以下将结合附图对本发明作进一步的描述,需要说明的是,本实施例以本技术方案为前提,给出了详细的实施方式和具体的操作过程,但本发明的保护范围并不限于本实施例。
如图1所示,本发明使用的系统模型包括:一个发射机、一个接收机和一个侦听机。其中,发射机的天线数目为1,侦听机的天线数目为Nr,信道阶数为L,发送信号采用4QAM调制方式,接收窗长为N,N≥L。
本实施例中,侦听机的天线数目Nr=8,信道阶数L=3,接收窗长N=5。
参照图2,本发明方法实现的步骤如下:
S1构造接收信号矩阵和信道矩阵:
本发明研究的为1×Nr的SIMO系统,每个子信道可表示为一个FIR滤波器。记第i根接收天线在时刻n接收到的观测输出为yi(n),其中
s(n)为时刻n的发送信号,发送信号之间相互独立;hi=[hi0(n),hi1(n),…,hiL(n)]表示第i条子信道对应的传输系数;vi(n)为均值为零、方差为σ2的加性高斯白噪声,与接收信号相互独立。
将接收窗长内N个连续的观测输出向量堆叠在一起,即得到YN(n)=[y(n),y(n-1),…,y(n-N-1)]Τ。建立的最终系统模型为:
上式中, 矩阵是一个NrN×(N+L)维的Toeplitz矩阵,定义为:
其中,
S2计算接收信号协方差矩阵R=E{YNYN Η}。其中
其中,RS表示发射信号的自协方差矩阵,是列满秩矩阵。
S3子空间分解
3.1)对接收信号的自协方差矩阵R进行特征值分解:
R=UΛUΗ;
其中,U为NNr×NNr的酉矩阵,Λ是R的特征值构成的对角矩阵。
3.2)将R的特征值从大到小排列:
l1>l2>…>lN+L≥lN+L+1≥…≥lNNr
R的特征值可以被分为信号子空间Ls={l1,l2,...,lN+L,}和噪声子空间Ls为的特征值按降序排列的集合,Ln为σ2I的特征值按降序排列的集合。
S4构造检测统计量:
本实施例中,噪声V为高斯白噪声,则σ2I=E{VVH}是一个Wishart矩阵。根据随机矩阵理论,利用Wishart矩阵的最大特征值和迹的平均值,可以构造检验统计量如下:
当Nr,N→∞,时,特征值lN+L+1与纯噪声Wishart矩阵的最大特征值渐进服从同一Tracy-Widom分布。因此,UN+L+1的分布函数可以表示为
其中,参数和是Nr和N的函数。
当β=1时代表实数噪声:
当β=2时代表复数噪声:
步骤5,求解判决门限:
其中,表示第一类或第二类Tracy-Widom分布的逆累积分布函数,Pfa表示预先设定的虚警概率。
步骤6,确定检测统计量属性:
如果检测统计量Uk大于判决门限γk,则将检测统计量Uk判定为信号统计量;
如果检测统计量Uk小于或等于判决门限γk,则将对应检测统计量Uk判定为噪声统计量;
步骤7,确定SIMO系统的信道阶数的估计值:
确定信号统计量的数目K,并将该数目与接收窗长的差值作为SIMO-FIR系统信道的估计值。
本发明的方法所带来的有益效果可以通过以下仿真来进一步进行说明。
一、仿真条件
1个具有单天线发射机,1个具有8天线的侦听机,发射机和侦听机之间的信道为频率选择性信道,信道阶数为5,发射信号调制方式为4QAM,虚警概率Pfa=0.1。
二、仿真内容及仿真结果
仿真1,信噪比为5dB,分别采用本发明的方法(简称WME)、基于模式选择的MDL方法和AIC方法、LC方法,以及MMRR方法在不同采样点下进行正确识别概率的仿真对比。
仿真结果:如图3所示,在采样点少的情况下,本发明的方法识别性能优于其它方法,可以确保满足通信信号盲识别的工程要求
仿真2,采样点数为64,分别采用本发明的方法、基于模式选择的MDL方法和AIC方法、LC方法,以及MMRR方法在不同信噪比条件下进行正确识别概率的仿真对比。
仿真结果:如图4所示,在小信噪比情况下,本发明的方法识别性能优于其它方法,可以确保满足通信信号盲识别的工程要求。
由此可见,本发明的方法能够在样本规模较小或者信噪比较低的条件下,实现高识别率的信道阶数估计,解决了低信噪比环境下信道阶数识别率低的问题。
本发明的方法可用于通信信号盲识别和认知无线电,可以保证在军事通信对抗和SIMO-FIR信道中的适用性。
对于本领域的技术人员来说,可以根据以上的技术方案和构思,作出各种相应的改变和变形,而所有的这些改变和变形都应该包括在本发明权利要求的保护范围之内。
Claims (2)
1.一种SIMO-FIR信道阶数估计方法,所述方法使用的系统模型为1×Nr的SIMO系统,包括一个发射机、一个接收机和一个侦听机,每个子信道表示为一个FIR滤波器,其中,发射机的天线数目为1,侦听机的天线数目为Nr,信道的实际阶数为L,发送信号采用4QAM调制方式,接收窗长为N,N≥L;其特征在于,所述方法包括如下步骤:
S1构造接收信号矩阵和信道矩阵:
侦听机第i根接收天线在时刻n接收到的观测输出为yi(n),将接收窗长内N个连续的观测输出向量堆叠在一起,构成NrN×1维的接收矩阵YN(n);
S2计算接收信号自协方差矩阵R:
R=E{YNYN Η};
S3对R进行特征值分解:R=UΛUΗ,然后将R的特征值从大到小排列;其中,U为NNr×NNr的酉矩阵,Λ是R的特征值构成的对角矩阵;
S4构造检验统计量:
lk表示步骤S3中经过排列的R的第k个特征值;
S5求解判决门限
其中,表示第一类或第二类Tracy-Widom分布的逆累积分布函数,Pfa表示预先设定的虚警概率;参数和是Nr和N的函数;β为的下标,当β=1时,表示第一类Tracy-Widom分布的逆累积分布函数,当β=2时,表示第二类Tracy-Widom分布的逆累积分布函数;
S6确定检测统计量属性:
如果检测统计量Uk大于判决门限γk,则将检测统计量Uk判定为信号统计量;
如果检测统计量Uk小于或等于判决门限γk,则将对应检测统计量Uk判定为噪声统计量;
统计得到信号统计量的总数K;
S7确定SIMO系统信道阶数估计值:
2.根据权利要求1所述的SIMO-FIR信道阶数估计方法,其特征在于,步骤S5中,参数和通过下式确定:
当β=1时代表实数噪声:
当β=2时代表复数噪声:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710126715.5A CN107046512B (zh) | 2017-03-06 | 2017-03-06 | 一种simo-fir信道阶数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710126715.5A CN107046512B (zh) | 2017-03-06 | 2017-03-06 | 一种simo-fir信道阶数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107046512A CN107046512A (zh) | 2017-08-15 |
CN107046512B true CN107046512B (zh) | 2019-12-06 |
Family
ID=59544628
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710126715.5A Active CN107046512B (zh) | 2017-03-06 | 2017-03-06 | 一种simo-fir信道阶数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107046512B (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103532649A (zh) * | 2013-10-22 | 2014-01-22 | 北京邮电大学 | 一种适用于空天信息网络的环境认知技术及设备 |
CN104301281A (zh) * | 2014-10-27 | 2015-01-21 | 西安电子科技大学 | 频率选择性衰落信道下mimo-ofdm系统的发射天线数目估计方法 |
CN104601264B (zh) * | 2015-02-27 | 2018-06-05 | 吉首大学 | 一种适用于高维有限样本条件的多天线频谱感知方法 |
CN106254002B (zh) * | 2016-09-22 | 2018-12-14 | 哈尔滨工业大学 | 认知网络中加权的基于信号相关特性的频谱检测方法 |
-
2017
- 2017-03-06 CN CN201710126715.5A patent/CN107046512B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107046512A (zh) | 2017-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102404063B (zh) | 一种基于过采样的glrt检测方法 | |
US20050254589A1 (en) | Channel estimation device, channel estimation method, and wireless receiver | |
US10972141B2 (en) | Method for estimating arrival time based on noise cancellation | |
CN103795479A (zh) | 一种基于特征值的协作频谱感知方法 | |
CN107276933B (zh) | 用于上行多用户mimo系统中基于二阶统计量的信道估计方法 | |
CN106713190B (zh) | 基于随机矩阵理论和特征阈值估计的mimo发射天线数目盲估计算法 | |
CN102412918A (zh) | 一种基于过采样的空时相关的glrt检测方法 | |
CN100571098C (zh) | 通信系统中低复杂度的极大似然检测方法及装置 | |
CN104780006A (zh) | 基于最小错误概率准则的频谱检测器软融合方法 | |
EP2596605B1 (en) | Method and system for wireless communication channel estimation | |
KR20090003006A (ko) | 다중 입출력 무선통신 시스템에서 다변 다항식을 이용한신호검출 장치 및 방법 | |
US9094080B1 (en) | Method and apparatus for estimating statistics in wireless systems | |
Kortun et al. | Exact performance analysis of blindly combined energy detection for spectrum sensing | |
Ichikawa et al. | Radio environment map construction using Hidden Markov Model in multiple primary user environment | |
CN107046512B (zh) | 一种simo-fir信道阶数估计方法 | |
CN107465472B (zh) | 一种基于路径合成的多径时延估计方法 | |
CN110190908B (zh) | 一种最小化isi信道中非相干大规模simo系统误码率的星座图设计方法 | |
EP1585245B1 (en) | Estimation method for noise space correlated characteristics in time slot cdma system | |
Tseng | Optimum diversity combining and equalization over interference-limited cellular radio channel | |
CN110138478B (zh) | 一种用于非圆信号的多天线频谱感知方法 | |
CN112910518A (zh) | 无人机通信中非高斯噪声下mimo系统发射天线数估计方法 | |
CN109525521B (zh) | 一种信道估计方法 | |
Kaplan et al. | A message transmission scheme for linear time-varying multipath channels | |
CN114422054B (zh) | 一种基于角度信息的无线信道准平稳区间计算装置及方法 | |
Beygi et al. | Time-varying narrowband channel estimation: Exploiting low-rank and sparsity via bilinear representation |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |