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('加窗后的信号频谱图')
msgbox('MATLAB编程答疑,请加QQ: 1530497909','MATLAB答疑','help')
web http://url.cn/TKcdXk -browser