CN106236147A - 一种超声回波数字解调方法及装置 - Google Patents

一种超声回波数字解调方法及装置 Download PDF

Info

Publication number
CN106236147A
CN106236147A CN201610905267.4A CN201610905267A CN106236147A CN 106236147 A CN106236147 A CN 106236147A CN 201610905267 A CN201610905267 A CN 201610905267A CN 106236147 A CN106236147 A CN 106236147A
Authority
CN
China
Prior art keywords
signal
data
window function
ultrasound echo
function sequence
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
CN201610905267.4A
Other languages
English (en)
Other versions
CN106236147B (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.)
Optically Focused Electricity Development In Science And Technology Co Ltd In Hunan
Original Assignee
Optically Focused Electricity Development In Science And Technology Co Ltd In Hunan
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 Optically Focused Electricity Development In Science And Technology Co Ltd In Hunan filed Critical Optically Focused Electricity Development In Science And Technology Co Ltd In Hunan
Priority to CN201610905267.4A priority Critical patent/CN106236147B/zh
Publication of CN106236147A publication Critical patent/CN106236147A/zh
Application granted granted Critical
Publication of CN106236147B publication Critical patent/CN106236147B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/02—Measuring pulse or heart rate
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0866—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving foetal diagnosis; pre-natal or peri-natal diagnosis of the baby
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/488—Diagnostic techniques involving Doppler signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Gynecology & Obstetrics (AREA)
  • Pregnancy & Childbirth (AREA)
  • Cardiology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种超声回波数字解调方法及装置。该方法包括以下步骤：确定待解调的超声回波信号；采用预先确定的I信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据；采用预先确定的Q信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据。应用本发明实施例所提供的技术方案，采用I信号窗函数序列对超声回波信号进行加窗处理的过程或采用Q信号窗函数序列对超声回波信号进行加窗处理的过程只需要单个乘法器，对加窗之后的数据进行累加，也只需要少量加法器即可实现，节约了硬件成本，并使得在处理器芯片上的选型的自由度更大。

Description

