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

具体内容：利用FFT分析连续正弦信号的频谱：改变时间截断长度研究其对频谱泄漏的影响，并研究采用不同窗函数对信号进行时域截断时，对频谱泄漏的影响；

信号周期：0.04s；

时间长度分别为0.08s、0.044s；

0.08s为周期的整数倍，不会泄露；0.044s不是周期的整数倍，会发生泄露；

%a时间截断长度对泄露的影响
clc;
clear;
close all;

TT=0.08;%总时长
Freq=25;%信号频率 0.04s为周期
fs=1000;%采样频率
Ts=1/fs;%时间间隔

t1=0:Ts:(TT-Ts);
N1=length(t1);

x1=sin(2*pi*Freq*t1);
y1=fft(x1)/N1;
Z1=abs(y1);%振幅
Z2=angle(y1);%相位

for i=1:N1
    if Z1(i)<1E-10;
        Z2(i)=0;
    end
end

k1=-N1/2:(N1/2-1);
f1=k1*fs/N1;

figure(1);
subplot(3,1,1);
plot(t1,x1);
title(sprintf('信号x(t)=sin(%0.0fpit);采样频率：%0.0fHz',Freq*2,fs));
subplot(3,1,2);
stem(f1,fftshift(Z1));
subplot(3,1,3);
stem(f1,fftshift(Z2));

TT=0.044;%总时长
Freq=25;%信号频率 0.04s为周期
Ts=0.001;%时间间隔
fs=1000;%采样频率

t1=0:Ts:(TT-Ts);
N1=length(t1);

x1=sin(2*pi*Freq*t1);
y1=fft(x1)/N1;
Z1=abs(y1);%振幅
Z2=angle(y1);%相位

k1=-N1/2:(N1/2-1);
f1=k1*fs/N1;

figure(2);
subplot(3,1,1);
plot(t1,x1);
title(sprintf('信号x(t)=sin(%0.0fpit);采样频率：%0.0fHz',Freq*2,fs));
subplot(3,1,2);
stem(f1,fftshift(Z1));
subplot(3,1,3);
stem(f1,fftshift(Z2));

时间截断长度为0.08s，不泄露： 

 

 时间截断长度为0.044s，泄露： 

从上到下：信号时域图，信号幅频图，信号相位图