[发明专利]基于IpDFT的质子磁力仪自由感应衰减信号的频率估计方法

权利要求书

1.基于IpDFT的质子磁力仪自由感应衰减信号的频率估计方法，其特征在于：包括以下步骤：

步骤1：采集N点离散时间FID信号s(n)；

步骤2：获取信号s(n)的N点离散傅里叶变换序列，记做S(k)，其中因子因子是其共轭，因子λ＝e^-βe^jω，λ^*＝e^-βe^-jω是其共轭，；

步骤3：获取到频谱中k₁，k₂，k₃处三个值最大的频谱值S(k₁)、S(k₂)和S(k₃)，其中k₁，k₂，k₃的值可以通过粗估计获得，k₂＝k₀，是l₀整数部分，可以通过寻找信号频谱的谱峰所在位置，即得到k₂的值，k₁＝k₀-1，k₃＝k₀+1；

步骤4：根据插值傅里叶变换，利用S(k₁)、S(k₂)和S(k₃)三点频谱值得到信号频率和阻尼因子的估计值。

2.根据权利要求1所述的基于IpDFT的质子磁力仪自由感应衰减信号的频率估计方法，其特征在于：所述步骤4利用三点频谱值对信号频率和阻尼因子进行估计的方法具体包括以下步骤：

步骤41：因为有关系式(λ+λ^*)＝2e^-βcos(ω₀)，λλ^*＝e^-2β，所以解出λ的值就可以得到信号阻尼因子β和频率ω₀的估计值，令(λ+λ^*)＝x，λ·λ^*＝y，通过对S(k₁)、S(k₂)和S(k₃)表达式进行一系列的化简运算可以得到一元二次方程组

其中[]^*表示取共轭，公式(1)中各系数的值可以通过下面的公式求得

步骤42：求解公式(1)所示一元二次方程组，可以得到

其中代表x的估计值，代表y的估计值；

步骤43：根据关系式x＝2e^-βcos(ω₀)，y＝e^-2β，通过公式(3)最终可以得到信号频率ω₀和阻尼因子β的估计值，为

代表ω₀的估计值，代表β的估计值。