[发明专利]一种短群延时数字滤波方法

说明书

技术领域

本发明涉及一种数字滤波方法，具体涉及一种短群延时数字滤波方法，适用于核磁共振仪器以及其它需要进行数字滤波的信号处理系统。

背景技术

现代核磁共振仪器系统为了消除模拟滤波器所带来的失真，普遍采用了基于过采样的数字抽取滤波技术。现有的数字抽取滤波技术主要由CIC(Cascaded Integrator-Comb)数字滤波器与FIR(Finite Impulse Response)数字滤波器级联而成。CIC滤波器由于不需要乘法操作，具有硬件资源消耗小、执行效率高、群延时小的优点，但是其不平坦的通带幅频响应使得谱峰的相对幅度严重失真；而FIR滤波器具有可变的幅频特性和严格的线性相频特性，现代核磁共振仪器通常在CIC滤波器后级联高阶FIR作为补偿滤波器，得到了较好的通带平坦度、带外抑制度和过渡带陡峭度。然而，采用FIR作为补偿滤波器也存在一定的不足之处，主要是FIR滤波器带来较长的群延时。为了同时得到好的通带平坦度、带外抑制度和较窄的过渡带，FIR滤波器需要较大的阶数，阶数越大，群延时越大，信号向后推迟的程度就越大，核磁信号有效数据点的损失也就越大，从而一定程度上造成了信号的失真，这种失真在信号采集点数较少时则更加明显。

针对核磁共振仪器等系统中数字滤波器带来的群延时问题，科研工作者们提出了一些改进方法。Jon G.Wurl等(美国专利“DIGITAL FILTER PRE-CHARGING”US005652518)基于改变系统初始条件来优化滤波器瞬态响应的思想，提出核磁仪器接收部分在核磁信号前面根据第一个点的值增加一组平滑数据，经过数字滤波后，可实时的得到一组不含群延时的核磁共振FID(Free Induction Decay)信号。这种方法对核磁信号的第一个点的依赖性较大，而核磁信号的第一个点比较容易出现偏差，故该方法在核磁信号第一个点出现偏差时会产生谱图的基线扭曲；Petr SADOVSKY等从理论上对这种改变系统初始条件来优化滤波器瞬态响应的方法进行了分析与探讨(OPTIMISATION OF THE TRANSIENT RESPONSE OF A DIGITAL FILTER.Radioengineering 2000,9(2))，但第一个点出现偏差时谱图基线扭曲的问题并没有得到解决。

本专利为改善以上数字滤波器群延时带来的相关问题，提出一种短群延时数字滤波方法。具体内容为：先对过采样数据采用CIC实时数字滤波得到FID数字信号，然后在频域进行CIC的频域补偿并反傅立叶变换得到补偿后的FID信号。从而在不损失数字滤波器性能的前提下，满足了核磁共振实验的需求。

发明内容

本发明的目的是针对现有技术存在的问题，提出的一种短群延时数字滤波方法。

为了实现上述目的，本发明采用了以下技术措施：

一种短群延时数字滤波方法，包含以下步骤：

步骤1、核磁共振仪器载入设定的脉冲序列和采样参数，并按照设定的过采样频率对FID信号进行采样；

步骤2、将步骤1过采样得到的数字信号通过CIC数字滤波器进行CIC数字滤波；

步骤3、根据步骤2中CIC数字滤波器的幅频响应曲线H_cic(f)，计算CIC数字滤波器的频域幅度补偿曲线H_comp(f)，将步骤2中CIC数字滤波后的数据进行傅立叶变换得到谱数据S_cic，并将谱数据S_cic与频域幅度补偿曲线H_comp(f)逐点相乘，得到补偿后的谱数据S_comp；

步骤4、记步骤3中CIC补偿后的谱数据S_comp长度为L，截取S_comp居中的L/P个数据点，P>1，记为谱数据S_1/P；

步骤5、对步骤4中得到的谱数据S_1/P进行反傅立叶变换，得到经过CIC补偿后的FID信号。

如上所述的步骤3包括以下步骤：

步骤3.1、计算CIC数字滤波器的总幅频响应H_cic(f)：

其中，m表示级联CIC数字滤波器的个数，m>＝1，M_m，R_m，N_m，f_m分别表示第m级CIC数字滤波器的延时因子、抽取率、阶数、频率与主瓣宽度的归一化值；