一种超声回波数字解调方法及装置
技术领域
本发明涉及超声应用技术领域，特别是涉及一种超声回波数字解调方法及装置。
背景技术
在超声应用技术领域中，基于超声回波进行解调处理，得到被检测物体的组织幅度信息，是进行超声成像与多普勒信号提取的基础。目前，在多普勒胎心监测和脉冲多普勒计算中，对超声回波进行解调处理通常需经过五个步骤：步骤1：距离选通、步骤2：信号加窗、步骤3：调制频率信号相乘、步骤4：低通滤波和步骤5：信号累加之后，获得I(In-phase，同相)数据和Q(Quadrature，正交)数据，具体如图1所示，这也称为标准解调流程。
在图1中，对输入的信号进行步骤1距离选通处理，得到特定深度的超声回波信号。对该超声回波信号进行步骤2信号加窗处理，得到加窗处理信号，通过步骤2对超声回波信号进行加窗处理，使得加窗处理信号在时域上的截断不至于在频域上出现吉伯斯抖动。之后，在步骤3，将加窗处理信号分为两路，一路与载波同频率、同相位的三角函数信号相乘，另一路与相位偏差90°的三角函数信号相乘，分别得到第一步解调信号。每路第一步解调信号经过步骤4低通滤波之后，可以将第一步解调信号中的二倍频信号滤除，得到超声回波的同相信号(I信号)和正交信号(Q信号)。在超声脉冲模式下，将步骤4计算得到的I信号和Q信号分别进行累加，得到每个脉冲重复周期内的I数据和Q数据，该两个数据即为进行速度计算或者频谱计算的数据基础。
在步骤4和步骤5这两个过程中，步骤4的低通滤波过程一般需要采用多个乘法器，将特定长度的FIR滤波器参数卷积步骤3处理得到的第一步解调信号，然后将所有乘法器得到的数据进行相加，得到滤波数据，最后再将滤波数据进行累加，分别得到I信号和Q信号，具体过程参见图2所示。
在图2中，X输入为图1中步骤3处理之后的数据，即第一步解调信号，X(N)至X(1)为N个寄存器，N为FIR滤波器长度，FIR_(1)至FIR_(N)为FIR滤波器参数，Mult_1至Mult_N为N个乘法器，中间的Add器件为N阶加法器，可以将N个乘法器相乘得到的数据进行相加。至此，低通滤波器的运行流程结束。Y为单个寄存器，将右侧Add加法器得到的数据进行存储并再次作为加数送入该Add加法器，实现累加的功能，最终得到I信号或者Q信号。
由图1和图2可以看出，目前的信号解调过程虽然比较成熟，但是运行过程中涉及到较多的乘法器，占用了较多的硬件资源，将导致在处理器芯片的选型上自由度较小，需要在FPGA级别的芯片上才能进行实时处理运算，而且，FPGA芯片价格较高，导致硬件成本较高。
发明内容
本发明的目的是提供一种超声回波数字解调方法及装置，以提高在处理器芯片选型上的自由度，节省硬件成本。
为解决上述技术问题，本发明提供如下技术方案：
一种超声回波数字解调方法，包括：
确定待解调的超声回波信号；
采用预先确定的I信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据；
采用预先确定的Q信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据；
其中，所述I信号窗函数序列和所述Q信号窗函数序列分别为根据窗函数矩阵、三角函数矩阵、滤波器组和单位矩阵确定的。
在本发明的一种具体实施方式中，所述采用预先确定的I信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据，包括：
从所述超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至所述超声回波信号和预先确定的I信号窗函数序列全部送入运算完毕，获得I数据：
将所述超声回波信号中的元素和所述I信号窗函数序列中相应位置的元素送入乘法器；
使用I信号加法器将所述乘法器的数据与第一寄存器中的数据相加。
在本发明的一种具体实施方式中，所述采用预先确定的Q信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据，包括：
从所述超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至所述超声回波信号和预先确定的Q信号窗函数序列全部送入运算完毕，获得Q数据：
将所述超声回波信号中的元素和所述Q信号窗函数序列中相应位置的元素送入所述乘法器；
使用Q信号加法器将所述乘法器的数据与第二寄存器中的数据相加。
在本发明的一种具体实施方式中，所述I信号窗函数序列ZI为：
ZI＝G COS F E＝{zI(1),zI(2),zI(3),...,zI(n)}；
其中，所述G为窗函数矩阵，所述COS为三角函数矩阵，所述F为滤波器组，所述E为单位矩阵。
在本发明的一种具体实施方式中，所述Q信号窗函数序列ZQ为：
ZQ＝G SIN F E＝{zQ(1),zQ(2),zQ(3),...,zQ(n)}；
其中，所述SIN为与所述COS相位偏差90°的三角函数矩阵。
一种超声回波数字解调装置，包括：
超声回波信号确定模块，用于确定待解调的超声回波信号；
I数据获得模块，用于采用预先确定的I信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据；
Q数据获得模块，用于采用预先确定的Q信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据；
其中，所述I信号窗函数序列和所述Q信号窗函数序列分别为根据窗函数矩阵、三角函数矩阵、滤波器组和单位矩阵确定的。
在本发明的一种具体实施方式中，所述I数据获得模块，具体用于：
从所述超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至所述超声回波信号和预先确定的I信号窗函数序列全部送入运算完毕，获得I数据：
将所述超声回波信号中的元素和所述I信号窗函数序列中相应位置的元素送入乘法器；
使用I信号加法器将所述乘法器的数据与第一寄存器中的数据相加。
在本发明的一种具体实施方式中，所述Q数据获得模块，具体用于：
从所述超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至所述超声回波信号和预先确定的Q信号窗函数序列全部送入运算完毕，获得Q数据：
将所述超声回波信号中的元素和所述Q信号窗函数序列中相应位置的元素送入所述乘法器；
使用Q信号加法器将所述乘法器的数据与第二寄存器中的数据相加。
在本发明的一种具体实施方式中，所述I信号窗函数序列ZI为：
ZI＝G COS F E＝{zI(1),zI(2),zI(3),...,zI(n)}；
其中，所述G为窗函数矩阵，所述COS为三角函数矩阵，所述F为滤波器组，所述E为单位矩阵。
在本发明的一种具体实施方式中，所述Q信号窗函数序列ZQ为：
ZQ＝G SIN F E＝{zQ(1),zQ(2),zQ(3),...,zQ(n)}；
其中，所述SIN为与所述COS相位偏差90°的三角函数矩阵。
应用本发明实施例所提供的技术方案，确定待解调的超声回波信号后，采用预先确定的I信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据，同时采用预先确定的Q信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据。采用I信号窗函数序列对超声回波信号进行加窗处理的过程或采用Q信号窗函数序列对超声回波信号进行加窗处理的过程只需要单个乘法器，对加窗之后的数据进行累加，也只需要少量加法器即可实现，节约了硬件成本，并使得在处理器芯片上的选型的自由度更大。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为现有技术中信号解调处理流程的示意图；
图2为现有技术中FIR滤波和累加过程的示意图；
图3为本发明实施例中一种超声回波数字解调方法的实施流程图；
图4为本发明实施例中加窗累加流程的示意图；
图5为本发明实施例中超声回波数字解调方法对应的一种电路逻辑图；
图6为本发明实施例中一种超声回波数字解调装置的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面结合附图和具体实施方式对本发明作进一步的详细说明。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
参见图3所示，为本发明实施例所提供的一种超声回波数字解调方法的实施流程图，该方法可以包括以下步骤：
S110：确定待解调的超声回波信号。
在本发明实施例中，待解调的超声回波信号可以是经过距离选通处理后的信号或者为其他方式输入的特定信号。从图1中可知，将输入的信号进行距离选通处理后，可以得到特定深度的超声回波信号。
S120：采用预先确定的I信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据。
S130：采用预先确定的Q信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据。
其中，所述I信号窗函数序列和所述Q信号窗函数序列分别为根据窗函数矩阵、三角函数矩阵、滤波器组和单位矩阵确定的。
为便于描述，将上述两个步骤结合起来进行说明。
在步骤S110，确定待解调的超声回波信号后，可以采用预先确定的I信号窗函数序列对该超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据，同时，可以采用预先确定的Q信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据。具体参见图4所示。
在本发明实施例中，可以通过以下推导过程预先确定I信号窗函数序列和Q信号窗函数序列：
以图1所示的标准数字解调流程为推导基础，设经过步骤1距离选通处理之后的超声回波信号的数据矩阵为X，长度为n，设窗函数矩阵为G，g(n)为窗函数系数，经过步骤2信号加窗处理之后的数据P为：
P＝XG；
其中，
G = g ( 1 ) 0 0 0 .. 0 0 g ( 2 ) 0 0 .. 0 0 0 g ( 3 ) 0 .. 0 0 0 0 g ( 4 ) .. 0 .. .. .. .. .. .. .. .. .. .. .. .. 0 0 0 0 .. g ( n ) ( n , n ) .
以I信号计算为例，设COS为三角函数矩阵，经过图1所示步骤3的三角函数序列相乘，其矩阵相乘的结果D为：
D＝PCOS；
其中：
C O S = cos ( w t ) 0 0 0 .. 0 0 cos ( w ( t + Δ t ) ) 0 0 .. 0 0 0 cos ( w ( t + 2 Δ t ) ) 0 .. 0 0 0 0 cos ( w ( t + 3 Δ t ) ) .. 0 .. .. .. .. .. .. .. .. .. .. .. .. 0 0 0 0 .. cos ( w ( t + n Δ t ) ) ( n , n )
其中，f为载波信号频率，fsample为对信号的采样频率。
图1所示步骤4为低通滤波，设Y为低通滤波之后得到的序列，由y(1)至y(m)构成，F为滤波器组，其长度为n，由f(1)至f(n)组成，d(1)至d(m+l)组成D矩阵，则可以得到：
y(1)＝d(1)*f(1)+d(2)*f(2)+d(3)*f(3)+...+d(l)*f(l)；
y(2)＝d(2)*f(1)+d(3)*f(2)+d(4)*f(3)+...+d(l+1)*f(l)；
y(3)＝d(3)*f(1)+d(4)*f(2)+d(5)*f(3)+...+d(l+2)*f(l)；
……
y(m)＝d(m)*f(1)+d(m+1)*f(2)+d(m+2)*f(3)+...+d(m+l)*f(l)。
以如下矩阵方式表示F，其中m+l等于超声回波数据矩阵的长度n：
F = f ( 1 ) 0 0 0 .. 0 f ( 2 ) f ( 1 ) 0 0 .. 0 f ( 3 ) f ( 2 ) f ( 1 ) 0 .. 0 .. .. .. .. .. .. .. .. .. .. .. .. f ( l ) f ( l - 1 ) f ( l - 2 ) f ( l - 3 ) .. f ( 1 ) 0 f ( l ) f ( l - 1 ) f ( l - 2 ) .. f ( 2 ) 0 0 f ( l ) f ( l - 1 ) .. f ( 3 ) .. .. .. .. .. .. .. .. .. .. .. .. 0 0 0 0 .. f ( l - 1 ) 0 0 0 0 .. f ( l ) ( n , m ) ;
则：Y＝DF。
经过步骤5的累加过程，累加后的数据矩阵表示为：
I＝YE；
其中，单位矩阵E＝[1 1 1 1 1 1 .. 1]'(1,m)；
从以上分析和描述可以看出，整个解调过程可以表示为：
I＝X G COS F E。
根据矩阵运算的性质，可将I进一步表示为：
I＝XZI，
其中，ZI＝G COS F E。
同理，对于正交Q数据，同样有：
Q＝XZQ，
其中，ZQ＝G SIN F E。
从以上分析可以看出，整个I数据和Q数据的运算结果均可以看做为两个矩阵的相乘结果。在实际应用中，G、COS、SIN、F、E都是已知的。
所以，可以预先确定I信号窗函数序列ZI为：
ZI＝G COS F E＝{zI(1),zI(2),zI(3),...,zI(n)}；
其中，G为窗函数矩阵，COS为三角函数矩阵，F为滤波器组，E为单位矩阵，即根据窗函数矩阵、三角函数矩阵、滤波器组和单位矩阵可以确定I信号窗函数序列。
可以预先确定Q信号窗函数序列ZQ为：
ZQ＝G SIN F E＝{zQ(1),zQ(2),zQ(3),...,zQ(n)}；
其中，SIN为与COS相位偏差90°的三角函数矩阵，即根据窗函数矩阵、三角函数矩阵、滤波器组和单位矩阵可以确定Q信号窗函数序列。。
通过对图1所示的解调过程的分析，可将图1所示的解调过程简化为一个加窗累加的过程，相应的电路逻辑图如图5所示。
在该电路逻辑图中，包含一个乘法器Mult，两个加法器Add_I和Add_Q，两个寄存器Y1和Y2。其中，乘法器Mult分别与加法器Add_I和加法器Add_Q相连，乘法器Mult的数据作为加法器Add_I和加法器Add_Q的一个输入，寄存器Y1作为加法器Add_I的另一个输入，寄存器Y2作为加法器Add_Q的另一个输入，加法器Add_I的输出与寄存器Y1相连，加法器Add_Q的输出与寄存器Y2相连。
在本发明的一种具体实施方式中，步骤S120可以包括以下步骤：
从超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至超声回波信号和预先确定的I信号窗函数序列全部送入运算完毕，获得I数据：
将超声回波信号中的元素和I信号窗函数序列中相应位置的元素送入乘法器；
使用I信号加法器将乘法器的数据与第一寄存器中的数据相加。
在本发明的另一种具体实施方式中，步骤S130可以包括以下步骤：
从超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至超声回波信号和预先确定的Q信号窗函数序列全部送入运算完毕，获得Q数据：
将超声回波信号中的元素和Q信号窗函数序列中相应位置的元素送入乘法器；
使用Q信号加法器将乘法器的数据与第二寄存器中的数据相加。
为便于理解，以图5为例，将上述解调过程进行详细说明。
第一步：将超声回波信号的第一个元素x(1)和I信号窗函数序列的第一个元素zI(1)送入乘法器Mult。
第二步：将加法器Add_I打开，加法器Add_Q关闭，将乘法器Mult的数据与第一寄存器即寄存器Y1中的数据相加。
第三步：将Q信号窗函数序列的第一个元素zQ(1)送入乘法器，x(1)不变。
第四步：将加法器Add_Q打开，加法器Add_I关闭，将乘法器Mult的数据与第二寄存器即寄存器Y2中的数据相加。
第五步：将超声回波信号的第二个元素x(2)和I信号窗函数序列的第二个元素zI(2)送入乘法器。
第六步：将加法器Add_I打开，加法器Add_Q关闭，将乘法器Mult的数据与寄存器Y1中的数据相加。
第七步：将Q信号窗函数序列的第二个元素zQ(2)送入乘法器，x(2)不变。
第八步：将加法器Add_Q打开，加法器Add_I关闭，将乘法器Mult的数据与寄存器Y2中的数据相加。
……
按照上述方法循环送入数据，将所有的超声回波信号和窗函数序列全部送入运算完毕，得到最后的I数据和Q数据。
需要说明的是，本发明实施例对步骤S120和步骤S130的执行顺序不做限制。
应用本发明实施例所提供的方法，确定待解调的超声回波信号后，采用预先确定的I信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据，同时采用预先确定的Q信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据。采用I信号窗函数序列对超声回波信号进行加窗处理的过程或采用Q信号窗函数序列对超声回波信号进行加窗处理的过程只需要单个乘法器，对加窗之后的数据进行累加，也只需要少量加法器即可实现，节约了硬件成本，并使得在处理器芯片上的选型的自由度更大。
相应于上面的方法实施例，本发明实施例还提供了一种超声回波数字解调装置，下文描述的一种超声回波数字解调装置与上文描述的一种超声回波数字解调方法可相互对应参照。
参见图6所示，该装置包括以下模块：
超声回波信号确定模块210，用于确定待解调的超声回波信号；
I数据获得模块220，用于采用预先确定的I信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据；
Q数据获得模块230，用于采用预先确定的Q信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据；
其中，所述I信号窗函数序列和所述Q信号窗函数序列分别为根据窗函数矩阵、三角函数矩阵、滤波器组和单位矩阵确定的。
应用本发明实施例所提供的装置，确定待解调的超声回波信号后，采用预先确定的I信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据，同时采用预先确定的Q信号窗函数序列对超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据。采用I信号窗函数序列对超声回波信号进行加窗处理的过程或采用Q信号窗函数序列对超声回波信号进行加窗处理的过程只需要单个乘法器，对加窗之后的数据进行累加，也只需要少量加法器即可实现，节约了硬件成本，并使得在处理器芯片上的选型的自由度更大。
在本发明的一种具体实施方式中，I数据获得模块220，具体用于：
从超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至超声回波信号和预先确定的I信号窗函数序列全部送入运算完毕，获得I数据：
将超声回波信号中的元素和I信号窗函数序列中相应位置的元素送入乘法器；
使用I信号加法器将乘法器的数据与第一寄存器中的数据相加。
在本发明的一种具体实施方式中，Q数据获得模块230，具体用于：
从超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至超声回波信号和预先确定的Q信号窗函数序列全部送入运算完毕，获得Q数据：
将超声回波信号中的元素和Q信号窗函数序列中相应位置的元素送入乘法器；
使用Q信号加法器将乘法器的数据与第二寄存器中的数据相加。
在本发明的一种具体实施方式中，I信号窗函数序列ZI为：
ZI＝G COS F E＝{zI(1),zI(2),zI(3),...,zI(n)}；
其中，G为窗函数矩阵，COS为三角函数矩阵，F为滤波器组，E为单位矩阵。
在本发明的一种具体实施方式中，Q信号窗函数序列ZQ为：
ZQ＝G SIN F E＝{zQ(1),zQ(2),zQ(3),...,zQ(n)}；
其中，SIN为与COS相位偏差90°的三角函数矩阵。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其它实施例的不同之处，各个实施例之间相同或相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的技术方案及其核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。

