发明名称 基于时域重排的超声导波频散补偿与多模式分离方法
摘要 本发明属于超声技术领域,具体为一种基于时域重排的超声导波频散补偿与多模式分离算法。本发明主要包括三个部分,第一部分为超声导波频散的正向仿真算法,在已知波导介质理论频散曲线的情况下,可求解不同频率成分的时延参数,从而获得任意激励下的多模式超声导波时域信号;第二部分为逆向频散补偿算法,可利用频散关系对单一模式导波各频率成分实现不同程度的时延修正,从而将原持续时间较长的频散模式能量压缩为瞬时冲击,使导波波形简单,易于分析。第三部分为基于逆向频散补偿算法,可实现超声导波的多模式分离。
申请公布号 CN101982943B 申请公布日期 2013.08.21
申请号 CN201010502661.6 申请日期 2010.10.11
申请人 复旦大学 发明人 他得安;许凯亮;王威琪
分类号 H04B11/00(2006.01)I;H04L25/02(2006.01)I;A61B8/00(2006.01)I;G06F19/00(2011.01)I 主分类号 H04B11/00(2006.01)I
代理机构 上海正旦专利代理有限公司 31200 代理人 陆飞;盛志范
主权项 1.一种基于时域重排的超声导波频散补偿与多模式分离方法,其特征在于具体步骤为:(1)计算群延迟:对某一次激励,由式(1)计算某导波模式各谐波成份所对应群延迟:<img file="358929DEST_PATH_IMAGE001.GIF" wi="181" he="64" />(1)其中,<img file="807228DEST_PATH_IMAGE002.GIF" wi="23" he="36" />为导波激励时刻,d为发射换能器和接收换能器距离,<img file="586965DEST_PATH_IMAGE003.GIF" wi="56" he="35" />为某一导波模式群速度与频率的关系;(2)求频散传递函数:根据群延迟,求得对应导波模式的频散传递函数,其表达式为:<img file="982174DEST_PATH_IMAGE004.GIF" wi="397" he="71" />(2)其中信号<img file="999809DEST_PATH_IMAGE005.GIF" wi="67" he="36" />为积分初值;(3)超声导波频散的正向仿真:根据频散传递函数,求解激励信号通过此频散系统的输出,即得某一模式对应该激励下的仿真信号,若激励信号频率为<img file="822271DEST_PATH_IMAGE006.GIF" wi="47" he="26" />,则接收信号频谱<img file="89305DEST_PATH_IMAGE007.GIF" wi="54" he="30" />算式为:<img file="943997DEST_PATH_IMAGE008.GIF" wi="205" he="30" />(3)(4)多导波模式叠加:对不同导波模式用步骤(1)-(3)的方法分别求解,并叠加得多模式信号;(5)求频散传递函数反函数:为实现对应模式的频散补偿,求解上述步骤(2)所得传递函数的<img file="612876DEST_PATH_IMAGE009.GIF" wi="42" he="22" />反函数,其表达式为:<img file="871819DEST_PATH_IMAGE010.GIF" wi="469" he="75" />(4)(6)导波模式频散补偿:用上式(4)对某一模式实现频散补偿,若频散信号频谱为<img file="626148DEST_PATH_IMAGE011.GIF" wi="44" he="22" />,则频散补偿算式为:<img file="363160DEST_PATH_IMAGE012.GIF" wi="217" he="33" />(5)(7)补偿模式能量分离:对频散补偿结果分离出补偿模式瞬时激励以及欠补偿模式的混合信号;能量分离的实现方法为:1)若激励可预先测定,直接去除激励即得剩余模式信号;2)若激励不可预先测定,分离出补偿结果在零时刻附近的所有信号能量算为补偿模式能量,剩余部分为欠补偿模式能量;(8)频散反补偿:对步骤(7)分离结果反补偿还原信号模式;反补偿具体算式为:<img file="293070DEST_PATH_IMAGE013.GIF" wi="230" he="33" />(6)其中<img file="722914DEST_PATH_IMAGE014.GIF" wi="46" he="22" />为步骤(7)分离结果,<img file="43168DEST_PATH_IMAGE015.GIF" wi="48" he="22" />为某一模式导波频散信号,<img file="318292DEST_PATH_IMAGE009.GIF" wi="42" he="22" />为对应导波模式的传递函数表达式;(9)对剩余导波模式重复步骤(6)-(8),至分离所有模式。
地址 200433 上海市杨浦区邯郸路220号