基于傅里叶变换法和Rife-Jane法计算信号的周期

1 傅里叶变换法

傅里叶变换法估计信号频率的原理是首先通过对信号进行傅里叶变换，然后找到频谱中幅值最大值所对应的频率，再求出其倒数即为信号的周期，使用该方法不需要对加速度信号进行预处理，但该方法容易受傅立叶变换时的频率分辨率的影响。

对加速度信号进行傅里叶变化，如下图所示。

傅里叶变换法求周期的程序可参考Matlab的库函数fft。

2 Rife-Jane方法

Rife算法是对傅里叶变换法的改进，是为了解决傅里叶变换法所存在的量化误差问题，其通过对频谱的最大谱线以及与其相邻的较大谱线进行插值，从而减小量化误差，提高了信号频率的估计精度。

对加速度信号进行傅里叶变化，如下图所示。

3 Matlab程序

%傅立叶变换法
N=756;
x=a;
fs=1/0.049*1.024;
delta_f=fs/N ; 
n=0:N-1; 
Xk=fft(x) ; 
P=abs(Xk) ; 
P(1)=0;
[~, k1] = max(P);
k0=k1-1 ; 
f0=delta_f * (k0);
x=fs/f0;
temp(1,it)=x;
if x==378
   temp(1,it)=314;
else
   temp(1,it)=296;
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Rife-Jane
N = 756;
x = a4;
fs=1/0.049*1.024;
delta_f = fs / N ; 
n = 0 : N-1 ; 
Xk = fft(x) ; 
P = abs(Xk) ; 
P(1) = 0;
[~, k1] = max(P);
k0 = k1 - 1 ; 
f0 = delta_f * (k0);
if P(k1 + 1) >= P(k1 - 1)
    r = 1;
else
    r = -1;
end
temp_f = (k0 + r * P(k1+r) / (P(k1) + P(k1 + r))) * fs / N;
tempp(1,it)=fs/temp_f;