Claims (10)

1.一种超声回波数字解调方法，其特征在于，包括：
确定待解调的超声回波信号；
采用预先确定的I信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据；
采用预先确定的Q信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据；
其中，所述I信号窗函数序列和所述Q信号窗函数序列分别为根据窗函数矩阵、三角函数矩阵、滤波器组和单位矩阵确定的。
2.根据权利要求1所述的超声回波数字解调方法，其特征在于，所述采用预先确定的I信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据，包括：
从所述超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至所述超声回波信号和预先确定的I信号窗函数序列全部送入运算完毕，获得I数据：
将所述超声回波信号中的元素和所述I信号窗函数序列中相应位置的元素送入乘法器；
使用I信号加法器将所述乘法器的数据与第一寄存器中的数据相加。
3.根据权利要求2所述的超声回波数字解调方法，其特征在于，所述采用预先确定的Q信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据，包括：
从所述超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至所述超声回波信号和预先确定的Q信号窗函数序列全部送入运算完毕，获得Q数据：
将所述超声回波信号中的元素和所述Q信号窗函数序列中相应位置的元素送入所述乘法器；
使用Q信号加法器将所述乘法器的数据与第二寄存器中的数据相加。
4.根据权利要求1至3任一项所述的超声回波数字解调方法，其特征在于，所述I信号窗函数序列ZI为：
ZI＝G COS F E＝{zI(1),zI(2),zI(3),...,zI(n)}；
其中，所述G为窗函数矩阵，所述COS为三角函数矩阵，所述F为滤波器组，所述E为单位矩阵。
5.根据权利要求4所述的超声回波数字解调方法，其特征在于，所述Q信号窗函数序列ZQ为：
ZQ＝G SIN F E＝{zQ(1),zQ(2),zQ(3),...,zQ(n)}；
其中，所述SIN为与所述COS相位偏差90°的三角函数矩阵。
6.一种超声回波数字解调装置，其特征在于，包括：
超声回波信号确定模块，用于确定待解调的超声回波信号；
I数据获得模块，用于采用预先确定的I信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第一数据进行累加，得到I数据；
Q数据获得模块，用于采用预先确定的Q信号窗函数序列对所述超声回波信号进行加窗处理，并将加窗之后的第二数据进行累加，得到Q数据；
其中，所述I信号窗函数序列和所述Q信号窗函数序列分别为根据窗函数矩阵、三角函数矩阵、滤波器组和单位矩阵确定的。
7.根据权利要求6所述的超声回波数字解调装置，其特征在于，所述I数据获得模块，具体用于：
从所述超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至所述超声回波信号和预先确定的I信号窗函数序列全部送入运算完毕，获得I数据：
将所述超声回波信号中的元素和所述I信号窗函数序列中相应位置的元素送入乘法器；
使用I信号加法器将所述乘法器的数据与第一寄存器中的数据相加。
8.根据权利要求7所述的超声回波数字解调装置，其特征在于，所述Q数据获得模块，具体用于：
从所述超声回波信号中的第一个元素开始，按顺序重复执行以下步骤，直至所述超声回波信号和预先确定的Q信号窗函数序列全部送入运算完毕，获得Q数据：
将所述超声回波信号中的元素和所述Q信号窗函数序列中相应位置的元素送入所述乘法器；
使用Q信号加法器将所述乘法器的数据与第二寄存器中的数据相加。
9.根据权利要求6至8任一项所述的超声回波数字解调装置，其特征在于，所述I信号窗函数序列ZI为：
ZI＝G COS F E＝{zI(1),zI(2),zI(3),...,zI(n)}；
其中，所述G为窗函数矩阵，所述COS为三角函数矩阵，所述F为滤波器组，所述E为单位矩阵。
10.根据权利要求9所述的超声回波数字解调装置，其特征在于，所述Q信号窗函数序列ZQ为：
ZQ＝G SIN F E＝{zQ(1),zQ(2),zQ(3),...,zQ(n)}；
其中，所述SIN为与所述COS相位偏差90°的三角函数矩阵。
CN201610905267.4A 2016-10-18 2016-10-18 一种超声回波数字解调方法及装置 Active CN106236147B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610905267.4A CN106236147B (zh) 2016-10-18 2016-10-18 一种超声回波数字解调方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610905267.4A CN106236147B (zh) 2016-10-18 2016-10-18 一种超声回波数字解调方法及装置

