CN104758005A - 一种单通道乘法器解调的超声多普勒血流速度测量方法 - Google Patents
一种单通道乘法器解调的超声多普勒血流速度测量方法 Download PDFInfo
- Publication number
- CN104758005A CN104758005A CN201510169602.4A CN201510169602A CN104758005A CN 104758005 A CN104758005 A CN 104758005A CN 201510169602 A CN201510169602 A CN 201510169602A CN 104758005 A CN104758005 A CN 104758005A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- doppler
- channel multiplier
- single channel
- 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
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/06—Measuring blood flow
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Data Mining & Analysis (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Theoretical Computer Science (AREA)
- Surgery (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pure & Applied Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Computational Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种单通道乘法器解调的超声多普勒血流速度测量方法,现有的多普勒血流测量系统多采用一双通道乘法器为核心的正交解调方案对超声回波信号中的多普勒信号进行解调,测量电路比较复杂。在本申请发明中,采用简单的单通道乘法器,通过选取合适的正弦波解调信号和频率轴的坐标变换,就可以提取多普勒血流信号的频谱;这种方案中对电路要求低,可以有效降低电路的复杂度,提高数据的测量速度,并具有较高的测量精度。
Description
技术领域
本发明属于医学超声应用领域,利用单通道乘法器和频差概念,实现多普勒的人体血管中血流速度测量。
背景技术
当声波波源运动时,测量到的频率较波源实际发出的频率有一定的变化,波源运动方向及速度决定频移的大小,这种现象叫做多普勒效应。利用多普勒效应可以测量运动物体的速度,而超声多普勒血流测量就是利用多普勒效应测量血管中红细胞的运动速度。
现有的多普勒血流测量系统多采用正交解调的方式得出超声回波信号中和速度有关的多普勒信号,然后经过计算机处理提取血流的运动信息。为了实现正交解调,需要双通道或两路在电路性能上完全一致的模拟乘法器,以及两个幅度相同相位成九十度的频率为信号载波频率的正弦波信号,否则将严重影响多普勒信号的提取,降低血流速度测量精度。
现有的多普勒血流测量处理过程如下:首先频率为fc的超声波入射到运动的物体(如血液中的红细胞)上,由于多普勒效应,回波信号中含有频率为fc+fd的多普勒信号,当然由于红细胞的速度有一定的分布范围,而且每个红细胞的运动方向不同,fd不是一个数值,其分布有一定的范围,比如分布在-fd1≤fd≤fd2。
含有多普勒信号的回波信号由超声接收探头接收,经过放大后输入到正交解调模块,正交解调模块的调制输入信号是频率为fc的正弦和余弦波信号,经滤波后的解调输出为两路信号,一路称为正交分量,另一路称为同相分量。再用两路AD将正交分量和同相分量采样合并为一个复数函数,通过傅里叶变换及其他方法处理就可以得到多普勒频移fd及血流速度。
从傅里叶变换的性质上可以看出,上述变换实际上利用频移特性,把信号的频谱从fc平移到了原点附近。由于实信号的频谱是关于原点对称,无法反映小于fc的多普勒信号,因此必须采用复信号的处理方式。
发明内容
本发明是为避免上述现有技术所存在的不足之处,提供一种单通道乘法器解调的超声多普勒血流速度测量方法,旨在简化现有测量电路的同时提高测量精度和测量速度。
本发明解决技术问题,采用如下技术方案:
本发明单通道乘法器解调的超声多普勒血流速度测量方法,其特点在于:以信号发生器产生中心频率为fc的正弦波驱动超声发射探头向运动物体发射超声波信号,并被所述运动物体反射,获得含有多普勒频移、且频率为fc+fd的回波信号,其中fd为多普勒频移;所述回波信号被超声接收探头接收,然后作为乘法器的一个输入进入到单通道乘法器;由所述信号发生器产生的中心频率为fc的正弦波进入频率合成器,产生频率为fc1的正弦波,作为所述单通道乘法器的另一个输入进入单通道乘法器,其中fc-fc1≥fd,且fc>fc1;
所述单通道乘法器的输出信号中包含一个频率为fc+fc1+fd的高频分量、一个频率为fc-fc1-fd的低频分量和一个频率为fc-fc1的低频分量;所述单通道乘法器的输出信号通过低通滤波器去除频率为fc+fc1+fd的高频分量,剩余的低频分量送入中心频率为fc-fc1的陷波器,通过所述陷波器去除频率为fc-fc1的低频分量,所述陷波器的输出信号通过模数转换单元转换为数字信号;
数据处理单元对所述数字信号进行傅里叶变换得到所述数字信号的频谱,然后进行频率坐标轴坐标平移,把频率点fc-fc1移到坐标原点,即获得与所述运动物体的运动速度相对应的多普勒信号频谱,从所述多普勒信号频谱中获得fd。
为满足fc-fc1≥fd,由频率合成器产生的正弦波的频率fc1可按如下方式确定:
设仪器所能测量的最大速度为vmax,其方向为靠近发射者,c为超声波的速度,由多普勒原理可知,此时所能得到的最大频偏fdmax为:
fc1由下式确定:
fc1≤fc-fdmax。
在本申请中采用简单的单通道或一路乘法器代替现在所用的正交解调模块。乘法器的一个乘数输入端输入经放大后含有多普勒频移的回波信号,另一个乘数输入端输入简单的正弦信号,其频率为fc1,取值小于fc,在考虑到最大的多普勒频移后,乘法器的输出信号通过低通滤波器后的频谱是以fc-fc1为中心的位置,因为fc-fc1≥fd,这样以fc-fc1为中心的位置距离原点距离是fd,如果fd取最大速度对应的多普勒平移,这样整个多普勒平移都落在了正频率范围内。然后对此信号进行模数转换和傅里叶变换得到信号的频谱,再通过坐标轴平移就实现了多普勒频谱的检测和速度的测量。
因此,本发明的有益效果体现在:
1、由于本发明仅采用一路乘法器就可以实现多普勒信号的提取,不需要两路在电路性能上完全一致的模拟乘法器,以及两路幅度相同、相位成九十度的频率为信号载波频率的正弦波信号,不仅简化了电路,降低了多普勒信号解调对电路的要求,而且可以在一定程度上提高多普勒信号提取的精度。同时本发明在将模拟信号转换为数字信号的部分,使原来需要对两路信号进行模数转换降低了一路,进一步简化了电路。
2、由于本发明中的数据处理单元是在实数域对信号进行处理,不需要原来的复数处理,这样对信号做傅里叶变换的速度提高了一倍,因此本发明的方法可以提高超声多普勒血流速度测量速度。
附图说明
图1为本发明超声多普勒血流速度测量方法所用系统的结构框图;
图2为本发明中数据处理单元的数据处理过程示意图;
图中标号:1信号发生器、2超声发射探头、3运动物体、4超声接收探头、5单通道乘法器、6频率合成器、7低通滤波器、8陷波器、9模数转换单元、10数据处理单元。
具体实施方式
本发明单通道乘法器解调的超声多普勒血流速度测量方法为:以信号发生器1产生中心频率为fc的正弦波驱动超声发射探头2向运动物体3发射超声波信号,并被运动物体3反射,获得含有多普勒频移、且频率为fc+fd的回波信号,其中fd为多普勒频移;回波信号被超声接收探头4接收,然后作为乘法器的一个输入进入到单通道乘法器5;由信号发生器1产生的中心频率为fc的正弦波进入频率合成器6,产生频率为fc1的正弦波,作为单通道乘法器5的另一个输入进入单通道乘法器5,其中fc-fc1≥fd,且fc>fc1;
单通道乘法器5的输出信号中包含一个频率为fc+fc1+fd的高频分量、一个频率为fc-fc1-fd的低频分量和一个频率为fc-fc1的低频分量;单通道乘法器5的输出信号通过低通滤波器7去除频率为fc+fc1+fd的高频分量,剩余的低频分量送入中心频率为fc-fc1的陷波器8,通过陷波器8去除频率为fc-fc1的低频分量,陷波器8的输出信号通过模数转换单元9转换为数字信号;
数据处理单元10对数字信号进行傅里叶变换得到数字信号的频谱,然后进行频率坐标轴坐标平移,把频率点fc-fc1移到坐标原点,即获得与运动物体的运动速度相对应的多普勒信号频谱,从所述多普勒信号频谱中获得fd。
下面结合附图对本发明具体实施方式作一步说明。
本发明测量速度的具体实施方式是:
如图1所示,信号发生器1驱动超声发射探头2发射中心频率为fc、且持续一定时间的超声信号x(t):
x(t)=cos(2πfct) (1);
超声信号x(t)入射到例如血细胞的运动物体3并被反射形成回波信号r(t),根据多普勒效应,回波信号中含有多普勒频移fd,即:
超声接收探头4接收到回波信号并输入到单通道乘法器5进行处理,获得输出信号y(t):
y(t)=r(t)×cos(2πfc1t) (3);
(3)式表明单通道乘法器的两个乘数分别是r(t)和cos(2πfc1t),后者cos(2πfc1t)由频率合成器6产生。由于血流的方向不同,多普勒频移fd可正可负。对于(3)式,fc1取值得条件是:
fc1≤fc-fdmax (4);
其中fdmax为多普勒频偏的最大绝对值。
因为多普勒频移远小于fc,所以fc1接近fc。(3)式所乘结果将包含三个频率分量,一个是频率为fc+fc1+fd的高频分量,通过低通滤波器7就可以滤去。另两个分量是频率为fc-fc1-fd的低频分量和一个频率为fc-fc1的低频分量。显然,由于(4)式的关系,fc-fc1-fd1≥0,这说明通过(3)式处理,信号的频移分量全部在正频段,可以直接通过实函数的傅里叶变换就可以得出多普勒信号的频谱。
低通滤波器6提取的频率为fc-fc1-fd的低频分量和频率为fc-fc1的低频分量,送入陷波器8,陷波器的中心频率点为fc-fc1,其目的是滤除频率为fc-fc1的低频分量,因为这一部分实际是静止物体信号和缓慢变化的血管壁信号,而非多普勒血流信号。
陷波器的输出信号经过模数转换单元9转换为数字信号。
如图2所示,数据处理单元10对数字信号进行傅里叶变换得到信号的频谱,然后进行频率坐标轴坐标平移,把频率点fc-fc1移到坐标原点,即可得到和运动物体的运动速度对应的正负频率的多普勒信号频谱,在所述多普勒信号频谱中的最大值所对应的频率值,即是多普勒频移fd。
Claims (1)
1.一种单通道乘法器解调的超声多普勒血流速度测量方法,其特征在于:以信号发生器(1)产生中心频率为fc的正弦波驱动超声发射探头(2)向运动物体(3)发射超声波信号,并被所述运动物体(3)反射,获得含有多普勒频移、且频率为fc+fd的回波信号,其中fd为多普勒频移;所述回波信号被超声接收探头(4)接收,然后作为乘法器的一个输入进入到单通道乘法器(5);由所述信号发生器(1)产生的中心频率为fc的正弦波进入频率合成器(6),产生频率为fc1的正弦波,作为所述单通道乘法器(5)的另一个输入进入单通道乘法器(5),其中fc-fc1≥fd,且fc>fc1;
所述单通道乘法器(5)的输出信号中包含一个频率为fc+fc1+fd的高频分量、一个频率为fc-fc1-fd的低频分量和一个频率为fc-fc1的低频分量;所述单通道乘法器(5)的输出信号通过低通滤波器(7)去除频率为fc+fc1+fd的高频分量,剩余的低频分量送入中心频率为fc-fc1的陷波器(8),通过所述陷波器(8)去除频率为fc-fc1的低频分量,所述陷波器(8)的输出信号通过模数转换单元(9)转换为数字信号;
数据处理单元(10)对所述数字信号进行傅里叶变换得到所述数字信号的频谱,然后进行频率坐标轴坐标平移,把频率点fc-fc1移到坐标原点,即获得与所述运动物体的运动速度相对应的多普勒信号频谱,从所述多普勒信号频谱中获得fd。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510169602.4A CN104758005B (zh) | 2015-04-10 | 2015-04-10 | 一种单通道乘法器解调的超声多普勒血流速度测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510169602.4A CN104758005B (zh) | 2015-04-10 | 2015-04-10 | 一种单通道乘法器解调的超声多普勒血流速度测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104758005A true CN104758005A (zh) | 2015-07-08 |
CN104758005B CN104758005B (zh) | 2017-06-06 |
Family
ID=53640357
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510169602.4A Active CN104758005B (zh) | 2015-04-10 | 2015-04-10 | 一种单通道乘法器解调的超声多普勒血流速度测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104758005B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105574861A (zh) * | 2015-12-14 | 2016-05-11 | 上海交通大学 | 无标记的血流成像方法及系统 |
CN106236147A (zh) * | 2016-10-18 | 2016-12-21 | 湖南中聚光电科技发展有限公司 | 一种超声回波数字解调方法及装置 |
CN107370469A (zh) * | 2017-07-28 | 2017-11-21 | 电子科技大学 | 一种用于气侵检测的超声波接收电路 |
JP2018051162A (ja) * | 2016-09-30 | 2018-04-05 | 学校法人慶應義塾 | 生体信号検出システム、生体信号検出方法 |
CN115778434A (zh) * | 2022-12-02 | 2023-03-14 | 业成科技(成都)有限公司 | 讯号量测方法与都卜勒超音波系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4114156A (en) * | 1962-02-28 | 1978-09-12 | The United States Of America As Represented By The Secretary Of The Army | Frequency shift rate detector |
CN101416886A (zh) * | 2008-12-02 | 2009-04-29 | 深圳市蓝韵实业有限公司 | 一种连续多普勒信号处理装置和处理方法 |
CN102164544A (zh) * | 2008-09-23 | 2011-08-24 | 亨特来夫工业技术有限公司 | 胎儿心脏监视 |
JP2011217898A (ja) * | 2010-04-07 | 2011-11-04 | Hitachi Aloka Medical Ltd | 超音波診断装置 |
-
2015
- 2015-04-10 CN CN201510169602.4A patent/CN104758005B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4114156A (en) * | 1962-02-28 | 1978-09-12 | The United States Of America As Represented By The Secretary Of The Army | Frequency shift rate detector |
CN102164544A (zh) * | 2008-09-23 | 2011-08-24 | 亨特来夫工业技术有限公司 | 胎儿心脏监视 |
CN101416886A (zh) * | 2008-12-02 | 2009-04-29 | 深圳市蓝韵实业有限公司 | 一种连续多普勒信号处理装置和处理方法 |
JP2011217898A (ja) * | 2010-04-07 | 2011-11-04 | Hitachi Aloka Medical Ltd | 超音波診断装置 |
Non-Patent Citations (1)
Title |
---|
彭京思 等: "基于正弦调频连续波的超声多普勒血流速度的测量", 《生物医学工程研究》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105574861A (zh) * | 2015-12-14 | 2016-05-11 | 上海交通大学 | 无标记的血流成像方法及系统 |
CN105574861B (zh) * | 2015-12-14 | 2018-05-08 | 上海交通大学 | 无标记的血流成像方法及系统 |
JP2018051162A (ja) * | 2016-09-30 | 2018-04-05 | 学校法人慶應義塾 | 生体信号検出システム、生体信号検出方法 |
CN106236147A (zh) * | 2016-10-18 | 2016-12-21 | 湖南中聚光电科技发展有限公司 | 一种超声回波数字解调方法及装置 |
CN107370469A (zh) * | 2017-07-28 | 2017-11-21 | 电子科技大学 | 一种用于气侵检测的超声波接收电路 |
CN115778434A (zh) * | 2022-12-02 | 2023-03-14 | 业成科技(成都)有限公司 | 讯号量测方法与都卜勒超音波系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104758005B (zh) | 2017-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104758005A (zh) | 一种单通道乘法器解调的超声多普勒血流速度测量方法 | |
CN104490422B (zh) | 处理超声彩色血流成像的系统和方法 | |
CN106644030A (zh) | 一种基于多普勒雷达的非接触式振动测量方法 | |
US20120313813A1 (en) | Direct-to-Digital Software-Defined Radar | |
CN104316160B (zh) | 基于小波脊的水下声信号瞬时频率解调方法 | |
CN102508031A (zh) | 一种基于傅里叶级数的局部放电脉冲相角测量方法 | |
CN105212964B (zh) | 基于rf数据超声成像处理方法及系统 | |
CN101561501A (zh) | 基于dsp的雷达目标跟踪识别器 | |
CN102860843B (zh) | 一种胎心信号采集的方法及装置 | |
JP2006288974A (ja) | 超音波診断装置 | |
EP2248468B1 (en) | Ultrasound system and method for adaptively performing clutter filtering | |
JPS62117535A (ja) | 超音波ドプラ装置 | |
JP4369427B2 (ja) | ドプラ速度検出装置及びそれを用いた超音波診断装置 | |
CN101576618A (zh) | 基于小波变换的水声定位测量方法和测量系统 | |
KR100936173B1 (ko) | 표적 추적장치에서의 도플러 주파수를 추출하는 필터뱅크및 그 방법 | |
CN108195443B (zh) | 水位测量方法、系统及设备 | |
CN111750972A (zh) | 一种适用于双频连续激光雷达的目标振动测量方法 | |
CN111175724A (zh) | 一种基于质心计算的目标微多普勒分量提取方法 | |
CN201555950U (zh) | 基于dsp的雷达信号处理器 | |
CN105745550A (zh) | 模拟至信息转换器 | |
JPS6241644A (ja) | 超音波パルスドツプラ血流計 | |
CN205157648U (zh) | 一种无盲区数字相位计装置 | |
US20240063829A1 (en) | Smart fft measurement for reconfigurable sensor using a wideband digital receiver | |
JPH0318457B2 (zh) | ||
CN103134946A (zh) | 一种基于线性调频信号的空间介质流速测量方法与装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |