CN104702546B - 一种基于滤波的希尔伯特实时解调实现方法 - Google Patents
一种基于滤波的希尔伯特实时解调实现方法 Download PDFInfo
- Publication number
- CN104702546B CN104702546B CN201310670304.4A CN201310670304A CN104702546B CN 104702546 B CN104702546 B CN 104702546B CN 201310670304 A CN201310670304 A CN 201310670304A CN 104702546 B CN104702546 B CN 104702546B
- Authority
- CN
- China
- Prior art keywords
- output
- filter
- hilbert
- real
- coefficient
- 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
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明涉及一种基于滤波的希尔伯特实时解调实现方法,该方法包括以下步骤:1)利用数字滤波器实现一个带通的希尔伯特变换器和一个带通滤波器,所述的希尔伯特变换器和带通滤波器的响应一致,相位相差90°;2)将上述希尔伯特变换器和带通滤波器的输出作为复数滤波器的两组系数存入复数滤波器;3)根据复数滤波器的输出进行包络解调计算,得到包络信号。与现有技术相比,本发明具有计算简单、计算量小、可实现连续的实时包络等优点。
Description
技术领域
本发明涉及一种希尔伯特变换实现方法,尤其是涉及一种基于滤波的希尔伯特实时解调实现方法。
背景技术
在许多科学领域的信号处理中,都需要提取信号的包络信息,如在齿轮和滚动轴承等机械故障诊断中经常用到包络解调,用振动信号的包络实现诊断是一种行之有效的方法。包络解调是将高频部分的包络解调到低频,便于观察故障现象。目前最常用的信号包络提取方法是采用希尔伯特(Hilbert)变换实现。传统的Hilbert是利用卷积方法实现的,或利用FFT的方法计算,这种方法计算量大,计算时间长,不易实现连续的实时包络。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种计算简单、计算量小、可实现连续的实时包络的基于滤波的希尔伯特实时解调实现方法。
本发明的目的可以通过以下技术方案来实现:
一种基于滤波的希尔伯特实时解调实现方法,该方法包括以下步骤:
1)利用数字滤波器实现一个带通的希尔伯特变换器和一个带通滤波器,所述的希尔伯特变换器和带通滤波器的响应一致,相位相差90°;
2)将上述希尔伯特变换器和带通滤波器的输出作为复数滤波器的两组系数存入复数滤波器;
3)根据复数滤波器的输出进行包络解调计算,得到包络信号。
所述的复数滤波器包括输入寄存器、输出寄存器、输入数据存储器和系数存储器,所述的输入寄存器、输入数据存储器、系数存储器、输出寄存器依次连接,希尔伯特变换器和带通滤波器的输出存储在系数存储器中;原始信号通过输入寄存器分区保存至输入数据存储器中,根据原始信号及系数存储器中的系数获得复数滤波器的输出,并保存至输出寄存器中,所述的复数滤波器的输出F(n)通过以下公式获得:
式中,Re(H(i))为带通滤波器获得的系数,Im(H(i))为希尔伯特变换器获得的系数,Re(x(n-i))为原始信号x(n-i)的实部,Im(x(n-i))为原始信号x(n-i)的虚部,Re(F(n))为输出的实部,Im(F(n))为输出的虚部。
所述的包络信号具体计算公式为:
式中,y(n)为包络信号,x(n)=Re(F(n)),x′(n)=Im(F(n))。
所述的输入寄存器、输出寄存器、输入数据存储器和系数存储器间的数据传输采用DMA方式。
与现有技术相比,本发明具有以下优点:
1、本发明采用硬件复数滤波器来实现这个希尔伯特实时解调,计算量小,计算时间短;
2、输入连续采样AD信号,用DMA方法可以产生连续实时的包络信号,不需要CPU参与,将数据传送到显示屏上可以实时显示包络信号,数据传输快,有利于故障的诊断。
附图说明
图1为本发明复数滤波器的结构示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
一种基于滤波的希尔伯特实时解调实现方法,该方法包括以下步骤:
1)利用数字滤波器实现一个带通的希尔伯特变换器和一个带通滤波器,所述的希尔伯特变换器和带通滤波器的响应一致,相位相差90°,可采用DSP自带的数字滤波器实现上述功能。
2)将上述希尔伯特变换器和带通滤波器的输出作为复数滤波器的两组系数存入复数滤波器。
如图1所示,复数滤波器包括输入寄存器1、输出寄存器4、输入数据存储器2和系数存储器3,所述的输入寄存器1、输入数据存储器2、系数存储器3、输出寄存器4依次连接,希尔伯特变换器和带通滤波器的输出存储在系数存储器中;原始信号通过输入寄存器1分区保存至输入数据存储器2中,根据原始信号及系数存储器3中的系数(B0,B1,…,BN,N为系数个数)获得复数滤波器的输出,并保存至输出寄存器4中,所述的复数滤波器的输出F(n)通过以下公式获得:
式中,Re(H(i))为带通滤波器获得的系数,Im(H(i))为希尔伯特变换器获得的系数,Re(x(n-i))为原始信号x(n-i)的实部,Im(x(n-i))为原始信号x(n-i)的虚部,Re(F(n))为输出的实部,Im(F(n))为输出的虚部。
输入寄存器、输出寄存器、输入数据存储器和系数存储器间的数据传输采用DMA方式,不需要CPU参与。
3)根据复数滤波器的输出进行包络解调计算,得到包络信号,具体计算公式为:
式中,y(n)为包络信号,x(n)=Re(F(n)),x′(n)=Im(F(n))。
振动分析中一般是对低频的包络信号感兴趣,将上述包络数据经过低通滤波器滤出需要的频率,即可得到所需要的包络解调数据。
Hilbert变换可以提供90°的相位变化而不改变频谱分量的幅度,即对信号进行希尔伯特变换就相当于对该信号进行正交移相,构建解析信号,对解析信号求模即可得到包络信号。包络解调是将高频部分的包络解调到低频,便于观察故障现象,可用于齿轮和滚动轴承的故障诊断中。
Claims (2)
1.一种基于滤波的希尔伯特实时解调实现方法,其特征在于,该方法包括以下步骤:
1)利用数字滤波器实现一个带通的希尔伯特变换器和一个带通滤波器,所述的希尔伯特变换器和带通滤波器的响应一致,相位相差90°;
2)将上述希尔伯特变换器和带通滤波器的输出作为复数滤波器的两组系数存入复数滤波器;
3)根据复数滤波器的输出进行包络解调计算,得到包络信号;
所述的复数滤波器包括输入寄存器、输出寄存器、输入数据存储器和系数存储器,所述的输入寄存器、输入数据存储器、系数存储器、输出寄存器依次连接,希尔伯特变换器和带通滤波器的输出存储在系数存储器中;原始信号通过输入寄存器分区保存至输入数据存储器中,根据原始信号及系数存储器中的系数获得复数滤波器的输出,并保存至输出寄存器中,所述的复数滤波器的输出F(n)通过以下公式获得:
式中,Re(H(i))为带通滤波器获得的系数,Im(H(i))为希尔伯特变换器获得的系数,Re(x(n-i))为原始信号x(n-i)的实部,Im(x(n-i))为原始信号x(n-i)的虚部,Re(F(n))为输出的实部,Im(F(n))为输出的虚部;
所述的输入寄存器、输出寄存器、输入数据存储器和系数存储器间的数据传输采用DMA方式。
2.根据权利要求1所述的一种基于滤波的希尔伯特实时解调实现方法,其特征在于,所述的包络信号具体计算公式为:
式中,y(n)为包络信号,x(n)=Re(F(n)),x'(n)=Im(F(n))。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310670304.4A CN104702546B (zh) | 2013-12-10 | 2013-12-10 | 一种基于滤波的希尔伯特实时解调实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310670304.4A CN104702546B (zh) | 2013-12-10 | 2013-12-10 | 一种基于滤波的希尔伯特实时解调实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104702546A CN104702546A (zh) | 2015-06-10 |
CN104702546B true CN104702546B (zh) | 2018-07-20 |
Family
ID=53349323
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310670304.4A Active CN104702546B (zh) | 2013-12-10 | 2013-12-10 | 一种基于滤波的希尔伯特实时解调实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104702546B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107765235A (zh) * | 2016-08-16 | 2018-03-06 | 合肥科盛微电子科技有限公司 | 基于数字滤波、数字包络提取的超声波测距算法 |
CN110138699B (zh) * | 2019-06-14 | 2022-03-18 | 哈尔滨工程大学 | 一种基于复数域滤波器的基带2fsk信号非相干解调方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1459143A (zh) * | 2000-09-19 | 2003-11-26 | 汤姆·瑞雷 | 复值德耳塔西格马锁相环解调器 |
CN101335023A (zh) * | 2007-06-29 | 2008-12-31 | 索尼株式会社 | 跟踪误差信号检测设备和光盘装置 |
CN102967334A (zh) * | 2012-09-26 | 2013-03-13 | 朱作行 | 利用对信号包络线处理测量流体流量的系统及方法 |
CN103116032A (zh) * | 2013-01-18 | 2013-05-22 | 北京金风科创风电设备有限公司 | 用于获取风力发电机组的转速的方法和设备 |
CN103424258A (zh) * | 2013-08-06 | 2013-12-04 | 昆明理工大学 | 一种滚动轴承的故障诊断方法 |
-
2013
- 2013-12-10 CN CN201310670304.4A patent/CN104702546B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1459143A (zh) * | 2000-09-19 | 2003-11-26 | 汤姆·瑞雷 | 复值德耳塔西格马锁相环解调器 |
CN101335023A (zh) * | 2007-06-29 | 2008-12-31 | 索尼株式会社 | 跟踪误差信号检测设备和光盘装置 |
CN102967334A (zh) * | 2012-09-26 | 2013-03-13 | 朱作行 | 利用对信号包络线处理测量流体流量的系统及方法 |
CN103116032A (zh) * | 2013-01-18 | 2013-05-22 | 北京金风科创风电设备有限公司 | 用于获取风力发电机组的转速的方法和设备 |
CN103424258A (zh) * | 2013-08-06 | 2013-12-04 | 昆明理工大学 | 一种滚动轴承的故障诊断方法 |
Non-Patent Citations (1)
Title |
---|
采用复数波数字滤波器实现希尔伯特(Hilbert)变换器;黄大卫;《通信学报》;19841031(第4期);第42页第2段至第45页最后一段 * |
Also Published As
Publication number | Publication date |
---|---|
CN104702546A (zh) | 2015-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Bearing fault diagnosis method based on adaptive maximum cyclostationarity blind deconvolution | |
Qin et al. | Higher-density dyadic wavelet transform and its application | |
Chen et al. | Improved spectral kurtosis with adaptive redundant multiwavelet packet and its applications for rotating machinery fault detection | |
CN107276560B (zh) | 一种fri稀疏采样核函数构建方法及电路 | |
CN104301052A (zh) | 一种基于fpga的无缝采集及实时频谱监测的实现方法 | |
CN101425176A (zh) | 一种基于中值滤波的图像小波去噪方法 | |
CN104702546B (zh) | 一种基于滤波的希尔伯特实时解调实现方法 | |
CN104124976A (zh) | 有限新息率信号结构化亚奈奎斯特率采样方法 | |
CN105306405B (zh) | 一种无源互调信号时延、频率和相位估计装置与方法 | |
CN103823133A (zh) | 一种基于压缩感知的在线电能质量监测系统 | |
Chen et al. | Construction of adaptive redundant multiwavelet packet and its application to compound faults detection of rotating machinery | |
CN110108938B (zh) | 基于gpu数据流驱动的实时电磁频谱分析系统及方法 | |
Qin et al. | Dense framelets with two generators and their application in mechanical fault diagnosis | |
Zhang et al. | Compound fault extraction method via self-adaptively determining the number of decomposition layers of the variational mode decomposition | |
Zhang et al. | Gearbox composite fault diagnosis method based on minimum entropy deconvolution and improved dual-tree complex wavelet transform | |
CN102129073A (zh) | 一种合成孔径雷达方位预处理的两级滤波的方法 | |
Ben et al. | MEMS gyro denoising based on second generation wavelet transform | |
CN202982024U (zh) | 一种用于x射线机的信号跟随滤噪电路 | |
Li et al. | A wideband reconnaissance receiver design based on real-time spectrum analysis technology | |
CN105578083B (zh) | 低噪声ccd信号成像电子学系统 | |
Li et al. | Compensation method for the CIC filter in digital down converter | |
CN104393841A (zh) | 一种数字中频信号正交下变频的实现方法及模块 | |
CN109672489A (zh) | 一种基于互素分析的海上全景频谱感知分析方法及其装置 | |
CN105743505B (zh) | 信号自相关调制压缩采样与恢复方法 | |
CN203502108U (zh) | 一种声波采集电路 |
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 |