52matlab技术网站,matlab教程,matlab安装教程,matlab下载

 找回密码
 立即注册
搜索
热搜: 活动 交友 discuz
查看: 3576|回复: 0

强干扰下跳频信号的参数估计

[复制链接]

123

主题

207

帖子

2992

积分

版主

Rank: 7Rank: 7Rank: 7

积分
2992
发表于 2018-5-8 12:24:23 | 显示全部楼层 |阅读模式
close all
clear
clc


Fs = 10*10^3;      % 采样频率
Ts = 1/Fs;          % 采样周期
Fh = 20;           % 跳频频率
Th = 1/Fh;          % 跳频周期
SNR = 20;           % 信噪比
snr = 10^(SNR/10);  % 线性信噪比
f_set = [5 15 10 30 32 20 15 30 45 35 40 44]*10^2;        % 跳频频率集
noise_f = 24*10^2;                           %定频干扰频率
N = length(f_set);                                  % 跳频点数
t = 0:1/Fs:N*Th-1/Fs;                               % 时间
% 产生跳频信号,每个跳频周期对就一个频率
% 每个跳频周期的时间点为Th/Ts
for i = 1:N
    sig(1+(i-1)*fix(Th/Ts):i*fix(Th/Ts)) = sqrt(2)*cos(2*pi*f_set(i)*t(1+(i-1)*fix(Th/Ts):i*fix(Th/Ts)));
    noise_sig(1+(i-1)*fix(Th/Ts):i*fix(Th/Ts)) = sqrt(2)*cos(2*pi*noise_f*t(1+(i-1)*fix(Th/Ts):i*fix(Th/Ts)));%定频干扰
end

x = sig + noise_sig ;%加定频干扰
figure
subplot(211);
plot(t,x);
xlabel('时间 t/s');
ylabel('幅度');
title('跳频信号');
grid on;
subplot(212);
L = length(x);
NFFT = 2^nextpow2(L);
Y = abs(fft(x,NFFT));
f = (0:length(Y)-1)*Fs/length(Y);
plot(f(1:NFFT/2), 2*Y(1:NFFT/2));
xlabel('频率 f/Hz');
ylabel('|Y(f)|');
title('跳频信号频谱');
grid on;

x = x';
x = hilbert(x);
N = 2^nextpow2(length(x));
H = window(@hamming,1023);
[TFR,T,F]=tfrstft(x,1:length(x),N,H);


figure
contour(t,F(1:N/2),abs(TFR(1:N/2,: )));
axis([min(t), max(t), 0, 0.5])
ylabel('归一化频率 f/Hz');
xlabel('时间 t/s');
title('STFT变换时频图');
colormap
grid on;

f = F(1:N/2);
SQTFR = abs(TFR(1:N/2,: )).^2;

ff = sum(repmat(f,1,size(SQTFR,2)).*SQTFR,1)./sum(SQTFR,1);
figure
plot(t,ff)
xlabel('时间 t/ms');
ylabel('归一化频率 f/Hz');
title('时频重心曲线')

wf = wden(ff,'sqtwolog','s','mln',3,'haar');
figure
plot(t,wf)
xlabel('时间 t/ms');
ylabel('模值');
title('小波变换结果')



本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Archiver|手机版|小黑屋|52matlab技术网站 ( 粤ICP备14005920号-5 )

GMT+8, 2024-3-29 22:53 , Processed in 0.090471 second(s), 21 queries .

Powered by Discuz! X3.2 Licensed

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表