CN101190137B - 用于超声成像系统的实时数字正交解调方法与装置 - Google Patents

用于超声成像系统的实时数字正交解调方法与装置 Download PDF

Info

Publication number
CN101190137B
CN101190137B CN2006101469801A CN200610146980A CN101190137B CN 101190137 B CN101190137 B CN 101190137B CN 2006101469801 A CN2006101469801 A CN 2006101469801A CN 200610146980 A CN200610146980 A CN 200610146980A CN 101190137 B CN101190137 B CN 101190137B
Authority
CN
China
Prior art keywords
real
frequency
quadrature demodulation
digital quadrature
value
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.)
Active
Application number
CN2006101469801A
Other languages
English (en)
Other versions
CN101190137A (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.)
Shenzhen Mindray Scientific Co Ltd
Original Assignee
Shenzhen Mindray Bio Medical Electronics Co Ltd
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 Shenzhen Mindray Bio Medical Electronics Co Ltd filed Critical Shenzhen Mindray Bio Medical Electronics Co Ltd
Priority to CN2006101469801A priority Critical patent/CN101190137B/zh
Priority to US11/617,959 priority patent/US7706476B2/en
Publication of CN101190137A publication Critical patent/CN101190137A/zh
Priority to US12/767,642 priority patent/US8179997B2/en
Application granted granted Critical
Publication of CN101190137B publication Critical patent/CN101190137B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979Combined Doppler and pulse-echo imaging systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52025Details of receivers for pulse systems
    • G01S7/52026Extracting wanted echo signals
    • G01S7/52028Extracting wanted echo signals using digital techniques
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03DDEMODULATION OR TRANSFERENCE OF MODULATION FROM ONE CARRIER TO ANOTHER
    • H03D3/00Demodulation of angle-, frequency- or phase- modulated oscillations
    • H03D3/007Demodulation of angle-, frequency- or phase- modulated oscillations by converting the oscillations into two quadrature related signals

Abstract

本发明公开了一种用于超声成像系统的实时数字正交解调方法与装置。所述方法除包括乘法步骤和滤波步骤外,还包括正、余弦表生成步骤,用于实时生成正、余弦表,以及滤波参数生成步骤,用于实时生成相应的滤波参数以对来自乘法器的信号进行滤波。所述装置包括两个乘法器,两个滤波器,正、余弦表生成模块,滤波参数生成模块,以及两个参数存储器。按照本发明的用于超声成像系统的实时数字正交解调方法与装置,有效地节约了存储器资源,也更便于控制。

Description

