一种B超回波信号包络的快速提取方法.pdfVIP

Pt．2 第29卷第6期 声学技术 V01．29，No．6 2010年12月 TechnicalAcoustics Dec．，2010 一种B超回波信号包络的快速提取方法 蒲杰，孟晓辉，肖 灵，吴文焘 中国科学院声学研究所，北京，100190 A forfast the ofultrasonicecho approachextractingenvelop signal PU Wen-tao Jie，MENGXiao-hui，XIAOLing，WU InstituteofAcoustics，ChineseAcademyofSciences，Beijing10019nChina 1引 言 3 B超成像中，为满足回波波束形成的延时要求， 回波采样率通常较高，一般为40．60MHz，远大于回 波带宽，波束形成后的数据量较大，无法实现成像 出回波包络P f ：正丽瓦虿万。 数据的实时传输和处理，因此需要降采样后提取用 由 2 ， 3 式可知，可以直接在频域截取带通信 于成像的回波包络信息。本文设计了一种超声回波 号通带内频谱，求得正交解调信号的频谱。 信号包络的快速提取方法，对回波波束形成后信号 2．2频域降采样包络提取方法 后文简称回波信号 降采样的同时，完成回波包 根据上述包络提取原理，本文提出一种频域降 络的计算，具有运算量低和存储空间需求小的特点。 采样包络提取方法。 2原理及实现方法 例如，对于中心频率为3．5MHz的超声探头， 其回波信号带宽约为5MHz，若系统设计采样率为 本文利用回波信号为带通信号的特性，在快速 傅里叶变换过程中仅计算回波信号通带内的频谱， 可以很好地降低计算量，减小计算过程中对存储空 间需求。计算得到的通带内的频谱可视作频率搬移 采样的回波信号长度约为8192点。 后已进行抗混叠滤波的降采样信号的频谱，利用正 假定回波信号的快速傅里叶变换结果为只约， 交解调原理，对所得频谱进行适当处理后进行快速 傅里叶反变换，可直接获取降采样后回波信号的正 心频率3．5MHz，带宽5MHz ，设定 交解调信号，从而计算出回波包络，避免不必要数 据点的计算。 D以m 即对以约进行频谱搬移以及降采样。 利用傅里叶变换基本性质， 2．1包络提取原理 回波信号包络的提取是基于正交解调原理u1， 若f t 为实信号， 假定s 力为带通信号，09。为其中心调制频率，e O ／O ——塑竺!二专F co ，则F co F+ 一国 为带通信号包络，则有： 可推知：DF m 的傅里叶反变换An 的实部为正交 J f e t cos co。f+9 —塑!!!马s co 1 s， f s t cos coof 0．5[e t cosO+e t cos 2coof+臼 】 2 ·—』竺竺斗0．5×【S co+coo +S co一090 】 利用∞? 和Q@ 就可计算出降采样后的回波包 作者简介：蒲杰 1983一 ，男，重庆，汉，在读博士，研究方向为信号处理 络e n 。 Email：pujie@mail．ioa．ac．cn 由上面计算过程可知，计算中仅仅使用到回波 第6期 蒲杰等：一种B超回被信号包络的快速提取方法 信号通带内的正频率频谱DF m ，因此不需要计算 8192点回波数据的所有频谱点，本文拟采用快速傅复数加法。 里叶变换蝶形算法结构形式…，但仅计算所需要的 采用本文所提算法，需要进行 频率分量DF m ，其余频率分量不予计算，这样可 以减少计算量，同时降低计算过程中对高速存储空 次复数乘法与复数加法。由此可见本文提出的方法 间的需求。具体实现形式如图1所示： 较常规方法减少了约一半运算量。 3仿真结果 为了验证本方法的正确性，使用FleldII程序引 生成一条回波信号进行仿真，仿真结果如图2所示 取部分数据点显示 。 目#＆《m＆月口e ’ 。2 0 _o2 座机电话号码3200 座机电话号$女*月女#$白勺镕月 02 O’ 目1运#镕目i自目 g面 4面 秸0 500 s50 ※月$≈1白#女T※#m自镕*#镕* Fi91Calcelatt衄StrUCture 计算思路是，对回波信号进行bit-reverse排序 后顺序的分为8块，每块1024个点。先使用FFr 计算Al和A2块的结果，利用蝶形运算结构，计算 。 一 ：8 II级B1块中用于ⅡI级cl块计算所需的1024个点； 同理使用FFT计算A3和A4块的结果，然后计算 目2 M#a#口#镕*№较 一 2 FiB oftlwtwomethods II级B2块中用于III级cl块计算所需的1024个点； R∞uk—padson 最后利用B1和B2块中计算出的1024个点，采用由图2可