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