用于超声成像系统的实时数字正交解调方法与装置
技术领域
本发明涉及一种数字正交解调方法与装置,特别是涉及一种用于超声成像系统的实时数字正交解调方法与装置。
背景技术
数字正交解调装置是全数字超声成像系统,特别是全数字彩色多普勒超声成像仪中不可缺少的部分。
图1为一个典型的超声成像系统(发射部分由于和本发明无太大关系,未画出)。一个常规的B型成像流程为:探头发射出脉冲,参与接收的各阵元接收到回波,经过放大,A/D转换,在波束合成器中以不同的延时量相加得到射频数据。射频数据(即下面所述的x(t))进入正交解调环节后分为I,Q两路,根据成像模式的不同进入B信号处理模块、Color(或者Colorflow)信号处理模块或者Doppler信号处理模块。处理完毕后,经过DSC(数字扫描转换)变成便于理解的图像形式显示在屏幕上。主控CPU负责每个模块的参数更新。
图2为一个超声成像系统中常见的正交解调装置,图中的I1、Q1以及I、Q与下面的公式(1)和(2)的意义相同。波束合成器输出的RF信号被同步分到两个乘法器和正弦表及余弦表相乘。正弦表和余弦表数值由正弦表存储器和余弦表存储器中查表得出。乘法器的输出进入低通滤波器。低通滤波器根据进入信号所代表的深度,从滤波参数存储器中取出对应滤波参数。用该滤波参数对输入信号滤波,输出即正交解调结果。
医疗设备中的超声波具有以下特点:中心频率沿深度会发生变化,信号带宽沿深度会发生变化。与超声波的这种特点相对应,数字正交解调装置要求支持解调信号的中心频率可变,带宽可变。一般数字超声设备均是用硬件(FPGA)实现的数字正交解调,要满足中心频率沿深度(或者随时间)变化,带宽沿深度(或者随时间)变化,需储存大量参数,增大了储存器开销。
假定超声回波信号的表达式为:x(t)=A(t)cos(ω(t)*t+Φ(t)),式中,A(t)是一低频信号,代表回波沿不同深度的幅度变化,ω(t)代表回波的频率,随时间也发生变化,正交解调的过程分以下两步:
1) I 1 ( t ) = x ( t ) × cos ( ω ( t ) × t ) = A ( t ) 2 ( cos ( Φ ( t ) ) + cos ( 2 ω ( t ) × t + Φ ( t ) ) ) - - - ( 1 )
Q 1 ( t ) = x ( t ) × sin ( ω ( t ) × t ) = - A ( t ) 2 ( sin ( Φ ( t ) ) - sin ( 2 ω ( t ) × t + Φ ( t ) ) )
2)从上式I1(t)和Q1(t)的表达式可以看出I1(t)和Q1(t)均是由两部分信号组成,频率在0附近的低频信号和频率在2ω(t)附近的高频信号,将I1(t)和Q1(t)送入一个低通滤波器,滤去高频部分,假定滤波器的单位冲激响应为h(t),即:
I ( t ) = I 1 ( t ) ⊗ h ( t ) = A ( t ) / 2 × cos ( Φ ( t ) )
Q ( t ) = Q 1 ( t ) ⊗ h ( t ) = A ( t ) / 2 × sin ( Φ ( t ) ) - - - ( 2 )
I、Q即为正交解调的结果。即正交解调之后,信号被分成了相互正交的两路I和Q,I和Q代表的意思是同相(In-Phase)和正交(Quadrature),两路信号主要保留了原信号的低频部分。
对I和Q求模,得到的结果显然是A(t)/2,即原信号的幅度信息(包络),这是B型成像的基础,对I和Q进行别的处理,可以得到诊断对象内的血流信息。
在全数字超声设备中,正交解调的实现一般是用数字处理方法,将叙述上述原理时使用的公式换成数字表达即可。从上面叙述可以看出,正交解调涉及参数包括正弦表,余弦表,低通滤波参数。一般系统中均是将这些参数保存在存储器,这些参数要考虑到解调频率可变,带宽可变,参数一般比较多。
假定一个系统A D采样率为40M,要求探测深度为30cm,需要保存的正弦表和余弦表长度为约16000点,正交解调要求的量化精度一般比较高,假定为16bit,则正、余弦表的存储量为16000×2×16=512Kbit;用于解调的低通滤波器假定为100阶(在实际情况下可能需要更高),位宽12bit,每隔64个采样点切换一组滤波参数(以满足沿深度带宽可变),需要的存储量为(16000/64)×50×12=150Kbit。这样的设计方式一方面必须在系统中外加存储器芯片,增加成本,另外控制也比较麻烦。
发明内容
本发明的目的是提供一种用于超声成像系统的实时数字正交解调方法与装置,可大量节约存储器,同时满足超声成像系统的需要。
按照本发明的第一方面,提供一种用于超声成像系统的实时数字正交解调方法,包括:乘法步骤,用于将同步分到两个乘法器的信号分别与正弦表和余弦表数值相乘;滤波步骤,用于根据来自乘法器的信号所代表的深度,利用相应的滤波参数对来自乘法器的信号进行滤波,以输出正交解调结果;其中还包括:正、余弦表生成步骤,用于实时生成正、余弦表;以及滤波参数生成步骤,用于实时生成相应的滤波参数以对来自乘法器的信号进行滤波。
按照本发明第一方面的实时数字正交解调方法,其中所述正、余弦表生成步骤包括:计数器产生从存储器中读取数据的读地址;将从存储器中读取的数据与累加器中的数据相加,结果作为cordic模块的输入角度;以及由cordic模块产生正、余弦值。
按照本发明第一方面的实时数字正交解调方法,其中所述滤波参数生成步骤包括:用当前系数索引n对存储有Wei_Win(n)值的第一存储器进行查表,以得到窗函数相关的值Ha;用Sin查找表地址对存储有正弦查找表的第二存储器查表,得到正弦值Sn;以及用Ha和正弦值Sn相乘,得到当前索引对应的滤波参数。其中生成所述查找表地址包括步骤:用当前系数索引n变换后的值Tn=2n-N+1与滤波器截止频率k相乘,得到查表地址,其中N为滤波器的阶数,n=0,1,...,N-1。
可选地的是,其中所述截止频率k是固定的,对2pi/M量化后的值即为k,其中M为2的整次幂。
还可选地是,其中所述截止频率k是随时间变化的,生成随时间变化的截止频率k包括步骤:用计数器的值对存储有截止频率变化量的第三存储器查表,得到当前截止频率变化量Ac;以及用Ac的值和保存在寄存器Reg中的上一个滤波器截止频率KD-1相加,得到当前的截止频率KD
按照本发明的第二方面,提供一种用于超声成像系统的实时数字正交解调装置,包括:两个乘法器,用于将信号分别与正弦表和余弦表数值相乘;两个滤波器,用于根据来自所述乘法器的信号所代表的深度,利用相应的滤波参数对来自所述乘法器的信号进行滤波,以输出正交解调结果;其中还包括:正、余弦表生成模块,用于实时生成正、余弦表;滤波参数生成模块,用于实时生成相应的滤波参数以对来自所述乘法器的信号进行滤波;以及两个参数存储器,分别与正、余弦表生成模块和滤波参数生成模块相连接,用于分别存储实时生成正、余弦表和实时生成相应的滤波参数所需的输入参数。
按照本发明第二方面的实时数字正交解调装置,其中所述正、余弦表生成模块包括:存储器,用于存储角度数据;第一计数器,用于产生从存储器中读取数据的读地址;第一加法器,用于将从所述存储器读取的数据与累加器中的数据相加;以及CORDIC模块,用于将相加的结果作为输入角度来产生正、余弦值。
按照本发明第二方面的实时数字正交解调装置,其中所述滤波参数生成模块包括:第一存储器,用于存储Wei_Win(n)的值,以便用当前系数索引n对所述第一存储器查表从而得到窗函数相关的值Ha;Sin查找表地址生成模块,用于生成Sin查找表地址;第二存储器,用于存储正弦查找表,以便用Sin查找表地址对所述第二存储器查表从而得到正弦值Sn;以及第一乘法器,用于将Ha和正弦值Sn相乘,得到当前索引对应的滤波参数。
其中所述Sin查找表地址生成模块进一步包括:变换器,用于将当前系数索引n变换为Tn=2n-N+1,其中N为滤波器的阶数,n=0,1,...,N-1;以及第二乘法器,用于将Tn值与滤波器截止频率k相乘,以得Sin查找表地址。其中所述截止频率k是固定的,对2pi/M量化后的值即为k,其中M为2的整次幂。
可选地是,按照本发明第二方面的实时数字正交解调装置还包括:截止频率控制模块,用于生成随时间变化的截止频率k;其中所述截止频率控制模块进一步包括:第三存储器,用于存储截止频率的变化量;第二计数器,用于对所述第三存储器查表,得到当前截止频率变化量Ac;寄存器,用于保存上一个滤波器截止频率KD-1;以及第二加法器,用于将所述Ac值与所述上一个滤波器截止频率KD-1相加,以得到当前的截止频率KD
优选地是,所述计数器是均匀计数器;所述滤波器为低通滤波器;所述存储器、第一存储器和第三存储器为RAM;所述第二存储器为ROM。
采用本发明技术方案的用于超声成像系统的实时数字正交解调方法与装置,有效地节约了存储器资源,也更便于控制。
附图说明
图1为典型的超声成像系统结构框图;
图2为常见的正交解调装置结构框图;
图3为按照本发明的产生正、余弦表的流程图;
图4为按照本发明的滤波参数实时计算的流程图;
图5为按照本发明的产生随时间变化的截止频率的流程图;
图6为按照本发明的本发明的正交解调装置的结构框图;
图7为按照本发明的产生正、余弦表的结构框图;
图8为按照本发明的截止频率固定的滤波参数实时计算的结构框图;
图9为按照本发明的截止频率可变的滤波参数实时计算的硬件结构框图。
具体实施方式
1.实时数字正交解调方法
在本实施例中,用于超声成像系统的实时数字正交解调方法,包括:乘法步骤,滤波步骤,正、余弦表生成步骤,以及滤波参数生成步骤。下面具体说明实时生成正、余弦表的方法以及实时生成滤波参数的方法。
(1)正、余弦表实时生成方法
假定余弦表和正弦表用数组表示为cosLut[i],sinLut[i],其中i=0~L-1,L代表需要的最大长度,与最大可能的信号(也就是RF信号)长度相同。如果是40M采样率,30cm探测深度,则L约为16000。正、余弦表用硬件实时生成的原理为:
ω0=0
令                                (3)
ωi+1=ωi+Δθk
cos Lut[i]=cos(ωi)
则                                (4)
sin Lut[i]=sin(ωi)
用数字电路常用的Cordic(坐标旋转数字计算机)设计可以同时算出cos(ωi)及sin(ωi)。上式中,Δθk为供正、余弦表生成模块使用的相关参数。如果需要产生的是一个代表单频的正、余弦表,则Δθk是一个固定数,改变Δθk可以产生不同频率单频正、余弦表;如果需要产生的是频率随时间变化的正、余弦表,则只要求让Δθk的值和i(即时间)相关,每隔一段时间切换一次Δθk(一个典型实施例中,每64个采样点切换一次Δθk),这样产生出的正、余弦表频率将随时间变化。
图3为产生正、余弦表的流程图。计数器产生从存储Δθk的存储器中读取数据的读地址(步骤302);用所产生的读地址对存储器查表(步骤304);将读取的angle_ram数据与累加器中的数据相加(步骤306);将相加后的结果进行存储以供下次使用(步骤308),另外将相加后的结果作为cordic模块的输入角度,由cordic模块产生正、余弦值(步骤310)。采用这样的正、余弦表硬件实时生成方法,存储量可以降到直接存储正、余弦表的1/2D(D表示每D点换一次Δθk)。如果是用逻辑实现,只需要使用逻辑器件的内置RAM即可。
(2)滤波参数实时生成方法
假定低通抽样滤波器的阶数为N阶(N为偶数),滤波参数共有N个,设其为h(n),用窗函数法计算滤波参数的公式为:
h ( n ) = sin [ ( n - N - 1 2 ) ω c ] · Win ( n ) pi ( n - ( N - 1 ) / 2 ) , n = 0,1 . . . . . . . . N - 1 - - - ( 5 )
上式中,ωc为低通滤波器的截止频率(也就是截止角频率),Win(n)为窗函数,常见的有hamming,hanning,高斯窗等。假如保证ωc的取值为k×2pi/M,k=0~M-1,M取2的整次幂,如256,上式乘号左边的项可以通过查找一个长度为M的正弦查找表(保存由0到pi均匀分成M个角度的正弦值)得到。
假定ωc的取值为k×2pi/M(k相当于截止频率对2pi/M量化后的值),
sin [ ( n - N - 1 2 ) ω c ] = sin [ ( 2 n - N + 1 ) × k × pi / M ] - - - ( 6 )
实际查表时,首先计算(2n-N+1)*k,结果无论正负均保留低log2(M)位,将这log2(M)位中的低log2(M)-1位作为查表地址。最高位如果是0,直接将对应地址的数取出即是要求的正弦值;如果是1,需要将对应地址的数取相反数才是要求的正弦值。这样的设计利用了正弦函数的奇对称性。如果保存的值包含了0~2pi所有角度的值,也可以不这样做。
式(5)乘号右边的项记作Wei_Win(n),其中n=0,1,...,N-1,可以由软件提前算出存在RAM中,也可通过查表得到不同的n对应的值。主要通过两次查表和一次乘法就可以得到滤波器的一个参数。
图4为滤波参数实时计算的流程图。用一个存储器(ram)存放上面所述的Wei_Win(n),用另一个存储器(rom)存放正弦查找表。用当前系数索引n查表存储器ram得到索引n对应窗函数值Ha(即式(5)中的
Figure G061E6980120061212D000081
)(步骤402);用当前系数索引n变换(步骤404)后的值Tn(Tn=2n-N+1)和量化后的截止频率k相乘(步骤406),结果即是查表地址;用该查表地址对存储器rom进行查表(步骤408),得到正弦值Sn;用Ha和正弦值Sn相乘(步骤410),便得到当前索引n对应的系数Coe。
在上面的滤波参数生成过程中,滤波器截止频率是固定的(对2pi/M量化后值为k)。如上所述,超声设备中的正交解调装置要求带宽可变,即滤波器的截止频率ωc可变。可以把不同深度对应的ωc存在存储器中,滤波参数计算模块在对不同深度信号解调时,从存储器取出对应的ωc计算参数送入滤波器模块,这样就实现了带宽可变的正交解调装置。滤波参数实时计算硬件需要存储Wei_Win(n)和不同深度的ωc,相对于不同深度对应的N/2个滤波参数总量,这些存储量几乎可以忽略不计。
图5为产生随时间变化的截止频率的流程图。用一个存储器(Accum ram)保存截止频率的变化量(对于截止频率变化不快的情况,保存变化量比直接保存截止频率值节约储存空间),用一个寄存器保存上一个滤波器截止频率。使用均匀计数器产生的值(步骤502)对Accum ram查表(步骤504),得到当前截止频率变化量Ac;用Ac值和保存在寄存器Reg中的上一个滤波器截止频率KD-1(步骤506)相加(步骤508),得到当前的截止频率KD。使用该随时间变化的截止频率KD来生成滤波参数的方法与图4所示的完全一样,这样就实现了生成截止频率可以随时间变化的滤波参数。
本发明提供的低通滤波参数实时计算方法可以方便地扩展成计算带通、高通滤波器的实时计算方法。
2.正交解调装置
图6为本发明的正交解调装置结构框图。按照本发明的正交解调装置包括:两个乘法器602和604,用于将RF信号分别与正弦表和余弦表数值相乘;两个滤波器606和608,用于根据来自乘法器602和604的信号所代表的深度,利用相应的滤波参数对来自所述乘法器的信号进行滤波,以输出正交解调结果;正、余弦表生成模块610,用于实时生成正、余弦表;滤波参数生成模块612,用于实时生成相应的滤波参数以供滤波器对来自所述乘法器的信号进行滤波;以及两个参数存储器614和616,分别与正、余弦表生成模块610和滤波参数生成模块612相接,用于分别存储实时生成正、余弦表和实时生成相应的滤波参数所需的输入参数。在本发明的正交解调装置中,不需要正弦表存储器、余弦表存储器及低通滤波参数存储器,而用专用的正、余弦表生成模块及滤波参数生成模块(具体到一个实施例中,用数字逻辑电路实现)。这些模块也需要一些输入参数,放置在图6所示的参数存储器614和616中。由于输入的参数量很小,不需要专用的外置存储器。
下面详细说明本发明正交解调装置中的产生正、余弦表的结构和产生滤波参数的结构。
图7为产生正、余弦表的结构框图。正、余弦表生成模块包括:存储器702,第一计数器704,第一加法器706,累加器708,以及CORDIC模块710。第一计数器704产生从存储Δθk的存储器702中读取数据的读地址,由第一加法器706将从存储器702中读取的angle_ram数据与累加器708中的数据相加,结果作为cordic模块710的输入角度,由cordic模块710产生正、余弦值。采用这样的正、余弦表硬件实时生成方法,存储量可以降到直接存储正、余弦表的1/2D(D表示每D点换一次Δθk)。如果是用逻辑实现,只需要使用逻辑器件的内置RAM即可。
图8为截止频率固定的滤波参数实时计算的结构框图。滤波参数生成模块包括:第一存储器(ram)802,Sin查找表地址生成模块804,第二存储器(rom)806,以及第一乘法器808。第一存储器(ram)802存放了上文提到的Wei_Win(n),第二存储器(rom)806存放正弦查找表。用当前系数索引n查表第一存储器(ram)802得到窗函数相关的值Ha;Sin查找表地址生成模块804对第二存储器(rom)806查表,得到正弦值Sn;由第一乘法器808将Ha和正弦值Sn相乘,便得到当前索引i对应的系数Coe。其中Sin查找表地址生成模块804包括:变换器810,以及第二乘法器812。Sin查找表地址生成模块804的一种实现方式为:用当前系数索引n变换后的值Tn(Tn=2n-N+1)与量化后的截止频率k相乘,结果即是查表地址。
图9为截止频率可变的滤波参数实时计算的结构框图。与截止频率固定的滤波参数实时计算硬件结构相比,增加了一个截止频率控制模块900,以产生随时间变化的截止频率。该模块900的一种实现方式主要包括:第二计数器902,第三存储器(Accum ram)904,寄存器(Reg)906,以及加法器908。第三存储器(Accum ram)904保存了截止频率的变化量(对于截止频率变化不快的情况,保存变化量比直接保存截止频率值节约储存空间)。第二计数器902为均匀计数器,用该计数器的值对第三存储器(Accum ram)904查表,得到当前截止频率变化量Ac。用Ac的值和保存在寄存器(Reg)906中的上一个滤波器截止频率KD-1相加,得到当前的截止频率KD。用该截止频率KD计算滤波参数的硬件结构的其余部分与图8所示的完全一样。
在本实施例中,计数器选用均匀计数器,滤波器为低通滤波器。
以上通过具体的实例对本发明进行了说明,但本发明并不限于这些具体的实例。在本发明的说明书和权利要求书中所使用的一些术语,如“第一”、“第二”、“第三”等等,仅仅是为了便于进行描述,而不具有任何限制意义。

