跳频信号参数估计

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% project: 跳频信号参数估计
% standard:MSE
% methods:差分
% date:2014/6/3
% note:能达到很好的性能要求
clear,clc,close all;
warning('off');
% s = RandStream('mt19937ar','seed',2);
% RandStream.setGlobalStream(s);
Fs = 10e6;      % 采样频率
Ts = 1/Fs;          % 采样周期
Fh = 2000;          % 跳频速率
Th = 1/Fh;          % 跳频周期
Ht = Th/2;          % 跳时
SNR = -5:5:10;           % 信噪比
delta_f = 25e4;
Nu = 8;
f_set = [25 100 125 75 175 150 200 50]*10^4;        % 跳频频率集
t = 0:1/Fs:Nu*Th-1/Fs;                               % 时间
Monte = 100;
% 产生跳频信号，每个跳频周期对就一个频率
% 每个跳频周期的时间点为Th/Ts
for i = 1:Nu
    sig(1+(i-1)*Th/Ts:i*Th/Ts) = sqrt(2)*cos(2*pi*f_set(i)*t(1+(i-1)*Th/Ts:i*Th/Ts));
end
t = t(1:end-Ht*Fs);
%% 参数估计
hopping_rate = zeros(length(SNR),Monte);        % 跳速
ft = zeros(length(SNR),Monte);                  % 跳频周期