CN105548373B - 一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法 - Google Patents
一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法 Download PDFInfo
- Publication number
- CN105548373B CN105548373B CN201510934671.XA CN201510934671A CN105548373B CN 105548373 B CN105548373 B CN 105548373B CN 201510934671 A CN201510934671 A CN 201510934671A CN 105548373 B CN105548373 B CN 105548373B
- Authority
- CN
- China
- Prior art keywords
- golay
- code
- excited
- quasi
- codes
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/34—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor
- G01N29/341—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor with time characteristics
- G01N29/345—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor with time characteristics continuous waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/10—Number of transducers
- G01N2291/106—Number of transducers one or more transducer arrays
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法,该方法包括:由Golay A码进行编码激励,通过软件算法对接收回波实现Golay A→B码激励转换,从而得到Golay A、B两次激励回波(形成正交互补Golay编码激励条件),分别进行脉冲压缩,再进行矢量求和实现解码,实现准单次正交互补Golay编码激励技术,替代传统物理硬件上双次发射方案;本发明采用有效的噪声抑制‑反卷积算法,通过FPGA(Field Programmable Gate Array)实现序列Golay A码激励到Golay B码激励转换算法,达到准单次激励目的,改进了传统Golay A、B码分次激励的方案,提高相控阵系统50%运行效率及大幅度提高信噪比(与编码长度Lc相关),且在一些动态扫查过程中可避免由于两次发射导致位置变化使得接收波形不一致而影响解码效果问题,在实际工程中具有重要应用价值。
Description
技术领域
本发明涉及一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法。
背景技术
超声相控阵系统性能关键技术中,通过增强发射能量来提高仪器的穿透力、分辨力是工业无损探伤关键技术之一,如增大收发孔径、改善探头匹配电路、提高发射电压,但这使硬件系统变得复杂且可控性变差,直接影响到仪器的稳定、可靠性,而采用编码激励技术发射多个连续波间接增强发射能量(接收端对低电压的连续波进行能量聚合)能避免该问题产生。单次激励如Yannis S.Avrithis等(1998)采用CDMA(Code Division MultipleAccess)伪随机编码激励超声系统,较传统单脉冲激励方式有更高图像采集率、横向分辨力和对比度;韩国Jeong J S.学者(2013)使用Barker编码技术抑制高强度聚焦超声带来干扰达-40dB;美国Vanderbilt大学Byram B.等学者(2014)采用Chirp码调频模式,有效抑制超声多路径与波束形成的杂波;但上述解码效果不理想(产生水距离平旁瓣,不能还原成理想的回波信号)且发射电路复杂。双次激励如利福尼亚大学Jinhyoung Park(2010)采用Golay码激励技术与6dB带宽放大器,实现10~110MHz频带内纹波在4dB内,体现卓越性能;中国科学技术大学(2010)用Chirp信号调制Golay互补码激励,增加医学超声透射深度与抗干扰能力;中国医学科学院(2014)用长度为4的Golay互补序列获得了更高的信噪比;理论上虽达到理想解码效果,但需两次发射,效率低、可靠性差。
本文基于优良正交互补Golay码分析编码激励原理,提出一种准单次激励解码技术,能达到理想的解码效果,改善传统Golay编解码技术:两次发射实现一次接收解码以牺牲脉冲重复频率PRF(Pulse Recurrence Frequency)50%的缺点,最后通过构建8位正交互补Golay码实施超声相控阵仪器收发系统控制方案,通过ModelSim仿真其性能,回波幅度增大约16倍,由FPGA(Field Programmable Gate Array)验证其效果。
发明内容
为解决上述技术问题,本发明的目的是提供一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法,对传统A、B码分次激励的方案进行改进,提出准单次激励方案,该方法使得相控阵仪器基于Golay编码激励方式的扫查效率提高50%、避免两次发射、接收信号不一致带来的解码误差,且大幅度提高信噪比。
本发明目的通过以下技术方案来实现:
一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法,该方法包括:
a由Golay A码超声编码激励波形;
b回波通过噪声抑制-反卷积法对A码求反卷积、B码求卷积,同时得到B码激励回波,形成准单次激励;
c分别对A、B码激励回波进行脉冲压缩,再矢量求和实现解码。
与现有技术相比,本发明的一个或多个实施例可以具有如下优点:对传统GolayA、B码分次激励方案进行改进,提出准单次激励方案,通过Golay A码实际激励一次,由软件算法实现Golay A→B码激励转换因子,通过解反卷积方式求解转换因子,考虑到实际应用情况,接收回波信号都携带噪声,提出了基于噪声抑制-反卷积模型,有效抑制噪声的放大。采用FPGA作为实现工具,能快速现实两次激励到准单次激励的转换,避免了物理硬件上以牺牲效率的Golay(A,B)两次激励技术,在超声相控阵仪器前端收发系统中具有重要的实际意义,从而提高其扫查效率50%与信噪比。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例共同用于解释本发明,并不构成对本发明的限制。在附图中:
图1是一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法流程图;
图2是编码电路控制器组成原理框图;
图3是编码激励波形生产过程;
图4是Golay A→B激励转换离散差分方程方框图;
图5是准单次Golay码激励实际解码模型;
图6是码长Lc=8、5MHz探头Golay(A、B)编码激励与解码ModelSim仿真图;
图7是准单次编码激励FPGA解码效果图。
具体实施方式
容易理解,根据本发明的技术方案,在不变更本发明的实质精神下,本领域的一般技术人员可以提出本发明的多个结构方式和制作方法。因此以下具体实施方式以及附图仅是本发明的技术方案的具体说明,而不应当视为本发明的全部或者视为本发明技术方案的限定或限制。
下面结合实施例及附图对本发明作进一步详细的描述。
图1是一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法流程图,包括:
步骤10由Golay A码编码激励的方案,实现步骤如下:
Golay码是由简单的数字符号1、-1组成,不能直接激励探头,必须由经符号→波形转换产生激励波形。图2为编码电路控制器组成原理框图,其主要步骤如下:
①基波生成模块,产生标准的双极性方波BW(n),周期为探头振荡周期TP,与探头激励系统形成共振,能量输出最大化;
②码元经R倍内插保持器(R≥Δdc),Δdc为探头一次激励振荡的时间长度,内插后得到GR(n),保证一个码元对应探头振荡周期时间;
③卷积调相,GR(n)与基波进行卷积,得到激励探头波形PW(n)(PW(n)=BW(n)*GR(n))(*表卷积)对基波进行相位调制,实现相位变换产生编码激励信号,如图3编码激励波形生产过程。
步骤20回波通过噪声抑制-反卷积法对Golay A码求反卷积、Golay B码求卷积,同时得到Golay B码激励回波,形成准单次激励,实现方法如下:
①设整个硬件系统的传递函数为hd(n),单脉冲激励信号e(n),接收回波信号y(n)=e(n)*hd(n),采用二元两相调制编码激励函数F[e(n)],则输出y(n)=F[e(n)]*hd(n),若形成理想编码激励,则有:
F[e(n)]*hd(n)=Np·δ(n)*e(n)*hd(n),L∈Z+ (1)
②目前正交互补Golay码具有理想编码激励性质,设一对Lc长度的正交互补Golay(A、B)码其码元序列分别为二相序列an、bn,令A={an},an∈(+1,-1),{n∈N|n≤Lc-1}和B={bn},bn∈(+1,-1),{n∈N|n≤Lc-1},根据Golay码定义,则A、B序列满足:
{an}*{a(-n)}+{bn}*{b(-n)}=2Lc·δ(n) (2)
方便表示可用an代表A码编码函数,则其激励探头回波信号为:
③图4为Golay A→B激励转换离散差分方程方框图,由hd(n)*a(n)计算hd(n)*b(n),可得到双次激励回波,为正确解码提供保证。B码激励转换因子b(n)*-1a(n)(对A码求反卷积),即
式中*-1表反卷积运算,σa、σb∈{-1,1}为编码符号极性。根据实际应用对含噪回波信号加以噪声抑制处理以降低反卷积对噪声的放大作用,可通过添加降噪模块如FIR(Finite Impulse Response)或Wavelet,如图5准单次Golay码激励实际解码模型中①②③标签位置模块,对噪声进行抑制。
步骤30分别对Golay A、Golay B码进行脉冲压缩再矢量求和实现解码,实现步骤如下:
①脉冲压缩
由于发送的是一串连续波,则对应接收也为一串连续波,必须压缩成单脉冲回波,才能有效对超声回波分析、判断。
下面描述解码过程,由式(2)可用相同的编码序列与对应的回波做自相关运算,进行脉冲压缩,分别得:
其中表示相关运算。
②矢量求和
由于Golay A、B码激励回波,经脉冲压缩后,除了大部分主瓣能量外,含有少量的旁瓣能量,脉冲压缩后的两路波形主瓣能量相同相位一致,而旁瓣能量相同相位相反,根据这个性质,可进行矢量求和消除旁瓣同时聚合主瓣能量实现解码。
经矢量求和后得到解码输出:
可见,码长Lc的正交互补Golay编码激励可产生2Lc倍增益,达到高能单脉冲激励效果,实现理想编解码。
如8位正交互补Golay码A={1 1 1 -1 1 1 -1 1}、B={1 -1 1 1 1 -1 -1 -1}
Yout=Ya+Yb=hd(n)*(0,0,0,0,0,0,0,16,0,0,0,0,0,0,0)=16·δ(n-7)*hd(n)。
解码结果只是有些延时,可知是符合正交互补Golay码理想编码激励性质。
图6为码长Lc=8、5MHz探头Golay(A、B)编码激励与解码ModelSim仿真图,
图7为准单次编码激励FPGA解码效果图。
上述激励方法通过噪声抑制-反卷积法,实现A码激励回波到B码激励回波的转换,应用FPGA快速运算能力,快速实现求解反卷积与卷积运算,即求解A→B激励转换因子:
Golay A→B码激励转换是物理硬件上实现双次发射到单次发射转换基础,为避免正交互补Golay码两次激励(A,B)序列码方能进行一次有效解码的劣势,提出由Golay A码激励一次,通过软件算法,变换到Golay B码激励,变相实现正交序列对的两次激励,即实现Golay A→B码激励转换。实际应用中,超声回波信号会携带噪声,准单次Golay码激励技术单纯的解反卷积,会使得部分噪声放大,导致恢复的原始信号严重失真,可通过添加降噪模块如FIR或Wavelet对噪声进行抑制,获得较好效果。
准单次正交互补Golay(A,B)码超声相控阵编码激励方法,改变了基于传统Golay编码激励所需两次激励方能达到理想解码方案,通过Golay A→B码激励转换,同时形成Golay A、B码激励波形,由相关逆向滤波法分别进行脉冲压缩再进行矢量叠加实现解码,可提高50%效率与信噪比,所述步骤具体包括:
脉冲压缩方式解码过程是一个自相关过程即回波与激励码序列做相关运算,可把多个连续回波压缩成单个波形,结合Golay码相关互补性质,对A、B码激励回波做脉冲压缩,再矢量相加后得到解码输出,而不会产生附加噪声(距离旁瓣为0)。
正交互补Golay(A、B)编码激励技术能提高回波增益,而不会带来噪声,达到理想解码效果。
⑴研究准单次激励正交互补Golay编解码实现技术。在不修改硬件电路条件下,通过发射一系列编码的连续波形,再对接收的连续波进行脉冲压缩成单脉冲激励回波实现解码,间接增大发射能量,具有通用性强、移植性好、灵活性高等特点,有重要意义;
⑵本发明设计准单次激励正交互补Golay编解码总体方案,从互补Golay码性质入手,分析传统双次Golay激励模式,通过Golay A→B码激励转换因子、脉冲压缩与矢量合成等核心内容,提出准单次激励方案,使得相控阵仪器基于Golay编码激励的扫查效率提高50%,弥补了双次发射带来的误差与低效率的不足,大幅度提高了信噪比。
虽然本发明所揭露的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (3)
1.一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法,该方法包括:
a由Golay A码超声编码激励波形;
b回波通过噪声抑制-反卷积法对Golay A码求反卷积、Golay B码求卷积,同时得到Golay B码激励回波,形成准单次激励;
c分别对Golay A、Golay B码激励回波进行脉冲压缩,再矢量求和实现解码;
所述步骤a具体包括:
通过对一定码长Lc的Golay A码元进行二相编码调制成激励探头的波形,其中码元符号相位映射为1→0、-1→π,包括根据实际应用设计编码距离Δdc与对码元进行重采样进行长度扩张,再与基波信号卷积,形成可以驱动发射电路的编码序列串激励探头。
2.根据权利要求1所述的准单次正交互补Golay(A,B)码超声相控阵编码激励方法,其特征在于,所述步骤b具体包括:
接收回波通过噪声抑制法对Golay A码求反卷积、Golay B码求卷积得到Golay B码激励回波,即由Golay A码激励一次,可获得Golay A、Golay B两次激励回波,达到互补Golay(A、B)分次激励效果,实现准单次激励。
3.根据权利要求1所述的准单次正交互补Golay(A,B)码超声相控阵编码激励方法,其特征在于,所述步骤c具体包括:
依靠Golay A、Golay B码正交互补关系,分别对Golay A、Golay B码激励回波进行相关逆向滤波实现脉冲压缩,把接收到的多个连续波进行能量压缩,再矢量求和抵消旁瓣、增加主瓣,实现解码,增大回波幅度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510934671.XA CN105548373B (zh) | 2015-12-15 | 2015-12-15 | 一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510934671.XA CN105548373B (zh) | 2015-12-15 | 2015-12-15 | 一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105548373A CN105548373A (zh) | 2016-05-04 |
CN105548373B true CN105548373B (zh) | 2018-03-27 |
Family
ID=55827706
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510934671.XA Active CN105548373B (zh) | 2015-12-15 | 2015-12-15 | 一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105548373B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105866257B (zh) * | 2016-06-12 | 2018-05-08 | 华南理工大学 | 一种基于Golay互补卷积码的多元调制超声编码单次激励方法 |
CN111123272B (zh) * | 2018-10-31 | 2022-02-22 | 无锡祥生医疗科技股份有限公司 | 单极系统的戈莱码编码激励方法和解码方法 |
CN110575201B (zh) * | 2019-10-09 | 2021-12-24 | 珠海医凯电子科技有限公司 | 基于反相Golay码的超声微泡空化成像方法及装置 |
CN111307956B (zh) * | 2019-12-19 | 2023-08-01 | 中国航空工业集团公司北京长城航空测控技术研究所 | 一种基于线性调频信号的导波信号激励电路 |
CN113203374B (zh) * | 2021-06-01 | 2022-05-03 | 哈尔滨工业大学 | 一种基于脉冲压缩的电磁超声测厚装置及方法 |
CN113849360B (zh) * | 2021-09-24 | 2024-01-26 | 北京润科通用技术有限公司 | 一种波形激励方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1266989A (zh) * | 1999-03-16 | 2000-09-20 | 通用电气公司 | 有使用单极或者双极编码激励的波束形成的超声成象系统 |
CN101149433A (zh) * | 2007-11-01 | 2008-03-26 | 绵阳索尼克电子有限责任公司 | 一种使用反向Golay码译码的超声运动检测方法 |
CN104765027A (zh) * | 2014-12-29 | 2015-07-08 | 中国科学院苏州生物医学工程技术研究所 | 高频超声编码激励发射系统及脉冲产生方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7523640B2 (en) * | 2005-08-01 | 2009-04-28 | Baker Hughes Incorporated | Acoustic fluid analyzer |
TWI507709B (zh) * | 2013-03-20 | 2015-11-11 | Univ Nat Taiwan | 藉由格雷編碼激發之超音波都卜勒偵測方法 |
-
2015
- 2015-12-15 CN CN201510934671.XA patent/CN105548373B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1266989A (zh) * | 1999-03-16 | 2000-09-20 | 通用电气公司 | 有使用单极或者双极编码激励的波束形成的超声成象系统 |
CN101149433A (zh) * | 2007-11-01 | 2008-03-26 | 绵阳索尼克电子有限责任公司 | 一种使用反向Golay码译码的超声运动检测方法 |
CN104765027A (zh) * | 2014-12-29 | 2015-07-08 | 中国科学院苏州生物医学工程技术研究所 | 高频超声编码激励发射系统及脉冲产生方法 |
Non-Patent Citations (2)
Title |
---|
Design and Performance of Huffman Sequences in Medical Ultrasound Coded Excitation;Polpetta Alessandro et al.;《IEEE TRANSACTIONS ON ULTRASONICS FERROELECTRICS AND FREQUENCY CONTROL》;20120430;第59卷(第4期);第630-647页 * |
基于相位编码技术的超大功率超声收发系统研究;邓勇;《电子设计工程》;20150630;第23卷(第11期);第122-124页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105548373A (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105548373B (zh) | 一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法 | |
US6558328B2 (en) | Method and apparatus using coded excitation in ultrasonic perfusion imaging | |
CN100475148C (zh) | 超声成像方法与超声成像设备 | |
CN110998361B (zh) | 用于编码的多脉冲超声造影成像的方法 | |
CN105120759B (zh) | 用于从旋转血管内超声设备产生图像的系统和方法 | |
CN105866257B (zh) | 一种基于Golay互补卷积码的多元调制超声编码单次激励方法 | |
US6050947A (en) | Method and apparatus for harmonic tissue imaging and contrast imaging using coded transmission | |
Fu et al. | Barker coded excitation with linear frequency modulated carrier for ultrasonic imaging | |
Kim et al. | Barker-sequence-modulated golay coded excitation technique for ultrasound imaging | |
Wang et al. | Pseudo Chirp-Barker-Golay coded excitation in ultrasound imaging | |
JP2003190157A (ja) | 超音波診断装置 | |
Zhao et al. | Coded excitation for diverging wave cardiac imaging: A feasibility study | |
US12004908B2 (en) | Ultrasound imaging system for high resolution wideband harmonic imaging | |
CN106814352B (zh) | 一种基于格雷互补波形的多目标检测方法 | |
US8388536B2 (en) | Ultrasonic imaging apparatus | |
CN104793211A (zh) | 一种风廓线雷达干扰抑制方法 | |
Fu et al. | Barker coded excitation using LFM carrier for improving axial resolution in ultrasound imaging | |
Polpetta et al. | Design and performance of Huffman sequences in medical ultrasound coded excitation | |
Lahav et al. | FoCUS: Fourier-based coded ultrasound | |
CN101642378B (zh) | 控制谐波成像的发射脉冲频谱的方法和装置 | |
Tamraoui et al. | Complete Complementary Coded Excitation Scheme for SNR Improvement of 2D Sparse Array Ultrasound Imaging | |
JP5612984B2 (ja) | 超音波診断装置 | |
Rosli et al. | Design of amplitude and phase modulated pulse trains with good auttocorrelation properties for radar communications | |
Yang et al. | Design of orthogonal coded excitation for synthetic aperture imaging in ultrasound systems | |
US20230258788A1 (en) | Method and system using coded excitation with invertible skew-symmetric binary sequences and inverse filtering |
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 |