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

权利要求书

1.一种短群延时数字滤波方法，其特征在于，包含以下步骤：

步骤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信号。

2.根据权利要求1所述的一种短群延时数字滤波方法，其特征在于，所述的步骤3包括以下步骤：

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

Hcic(f)=|sin(πM1f1)sin(πf1/R1)|N1*|sin(πM2f2)sin(πf2/R2)|N2*......*|sin(πMmfm)sin(πfm/Rm)|Nm]]>

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

步骤3.2、根据步骤3.1获得的CIC数字滤波器的总幅频响应H_cic(f)计算频域幅度补偿曲线H_comp(f)：

Hcomp(f)=1Hcic(f);]]>

步骤3.3、将步骤2中CIC数字滤波得到的数据进行傅立叶变换得到谱数据S_cic；

步骤3.4、将步骤3.3中的谱数据S_cic与步骤3.2中的频域幅度补偿曲线H_comp(f)逐点相乘，得到补偿后的谱数据S_comp。