【matlab】信号频谱分析之时间截断长度对信号泄露的影响

具体内容;利用FFT，分析连续周期方波信号的频谱：通过改变采样率与时间截断长度,分别研究频谱的与泄漏;

信号周期：0.2s

时间截断长度分别分：1s（5个周期），0.5s（2.5个周期）

泄露;整数个周期不会发生泄露，非整数个周期会发生泄露

原信号：方波信号

%b泄露
clc;
clear;
close all;
fs=100;
ts=1/fs;
Freq=5;
TT=1;%=0.5
t1=0:ts:(TT-ts);
x1=(square(2*pi*Freq*t1)+1)/2;
N=length(t1);
if mod(N,2)==0
    f=(-N/2:N/2-1)*fs/N;
else
    f=(-(N-1)/2:(N-1)/2)*fs/N;
end
y1=fft(x1)/N;
Z1=abs(y1);
figure(1);
subplot(2,1,1);
plot(t1,x1);
ylim([-0.5,1.5]);
title('方波信号 5个周期');
subplot(2,1,2);
stem(f,fftshift(Z1));
axis([-50,50,0,1]);

%%%%%%%%%%%%%%%%
fs=100;
ts=1/fs;
TT=0.5;%=0.5
t1=0:ts:(TT-ts);
x1=(square(2*pi*Freq*t1)+1)/2;
N=length(t1);
if mod(N,2)==0
    f=(-N/2:N/2-1)*fs/N;
else
    f=(-(N-1)/2:(N-1)/2)*fs/N;
end
y1=fft(x1)/N;
Z1=abs(y1);
figure(2);
subplot(2,1,1);
plot(t1,x1);
ylim([-0.5,1.5]);
title('方波信号2.5个周期');
subplot(2,1,2);
stem(f,fftshift(Z1));
axis([-50,50,0,1]);

 

 

 

 

上：时域图

下：幅频图