fftplot()函数实现，在matlab中查看离散信号频率

ZDA2022

已于 2023-08-01 16:15:42 修改

阅读量1.2k

点赞数

文章标签： matlab 开发语言

于 2023-03-23 15:26:09 首次发布

本文链接：https://blog.csdn.net/qq_43686057/article/details/129730658

版权

本文展示了如何使用快速傅里叶变换(FFT)和短时傅里叶变换(STFT)来分析信号的频率内容。通过示例代码，解释了如何处理实信号和复信号，特别是针对线性调频波形(LFM)进行频谱分析。同时，讨论了不同参数对信号观察的影响，包括窗口长度、重叠数量和FFT长度等。

摘要由CSDN通过智能技术生成

利用fft观察频率

%实信号
fs=50;%采样频率50Hz
f1=15;
f2=20;
t = 0:1/Fs:10-1/Fs; % 0-9.98s 一共500个点
x = sin(2*pi*f1*t) + sin(2*pi*f2*t);%原始信号
fftplot(x,fs);

%%
%复信号
fs = 4e9;
bw = 20e6;
pw = 2e-6;
pri = 200e-6;

LFM = phased.LinearFMWaveform('SampleRate',fs,...
    'SweepBandwidth',200e6,...
    'PulseWidth',pw,'PRF',1/pri,'FrequencyOffset',1.2e9);
x= LFM();
fftplot(x,fs);


function fftplot(signal, fs)

% signal: 离散信号 % time: 离散时间

N = length(signal); 
X = abs(fftshift(fft(signal))); 
f = linspace(-fs/2, fs/2, N);

figure; 
plot(f,X); 
xlabel('Frequency (Hz)'); 
ylabel('Amplitude'); 
title('FFT plot'); 
grid on;

end

利用stft观察频率

对于一般的只是向看一下频谱特点时则只需要传入离散信号与采样率即可，stft函数会自动绘图。对于其他有细节操作限制时，则需要另外规定。

% 创建LinearFMWaveform对象
LFM = phased.LinearFMWaveform('SampleRate', 4e9, 'PulseWidth', 2e-6, ...
    'PRF', 1/4e-6, 'OutputFormat', 'Pulses', ...
    'SweepBandwidth', 200e6, 'SweepDirection', 'Up', 'FrequencyOffset', 1.2e9);
signal = LFM();
figure();
plot(abs(signal));
figure();

%仅传入信号与采样率
stft(signal,LFM.SampleRate);

对图像三维可视化

对信号操作有特殊要求

% 创建LinearFMWaveform对象
LFM = phased.LinearFMWaveform('SampleRate', 4e9, 'PulseWidth', 2e-6, ...
    'PRF', 1/4e-6, 'OutputFormat', 'Pulses', ...
    'SweepBandwidth', 200e6, 'SweepDirection', 'Up', 'FrequencyOffset', 1.2e9);

% 发送信号
signal = LFM();

window_len = round(length(signal)/10); % 窗口长度
num_overlap = round(window_len/2); % 窗口重叠数量
nfft = 2^nextpow2(window_len); % FFT长度
stft(signal, LFM.SampleRate, 'Window', hamming(window_len), ...
    'OverlapLength', num_overlap, 'FFTLength', nfft);