Claims (19)

1.一种用于超声成像系统的实时数字正交解调方法,包括:
乘法步骤,用于将同步分到两个乘法器的信号分别与正弦表和余弦表数值相乘;
滤波步骤,用于根据来自乘法器的信号所代表的深度,利用相应的滤波参数对来自乘法器的信号进行滤波,以输出正交解调结果;
其特征在于,还包括:
正、余弦表生成步骤,用于实时生成正、余弦表数据;以及
滤波参数生成步骤,用于实时生成相应的滤波参数以供对来自乘法器的信号进行滤波。
2.如权利要求1所述的用于超声成像系统的实时数字正交解调方法,其特征在于,所述正、余弦表生成步骤包括:
计数器产生从存储器中读取数据的读地址;
将从存储器中读取的数据与累加器中的数据相加,结果作为cordic模块的输入角度;
由cordic模块产生正、余弦值。
3.如权利要求1所述的用于超声成像系统的实时数字正交解调方法,其特征在于,所述滤波参数生成步骤包括:
用当前系数索引n对存储有窗函数Win(n)相关的值Ha的第一存储器进行查表,以得到窗函数Win(n)相关的值Ha;其中n=1,2,...,N-1,N为滤波器阶数,
Figure FSB00000387591000011
用Sin查找表地址对存储有正弦查找表的第二存储器查表,得到正弦值Sn
将Ha和正弦值Sn相乘,得到当前索引对应的滤波参数。
4.如权利要求3所述的用于超声成像系统的实时数字正交解调方法,其特征在于,生成所述查找表地址包括步骤:
用当前系数索引n变换后的值Tn=2n-N+1与滤波器截止频率k相乘,得到查表地址,其中N为滤波器的阶数,n=0,1,...,N-1。
5.如权利要求4所述的用于超声成像系统的实时数字正交解调方法,其特征在于:
所述截止频率k是固定的,对2π/M量化后的值即为k,其中M为2的整次幂。
6.如权利要求4所述的用于超声成像系统的实时数字正交解调方法,其特征在于,所述截止频率k是随时间变化的,生成随时间变化的截止频率k包括步骤:
用计数器的值对存储有截止频率变化量的第三存储器查表,得到当前截止频率变化量Ac
用Ac的值和保存在寄存器Reg中的上一个滤波器截止频率KD-1相加,得到当前的截止频率KD
7.一种用于超声成像系统的实时数字正交解调装置,包括:
两个乘法器,用于将RF信号分别与正弦表和余弦表数值相乘;
两个滤波器,用于根据来自所述乘法器的信号所代表的深度,利用相应的滤波参数对来自所述乘法器的信号进行滤波,以输出正交解调结果;
其特征在于还包括:
正、余弦表生成模块,用于实时生成正、余弦表;
滤波参数生成模块,用于实时生成相应的滤波参数以对来自所述乘法器的信号进行滤波;以及
两个参数存储器,分别与正、余弦表生成模块和滤波参数生成模块相连接,用于分别存储实时生成正、余弦表和实时生成相应的滤波参数所需的输入参数。
8.如权利要求7所述的用于超声成像系统的实时数字正交解调装置,其特征在于,所述正、余弦表生成模块包括:
存储器,用于存储角度数据;
第一计数器,用于产生从存储角度数据的存储器中读取数据的读地址;
累加器;
第一加法器,用于将从所述存储角度数据的存储器读取的数据与所述累加器中的数据相加;以及
CORDIC模块,用于将相加的结果作为输入角度以产生正、余弦值。
9.如权利要求7所述的用于超声成像系统的实时数字正交解调装置,其特征在于,所述滤波参数生成模块包括:
第一存储器,用于存储窗函数Win(n)相关的值Ha,以便用当前系数索引n对所述第一存储器查表从而得到窗函数Win(n)相关的值Ha;其中n=1,2,...,N-1,N为滤波器阶数,
Figure FSB00000387591000031
S in查找表地址生成模块,用于生成Sin查找表地址;
第二存储器,用于存储正弦查找表,以便用Sin查找表地址对所述第二存储器查表从而得到正弦值Sn
第一乘法器,用于将Ha和正弦值Sn相乘,得到当前索引对应的滤波参数。
10.如权利要求9所述的用于超声成像系统的实时数字正交解调装置,其特征在于,所述Sin查找表地址生成模块包括:
变换器,用于将当前系数索引n变换为Tn=2n-N+1,其中N为滤波器的阶数,n=0,1,...,N-1;以及
第二乘法器,用于将Tn值与滤波器截止频率k相乘,以得到Sin查找表地址。
11.如权利要求10所述的用于超声成像系统的实时数字正交解调装置,其特征在于:
所述截止频率k是固定的,对2π/M量化后的值即为k,其中M为2的整次幂。
12.如权利要求10所述的用于超声成像系统的实时数字正交解调装置,其特征在于,还包括:
截止频率控制模块,用于生成随时间变化的截止频率k;
其中所述截止频率控制模块进一步包括:
第三存储器,用于存储截止频率的变化量;
第二计数器,用于对所述第三存储器查表,得到当前截止频率变化量Ac
寄存器,用于保存上一个滤波器截止频率KD-1
第二加法器,用于将所述Ac值与所述上一个滤波器截止频率KD-1相加,以得到当前的截止频率KD
13.如权利要求8所述的用于超声成像系统的实时数字正交解调装置,其特征在于:所述第一计数器是均匀计数器。
14.如权利要求12所述的用于超声成像系统的实时数字正交解调装置,其特征在于:所述第二计数器是均匀计数器。
15.如权利要求7至12中的任何一项所述的用于超声成像系统的实时数字正交解调装置,其特征在于:所述两个滤波器为低通滤波器。
16.如权利要求7所述的用于超声成像系统的实时数字正交解调装置,其特征在于:所述参数存储器为RAM。
17.如权利要求8所述的用于超声成像系统的实时数字正交解调装置,其特征在于:所述存储角度数据的存储器为RAM。
18.如权利要求9所述的用于超声成像系统的实时数字正交解调装置,其特征在于:所述第一存储器为RAM;所述第二存储器为ROM。
19.如权利要求12所述的用于超声成像系统的实时数字正交解调装置,其特征在于:所述第三存储器为RAM。
CN2006101469801A 2006-11-27 2006-11-27 用于超声成像系统的实时数字正交解调方法与装置 Active CN101190137B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN2006101469801A CN101190137B (zh) 2006-11-27 2006-11-27 用于超声成像系统的实时数字正交解调方法与装置
US11/617,959 US7706476B2 (en) 2006-11-27 2006-12-29 Real-time digital quadrature demodulation method and device for ultrasonic imaging system
US12/767,642 US8179997B2 (en) 2006-11-27 2010-04-26 Real-time digital quadrature demodulation method and device for ultrasonic imaging system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2006101469801A CN101190137B (zh) 2006-11-27 2006-11-27 用于超声成像系统的实时数字正交解调方法与装置

