clear;clc;close all
Fs=600;
N=600;
t=(0:N-1)/Fs; % 采样时间序列
x=sin(2*pi*50*t)+cos(2*pi*80*t)+0.6; %采样信号
%% 绘制原始采样信号
subplot(3,1,1)
plot(t,x);
title('原始信号')
N=length(x);
wn=hanning(N); %汉宁窗
x=x.*wn'; % 原始信号时域加窗
%% 绘制加窗后的信号
subplot(3,1,2)
plot(t,x);
title('加窗后的信号')
%% 对加窗信号做FFT
X=fft(x,N);
f=(0:N-1)*Fs/N;
%% 绘制频谱图
subplot(3,1,3)
plot(f,abs(X));
title('加窗后的信号频谱图')
- 1
- 2
- 3
- 4
- 5
- 6
前往页