Publications (2)

Publication Number Publication Date
CN106236147A true CN106236147A (zh) 2016-12-21
CN106236147B CN106236147B (zh) 2017-10-13

Family

ID=57611754

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610905267.4A Active CN106236147B (zh) 2016-10-18 2016-10-18 一种超声回波数字解调方法及装置

Country Status (1)

Country Link
CN (1) CN106236147B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108042155A (zh) * 2017-12-20 2018-05-18 飞依诺科技(苏州)有限公司 超声回波信号自动时间频率解调方法及处理系统
CN109632970A (zh) * 2018-11-16 2019-04-16 中聚科技股份有限公司 一种用于编码超声的单通道超声脉冲发射与回波放大系统
WO2024027529A1 (zh) * 2022-08-05 2024-02-08 飞依诺科技股份有限公司 一种抑制杂波的超声波束合成方法、系统及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5482044A (en) * 1992-01-14 1996-01-09 Diasonics Ultrasound, Inc. Direct demodulation in ultrasound instruments
CN101912277A (zh) * 2010-08-06 2010-12-15 深圳市蓝韵实业有限公司 一种基于流水线设计的实时数字正交解调的方法与装置
CN201847706U (zh) * 2010-11-11 2011-06-01 深圳市信步科技有限公司 一种多普勒诊断仪解调电路
CN104758005A (zh) * 2015-04-10 2015-07-08 合肥工业大学 一种单通道乘法器解调的超声多普勒血流速度测量方法
CN105662464A (zh) * 2016-01-08 2016-06-15 深圳大学 一种超声波胎心监测仪及其回波信号的数字解调方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5482044A (en) * 1992-01-14 1996-01-09 Diasonics Ultrasound, Inc. Direct demodulation in ultrasound instruments
CN101912277A (zh) * 2010-08-06 2010-12-15 深圳市蓝韵实业有限公司 一种基于流水线设计的实时数字正交解调的方法与装置
CN201847706U (zh) * 2010-11-11 2011-06-01 深圳市信步科技有限公司 一种多普勒诊断仪解调电路
CN104758005A (zh) * 2015-04-10 2015-07-08 合肥工业大学 一种单通道乘法器解调的超声多普勒血流速度测量方法
CN105662464A (zh) * 2016-01-08 2016-06-15 深圳大学 一种超声波胎心监测仪及其回波信号的数字解调方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108042155A (zh) * 2017-12-20 2018-05-18 飞依诺科技(苏州)有限公司 超声回波信号自动时间频率解调方法及处理系统
CN109632970A (zh) * 2018-11-16 2019-04-16 中聚科技股份有限公司 一种用于编码超声的单通道超声脉冲发射与回波放大系统
WO2024027529A1 (zh) * 2022-08-05 2024-02-08 飞依诺科技股份有限公司 一种抑制杂波的超声波束合成方法、系统及存储介质

