CN1819572A - 基于希尔伯特-黄变换的二进制频移键控***解调方法 - Google Patents
基于希尔伯特-黄变换的二进制频移键控***解调方法 Download PDFInfo
- Publication number
- CN1819572A CN1819572A CN 200610024974 CN200610024974A CN1819572A CN 1819572 A CN1819572 A CN 1819572A CN 200610024974 CN200610024974 CN 200610024974 CN 200610024974 A CN200610024974 A CN 200610024974A CN 1819572 A CN1819572 A CN 1819572A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- code element
- received signal
- time
- 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
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
一种基于希尔伯特-黄变换的二进制频移键控***解调方法,用于无线电技术领域。方法如下:1)将任意有限长度的二进制频移键控接收信号波形,进行必要的带通滤波处理;2)处理后的信号波形进行HHT处理,得到该接收信号波形的瞬时频率分布情况;3)在时频分布域内,根据给定的判别准则,确定接收信号的极性,完成解调。本发明立足于调制信号本身变化特征,直接使用HHT方法对调制信号进行处理;在判别码元极性时,根据信号最本质的瞬时频率分布情况给出了判决准则,计算复杂度低,物理实现简单。实验结果反映了本发明对于高比特传输率2FSK信号解调的有效性。
Description
技术领域
本发明涉及一种无线电技术领域的方法,具体涉及一种基于希尔伯特-黄变换(HHT)的二进制频移键控(2FSK)***解调方法。
背景技术
HHT是一种自适应的非平稳信号时频分析方法,它包括经验模式分解(EMD)和希尔伯特(Hilbert)谱分析两部分。该方法基于信号本身的变化特征,在时频域自适应分析信号的瞬时频率分布情况。目前,HHT方法主要应用于自然界的一些非平稳信号(如地震波,潮汐,声纳回波等)分析与处理。Norden E.Huang和Steven R.Long等人发表的“The Empirical Mode Decomposition and theHilbert Spectrum for Nonlinear and Non-stationary Time Series Analysis”(非线性和非平稳时间序列分析中的经验模式分解和希尔伯特谱,Proc.R.Soc.Lond.A(1998)454,903-995),论述了HHT的一般过程。在通信***中,2FSK是一种典型的调制***。目前对于2FSK***的解调,基本模式为非相干解调与相干解调,两种解调方法的不足表现为:前者在解调中由于包络检波过程的非线性特性会将有用信号扰乱成噪声,所以抗噪声性能低;后者虽然在理论上能实现信号与噪声单独解调,抗噪声性能优,但由于在解调过程中必需构造本地同步载波,故解调***复杂。另外,两种方法的解调过程与接收信号本身变动特征无关,有效性降低。
经对现有的技术文献检索发现,孙晖、朱善安在《浙江大学学报》(工学版)2005年第12期第1998页上发表的“基于时延自相关预处理的Hilbert-Huang变换解调”,该文提出对加噪调频信号进行时延自相关预处理后,对自相关函数进行HHT处理,根据幅度谱识别信号的调制频率。该方法在信号检测方面效果明显,但其不足在于:如果通过自相关函数获取幅度谱,则在单位码元宽度内需要对调制信号采样多次(理论上采样次数应趋于无穷),对于通信***中高比特率数据传输,在物理实现上将存在困难;从幅度谱的分布情况虽然能识别信号的调制信息,但文献中没有给出码元极性判别依据。
发明内容
本发明的目的在于克服现有技术中存在的不足和缺陷,提出一种基于希尔伯特-黄变换的二进制频移键控***解调方法。使其立足于信号本身的变化特征,具有自适应性,解调过程容易通过软、硬件实现,并且抗噪声性能佳。
本发明是通过以下技术方案实现的,本发明是借用非平稳信号处理方法分析平稳信号,从接收信号的本身特征出发,通过HHT得到信号的瞬时频率分布情况,在时频域内结合一定的判定准则对接收信号的极性作出判断。具体步骤如下:
1)取任意有限长度的二进制频移键控接收信号波形,进行带通滤波处理。
2)处理后的信号波形进行HHT处理,得到该信号波形的瞬时频率分布情况。
3)在时频分布域内,根据给定的判别准则,确定接收信号的极性,完成解调。
所述步骤2),具体实现如下:
以X(t)表示接收信号序列,长度可以取任意个码元宽度,则信号经过EMD处理后可以表示为:
式中cj(t)表示信号的第j个特征模式函数(IMF)分量,rn(t)表示信号经过n次分解后的残余函数分量,IMF分量满足Hilbert变换的前提条件。对分解结果中的各IMF分量进行Hilbert变换,并基于该变换计算出信号的瞬时频率,此过程可以表示为:
式中ai(t)和φi(t)分别表示第i个IMF分量的时变幅度和相位。经过数据综合,得到接收信号的瞬时频率分布(即f-t分布)。
所述步骤3),具体实现如下:
按照如下方法在时频域对接收信号序列进行极性判别:在时频域每个码元宽度内,估计出信号频率的平均值
f,根据2FSK信号的发送频率(当发送码元“0”时载频为f1,当发送码元“1”时载频为f2),以载频的均值f0为判决门限,如果某码元宽度内的频率平均值小于判决门限,该码元判为“0”,如果某码元宽度内的频率平均值大于判决门限,该码元判为“1”,该过程表示为:
f0=(f1+f2)/2;
当
f<f0时,码元判为“0”;
当
f>f0时,码元判为“1”。
本发明立足于调制信号本身变化特征,直接使用HHT方法对调制信号进行处理;在判别码元极性时,根据信号最本质的瞬时频率分布情况给出了判决准则,计算复杂度低,物理实现简单。实验结果反映了本发明对于高比特传输率2FSK信号解调的有效性。
本发明和现有方法相对照,其效果是积极和明显的。本发明具有如下优点:1)基于信号的时频分布特征进行自适应处理,充分利用了信号的有用成分,抗噪声性能佳。2)在解调过程中不需要构造本地载波,容易通过软硬件实现。3)灵活处理信号同步,时窗宽度选取等方面的问题。
附图说明
图1-a二进制频移键控***非相干解调***模块图。
图1-b二进制频移键控***相干解调***模块图。
图2为基于希尔伯特-黄变换的二进制频移键控解调***模块图。
图3为某接收信号波形的瞬时频率分布图。
图4为采用本方法和采用图1所示的非相干、相干方法的仿真结果比较。
具体实施方式
如图2所示,为更好地理解本发明,以下结合附图和具体实施例对方法的实现作进一步描述,实例参数设置如下:设二进制频移键控发送序列为1010 10101000 0110,每个码元宽度为Ts,当发送码元“0”时载频为f1=6000Hz,当发送码元“1”时载频为f2=9200Hz,信噪比为SNR=3dB,接收信号序列的长度为16Ts。本实例的具体实现过程如下:
1)首先对信号序列进行带通滤波处理。
2)对接收信号序列进行HHT处理,得到如图3所示的信号瞬时频率分布。
3)在时频域内,判决门限频率取f0=7600Hz。在某码元宽度内,估计出频率的平均值
f。判决准则为:在某码元宽度内如果信号的平均频率
f小于判决门限f0,该码元判为“0”;如果信号的平均频率
f大于判决门限f0,该码元判为“1”。根据此判决准则,可以确定接收信号码元序列判为1010 1010 10000110。
2FSK是在通信***中典型的调制***,对于2FSK***的解调,基本模式为图1所示的非相干解调(图1-a)与相干解调(图1-b)。图4给出了采用本方法和采用图1所示的非相干、相干方法得到的误码率Monte Carlo仿真结果,仿真次数大于1000次。从图中可以看到,HHT解调与非相干解调虽然都涉及到对信号波形的包络处理,但由于HHT方法基于信号本身的时频分布特征与分析,所以抗噪声性能优于非相干解调。相干解调方法抗噪声性能优,但由于解调过程中需要构造本地同步载波,实现过程比HHT方法复杂。
Claims (3)
1、一种基于希尔伯特-黄变换的二进制频移键控***解调方法,其特征在于,具体步骤如下:
1)取任意有限长度的二进制频移键控接收信号波形,进行带通滤波处理;
2)处理后的信号波形进行HHT处理,得到该信号波形的瞬时频率分布情况;
3)在时频分布域内,根据给定的判别准则,确定接收信号的极性,完成解调。
2、根据权利要求1所述的基于希尔伯特-黄变换的二进制频移键控***解调方法,其特征是,所述步骤2),具体实现如下:
以X(t)表示接收信号序列,长度取任意个码元长度,则信号经过EMD处理后表示为:
式中cj(t)表示信号的第j个特征模式函数IMF分量,rn(t)表示信号经过n次分解后的残余函数分量,IMF分量满足Hilbert变换的前提条件,对分解结果中的各IMF分量进行Hilbert变换,并基于该变换计算出信号的瞬时频率,此过程表示为:
式中αi(t)和φi(t)分别表示第i个IMF分量的时变幅度和相位,经过数据综合,得到接收信号的瞬时频率分布即f-t分布。
3、根据权利要求1所述的基于希尔伯特-黄变换的二进制频移键控***解调方法,其特征是,所述步骤3),具体实现如下:
按照如下方法在时频域对接收信号序列进行极性判别:在时频域每个码元宽度内,估计出信号频率的平均值
f,根据2FSK信号的发送频率,当发送码元“0”时载频为f1,当发送码元“1”时载频为f2,以载频的均值f0为判决门限,如果某码元宽度内的频率平均值小于判决门限,该码元判为“0”,如果某码元宽度内的频率平均值大于判决门限,该码元判为“1”,该过程表示为:
f0=(f1+f2)/2;
当
f<f0时,码元判为“0”;
当
f>f0时,码元判为“1”。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100249749A CN100496033C (zh) | 2006-03-23 | 2006-03-23 | 基于希尔伯特-黄变换的二进制频移键控***解调方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100249749A CN100496033C (zh) | 2006-03-23 | 2006-03-23 | 基于希尔伯特-黄变换的二进制频移键控***解调方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1819572A true CN1819572A (zh) | 2006-08-16 |
CN100496033C CN100496033C (zh) | 2009-06-03 |
Family
ID=36919260
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100249749A Expired - Fee Related CN100496033C (zh) | 2006-03-23 | 2006-03-23 | 基于希尔伯特-黄变换的二进制频移键控***解调方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100496033C (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102086784A (zh) * | 2010-12-16 | 2011-06-08 | 浙江大学 | 大型汽轮发电机组分布式远程振动监测与故障诊断*** |
CN102201824A (zh) * | 2010-03-23 | 2011-09-28 | 上海华虹集成电路有限责任公司 | 一种联合数字下变频和频偏校正电路及其实现方法 |
CN102346064A (zh) * | 2011-09-08 | 2012-02-08 | 浙江工业大学 | 装载机称重装置 |
CN106464624A (zh) * | 2016-07-12 | 2017-02-22 | 深圳市汇顶科技股份有限公司 | 一种应用于封闭通信***的信号解调装置及方法 |
CN107465639A (zh) * | 2017-08-23 | 2017-12-12 | 国家电网公司 | 一种基于短时离散傅立叶变换的多路延时同步判决解调方法 |
CN108614926A (zh) * | 2018-04-12 | 2018-10-02 | 西安交通大学 | 一种基于流形学习与希尔伯特-黄变换相结合的结构模态参数辨识方法 |
CN109359633A (zh) * | 2018-12-10 | 2019-02-19 | 西北工业大学 | 基于希尔伯特-黄变换和小波脊线的信号联合分类方法 |
CN110138699A (zh) * | 2019-06-14 | 2019-08-16 | 哈尔滨工程大学 | 一种基于复数域滤波器的基带2fsk信号非相干解调方法 |
CN110622043A (zh) * | 2017-05-12 | 2019-12-27 | 株式会社普利司通 | 路面状态判别方法和路面状态判别装置 |
CN111487318A (zh) * | 2020-05-29 | 2020-08-04 | 福建农林大学 | 一种时变结构瞬时频率提取方法 |
CN111935044B (zh) * | 2020-08-20 | 2021-03-09 | 金陵科技学院 | 基于emd分解的psk及qam类信号调制识别方法 |
-
2006
- 2006-03-23 CN CNB2006100249749A patent/CN100496033C/zh not_active Expired - Fee Related
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102201824A (zh) * | 2010-03-23 | 2011-09-28 | 上海华虹集成电路有限责任公司 | 一种联合数字下变频和频偏校正电路及其实现方法 |
CN102201824B (zh) * | 2010-03-23 | 2013-12-04 | 上海华虹集成电路有限责任公司 | 一种联合数字下变频和频偏校正电路的实现方法 |
CN102086784A (zh) * | 2010-12-16 | 2011-06-08 | 浙江大学 | 大型汽轮发电机组分布式远程振动监测与故障诊断*** |
CN102346064A (zh) * | 2011-09-08 | 2012-02-08 | 浙江工业大学 | 装载机称重装置 |
US10291388B2 (en) | 2016-07-12 | 2019-05-14 | Shenzhen GOODIX Technology Co., Ltd. | Signal demodulation apparatus and method in closed communication system |
CN106464624A (zh) * | 2016-07-12 | 2017-02-22 | 深圳市汇顶科技股份有限公司 | 一种应用于封闭通信***的信号解调装置及方法 |
CN110622043B (zh) * | 2017-05-12 | 2021-11-09 | 株式会社普利司通 | 路面状态判别方法和路面状态判别装置 |
CN110622043A (zh) * | 2017-05-12 | 2019-12-27 | 株式会社普利司通 | 路面状态判别方法和路面状态判别装置 |
US11486702B2 (en) | 2017-05-12 | 2022-11-01 | Bridgestone Corporation | Road surface condition determination method and road surface condition determination apparatus |
CN107465639B (zh) * | 2017-08-23 | 2020-05-19 | 国家电网公司 | 一种基于短时离散傅立叶变换的多路延时同步判决解调方法 |
CN107465639A (zh) * | 2017-08-23 | 2017-12-12 | 国家电网公司 | 一种基于短时离散傅立叶变换的多路延时同步判决解调方法 |
CN108614926A (zh) * | 2018-04-12 | 2018-10-02 | 西安交通大学 | 一种基于流形学习与希尔伯特-黄变换相结合的结构模态参数辨识方法 |
CN109359633A (zh) * | 2018-12-10 | 2019-02-19 | 西北工业大学 | 基于希尔伯特-黄变换和小波脊线的信号联合分类方法 |
CN109359633B (zh) * | 2018-12-10 | 2021-12-07 | 西北工业大学 | 基于希尔伯特-黄变换和小波脊线的信号联合分类方法 |
CN110138699A (zh) * | 2019-06-14 | 2019-08-16 | 哈尔滨工程大学 | 一种基于复数域滤波器的基带2fsk信号非相干解调方法 |
CN110138699B (zh) * | 2019-06-14 | 2022-03-18 | 哈尔滨工程大学 | 一种基于复数域滤波器的基带2fsk信号非相干解调方法 |
CN111487318A (zh) * | 2020-05-29 | 2020-08-04 | 福建农林大学 | 一种时变结构瞬时频率提取方法 |
CN111487318B (zh) * | 2020-05-29 | 2023-03-24 | 福建农林大学 | 一种时变结构瞬时频率提取方法 |
CN111935044B (zh) * | 2020-08-20 | 2021-03-09 | 金陵科技学院 | 基于emd分解的psk及qam类信号调制识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100496033C (zh) | 2009-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1819572A (zh) | 基于希尔伯特-黄变换的二进制频移键控***解调方法 | |
Li et al. | A novel noise reduction technique for underwater acoustic signals based on complete ensemble empirical mode decomposition with adaptive noise, minimum mean square variance criterion and least mean square adaptive filter | |
Liu et al. | Property of self-similarity between baseband and modulated signals | |
CN101848177B (zh) | 基于变频的双稳态最优随机共振单频微弱信号检测方法 | |
WO2004075502A3 (en) | Receiving apparatus and method for digital multi-carrier transmission | |
CN101217378B (zh) | 一种基于流量统计的小波分析边界处理方法 | |
EP3748897B1 (en) | Quantum chaotic wave packet digital signal generating method | |
CN101988935A (zh) | 基于数字下变频—希尔伯特黄变换的瞬时频率测量方法 | |
CN106209701A (zh) | Alpha稳定分布噪声环境下MFSK信号符号速率估计方法及装置 | |
CN109359633B (zh) | 基于希尔伯特-黄变换和小波脊线的信号联合分类方法 | |
Jiang et al. | Modulation recognition of non-cooperation underwater acoustic communication signals using principal component analysis | |
CN102546061A (zh) | 基于小波变换的自适应时频空穴检测方法 | |
CN109379310A (zh) | 一种基于Rife-Quinn综合的MPSK信号载频估计方法 | |
Sisi et al. | Carrier frequency and symbol rate estimation based on cyclic spectrum | |
CN107688167B (zh) | 一种多时宽线性调频脉冲压缩信号幅度包络曲线生成方法 | |
CN105162434A (zh) | 基于时频分析的时变滤波参数产生与实现***及方法 | |
CN104502820A (zh) | 基于二维小波变换的变压器局部放电图谱去噪方法 | |
Geng et al. | Quantitative identification of pulse-like ground motions based on intrinsic time-scale decomposition technique | |
CN110113273B (zh) | 一种基于频谱对称的时频重叠信号盲分离方法 | |
CN106842141A (zh) | 一种高阶多次共轭滞后积脉内分析方法 | |
Kennedy et al. | A novel method for topological embedding of time-series data | |
Shan et al. | Nonlinear time-varying spectral analysis: HHT and MODWPT. | |
CN104753842B (zh) | 基于峰值位置判别的信号调制方式识别方法 | |
CN103166887B (zh) | 一种边带调制信号分类的方法及装置 | |
Zhu et al. | Incomplete phase-space method to reveal time delay from scalar time series |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090603 Termination date: 20190323 |