CN107884751B - 一种利用单通道接收信号作信源数目估计的方法 - Google Patents
一种利用单通道接收信号作信源数目估计的方法 Download PDFInfo
- Publication number
- CN107884751B CN107884751B CN201710904039.XA CN201710904039A CN107884751B CN 107884751 B CN107884751 B CN 107884751B CN 201710904039 A CN201710904039 A CN 201710904039A CN 107884751 B CN107884751 B CN 107884751B
- Authority
- CN
- China
- Prior art keywords
- signal
- channel
- formula
- characteristic value
- subsample
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radio Transmission System (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种利用单通道接收信号作信源数目估计的方法,该方法的步骤包括:首先,单通道接收信号通过EMD分解扩维获得虚拟阵列信号,对虚拟阵列作信号多次不重复采样,得到多个子样本集,再对每个子样本集构造自相关系数矩阵求其特征值,利用特征值移动增量准则或特征值移动标准差准则检测该子样本集的信源数目。统计各个子样本集检测信源数目,以出现频率最高的信源数目为最终的信源数目。仿真结果表明,采用本发明提出的方法,在信噪比为‑2dB的白噪声环境下,利用单通道接收的信号能准确的检测出3个远场窄带信号源。本发明应用于阵列信号处理技术领域中低信噪比、欠定环境下的信源数目估计技术。
Description
技术领域
本发明涉及一种利用单通道接收信号作信源数目估计的方法,属于信源数目估计的技术领域。
背景技术
在空间谱估计技术中,例如像多重信号分类(MUSIC)和选择子空间不变(ESPRIT)等一些高分辨率测向算法在雷达,声呐和通信等领域中已经得到广泛的应用。然而,这些算法都是建立在准确得知信源数目的前提下实现,如果信源数目估计不准确,会导致这些测向算法性能受到极大影响甚至失效,因此解决信源数目估计问题的意义重大。
目前,比较具有代表性的信源数目估计方法有基于Akaike信息论(AIC)准则、最小描述长度(MDL) 准则、贝叶斯(BIC)准则及其改进法。这类算法是结合特征值分解、最大似然函数以及惩罚函数一起判别信源数目。然而以上信源数目估计算法及其改进算法的共同点是在包含多通道信号的阵列信号模型上进行,然而当只有单个通道的接收信号可用时,这些方法就显得无能为力。
相比多通道,单个通道不仅能避免不同通道信号之间的相互干扰的问题,而且系统结构简单,造价低廉。由于单通道接收信号模型在水声信号、生物医学信号、机械振动信号和一些单通道盲源分离邻域中的广泛应用,使得利用单通道接收信号准确估计信源数目成为后续信号处理的先决条件。利用单通道接收信号作信源数目估计,通常的做法是首先将单通道接收信号变成虚拟的阵列信号,即信号扩维;再用常规的针对阵列信号的信源数目估计算法完成对信源数目的估计。文献“H.Shao,X.H.Shi,and L.Li,“Power signalseparation in milling process based on wavelet transform and independentcomponent analysis,”International Journal of Machine Tools&Manufacture,vol.51,no.9,pp.701-710,Jan 2011.”利用小波变换对单通道信号分解得到伪多路信号,之后用快速独立成分分析(fast-ICA)分离伪多路信号得到信源数。然而不足之处是小波变换不满足自适应,需要人为选择小波基,且合理的小波基难以选择;文献“Z.Dong,J.P.Hu,B.L. Du,and Y.Z.He,“Improvement of Source Number Estimation Method for SingleChannel Signal,”Plos One,vol. 11,no.10,pp.1-12,Oct 2016.”对单通道信号作延时处理得到虚拟多路信号,然后用信息论准则判别信源数目。但缺点是需要大量快拍数完成延时扩维处理,且延时拍数选择同样存在人为因素;文献“Z.Zou,J.P.Hu, Z.P.Huang,andY.M.Zhang,“Estimation of Source Number for Single-Channel Received Signal,”IEEE International Conference on Electronic,Information and ComputerEngineering(ICEICE),vol.44,pp.1-4,Apr 2016.”利用EMD对单通道接收信号分解,虚拟扩展维数,然后用对角加载后的AIC、MDL准则估计信源数目,但是这种方法仅仅适合于高信噪比环境。由此可见,利用单通道接收信号作信源数目估计时,信号虚拟扩维的自适应性和针对虚拟阵列信号处理的有效性(虚拟阵列信号毕竟有别于真实阵列信号)成为该技术的关键问题。
发明内容
为了克服上述不足,本发明以“提升低信噪比下利用单通道接收的信号作信源数目估计性能”为依据,提供了一种利用单通道接收信号作信源数目估计的方法,该方法的特点在于结合了经验模式分解 (Empirical Mode Decomposition,EMD)、信号的重采样技术及构建自相关系数矩阵实现了利用单通道接收信号作信源数目估计。
一种利用单通道接收信号作信源数目估计的方法,其特征在于,具体步骤如下:
1、获取单通道接收信号;
2、对单通道接收信号作EMD分解,获得虚拟多通道阵列信号;
3、对虚拟多通道阵列信号作多次不重复采样,获得多个子样本集;
4、针对每个子样本集构造各自的自相关系数矩阵,并对各自相关系数矩阵作特征值分解得到特征值;
5、针对每个子样本集所对应的特征值,利用特征值移动增量准则或特征值移动标准差准则检测该子样本集的信源数目;
6、统计各个子样本集检测的信源数目,以出现频率最高的信源数目为最终的信源数目
与现有技术相比,本发明的有效效果在于:
EMD是一种自适应的分解方法,具有能将非平稳、非线性信号分解成多个具有某种物理意义分量的特点。利用EMD对单通道接收信号分解成多路信号,既能虚拟扩展维数,也可以尽可能地深度提取原始数据的信息;重采样技术则是通过每次从总样本集中不重复采样得到多个子样本集数据,利用这些子样本集计算更精确的估计值,能表现出更强的鲁棒性且减少了样本估计偏差;而利用各个子样本集构造自相关系数矩阵而不是构造传统的自协方差矩阵主要是为了保证:当作特征值分解时,信号特征值和噪声特征值之间区别更大,最后利用特征值移动增量准则或特征值移动标准差准则检测该子样本集的信源数目时,计算复杂度相比于AIC、MDL准则更低。
附图说明
图1是本发明的算法流程图;
图2是含3个信源的单通道接收信号EMD分解图;
图3是本发明方法与其他方法的信源数目估计的性能比较图;
图4不同信噪比下本发明方法最多能估计的信号源数目;
具体实施方式
为了使本发明的技术手段、创作特征、工作流程、使用方法达成目的与功效易于明白了解,下面结合附图和具体实施例,进一步阐述本发明具体实施方式,本领域技术人员可以在不违背本发明内涵的情况下做类似的推广,因此本发明不受下面公开的具体实施例的限制。
为实现上述的发明目的,本发明采用下述的方案技术,如图1所示:
1、获取单通道接收信号x(t);
假设远场有p个互相独立窄带信源S(t)=[s1(t),s2(t),...,sp(t)]T在t时刻,分别以入射角 {θi,i=1,2,...,p}入射到1根接收天线上,则单通道接收信号为:
x(t)=AS(t)+N(t) (1)
式中A=[a(θ1),a(θ2),...,a(θp)]为1*p维阵列导向矩阵,N(t)为高斯白噪声。
2、对单通道接收信号x(t)作EMD分解,获得虚拟多通道阵列信号Y(t);
将单通道接收信号x(t)经过EMD分解,可表示为一系列的固有模态函数(Intrinsic Mode Functions, IMF)和残余分量之和,为。如式(2)所示
式中ci(t)为第i个IMF分量,r(t)为残余分量。由此构造包含n+1个虚拟阵元的多通道阵列信号 y(t)=[y1(t),y2(t),...,yn+1(t)]T,其中,,y2(t)=c1(t),...,yn(t)=cn-1(t), yn+1(t)=r(t)
3、对虚拟多通道阵列信号Y(t)作多次不重复采样,获得多个子样本集Yr Z;
假设对单通道信号x(t)包含有L个快拍数,则虚拟多通道阵列信号也包含L个快拍数,即
Y=[yi(1),yi(2),...yi(L)] (3)
式中i=1,2,...,n+1;单次重采样的方法如下:从Y中随机选取Lr列,其中Lr=r×L,r为重采样比例,且0.5<r<1;并保证任意两次重采样后的子样本集不完全相同。设第z次重采样获得的子样本集为Yr z,如式(4)所示
a)利用特征值移动增量准则检测信源数目
b)利用特征值移动增量准则检测信源数目
为证明本发明方法的优势,我们在基于MATLAB2015a环境下做了一些仿真实验:
实验一、本发明的方法与其他方法在信号源数目估计性能上的比较
为了简单起见,且不失一般性,假设空间中有3个独立的等功率高斯信号源,分别以20°,30°,45°的入射角入射到单个天线上,空间噪声为白噪声,快拍数为L=500,重采样采样比例r=0.8,重采样次数Z=20,信噪比SNR从-15dB开始以步长1dB变化到15dB。将本发明的方法(EMD+重采样+自相关系数矩阵+特征值移动增量准则(EJAMIC)或EMD+重采样+自相关系数矩阵+特征值移动标准差准则(EJAMSTD)) 与文献“Z.Zou,J.P.Hu,Z.P.Huang,andY.M.Zhang,“Estimation of Source Number for Single-Channel Received Signal,”IEEE International Conference on Electronic,Information and ComputerEngineering (ICEICE),vol.44,pp.1-4.Apr 2016.”中的方法(EMD+自协方差矩阵+对角加载+角加载+AIC(EADLAIC) 和EMD+自协方差矩阵+对角加载+MDL(EADLMDL))和文献“Z.Dong,J.P.Hu,B.L.Du,and Y.Z.He, “Improvement of Source Number EstimationMethod for Single Channel Signal”.Plos One,vol.11,no.10,pp. 1-12.Oct 2016.”中的方法(延时处理+自协方差矩阵+MDL(DAMDL)和延时处理+自协方差矩阵+对角加载+MDL(DADLMDL))进行比较。以上每个实验都重复进行Monte-Carlo实验100次,正确检测概率定义为
其中F为试验总次数,Fk为正确检测出信源数的次数。
对单通道的EMD分解结果如图2所示,本发明的方法与其他方法的性能比较结果如图3所示。由图 3可知,本发明的方法(EJAMIC和EJAMSTD)在SNR=-2dB时仍然有90%以上的正确检测概率,而其他参与比较的方法中,性能最好的方法是EADLAIC和EADLMDL,这两种方法的正确检测率只有在 SNR≥8dB时才能达到90%以上,而另外的两种方法DAMDL和DADLMDL,由于依赖数量较大的节拍数L,所以几乎是失效的。
实验二、本发明的方法最多能估计的信号源数目
为了简单起见,且不失一般性,假设空间中有8个独立的等功率高斯信号源,分别以10°,20°,30°, 40°,50°,60°,70°和80°的入射角入射到单个天线上,空间噪声为白噪声,快拍数为L=500,重采样采样比例r=0.8,重采样次数Z=20,每个高斯信号源都有开启和关闭两种状态,随机选择P个信号源开启, P从1逐次增大到8。当每次P个信号源同时以各自的入射角入射到单个天线时,则分别计算 SNR=-5dB,SNR=0dB,SNR=5dB和SNR=10dB时的信号源估计数目,而每次实验都重复进行Monte-Carlo实验100次,实验结果如图4所示。
由图4可知,当SNR=-5dB时,由于信噪比太低,本发明的方法无法进行欠定情况下的信号源数目估计,本发明方法失效;当SNR=0dB时,本发明的方法能够准确的估计出3个信号源;而在 SNR=5dB和SNR=10dB两种情况时,本发明的方法均能够准确的估计出4个信号源。
Claims (2)
1.一种利用单通道接收信号作信源数目估计的方法,其特征在于,包括以下几个步骤:
步骤一:获取单通道接收信号x(t);
步骤二:对单通道接收信号x(t)作EMD分解,获得虚拟多通道阵列信号Y(t);
获取单通道接收信号x(t):
假设远场有p个互相独立窄带信源S(t)=[s1(t),s2(t),...,sp(t)]T在t时刻,分别以入射角{θi,i=1,2,...,p}入射到1根接收天线上,则单通道接收信号为:
x(t)=AS(t)+N(t) (1)
式中A=[a(θ1),a(θ2),...,a(θp)]为1*p维阵列导向矩阵,N(t)为高斯白噪声;
对单通道接收信号x(t)作EMD分解,获得虚拟多通道阵列信号Y(t):
将单通道接收信号x(t)经过EMD分解,可表示为一系列的固有模态函数(IntrinsicMode Functions,IMF)和残余分量之和,为如式(2)所示
式中ci(t)为第i个IMF分量,r(t)为残余分量;由此构造包含n+1个虚拟阵元的多通道阵列信号y(t)=[y1(t),y2(t),...,yn+1(t)]T,其中,y1(t)=x(t),y2(t)=c1(t),…,yn(t)=cn-1(t),yn+1(t)=r(t);
假设对单通道信号x(t)包含有L个快拍数,则虚拟多通道阵列信号也包含L个快拍数,即
Y=[yi(1),yi(2),...yi(L)] (3)
式中Y∈C(n+1)×L,i=1,2,...,n+1;单次重采样的方法如下:从Y中随机选取Lr列,其中Lr=r×L,r为重采样比例,且0.5<r<1;并保证任意两次重采样后的子样本集不完全相同,设第z次重采样获得的子样本集为如式(4)所示
a)利用特征值移动增量准则检测信源数目
b)利用特征值移动增量准则检测信源数目
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710904039.XA CN107884751B (zh) | 2017-09-29 | 2017-09-29 | 一种利用单通道接收信号作信源数目估计的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710904039.XA CN107884751B (zh) | 2017-09-29 | 2017-09-29 | 一种利用单通道接收信号作信源数目估计的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107884751A CN107884751A (zh) | 2018-04-06 |
CN107884751B true CN107884751B (zh) | 2021-03-12 |
Family
ID=61780986
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710904039.XA Active CN107884751B (zh) | 2017-09-29 | 2017-09-29 | 一种利用单通道接收信号作信源数目估计的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107884751B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108647635A (zh) * | 2018-05-09 | 2018-10-12 | 黑龙江大学 | 基于均分循环补足接收模型的单通道图像盲提取方法 |
CN109406084B (zh) * | 2018-09-13 | 2020-11-06 | 温州大学 | 一种对复合材料产品实现全频域快速冲击监测的方法 |
CN111596285B (zh) * | 2019-11-21 | 2023-08-25 | 中国人民解放军63892部队 | 基于特征值对角加载和构造二阶统计量的信源数估计方法 |
CN111723686B (zh) * | 2020-06-01 | 2023-09-22 | 中国电子科技集团公司第三十六研究所 | 一种信号个数估计方法、装置、电子设备和存储介质 |
CN112133321B (zh) * | 2020-09-23 | 2021-05-14 | 青岛科技大学 | 一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法 |
CN113158741B (zh) * | 2021-01-29 | 2023-04-11 | 中国人民解放军63892部队 | 一种基于特征值对角加载的信源数估计方法 |
CN112861066B (zh) * | 2021-02-15 | 2022-05-17 | 青岛科技大学 | 基于机器学习和fft的盲源分离信源数目并行估计方法 |
CN113341369A (zh) * | 2021-05-27 | 2021-09-03 | 国家无线电监测中心 | 快速单通道空间谱测向方法 |
CN114520755A (zh) * | 2022-01-12 | 2022-05-20 | 中国人民解放军63892部队 | 一种改进式特征值对角加载信源数估计方法 |
CN117331027B (zh) * | 2023-09-27 | 2024-06-04 | 青岛哈尔滨工程大学创新发展中心 | 一种基于子空间匹配度量的声源数估计方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101272168A (zh) * | 2007-03-23 | 2008-09-24 | 中国科学院声学研究所 | 一种信源数估计方法及其波达方向估计方法 |
CN105072062A (zh) * | 2015-08-26 | 2015-11-18 | 天津大学 | 一种基于重采样的信源数目估计方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288285B (zh) * | 2011-05-24 | 2012-11-28 | 南京航空航天大学 | 一种单通道振动信号的盲源分离方法 |
US9549253B2 (en) * | 2012-09-26 | 2017-01-17 | Foundation for Research and Technology—Hellas (FORTH) Institute of Computer Science (ICS) | Sound source localization and isolation apparatuses, methods and systems |
CN104935292B (zh) * | 2014-03-17 | 2017-11-17 | 西南科技大学 | 一种基于源数估计的表面肌电信号自适应滤波方法 |
CN106371070B (zh) * | 2016-08-30 | 2018-09-28 | 电子信息系统复杂电磁环境效应国家重点实验室 | 一种改进的基于小波分析的欠定盲源分离源数估计方法 |
-
2017
- 2017-09-29 CN CN201710904039.XA patent/CN107884751B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101272168A (zh) * | 2007-03-23 | 2008-09-24 | 中国科学院声学研究所 | 一种信源数估计方法及其波达方向估计方法 |
CN105072062A (zh) * | 2015-08-26 | 2015-11-18 | 天津大学 | 一种基于重采样的信源数目估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107884751A (zh) | 2018-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107884751B (zh) | 一种利用单通道接收信号作信源数目估计的方法 | |
Shi et al. | On parameter identifiability of diversity-smoothing-based MIMO radar | |
CN109444820B (zh) | 杂波和干扰共存时多通道雷达先干扰抑制后目标检测方法 | |
CN111239678B (zh) | 一种基于l型阵列的二维doa估计方法 | |
CN110596687B (zh) | 一种基于黎曼流形的单基地mimo雷达目标检测方法 | |
CN106291449B (zh) | 对称稳定分布噪声下波达方向角估计方法 | |
CN106405543B (zh) | 一种认知型盲源分离辐射源提取方法及其评价方法 | |
CN111007488A (zh) | 基于盖尔圆变换和修正Rao得分检验的信源数估计方法 | |
US8799345B1 (en) | Low order multiple signal classification (MUSIC) method for high spectral resolution signal detection | |
CN113253235A (zh) | 一种严重非均匀环境中的自适应信号检测方法与系统 | |
CN112034416A (zh) | 一种自动探测测向无人机的方法 | |
CN111693937A (zh) | 一种基于稀疏重构的无需网格化的近场信号源定位方法 | |
CN112327292B (zh) | 一种二维稀疏阵列doa估计方法 | |
Santamaria et al. | Subspace averaging for source enumeration in large arrays | |
Alexiev | Implementation of hough transform as track detector | |
CN108594165B (zh) | 一种基于期望最大化算法的窄带信号波达方向估计方法 | |
Volkov et al. | Blind Source Separation Model Research Based on Complete Ensamble Empirical Mode Decomposition with Assistant Noise | |
Zhou-hua | The fourth order cumulants based modified MUSIC algorithm for DOA in colored noise | |
CN110888106B (zh) | 一种角度与频率联合估计的增广doa矩阵方法 | |
CN109861770B (zh) | 一种基于波束形成输出功率组合的宽带信号检测方法 | |
CN110208733B (zh) | 基于四阶累量的非圆信号阵列波达方向角估计方法 | |
Tzagkarakis et al. | Covariation-based subspace-augmented MUSIC for joint sparse support recovery in impulsive environments | |
CN112363106B (zh) | 基于量子粒子群的信号子空间波达方向检测方法及系统 | |
CN113514808B (zh) | 一种用于判定小型无人机目标个数的智能辨识方法 | |
CN113093097B (zh) | 一种使用互质阵列的概率假设密度doa跟踪的方法 |
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 |