Publications (2)

Publication Number Publication Date
CN101190137A CN101190137A (zh) 2008-06-04
CN101190137B true CN101190137B (zh) 2011-06-22

Family

ID=39465005

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2006101469801A Active CN101190137B (zh) 2006-11-27 2006-11-27 用于超声成像系统的实时数字正交解调方法与装置

Country Status (2)

Country Link
US (2) US7706476B2 (zh)
CN (1) CN101190137B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7043082B2 (en) * 2000-01-06 2006-05-09 Canon Kabushiki Kaisha Demodulation and phase estimation of two-dimensional patterns
US20100004542A1 (en) * 2008-07-03 2010-01-07 Texas Instruments Incorporated System and method for ultrasound color doppler imaging
CN101912277B (zh) * 2010-08-06 2013-03-27 深圳市蓝韵实业有限公司 一种基于流水线设计的实时数字正交解调的方法与装置
US9184771B2 (en) 2011-10-12 2015-11-10 Optis Cellular Technology, Llc Digital down conversion and demodulation
US9264282B2 (en) * 2013-03-15 2016-02-16 Innophase, Inc. Polar receiver signal processing apparatus and methods
US10357180B2 (en) * 2014-01-16 2019-07-23 D.T.R. Dermal Therapy Research Inc. Health monitoring system
KR101702260B1 (ko) * 2016-01-25 2017-02-13 임채은 초음파 진단방법 및 장치
CN110488228B (zh) * 2019-07-11 2021-12-07 中国科学院电子学研究所 线性调频信号生成方法、装置及存储介质

