基于Matlab模拟常见信号时频分析与模糊函数图

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法       神经网络预测       雷达通信       无线传感器        电力系统

信号处理              图像处理               路径规划       元胞自动机        无人机 

⛄ 内容介绍

时频分析的主要研究对象是非平稳或时变信号,它能清楚地揭示信号的时变频谱特性,在信号处理中有着非常重要的作用.根据线性调频相参脉冲信号模糊函数理论特点,对大时宽相参雷达脉冲序列模糊函数及投影图进行了仿真计算.

⛄ 部分代码

%常见雷达信号的时频分析

%LFM线性调频信号：

%参数设置

clc

close all

B=20e6;  %带宽20MHz

T=10e-6;  %脉宽10us

K=B/T;   %调频斜率

f0=30e6; %中频30MHz

A=10;    %信号幅值

NFFT=512;%FFT点数

fs=4*(f0);  %奈奎斯特采样率

Ts=1/fs;

N=T/Ts;

phi0=0;

%时域和频域波形分析

t=-T/2:Ts:T/2;

S2=A*cos(2*pi*(f0*t+K*t.^2/2.0)+phi0);

figure(1)

subplot(211)

plot(t*1e6, S2, 'b', 'LineWidth', 1.5);

axis on

grid on

xlabel('t/us');

ylabel('S(t)');

axis([2*(-T/2*1e6) 2*(T/2*1e6) -3*A 3*A])

title('LFM线性调频信号时域');

Y2=fft(S2,NFFT);

Y2=Y2/NFFT;

Y22=fftshift(Y2);

Amp=abs(Y22).^2;

F2=linspace(-0.5,0.5,NFFT)*fs;    %双边频率

subplot(212)

plot(F2, Amp,'b', 'LineWidth', 1.5);   %双边幅度

axis([-0.5*fs 0.5*fs 0 1.5*max(Amp)])

axis on

grid on

xlabel('frequence/Hz');

ylabel('Magnitude');

title('LFM线性调频信号频域');

%信号模糊函数分析

ff = hilbert(S2');

[naf,tau,xi]=ambifunb(ff);

figure(2)

subplot(211)

subplot(211)

surf(2*tau,xi,abs(naf).^2,'edgecolor','none')

xlabel('Delay'); 

ylabel('Doppler');

shading interp

title('Narrow-band ambiguity function');

subplot(212)

contour(tau,xi,abs(naf));

xlabel('Delay'); 

ylabel('Doppler');

shading interp

title('模糊函数俯视图');

⛄ 运行结果

⛄ 参考文献

[1] 姚帅, 刘昱含, 蒋宇轩,等. 一种基于区域核函数的多分量线性调频信号时频分析方法:, CN202110634684.0[P]. 2021.

[2] 石珂. 基于Matlab的线性调频信号的时频分析[J].  2006.

[3] 孙健. 基于MATLAB的语音信号时频分析研究[J]. 电脑知识与技术:学术交流, 2015.

[4] 郭徽东. 线性调频相参脉冲雷达信号模糊函数仿真分析[J]. 舰船电子对抗, 2021, 044(006):56-59,120.

⛳️ 代码获取关注我

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料