CN1731388A - 超声信号的数字信号处理方法 - Google Patents
超声信号的数字信号处理方法 Download PDFInfo
- Publication number
- CN1731388A CN1731388A CN 200410011403 CN200410011403A CN1731388A CN 1731388 A CN1731388 A CN 1731388A CN 200410011403 CN200410011403 CN 200410011403 CN 200410011403 A CN200410011403 A CN 200410011403A CN 1731388 A CN1731388 A CN 1731388A
- Authority
- CN
- China
- Prior art keywords
- signal
- ultrasonic
- time
- digital
- frequency analysis
- 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
Abstract
本发明涉及利用数字信号处理超声信号的方法。将数字超声信号通过提升小波包变换滤波、希尔伯特(Hilbert)变换、正交解调等处理,在合成超声成像的同时进行联合时频分析和特征提取。本发明提出的集成式的数字信号处理构架,可以实时地分析射频信号、解析信号以及正交解调信号。特别地,可以使用联合时频分析技术提取出较传统的快速傅立叶变换更加丰富的有关被测材料完整性与几何特征的重要信息。作为虚拟仪器开发的软件平台,Labwindows/CVI被用于构筑了基于计算机的上层测试与控制组件。这一系统可以实现基于小波的动态滤波、降噪、下采样、联合时频分析,以及专业的测试与控制。本发明可以应用于材料的无损检测与评估、声纳、距离与速度测量以及医学应用等领域。
Description
技术领域:本发明属于超声信号处理技术,涉及利用数字信号处理(DSP)超声信号的方法。
背景技术:图1为超声信号处理的传统模型。由于AD转换器的采样速度较低,所以在AD转换器前端必须加入抗混叠滤波器。通常这一抗混叠滤波器的设计非常烦琐,它的电路在印刷电路板上占据了相当大的面积,而且设计实现后无法改变,无法适应不同类型的超声应用。特征提取与噪声衰减方面采用傅立叶变换。
图2为目前人们采用的一种改进的模型。虽然AD转换器的采样速度较高,但是为了提高信噪比并降低抗混叠滤波器设计的复杂度,使用4倍甚至更高的过采样率,而提高的采样率带来更大的数据量。为了提取有用信息并降低数据量,人们采用复杂可编程逻辑器件(FPGA)来实现数字滤波和下采样。FPGA虽具有在线可编程的特点,但是它不可以实时修改,配置并不方便,而且数字滤波器的设计也仍然无法实现自动化。特征提取与噪声衰减方面采用傅立叶变换,缺点是单纯地进行时域或频域的分析对于缺陷特征的识别往往不够充分。
发明内容:为了解决上述背景技术中抗混叠滤波器设计实现后无法改变,无法适应不同类型超声应用的问题;可编程逻辑器件不能实时修改、配置不方便、设计不能自动化,对缺陷特征的识别不充分等问题,本发明将要提出一种集成式的对超声信号的数字信号处理(DSP)方法。
本发明超声信号的数字处理方法步骤如下:
a.首先将超声信号转换成数字超声信号,数字超声信号通过提升小波包变换滤波后产生超声射频信号;
b.利用希尔伯特(Hilbert)变换将步骤a产生的超声射频信号转换成具有其特征的超声解析信号;
c.将步骤b产生的超声解析信号通过正交解调产生正交解调信号;
d.将步骤c产生的正交解调信号再次通过提升小波包变换,产生压缩的正交解调信号;
e.将步骤d产生的压缩的正交解调信号在进行超声成像的同时进行联合时频分析,得到含有时域和频域的时频信息;
f.将步骤b产生的超声解析信号和步骤e时频信息进行特征提取,则完成超声信号的数字处理。
本发明的优点:
由于采用了提升小波包变换(WPT),使用正交与双正交小波从有用子带中提取超声信号来进行滤波,小波分解的级次和有用的子带根据采样频率、传感器的中心频率和带宽自动确定,就不再需要设计滤波器,解决了抗混叠滤波器设计实现后无法改变、无法适应不同类型超声应用的问题。小波包变换的最大优点是信号特征的提取非常容易,它在特征提取与噪声衰减方面的性能在许多领域中都胜过了传统的快速傅立叶变换。
由于本发明采用DSP构架,将超声信号进行处理,具有高度灵活性,可以实时修改,因而实时性佳,解决了采用可编程逻辑器件不能实时修改、配置不方便、设计不能自动化,对于缺陷特征的识别不充分等问题。特别地,这种架构使用联合时频分析技术,提取被测材料有关超声信号的完整性与几何特征的重要信息较传统的快速傅立叶变换更加丰富。
本发明可以应用于材料的无损检测与评估、声纳、距离与速度测量以及医学应用等领域。
附图说明:
图1是背景技术传统超声信号处理的结构示意图
图2是背景技术常用的超声信号处理的结构示意图
图3时本发明超声信号处理的结构示意图
图4是本发明超声测试系统中数字信号处理结构
具体实施方式:
实施例:如图3图4所示包括时间增益补偿电路、A/D转换器、数字信号处理电路、计算机。
时间增益补偿电路包括:可编程逻辑器件、A转换器、可编程增益放大器。可编程逻辑器件采用CPLD可编程逻辑器件;DA转换器采用MAX5110高速DA转换器;可编程增益放大器采用AD604可编程增益放大器。时间增益补偿电路由CPLD经高速DA转换器控制可编程增益放大器来实现。
A/D转换器:可采用40Mhz/60MHz的A/D转换器,如:MAX1183。
数字信号处理电路:采用TI的DSP开发板:TMS320C6701EVM,结构如图5所示。
计算机分析软件采用美国国家仪器公司的Labwindows/CVI虚拟仪器开发系统。
本发明中:
a.利用提升小波包变换将数字超声信号分解到二进制子频带,根据传感器超声频带,保留所需二进制子频带信号能量,屏蔽或压制其它频带信号能量,而后经小波反变换,综合出滤波后信号,即超声射频信号(1)。
b.以超声射频信号作为实部,超声射频信号经过Hilbert变换后的信号作虚部即构成超声解析信号(2),超声解析信号是复信号,只有正频率成分。
Hilbert变换就是将超声射频信号作90度相移,得到超声射频信号的正交信号。定义如下:
c.超声解析信号(2)的幅度即构成了超声射频信号(1)的复包络,即正交解调信号(3),所以将步骤b产生的超声解析信号(2)取模,即得到超声射频信号(1)的包络,即正交解调信号(3);
d.步骤c产生的正交解调信号(3)再次通过提升小波包变换,提取低频段信号能量即得到压缩的正交解调信号(4);
e.将压缩的正交解调信号(4)的信号幅度以色彩表示,每个扫描信号即构成一条色彩随幅度变换的线条,上千条扫描线即构成了彩色的超声图像。
利用Labwindows/CVI中相关的函数对压缩的正交解调信号(4)进行联合时频分析。联合时频分析将一维的时域信号转变成二维的能量对时间和频率的联合分布。可以将变换结果以3维或称之2.5维的方式表现出来(能力大小以光强或色彩值表现)。联合时频分析对应于一个信号,是信号能量在时间和频率方面的综合表达。它能准确的判定特征信号的出现时间和对应的频率及其强度。联合时频分析展现了信号的频率成分随时间变换的过程。通过时间和频率域的同时分析我们得到了更丰富的信息。
f.利用Labwindows/CVI中相关的函数提取信号的瞬时频率的方法,完成精确的目标定位,将步骤b产生的超声解析信号(2)和步骤e时频信息进行特征提取,则完成超声信号的数字处理。
Claims (1)
1、超声信号的数字处理方法,其特征在于其步骤如下:
a.首先将数字超声信号通过提升小波包变换滤波后产生超声射频信号;
b.利用希尔伯特(Hilbert)变换将步骤a产生的超声射频信号转换成具有其特征的超声解析信号;
c.将步骤b产生的超声解析信号通过正交解调产生正交解调信号;
d.步骤c产生的正交解调信号再次通过提升小波包变换,产生压缩的正交解调信号;
e.步骤d产生的压缩的正交解调信号在进行超声成像的同时进行联合时频分析,得到含有时域和频域的时频信息;
f.将步骤b产生的超声解析信号和步骤e时频信息进行特征提取,则完成超声信号的数字处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100114032A CN100410925C (zh) | 2004-12-30 | 2004-12-30 | 超声信号的数字信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100114032A CN100410925C (zh) | 2004-12-30 | 2004-12-30 | 超声信号的数字信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1731388A true CN1731388A (zh) | 2006-02-08 |
CN100410925C CN100410925C (zh) | 2008-08-13 |
Family
ID=35963742
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004100114032A Expired - Fee Related CN100410925C (zh) | 2004-12-30 | 2004-12-30 | 超声信号的数字信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100410925C (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100417191C (zh) * | 2006-09-01 | 2008-09-03 | 上海大学 | 一种图像组合降噪方法 |
WO2008125048A1 (fr) * | 2007-04-17 | 2008-10-23 | Huawei Technologies Co., Ltd. | Procédé et appareil de traitement de mesures |
CN101312529B (zh) * | 2007-05-24 | 2010-07-21 | 华为技术有限公司 | 生成上下采样滤波器及实现编码的方法、系统和装置 |
CN101473528B (zh) * | 2006-06-21 | 2011-05-11 | Nxp股份有限公司 | 对调制信号进行解调的方法、解调器和接收机 |
CN101190137B (zh) * | 2006-11-27 | 2011-06-22 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声成像系统的实时数字正交解调方法与装置 |
CN102176121A (zh) * | 2011-01-18 | 2011-09-07 | 河海大学 | 数字超声经颅多普勒数字解调和信号处理方法及装置 |
CN101467897B (zh) * | 2007-12-27 | 2013-01-16 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置、超声波图像处理方法 |
CN106488358A (zh) * | 2015-09-09 | 2017-03-08 | 上海其高电子科技有限公司 | 优化声场成像定位方法及系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1077532A1 (en) * | 1999-08-17 | 2001-02-21 | BRITISH TELECOMMUNICATIONS public limited company | Spread Spectrum Signal Generator and Decoder for Single Sideband Transmission |
WO2004048903A1 (en) * | 2002-11-25 | 2004-06-10 | Instromet Ultrasonics B.V. | Ultrasonic signal processing method and applications thereof |
-
2004
- 2004-12-30 CN CNB2004100114032A patent/CN100410925C/zh not_active Expired - Fee Related
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101473528B (zh) * | 2006-06-21 | 2011-05-11 | Nxp股份有限公司 | 对调制信号进行解调的方法、解调器和接收机 |
CN100417191C (zh) * | 2006-09-01 | 2008-09-03 | 上海大学 | 一种图像组合降噪方法 |
CN101190137B (zh) * | 2006-11-27 | 2011-06-22 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声成像系统的实时数字正交解调方法与装置 |
WO2008125048A1 (fr) * | 2007-04-17 | 2008-10-23 | Huawei Technologies Co., Ltd. | Procédé et appareil de traitement de mesures |
CN101291517B (zh) * | 2007-04-17 | 2012-01-04 | 华为技术有限公司 | 一种测量处理方法及装置 |
US8432991B2 (en) | 2007-04-17 | 2013-04-30 | Huawei Technologies Co., Ltd. | Measurement processing method and apparatus |
CN101312529B (zh) * | 2007-05-24 | 2010-07-21 | 华为技术有限公司 | 生成上下采样滤波器及实现编码的方法、系统和装置 |
CN101467897B (zh) * | 2007-12-27 | 2013-01-16 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置、超声波图像处理方法 |
CN102176121A (zh) * | 2011-01-18 | 2011-09-07 | 河海大学 | 数字超声经颅多普勒数字解调和信号处理方法及装置 |
CN106488358A (zh) * | 2015-09-09 | 2017-03-08 | 上海其高电子科技有限公司 | 优化声场成像定位方法及系统 |
CN106488358B (zh) * | 2015-09-09 | 2019-07-19 | 上海其高电子科技有限公司 | 优化声场成像定位方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN100410925C (zh) | 2008-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100498349C (zh) | 中频全数字化频谱仪 | |
CN105445022B (zh) | 一种基于双树复小波变换‑熵特征融合的行星齿轮故障诊断方法 | |
CN1731388A (zh) | 超声信号的数字信号处理方法 | |
CN109580787A (zh) | 用于变压器高压套管引线超声检测的超声回波去噪方法 | |
CN110559014A (zh) | 基于探头补偿的分数阶傅里叶变换回波成像方法与系统 | |
CN105675778B (zh) | 基于双树复小波变换的色谱重叠峰分解方法 | |
CN109884192A (zh) | 用于焊缝特征导波缺陷回波信号特性提取的稀疏表示方法 | |
Xu et al. | Advanced methods for time-of-flight estimation with application to Lamb wave structural health monitoring | |
CN109374748A (zh) | 声发射源时反定位成像方法 | |
Wang et al. | Wavelet packet analysis for the propagation of acoustic emission signals across turbine runners | |
Chen et al. | Ultrasonic signal identification by empirical mode decomposition and Hilbert transform | |
CN105611018B (zh) | 一种mipi lp信号测试系统及方法 | |
Kažys et al. | Ultrasonic detection of defects in strongly attenuating structures using the Hilbert–Huang transform | |
Zhang et al. | Advanced acoustic microimaging using sparse signal representation for the evaluation of microelectronic packages | |
Shi et al. | Signal matching wavelet for ultrasonic flaw detection in high background noise | |
CN105702260A (zh) | 基于matlab的多功能语音测试滤波系统 | |
Shoupeng et al. | Quadrature demodulation based circuit implementation of pulse stream for ultrasonic signal FRI sparse sampling | |
Xu et al. | Analysis of the time-frequency characteristics of internal combustion engine vibration signal based on Hilbert-Huang transform | |
CN102539542A (zh) | 复小波变换法提取陶瓷基芳纶纤维增强复合板中单模式超声导波信息 | |
CN109100687A (zh) | 一种雷达设备lfm脉冲信号pslr参数确定方法 | |
CN112326017B (zh) | 一种基于改进半经典信号分析的微弱信号检测方法 | |
CN100592667C (zh) | 时频紧支信号的小波去噪方法 | |
CN107101715A (zh) | 一种基于交叉小波的变压器振动信号幅频特征量提取方法 | |
CN110221262A (zh) | 一种雷达设备lfm信号主瓣降低确定平台及方法 | |
Lan et al. | Envelope Extraction Algorithm of the Vehicle Acoustic Signal Based on Step-Size Transformation |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080813 Termination date: 20100201 |