function?[y1y2y3]?=?FH_Parameter_Estimation（xF0Fshlengthdelta_f）
%?function:?跳頻信號(hào)參數(shù)估計(jì)
%?y1:?跳頻頻率
%?y2:?跳時(shí)
%?y3:?跳速
%?x:?接收信號(hào)
%?F0:?中頻
%?Fs:?采樣頻率
%?hlength:?短時(shí)傅里葉變換窗寬，須為奇數(shù)

%%?STFT
xx?=?x;
x?=?hilbert（x‘）;
H1?=?window（@hamminghlength）;
N?=?2^nextpow2（hlength）;
[TFR1TF]=tfrstft（x1:length（x）NH1）;
TFR?=?abs（TFR1（1:N/2:））;
F?=?F（1:N/2）;
figure
t?=?1:length（x）;
contour（tFabs（TFR））;
axis（[0?max（t）?0?0.5]）;
xlabel（‘時(shí)間?t/s‘）;
ylabel（‘歸一化頻率?f/Hz‘）;
title（[‘STFT變換時(shí)頻圖（組合窗）‘]）;
colormap?gray;
set（gcf‘color‘‘white‘）;
grid?on;
%%?跳速估計(jì)
?%----------時(shí)頻脊線----------%
?[vloc]?=?max（TFR）;
loc?=?loc./size（TFR1）;
figure
plot（loc）;
%----------差分處理----------%
y?=?abs（diff（loc））;
%?figure
%?plot（y）;
%----------去噪--------------%
y?=?（y?>=?delta_f/Fs/2）;
%?figure
%?stem（y）;
%-------截取信號(hào)y中第一跳和最后一跳間的信號(hào)-----%
[vc]?=?find（y?>?0）;
yy?=?y（c（1）:c（end））;
%?figure
%?stem（yy）;
%-------跳速估計(jì)-----------------%
num1?=?sum（yy?==?1）;????%?跳頻信號(hào)突變個(gè)數(shù)
num2?=?sum（yy?==?0）;????%?跳頻周期點(diǎn)數(shù)和???
time_point?=?（num2+num1-1）/（num1-1）;
ft?=?time_point./Fs;????????????%?跳頻周期
hopping_rate?=?1/ft;??????????????????%?跳速
%%?跳時(shí)估計(jì)
tm?=?sum（c）/length（c）-（length（c）-1）/2*ft*Fs;??%?從信號(hào)突變位置的均值計(jì)算跳時(shí)點(diǎn)
%?tm?=?c（1）;
if?rem（tmft*Fs）?==?0
????tm?=?0;
end
hop_timing?=?tm/Fs;
%%?跳頻頻率估計(jì)
k?=?hop_timing*Fs+1;
if?hop_timing?k?=?1;
end
number?=?1;?????????????%?跳頻頻率集中的頻率個(gè)數(shù)
while?number?????tmp?=?sum（TFR（:k:min（[length（x）?k+Fs/hopping_rate-1]））2）;??
????k?=?k+Fs/hopping_rate;
????%?value為一個(gè)跳頻周期內(nèi)最大的STFT值，loc為對(duì)應(yīng)的頻率點(diǎn)
????[v?loc1]?=?max（tmp）;
????fre_set（number）?=?loc1*Fs/N;
????number?=?number+1;
end
number?=?num1-1;
%%?CZT方法對(duì)跳頻頻率估計(jì)精細(xì)化
x_czt?=?xx（round（hop_timing*Fs）+1:c（end））;
for?k?=?1:number
????fw?=?delta_f;?f0?=?fre_set（k）-delta_f/2;
????delta_f1?=?10;????????????%?分辨率
????M?=?fw/delta_f1;??????????%?觀察點(diǎn)數(shù)
????W?=?exp（-j*2*pi*fw/（M-1）/Fs）;
????A?=?exp（j*2*pi*f0/Fs）;
????Y?=?abs（czt（x_cztMWA））;
????f?=?（（0:length（Y）-1）*fw/length（Y））+f0;
????[vloc]?=?max（Y）;
????fre_set（k）?=?loc*fw/length（Y）+f0;
end
y1?=?fre_set+Fs*floor（F0/Fs）;
y2?=?hop_timing;
y3?=?hopping_rate;
end

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????2334??2014-06-29?10:59??FH_Parameter_Estimation.m

?????文件???????5178??2014-06-26?09:38??Parameters_Estimation_Diff.m

?????文件????????799??2014-06-03?10:36??STFT.m

?????文件???????3242??2014-05-22?14:22??Time_Frequency_Analysis.m

?????文件????????163??2014-06-30?11:28??新建文本文檔.txt

-----------?---------??----------?-----??----

????????????????11716????????????????????5