Also Published As

Publication number Publication date
CN106236147B (zh) 2017-10-13

Similar Documents

Publication Publication Date Title
CN106236147A (zh) 一种超声回波数字解调方法及装置
US11786216B2 (en) Ultrasound contrast enhanced imaging method and ultrasound imaging system
US8179997B2 (en) Real-time digital quadrature demodulation method and device for ultrasonic imaging system
CN101912277B (zh) 一种基于流水线设计的实时数字正交解调的方法与装置
CN105406836B (zh) 一种降采样解调滤波方法
CN107229075A (zh) 深度域地震子波的确定方法和装置
CN101741788B (zh) 最小频移键控信号的调制解调方法及专用的复解析带通滤波器设计方法
CN106384121A (zh) 基于标号空间匹配的标号随机集滤波器分布式融合方法
Pham et al. An adaptive computation of contour representations for mode decomposition
CN113435304A (zh) 提取扭振信号的扭振信息的方法、系统、设备及存储介质
Bailey et al. A digital matched filter for reverse time chaos
Li A class of negatively fractal dimensional Gaussian random functions
Kuzume et al. FPGA-based lifting wavelet processor for real-time signal detection
CN105581812B (zh) 用于超声成像设备的自动调节方法及系统
Maruo et al. A generalized harmonic analysis of ultrasound waves propagating in cancellous bone
CN105959245A (zh) 一种宽带捷变隐匿通信波形生成方法
CN114675253A (zh) 水面微幅波频率估计方法、装置、电子设备及存储介质
Pawula Dichotomous-noise-driven oscillators
CN110349088A (zh) 一种图像处理方法和系统
CN111413416A (zh) 基于本振信号的正交解调方法、装置、系统及超声系统
CN105266790A (zh) 一种心率测量的数字信号处理方法、系统及设备
CN111077566A (zh) 一种基于矩阵分解的双程波叠前深度偏移的方法
Abarbanel et al. Nonlinear analysis of time series data
Shoji et al. Interconnected Turing patterns in three dimensions
CN105212963A (zh) 实现超声彩色多普勒血流成像的壁滤波方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 410000 Hunan, Changsha high tech Development Zone, Yuelu Avenue West, No. 588, Changsha core city science and Technology Park, building 6, building 6

Applicant after: China Polytron Technologies Inc

Address before: 410000 Changsha City, Hunan province high tech Development Zone, Yuelu Avenue West, No. 588, Changsha core city science and Technology Park, building 6, building 6

Applicant before: Optically focused electricity development in science and technology company limited in Hunan

CB02 Change of applicant information
CB03 Change of inventor or designer information

Inventor after: Zhan Kai

Inventor after: Liang Yao

Inventor after: Zhang Ruiqin

Inventor after: Li Wenyu

Inventor after: Pan Youhua

Inventor after: Wang Hongbo

Inventor after: Liang Lei

Inventor before: Zhan Kai

Inventor before: Li Wenyu

Inventor before: Wang Hongbo

Inventor before: Pan Youhua

Inventor before: Liang Lei

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant