CN102613989A - 侦测系统与其信号处理方法 - Google Patents

侦测系统与其信号处理方法 Download PDF

Info

Publication number
CN102613989A
CN102613989A CN201110123291XA CN201110123291A CN102613989A CN 102613989 A CN102613989 A CN 102613989A CN 201110123291X A CN201110123291X A CN 201110123291XA CN 201110123291 A CN201110123291 A CN 201110123291A CN 102613989 A CN102613989 A CN 102613989A
Authority
CN
China
Prior art keywords
vector matrix
matrix
data
weights
matrix data
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
Application number
CN201110123291XA
Other languages
English (en)
Other versions
CN102613989B (zh
Inventor
蔡铭嘉
吴安宇
李百祺
詹承洲
陈郁豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Publication of CN102613989A publication Critical patent/CN102613989A/zh
Application granted granted Critical
Publication of CN102613989B publication Critical patent/CN102613989B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52046—Techniques for image enhancement involving transmitter or receiver
    • G01S7/52047—Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
    • 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
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明涉及一种信号处理方法，适于处理多个向量矩阵数据，以侦测默认范围内的影像，而且这些向量矩阵数据是依据多个超音波波束在默认范围被反射所产生。本发明的信号处理方法包括将预设时间区间内所有的向量矩阵数据相加，而产生总向量矩阵。另外，将此总向量矩阵乘上被转置的总向量矩阵，而获得自相关矩阵，并且依据自相关矩阵的反矩阵而获得权值。接着，依据权值而对在预设时间区间内所获得的向量矩阵数据进行加权运算，而获得加权运算结果，以进行影像合成作业。本发明侦测系统与其信号处理方法可以侦测默认范围内的影像信息。

Description

侦测系统与其信号处理方法
技术领域
本发明涉及一种信号处理的方法，尤其涉及一种用于超音波成像系统中的信号处理的方法。
背景技术
超音波主要是由压晶体管在电场作用下产生的机械振动波，通常频率超过20 kHz即被认定为超音波。目前的应用以超音波为工具来检验、测量或控制，例如测量厚度、测量距离、医学治疗、医学诊断或超音波成像等。另外，也可以利用超音波处理物质进而改变或加速改变物质的一些物理、化学、生物特性或状态，例如利用超音波在液体中的「空穴效应」来完成加工、清洗、焊接、乳化、粉碎、脱气、促进化学反应或医疗等目的。
在现有的超音波成像系统中，当收到由反射的超音波波束所产生的向量矩阵时，会将此向量矩阵乘以转置后的向量矩阵，而产生自相关矩阵。接着，将预设时间内所获得的所有自相关矩阵相加，而产生总自相关矩阵。此时，再将总自相关矩阵进行反矩阵运算，以计算权值，作为后续影像合成的参数。
由于每次取得向量矩阵时，都需要将此向量矩阵进行转置运算，以获得自相关矩阵，因此增加了处理的时间以及运算复杂度。另外，由于总自相关矩阵非常庞大，因此在进行反矩阵运算时，也会增加运算的复杂度。如此一来，就导致了整个系统的复杂度上升。
发明内容
因此，本发明提供一种侦测系统，可以侦测默认范围内的影像信息。
另外，本发明也提供一种信号处理的方法，可以适用于超音波成像系统，并且简化系统运算的复杂度。
为达上述优点，本发明提出一种侦测系统，包括超音波模块、多个接收单元、多个模拟数字转换器、处理模块和影像合成单元。超音波模块包括以阵列方式排列的多个超音波单元，并能朝预设范围连续地发射多个超音波波束。当超音波波述在默认范围内被反射，而分别由接收单元所接收时，这些接收单元就会分别产生多个频道信号。而每一频道信号会由分别由对应的模拟数字转换器转换为数字数据，而产生向量矩阵数据。此时，处理模块会将在默认时间区间内所接收到的向量矩阵数据相加，而产生总向量矩阵数据，再将总向量矩阵数据和转置后的总向量矩阵数据相乘，而获得自相关矩阵。接着，处理模块会将此自相关矩阵进行反矩阵运算，并且依据反矩阵运算后的自相关矩阵而获得权值，以对在预设时间区间内所获得的向量矩阵数据进行加权运算，而产生加权运算结果。
在本发明的一个实施例中，上述的处理模块包括权值运算单元、参数运算单元和乘法器。权值运算单元是依据向量矩阵数据而产生自相关矩阵和权值；而参数运算单元则是依据向量矩阵数据而产生相关参数函数。另外，乘法器则是耦接权值运算单元和参数运算单元，以将相关参数函数乘上权值，而完成将向量矩阵数据进行加权运算，而产生加权运算结果。
本发明另提供一种信号处理方法，适于处理多个向量矩阵数据，以侦测默认范围内的影像，而且这些向量矩阵数据是依据多个超音波波束在默认范围被反射所产生。本发明的信号处理方法包括将预设时间区间内所有的向量矩阵数据相加，而产生总向量矩阵。另外，将此总向量矩阵乘上被转置的总向量矩阵，而获得自相关矩阵，并且依据自相关矩阵的反矩阵而获得权值。接着，依据权值而对在预设时间区间内所获得的向量矩阵数据进行加权运算，而获得加权运算结果，以进行影像合成作业。
由于在本发明侦测系统与其信号处理方法中，处理模块是先获得总向量矩阵，然后再计算自相关矩阵，并且进行自相关矩阵的反矩阵运算。因此，本发明可以有效地降低系统运算的复杂度。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂，以下特举较佳实施例，并配合附图，详细说明如下。
附图说明
图1为依照本发明的一个较佳实施例的一种侦测系统的方块图。
图2为依照本发明的一个较佳实施例的一种处理模块的架构图。
【主要元件符号说明】
100：侦测系统
102：超音波模块
104、106、108、110：超音波单元
120：信号接收级
122[0:N]：接收单元
124[0:N]：放大器
126[0:N]：模拟数字转换器(ADC)
130：信号处理级
132[0:N]：解调器
134[0:N]：缓冲器
136[0:N]：时间延迟相位旋转器
138：处理模块
140：后端影像合成级
142：缓冲器
144：低通滤波器(LPF)
146：影像合成单元
202：权值运算单元
204：参数运算单元
206：乘法器
CH[0:N]：频道信号 
DATA[0:N]：数字料信号
De_MOD[0:N]：解调信号
FCF(t)：弹性自相关参数函数
IMG：影像信息
M_Data：加权运算结果
x(t)：向量矩阵
WMVDR(t)：权值。
具体实施方式
为更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效，以下结合附图及较佳实施例，对依据本发明提出的侦测系统与其信号处理方法其具体实施方式、方法、步骤、结构、特征及功效，详细说明如后。
有关本发明的前述及其他技术内容、特点及功效，在以下配合参考图式的较佳实施例详细说明中将可清楚的呈现。通过具体实施方式的说明，可对本发明为达成预定目的所采取的技术手段及功效有一更加深入且具体的了解，然而所附图式仅是提供参考与说明之用，并非用来对本发明加以限制。
图1为依照本发明的一个较佳实施例的一种侦测系统的方块图。请参考图1，本实施例所提供的侦测系统100包括超音波模块102，其具有N个超音波单元，例如104、106、108和110，是以阵列方式排列，其中N为大于或等于1的正整数。在本实施例中，这些超音波单元104、106、108和110会朝预设范围连续地发出多个超音波波束。
请继续参考图1，侦测系统100还包括信号接收级120、信号处理级130和后端影像合成级140。信号接收级120包括多个接收单元122[0:N]、多个放大器124[0:N]和多个模拟数字转换器(ADC) 126[0:N]。接收单元122[0:N]可以分别接收在预设范围内被反射的超音波波束，并且产生多个频道信号CH[0:N]给放大器124[0:N]。接着，放大器124[0:N]会分别将所接收到的频道信号CH[0:N]进行放大，然后再传送至ADC 126[0:N]。此时，ADC 126[0:N]会将放大后的频道信号CH[0:N]转换为多笔数字数据信号DATA[0:N]给信号处理级130。
信号处理级130包括多个解调器132[0:N]、多个缓冲器134[0:N]、多个时间延迟相位旋转器136[0:N]和处理模块138。其中，解调器132[0:N]会分别耦接至ADC 126[0:N]，以接收数字数据DATA[0:N]，并且加以解调，而产生多个解调信号De_MOD[0:N]。这些解调信号De_MOD[0:N]会通过缓冲器134[0:N]，并且被送至时间延迟相位转换器136[0:N]，以进行时间延迟和相位旋转，并且产生向量矩阵数据x(t)。接着，此向量矩阵数据x(t)会被送至处理模块138进行处理。
特别的是，在本实施例中，当处理模块138收到向量矩阵x(t)后，并不是先进行自相关矩阵的运算，而是将预设时间内所有的向量矩阵数据x(t)相加，而产生总SU3向量矩阵。
图2为依照本发明的一个较佳实施例的一种处理模块的架构图。请参考图2，本实施例所提供的处理模块208包括权值运算单元202、参数运算单元204和乘法器206。权值运算单元202可以接收向量矩阵数据x(t)，并且将预设时间内所获得的向量矩阵数据x(t)相加，而获得总向量矩阵y(t)，而此其可以表示为：
。
其中，K为整数。
接着，当获得总向量矩阵后，可以将总向量矩阵乘以经过转置运算后的总向量矩阵，而获得自相关矩阵(
Figure 201110123291X100002DEST_PATH_IMAGE002
)，以上的叙述可以利用下式来表示：
Figure 209633DEST_PATH_IMAGE003
其中δ为常数，而I则为单位矩阵。
接着，权值运算单元202可以依据下式，而将自相关矩阵
Figure 989370DEST_PATH_IMAGE002
进行反矩阵运算：
由于在上式等号右侧的第二个操作数中的分母是常数，因此会使得整个表达式的计算变得简单。
另外，权值运算单元202还会依据自相关矩阵的反矩阵，来计算权值(WMVDR(t))，其表示如下：
Figure 201110123291X100002DEST_PATH_IMAGE006
其中a为单位向量。
请继续参考图2，另一方面，参数运算单元204也会接收向量矩阵数据x(t)，并且依据此向量矩阵数据x(t)而计算弹性自相关参数函数(FCF(t))，如下所示：
Figure 323585DEST_PATH_IMAGE007
其中m建议为大于0而小于等于1的值。
另外，权值运算单元202和参数运算单元204的输出端都会耦接乘法器206。因此，乘法器206会将权值WMVDR(t)乘上弹性自相关参数函数FCF(t)，以进行加权运算，并且获得加权运算结果M_Data给后端影像合成级140，以进行影像合成作业。
后端影像合成级140包括缓冲器142、低通滤波器(LPF) 144和影像合成单元146。当加权运算结果M_Data被送至后端影像合成级140后，会先由缓冲器142接收，并且将其输出给低通滤波器144来进行低通滤波，以滤除噪声。接着，进行完低通滤波后的加权运算结果M_Data会被送至影像合成单元146。藉此，影像合成单元146可以依据加权运算结果M_Data，而得到关于默认范围内的影像信息IMG。
综上所述，由于在本发明侦测系统与其信号处理方法中，处理模块是先将默认时间内所有的向量矩阵相加，再计算自相关矩阵，因此可以简化运算复杂度。另外，利用以上方式所获得的自相关矩阵的反矩阵运算也较为单纯，因此可以更进一步的简化系统运算的复杂度。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。

Claims (9)

1.一种侦测系统，其特征是，包括：
超音波模块，包括以阵列方式排列的多个超音波单元，并能朝预设范围连续地发射多个超音波波束，所述超音波波束在预设范围内被反射；
多个接收单元，分别接收被反射的超音波波束，并产生多个频道信号；
多个模拟数字转换器，分别将所述多个频道信号转换为数字数据，以产生向量矩阵数据；
处理模块，将在默认时间区间内所接收到的向量矩阵数据相加，而产生总向量矩阵数据，再将所述总向量矩阵数据和转置后的总向量矩阵数据相乘，而获得自相关矩阵，且所述处理模块还将所述自相关矩阵进行反矩阵运算，并依据反矩阵运算后的自相关矩阵而获得权值，以对在所述预设时间区间内所获得的向量矩阵数据进行加权运算，而产生加权运算结果；以及
影像合成单元，依据所述加权运算结果而合成影像数据。
2. 根据权利要求1所述的侦测系统，其特征是：所述处理模块还包括：
权值运算单元，用以依据所述向量矩阵数据而产生所述自相关矩阵和所述权值；
参数运算单元，依据所述向量矩阵数据而产生相关参数函数；以及
乘法器，耦接所述权值运算单元和所述参数运算单元，以将所述相关参数函数乘上所述权值，而完成将所述向量矩阵数据进行加权运算，而产生所述加权运算结果。
3. 根据权利要求1所述的侦测系统，其特征是：所述侦测系统还包括多个放大器，分别耦接所述这些接收单元，以将所述这些频道信号进行放大，并将放大后的频道信号送至所述模拟数字转换器。
4.根据权利要求1所述的侦测系统，其特征是：所述侦测系统还包括：
多个解调器，分别耦接所述模拟数字转换器，以将所述这些数字数据进行解调；
多个第一缓冲器，分别耦接所述这些解调器，以接收解调后的数字数据；以及
多个时间延迟相位旋转器，分别耦接所述这些第一缓冲器，以将所述这些解调后的数字数据进行时间延迟和相位旋转，而产生所述这些向量矩阵数据。
5. 根据权利要求1所述的侦测系统，其特征是：所述侦测系统还包括：
第二缓冲器，耦接所述处理模块，以接收所述加权运算结果；以及
低通滤波器，耦接所述第二缓冲器，用以将所述加权运算结果进行低通滤波处理来滤除噪声，并且输出给所述影像合成单元。
6. 一种信号处理方法，适于处理多个向量矩阵数据，以侦测默认范围内的影像，且所述这些向量矩阵数据是依据多个超音波波束在所述默认范围被反射所产生，其特征是，所述信号处理方法包括下列步骤：
将预设时间区间内所有的向量矩阵数据相加，而产生总向量矩阵；
将所述总向量矩阵乘上被转置的总向量矩阵，而获得自相关矩阵；
依据所述自相关矩阵的反矩阵而获得权值；以及
依据所述权值而对在所述预设时间区间内所获得的向量矩阵数据进行加权运算，而获得加权运算结果，以进行所述影像合成作业。
7. 根据权利要求6所述的信号处理方法，其特征是：所述产生所述自相关矩阵的反矩阵的步骤，包括执行下列运算：
Figure 201110123291X100001DEST_PATH_IMAGE002
其中y(t)是总向量矩阵、δ是常数、而I是单位矩阵。
8. 根据权利要求6所述的信号处理方法，其特征是：所述产生所述权值的步骤，包括执行下列运算：
其中
Figure 201110123291X100001DEST_PATH_IMAGE006
是所述自相关矩阵，而a则为单位向量。
9. 根据权利要求6所述的信号处理方法，其特征是：所述获得所述加权运算结果的步骤，是将所述权值乘以弹性自相关参数函数，而产生所述弹性自相关参数函数的步骤，包括执行下列运算：
Figure 201110123291X100001DEST_PATH_IMAGE008
其中xn(t)是每一所述这些反射的超音波束所对应的向量函数、N为所述这些超音波束的总数、而m则是大于0而小于等于1。
CN201110123291.XA 2011-01-27 2011-05-13 侦测系统与其信号处理方法 Expired - Fee Related CN102613989B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
TW100103157A TW201232476A (en) 2011-01-27 2011-01-27 Detection system and signal processing method thereof
TW100103157 2011-01-27

Publications (2)

Publication Number Publication Date
CN102613989A true CN102613989A (zh) 2012-08-01
CN102613989B CN102613989B (zh) 2014-07-02

Family

ID=46554454

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110123291.XA Expired - Fee Related CN102613989B (zh) 2011-01-27 2011-05-13 侦测系统与其信号处理方法

Country Status (3)

Country Link
US (1) US20120192650A1 (zh)
CN (1) CN102613989B (zh)
TW (1) TW201232476A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104956218A (zh) * 2013-02-07 2015-09-30 西门子公司 用于在不规则的测量时改善saft分析的方法和装置
CN113727240A (zh) * 2020-05-26 2021-11-30 南宁富桂精密工业有限公司 声音播放调整方法及可携式装置
CN114609245A (zh) * 2022-05-10 2022-06-10 西南交通大学 一种阵列式超声导波仪、钢结构损伤监测系统及方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5504675A (en) * 1994-12-22 1996-04-02 International Business Machines Corporation Method and apparatus for automatic selection and presentation of sales promotion programs
US6141393A (en) * 1999-03-03 2000-10-31 Motorola, Inc. Method and device for channel estimation, equalization, and interference suppression
CN1575772A (zh) * 2003-06-25 2005-02-09 阿洛卡株式会社 超声波诊断装置
US20090024685A1 (en) * 2007-07-19 2009-01-22 Itt Manufacturing Enterprises, Inc. High Speed and Efficient Matrix Multiplication Hardware Module
CN101904751A (zh) * 2009-06-02 2010-12-08 信飞系统公司 超声信号压缩
CN101919707A (zh) * 2009-06-09 2010-12-22 株式会社东芝 超声波诊断装置、医用图像处理装置以及图像处理方法
CN101933814A (zh) * 2009-06-29 2011-01-05 信飞系统公司 超声系统中的波束形成后压缩

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5504675A (en) * 1994-12-22 1996-04-02 International Business Machines Corporation Method and apparatus for automatic selection and presentation of sales promotion programs
US6141393A (en) * 1999-03-03 2000-10-31 Motorola, Inc. Method and device for channel estimation, equalization, and interference suppression
CN1575772A (zh) * 2003-06-25 2005-02-09 阿洛卡株式会社 超声波诊断装置
US20090024685A1 (en) * 2007-07-19 2009-01-22 Itt Manufacturing Enterprises, Inc. High Speed and Efficient Matrix Multiplication Hardware Module
CN101904751A (zh) * 2009-06-02 2010-12-08 信飞系统公司 超声信号压缩
CN101919707A (zh) * 2009-06-09 2010-12-22 株式会社东芝 超声波诊断装置、医用图像处理装置以及图像处理方法
CN101933814A (zh) * 2009-06-29 2011-01-05 信飞系统公司 超声系统中的波束形成后压缩

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104956218A (zh) * 2013-02-07 2015-09-30 西门子公司 用于在不规则的测量时改善saft分析的方法和装置
CN104956218B (zh) * 2013-02-07 2018-02-09 西门子公司 用于在不规则的测量时改善saft分析的方法和装置
US10222352B2 (en) 2013-02-07 2019-03-05 Siemens Aktiengesellschaft Method and device for improving the SAFT analysis when measuring irregularities
CN113727240A (zh) * 2020-05-26 2021-11-30 南宁富桂精密工业有限公司 声音播放调整方法及可携式装置
CN114609245A (zh) * 2022-05-10 2022-06-10 西南交通大学 一种阵列式超声导波仪、钢结构损伤监测系统及方法
CN114609245B (zh) * 2022-05-10 2022-08-23 西南交通大学 一种阵列式超声导波仪、钢结构损伤监测系统及方法

Also Published As

Publication number Publication date
TW201232476A (en) 2012-08-01
CN102613989B (zh) 2014-07-02
US20120192650A1 (en) 2012-08-02

Similar Documents

Publication Publication Date Title
KR20180006952A (ko) 자율적 초음파 프로브, 및 관련된 장치 및 방법
CN104414687A (zh) 超声波测定装置、超声波图像装置及超声波测定方法
US6530890B2 (en) Ultrasound diagnostic apparatus and method for measuring blood flow velocity using doppler effect
US9704474B2 (en) Ultrasound imaging probe with sigma-delta beamformer and apodization therein
CN102613989B (zh) 侦测系统与其信号处理方法
WO2015126446A1 (en) Acoustic signal phase rotation processor
CN101461720B (zh) 基于频谱多普勒的运动速度测量范围的调整方法及装置
US7002876B2 (en) Acoustic-propagation-time measuring apparatus
JP6198204B2 (ja) 超音波診断装置及び超音波診断装置の作動方法
JP2017227480A (ja) 音波監視装置および航走体
CN117192557B (zh) 一种基于厚膜电路的精准超声系统、方法、设备及介质
CN103181781B (zh) 超声波成像系统及其图像处理方法
JPH11347035A (ja) 超音波診断装置
Terez Fundamental frequency estimation using signal embedding in state space
Massa et al. A wideband moving coil electrodynamic transducer system for autonomous underwater vehicle-based geoacoustic inversion
Bogdan et al. Quiet micro boats: An inexpensive acoustic sensing suite
Mani et al. Rain parameter estimation using impact generated low‐frequency acoustic signals
CN103808809A (zh) 一种基于线性调频信号激励的超声回波信号实时处理装置及处理方法
Bochud et al. Toward real-time assessment of material properties using elastic guided waves
Weidner et al. Investigating bubble transport and fate in the watercolumn with calibrated broadband split-beam echosounder data
Cray et al. High‐frequency laser sonar system
Aranchuk et al. Fast acoustic landmine detection using multiple beam laser Doppler vibrometry
Dwyer Robust sequential detection using acoustical arrays in undefined noise
KR20130067015A (ko) 창 함수 비선형 에너지 연산자를 이용한 투과 초음파 단층 촬영기에서 최초 도달 펄스를 검출하기 위한 방법
Kim et al. Underwater acoustic intensity measurements using a pressure/particle‐acceleration intensity vector sensor

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: 20140702

Termination date: 20180513