Citations (3)

* 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
CN1440726A (zh) * 2002-10-01 2003-09-10 深圳迈瑞生物医疗电子股份有限公司 全数字超声频谱多普勒成像方法及装置
CN1731388A (zh) * 2004-12-30 2006-02-08 中国科学院长春光学精密机械与物理研究所 超声信号的数字信号处理方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3392028B2 (ja) * 1997-11-28 2003-03-31 株式会社ケンウッド 階層化伝送ディジタル復調器
CN2315567Y (zh) 1998-03-13 1999-04-21 清华大学 一种新型的超声多普勒血流测量仪
CN1712926A (zh) 2005-06-19 2005-12-28 中国海洋大学 单细胞藻流式分析微流控芯片
US7593459B1 (en) * 2005-09-14 2009-09-22 Redpine Signals, Inc. Wireless link simulator
US7511469B2 (en) * 2006-07-24 2009-03-31 California Institute Of Technology Phase measurement device using inphase and quadrature components for phase estimation

Patent Citations (3)

* 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
CN1440726A (zh) * 2002-10-01 2003-09-10 深圳迈瑞生物医疗电子股份有限公司 全数字超声频谱多普勒成像方法及装置
CN1731388A (zh) * 2004-12-30 2006-02-08 中国科学院长春光学精密机械与物理研究所 超声信号的数字信号处理方法

