function?[fe?k0]=mrife（xfsN）
%頻率估計計算-----------------------------------------%
%?x:輸入數據
%?fs:數據采樣率
%?N:實際FFT的點數???N?=?length（x）;
%?fe:估計出的頻率
%?k0:最長譜線的索引值
%采用mrife算法-------------------------------------------------%
%rife算法（幅度比值法）-----------------------------------------%
%n=0:1:N-1;%x（n）
delta_f?=?fs/N;
n=0:N-1;
Xk?=?fft（x）;?????????????????????????????????????????%對x進行fft變換
P?=?abs（Xk）;?????????????????????????????????????????%Xk的絕對值
%subplot（211）???stem（nP）?;grid?on????%顯示xn的頻譜
[Pm?k1]?=?max（P）;????????????????????????????????????%第k1條譜線有最大值。
k0=k1-1;?????????????????????????????????????????????%其實際索引值應為k1-1
f0?=?delta_f*（k0）;???????????????????????????????????%最高頻率處的實際頻率值
if?k1==1||P（k1+1）?>=?P（k1-1）
????r?=?1;
else
????r?=?-1;
end
fe?=?（k0+r*P（k1+r）/（P（k1）+P（k1+r）））*fs/N;????????????%rife算法頻率估計值
fe1=fe;
%M-rife（改進rife算法）————————————————————————————————%
if?f