CN105933259A - 一种基于压缩感知重建的循环相关熵谱的载频估计方法 - Google Patents
一种基于压缩感知重建的循环相关熵谱的载频估计方法 Download PDFInfo
- Publication number
- CN105933259A CN105933259A CN201610251321.8A CN201610251321A CN105933259A CN 105933259 A CN105933259 A CN 105933259A CN 201610251321 A CN201610251321 A CN 201610251321A CN 105933259 A CN105933259 A CN 105933259A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- correlation entropy
- circular correlation
- signal
- carrier frequency
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 29
- 125000004122 cyclic group Chemical group 0.000 title abstract 8
- 230000006835 compression Effects 0.000 claims abstract description 7
- 238000007906 compression Methods 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 abstract description 5
- 238000002474 experimental method Methods 0.000 abstract description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0024—Carrier regulation at the receiver end
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Image Analysis (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于循环平稳信号处理和压缩感知技术领域，利用压缩感知重建后的循环相关熵谱投影进行幅度调制信号载波频率估计的方法。本方法首先计算待估计信号的循环相关熵谱；然后，通过压缩感知方法对循环相关熵谱进行采样和压缩；其次，对根据压缩后的数据进行循环相关熵谱的重建；最后，根据调制方式和载波频率与循环频率之间的关系，通过重建后的循环相关熵谱投影峰值所对应的估计信号的载波频率。实验证明本发明算法性能良好，能够在非高斯噪声和同频带干扰并存的条件下对信号进行载频估计，并且具有数据量小，便于存储和传输的优点。
Description
技术领域
本发明属于循环平稳信号处理和压缩感知技术领域，涉及到调制信号载频估计方法，特别涉及一种利用压缩感知重建后的循环相关熵谱投影进行幅度调制信号载波频率估计的方法。
背景技术
传统的载波频率估计方法主要是基于平稳信号的频率谱或者功率谱来实现的，但是当存在同频带干扰时，这类方法就无法准确地实现载频估计的功能。为了解决该问题，研究工作者又提出了基于循环相关谱的载频估计方法，部分解决了同频带干扰条件下的载频估计问题，但是该类方法同样具有缺陷，即其严重依赖于高斯噪声的假设，当噪声模型不服从高斯分布时，此类方法的性能就会严重退化。之后，针对脉冲噪声和同频带干扰并存条件下的载频估计问题，研究工作者提又出了基于循环相关熵谱的方法，依据循环相关熵谱的特点，在载波频率不断升高的情况下，循环相关熵谱的数据量也会不断增加，导致存储和传输成本变大，导致在这种情况下，该类方法难以实际应用。
因此，本发明提出了一种基于压缩感知重建的循环相关熵谱的载频估计方法，为解决非高斯噪声与同频带干扰并存条件下，非平稳信号的载频估计提出了一种便于存储和传输的方法。
发明内容
针对现有技术的不足,本发明提供一种基于压缩感知重建的循环相关熵谱的载频估计方法,该方法具有抑制非高斯噪声和同频带干扰的特点，同时具有数据量小、便于存储和传输的优点。
为了达到上述目的,本发明的技术方案为:
一种基于压缩感知重建的循环相关熵谱的载频估计方法,包括以下步骤：
第一步,计算待估计信号的循环相关熵谱
(1)计算信号的循环相关熵，它是相关熵的傅里叶级数的系数。
(2)由循环相关熵计算信号的循环相关熵谱，它是循环相关熵的傅里叶变换。
第二步，通过压缩感知方法对循环相关熵谱进行采样
(1)将循环相关熵谱进行离散化。
(2)利用归一化的随机高斯矩阵对循环相关熵谱进行采样，得到观测向量。
第三步，根据压缩感知方法对循环相关熵谱进行重建；
(1)利用观测向量和对应的随机矩阵重建循环相关熵谱向量。
(2)利用重建的循环相关熵谱向量恢复循环相关熵谱。
第四步，计算重建后的循环相关熵谱在循环频率域的投影，并估计载频。
(1)将恢复的循环相关熵谱投影至循环频率域，获得信号的循环相关熵谱投影；
(2)由循环相关熵谱投影中的峰值，获得对应的循环频率ω；
(3)根据循环频率和载频的公式ω＝±afc，估计信号的载波频率为fc，其中，不同的调制方式常数a取值不同。
本发明首先计算信号的相关熵、循环相关熵和循环相关熵谱；然后利用压缩感知的方法，通过归一化的高斯随机矩阵，采样得到观测向量；其次根据观测向量和对应的随机矩阵重建循环相关熵谱；最后利用恢复的循环相关熵谱投影峰值所对应的循环频率与载频的关系，估计信号的载波频率。实验表明，该方法不但解决了脉冲噪声和同频带干扰并存条件下的非平稳信号的载波频率估计问题，还具有数据量小便于存储和传输的优点。
附图说明
图1是本发明载频估计方法的总体流程图；
图2是本发明以BPSK信号为例的循环相关熵谱；
图3是本发明以BPSK信号为例的截取的四分之一的循环相关熵谱；
图4是本发明以BPSK信号为例只利用极少量数据恢复的四分之一的循环相关熵谱；
图5是本发明以BPSK信号为例所恢复的四分之一循环相关熵谱投影。
具体实施方式
为使本发明实施例的目的、技术方案及其优点更加清楚，下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚完整的描述，整体算法流程图如图1所示：
第一步，计算待估计信号的循环相关熵谱
1.1)按公式Vx(t,τ)＝E[κσ(x(t)-x(t+τ))]计算信号的相关熵，其中，t表示时间自变量，x(t)表示信号，τ表示信号的时延，E表示数学期望操作符，кσ表示高斯核函数，它满足再按照公式计算信号的循环相关熵Vx(ω,τ)，公式为：其中，＜T＞表示被积分的区间长度为T，ω表示循环频率，〈·〉t表示求时间平均；
1.2)根据公式计算信号的循环相关熵谱Sx(ω,f)，循环相关熵谱是循环相关熵的傅里叶变换；如图2所示，是以BPSK信号为例的循环相关熵谱的示意图。
第二步，通过压缩感知方法对循环相关熵谱进行采样
2.1)利用公式Sx(nω,nf)＝Sx(nωΔω,nfΔf)将循环相关熵谱Sx(ω,f)进行离散化,得到图像矩阵Sxy＝[s1,s2,…,sN']，其中,Δω是循环频率采样间隔，Δf是频率采样间隔，(nω,nf)图像中元素的坐标位置；为了有效降低数据量，方便存储和传输，通常情况下，如图3所示，我们只需要截取四分之一的循环相关熵谱进行离散化。
2.2)利用归一化的随机高斯矩阵Φn对循环相关熵谱的图像矩阵Sxy＝[s1,s2,…,sN']的第n列向量sn∈RN×1(n＝1,2,…,N')进行采样，得到观测向量zn，其中，zn＝Φnsn；
第三步，根据压缩感知方法对循环相关熵谱进行重建
3.1)利用观测向量zn和对应的随机矩阵Φn，依据优化算法公式得到重建后的循环相关熵谱向量ε是误差容限，是一个无穷小量；
3.2)利用重建的循环相关熵谱向量恢复循环相关熵谱恢复后的四分之一的循环相关熵谱如图4所示；
第四步，计算恢复的循环相关熵谱在循环频率域的投影，并估计载频
4.1)将恢复的循环相关熵谱投影至循环频率域，如图5所示，获得信号的循环相关熵谱投影；
4.2)由循环相关熵谱投影中的峰值，获得对应的循环频率ω；
4.3)根据循环频率和载频的公式ω＝±2fc，估计信号的载波频率为fc，估计信号的载波频率为fc＝1000Hz。
Claims (1)
1.一种基于压缩感知重建的循环相关熵谱的载频估计方法,其特征在于,包括以下步骤:
第一步，计算待估计信号的循环相关熵谱
1.1)按公式Vx(t,τ)＝E[κσ(x(t)-x(t+τ))]计算信号的相关熵，其中，t表示时间自变量，x(t)表示信号，τ表示信号的时延，E表示数学期望操作符，кσ表示高斯核函数，它满足再按照公式计算信号的循环相关熵Vx(ω,τ)，公式为：其中，＜T＞表示被积分的区间长度为T，ω表示循环频率，〈·〉t表示求时间平均；
1.2)根据公式计算信号的循环相关熵谱Sx(ω,f)，循环相关熵谱是循环相关熵的傅里叶变换；
第二步，通过压缩感知方法对循环相关熵谱进行采样
2.1)利用公式Sx(nω,nf)＝Sx(nωΔω,nfΔf)将循环相关熵谱Sx(ω,f)进行离散化,得到图像矩阵Sxy＝[s1,s2,…,sN']，其中,Δω是循环频率采样间隔，Δf是频率采样间隔，(nω,nf)图像中元素的坐标位置；
2.2)利用归一化的随机高斯矩阵Φn对循环相关熵谱的图像矩阵，的第n列向量sn∈RN×1(n＝1,2,…,N')进行采样，得到观测向量zn，其中，zn＝Φnsn；
第三步，根据压缩感知方法对循环相关熵谱进行重建
3.1)利用观测向量zn和对应的随机矩阵Φn，依据优化算法公式得到重建后的循环相关熵谱向量ε是误差容限，是一个无穷小量；
3.2)利用重建的循环相关熵谱向量恢复循环相关熵谱
第四步，计算恢复的循环相关熵谱在循环频率域的投影，并估计载频
4.1)将恢复的循环相关熵谱投影至循环频率域，获得信号的循环相关熵谱投影；
4.2)由循环相关熵谱投影中的峰值，获得对应的循环频率ω；
4.3)根据循环频率和载频的公式ω＝±afc，估计信号的载波频率为fc，其中，不同的调制方式常数a取值不同。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610251321.8A CN105933259B (zh) | 2016-04-21 | 2016-04-21 | 一种基于压缩感知重建的循环相关熵谱的载频估计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610251321.8A CN105933259B (zh) | 2016-04-21 | 2016-04-21 | 一种基于压缩感知重建的循环相关熵谱的载频估计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105933259A true CN105933259A (zh) | 2016-09-07 |
| CN105933259B CN105933259B (zh) | 2018-10-16 |
Family
ID=56839613
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610251321.8A Expired - Fee Related CN105933259B (zh) | 2016-04-21 | 2016-04-21 | 一种基于压缩感知重建的循环相关熵谱的载频估计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105933259B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107529687A (zh) * | 2017-09-20 | 2018-01-02 | 大连理工大学 | 一种基于循环相关熵的时间延迟估计方法 |
| CN108120953A (zh) * | 2017-12-27 | 2018-06-05 | 大连理工大学 | 一种基于波达方向估计的无线电定位方法 |
| CN108282429A (zh) * | 2017-12-01 | 2018-07-13 | 西安电子科技大学 | 基于眼图和信息熵的载波参数估计方法 |
| CN109668732A (zh) * | 2018-12-11 | 2019-04-23 | 北京交通大学 | 基于循环相关熵的滚动轴承的故障诊断方法 |
| CN110850162A (zh) * | 2019-11-22 | 2020-02-28 | 西南交通大学 | 基于误差相关熵的三相电力系统的频率估计方法 |
| WO2021143126A1 (zh) * | 2020-01-13 | 2021-07-22 | 北京邮电大学 | 一种调制方式确定方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120014467A1 (en) * | 2009-03-06 | 2012-01-19 | Telefonaktiebolaget L M Ericsson (Publ) | Method of and Equipment For Compensating Carrier Frequency Offset in an Orthogonal Frequency Division Multiplexing Wireless Radio Transmission System |
| CN105302940A (zh) * | 2015-09-16 | 2016-02-03 | 大连理工大学 | 一种基于循环相关熵的载频估计方法 |
-
2016
- 2016-04-21 CN CN201610251321.8A patent/CN105933259B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120014467A1 (en) * | 2009-03-06 | 2012-01-19 | Telefonaktiebolaget L M Ericsson (Publ) | Method of and Equipment For Compensating Carrier Frequency Offset in an Orthogonal Frequency Division Multiplexing Wireless Radio Transmission System |
| CN105302940A (zh) * | 2015-09-16 | 2016-02-03 | 大连理工大学 | 一种基于循环相关熵的载频估计方法 |
Non-Patent Citations (2)
| Title |
|---|
| YANG LIU，TIANSHUANG QIU: "A Robust Signal Selective TDOA Estimation Algorithm for Cyclostationary Signals", 《 2012 IEEE 14TH INTERNATIONAL CONFERENCE ON COMMUNICATION TECHNOLOGY》 * |
| 刘锋，郑鹏，张鑫，陶然: "一种改进的循环谱估计算法", 《电路与系统学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107529687A (zh) * | 2017-09-20 | 2018-01-02 | 大连理工大学 | 一种基于循环相关熵的时间延迟估计方法 |
| CN107529687B (zh) * | 2017-09-20 | 2019-05-14 | 大连理工大学 | 一种基于循环相关熵的时间延迟估计方法 |
| CN108282429A (zh) * | 2017-12-01 | 2018-07-13 | 西安电子科技大学 | 基于眼图和信息熵的载波参数估计方法 |
| CN108282429B (zh) * | 2017-12-01 | 2020-11-20 | 西安电子科技大学 | 基于眼图和信息熵的载波参数估计方法 |
| CN108120953A (zh) * | 2017-12-27 | 2018-06-05 | 大连理工大学 | 一种基于波达方向估计的无线电定位方法 |
| CN109668732A (zh) * | 2018-12-11 | 2019-04-23 | 北京交通大学 | 基于循环相关熵的滚动轴承的故障诊断方法 |
| CN110850162A (zh) * | 2019-11-22 | 2020-02-28 | 西南交通大学 | 基于误差相关熵的三相电力系统的频率估计方法 |
| CN110850162B (zh) * | 2019-11-22 | 2020-09-29 | 西南交通大学 | 基于误差相关熵的三相电力系统的频率估计方法 |
| WO2021143126A1 (zh) * | 2020-01-13 | 2021-07-22 | 北京邮电大学 | 一种调制方式确定方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105933259B (zh) | 2018-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105933259A (zh) | 一种基于压缩感知重建的循环相关熵谱的载频估计方法 | |
| Dai et al. | Multipath mitigation via component analysis methods for GPS dynamic deformation monitoring | |
| CN102098074B (zh) | 一种用于直接序列扩频系统的高动态弱信号快速捕获方法 | |
| CN105974376B (zh) | 一种sar射频干扰抑制方法 | |
| CN105302940B (zh) | 一种基于循环相关熵的载频估计方法 | |
| CN103364783B (zh) | 基于单通道sar的运动目标径向速度无模糊估计方法 | |
| CN103454621B (zh) | 基于匹配追踪的宽带雷达目标复回波去噪方法 | |
| CN103400383A (zh) | 基于nsct和压缩投影的sar图像变化检测方法 | |
| CN101901476A (zh) | 基于NSCT域边缘检测和Bishrink模型的SAR图像去噪方法 | |
| CN103344948B (zh) | 一种利用稀疏傅里叶变换计算外辐射源雷达互模糊函数的方法 | |
| CN103235295A (zh) | 基于压缩卡尔曼滤波的小场景雷达目标距离像估计方法 | |
| CN104698459A (zh) | 一种应用于缺失数据的条带sar压缩感知成像方法 | |
| CN103364646A (zh) | 微波暗室天线远场快速测量方法 | |
| CN105891600A (zh) | 四相相移键控信号频谱估计方法 | |
| US8467347B2 (en) | Method and apparatus for processing signals | |
| Elango et al. | Weak signal acquisition enhancement in software GPS receivers–Pre-filtering combined post-correlation detection approach | |
| CN106772373B (zh) | 针对任意地面运动目标的sar成像方法 | |
| Cumming et al. | Adding sensitivity to the MLBF Doppler centroid estimator | |
| CN110492903B (zh) | 线性调频信号扩频因子获取方法及装置、可读存储介质 | |
| CN103605139A (zh) | 适用于gnss接收机的载波频率和相位估计方法及系统 | |
| Chen et al. | Ionograms denoising via curvelet transform | |
| CN105260992A (zh) | 基于鲁棒主成分分解和特征空间重构的交通图像去噪算法 | |
| CN104020452B (zh) | 频域空域极化域参数联合估计方法 | |
| Focsa et al. | Mixed Compressive Sensing Back-Projection for SAR Focusing on Geocoded Grid | |
| Zhao et al. | Improved wiener filter super-resolution algorithm for passive millimeter wave imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181016 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |