[发明专利]一种超声多普勒回波信号仿真方法

说明书

本发明公开了一种超声多普勒回波信号仿真方法，该仿真方法按如下步骤依次进行：步骤一、在程序上设置超声探头所发射的超声波为f(t)＝cos(2×π×f₀×t)，将K个周期内所采f(t)的值存入矩阵x中；步骤二、被测物体相对于超声探头以径向速度v移动，根据公式(3)计算得出T₂的值；步骤三、设在矩阵x中第i个周期内所需要抽取的点数为N_i，通过公式(4)计算获得N_si，并对N_si四舍五入得到正整数N_i；步骤四、在矩阵x中，针对第i个周期，在当前周期的采样点中平均抽取N_i个点，依次放入矩阵y中；步骤五、在一个新的坐标系里，将得到的矩阵y中的的P个f(t)值作为纵坐标，以P个t_x作为横坐标，得到P个坐标点,作为超声多普勒回波信号。本发明的优点：计算效率高、更直观。

技术领域

本发明涉及计算机信号仿真领域，尤其涉及的是一种超声多普勒回波信号仿真方法。

背景技术

自20世纪以来，多普勒技术在通信、医学等领域被越来越广泛的应用，超声多普勒信号分析方法的验证需要用到仿真的超声多普勒回波信号，现有仿真方法计算复杂，不能直观地真实的反应被测物体的运动状态。

发明内容

本发明的目的在于克服现有技术的不足，基于波阵面传播的物理概念，提供了一种计算效率高、更直观的超声多普勒回波信号仿真方法。

本发明是通过以下技术方案实现的：

一种超声多普勒回波信号仿真方法，所述仿真方法按如下步骤依次进行：

步骤一、在程序上设置超声探头所发射的超声波为f(t)＝cos(2×π×f₀×t)、超声波的中心频率为f₀、信号采样频率为f_s、采样时长为K个超声波的周期时长，其中，f(t)为超声波的振幅，则每个周期采样点数量为f_s/f₀，将K个周期内所采f(t)的值存入矩阵x中；

步骤二、被测物体相对于超声探头以径向速度v移动，且在每个发射超声波周期内，径向速度v是恒定的，设被测物体接收相邻两个超声波阵面的时间为T₁，设被测物体反射相邻两个超声波阵面时间为T₂，则有：

(c-v)×T₁＝λ＝c×T (1)

(c+v)×T₁＝λ₂＝c×T₂ (2)

其中，c为测试环境下的声速，T为发射超声波的周期时长，T＝1/f₀，λ为发射超声波的波长，λ₂为反射超声波的波长，

则由公式(1)和公式(2)得出：

根据公式(3)计算得出T₂的值；

步骤三、设在矩阵x中第i个周期内所需要抽取的点数为N_i，通过如下公式(4)计算获得N_si，

其中，i为正整数，且i＝(1,2，...K)，m为预设值且0m1；若N_si不是正整数，则对N_si四舍五入得到正整数N_i；