Also Published As

Publication number Publication date
US7706476B2 (en) 2010-04-27
US20100202567A1 (en) 2010-08-12
US8179997B2 (en) 2012-05-15
US20080126460A1 (en) 2008-05-29
CN101190137A (zh) 2008-06-04

Similar Documents

Publication Publication Date Title
CN101190137B (zh) 用于超声成像系统的实时数字正交解调方法与装置
US10371804B2 (en) Ultrasound signal processing circuitry and related apparatus and methods
EP3132281B1 (en) Ultrasonic imaging compression methods and apparatus
US9124244B2 (en) Method and device for generating a filter coefficient in real time
CN101112320B (zh) 波束合成的接收聚焦参数的实时计算方法及其装置
Ramírez et al. Fast RNS FPL-based communications receiver design and implementation
CN108614271A (zh) 一种带反馈校正的多通道超声波任意波形信号发生器
KR0156087B1 (ko) 다중채널 디지탈 수신 방법 및 장치와 초음파 진단 장치
CN102508251B (zh) 多波束图像声纳中扇形变换的快速实现方法
CN110488228B (zh) 线性调频信号生成方法、装置及存储介质
CN101912277B (zh) 一种基于流水线设计的实时数字正交解调的方法与装置
CN108784737B (zh) 一种超声成像的波束合成方法和装置
CN106772475A (zh) 一种基于fpga捕获微弱gps导航信号的电路及方法
US20210279298A1 (en) Fft engine having combined bit-reversal and memory transpose operations
CN102768358B (zh) 一种基于fpga的水下实时成像方法及系统
CN110633447B (zh) 基于fpga的球面距离定点计算方法及其计算装置
CN101690671A (zh) 一种多普勒超声成像系统动态解调装置
Camacho et al. A strict-time distributed architecture for digital beamforming of ultrasound signals
CN115438790A (zh) 量子态信息处理系统、量子测控系统、量子计算机
CN101461720A (zh) 基于频谱多普勒的运动速度测量范围的调整方法及装置
Leclère et al. Implementing super-efficient FFTs in Altera FPGAs
CN206387910U (zh) 基于fpga捕获微弱gps导航信号的电路
CN115173887B (zh) Chirp信号的生成方法、装置、终端及计算机可读存储介质
JP3944084B2 (ja) 超音波診断装置
CN108196358A (zh) Fpga电路板和望远镜数字运算板

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20180511

Address after: 518057 MINDRAY science and technology south twelve road MINDRAY high tech Industrial Park, Shenzhen, Guangdong

Co-patentee after: SHENZHEN MINDRAY SCIENTIFIC Co.,Ltd.

Patentee after: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd.

Address before: 518057 MINDRAY science and technology south twelve road MINDRAY high tech Industrial Park, Shenzhen, Guangdong

Patentee before: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd.

TR01 Transfer of patent right
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20080604

Assignee: Shenzhen Mindray Animal Medical Technology Co.,Ltd.

Assignor: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd.

Contract record no.: X2022440020009

Denomination of invention: Real-time digital quadrature demodulation method and device for ultrasonic imaging system

Granted publication date: 20110622

License type: Common License

Record date: 20220804

EE01 Entry into force of recordation of patent licensing contract