DFT 正弦曲线的基频,直观

信息处理 自由度 窗函数 基频
2022-02-06 06:45:44

我试图直观地理解正弦信号的基频从时域到频域的映射,使用DFT公式。给定时域样本{x[0],x[1],,x[N1]}; 我们从DFT公式:

(1)X(fk)=n=0N1x[n]exp(j2πfkn)

为了fk=k/Nk=0,1,,N1.

为此,我在 Matlab 及其窗口版本中使用了正弦信号。这是我的 Matlab 代码的摘录:

Fs = 1000; 
Ts = 1/Fs;
t = 0:Ts:(1000 - 1)*Ts; %1sec

A = 2;
F = 20; 
y1 = A*sin(2*pi*F*t); 
y2 = y1.*hann(length(y1))';

nfft = 2^nextpow2(length(y1));
Fy = (Fs/2)*linspace(0, 1, nfft/2+1);

Y1 = 2*abs(fft(y1, nfft));
Y1(1) = .5*Y1(1);
Y1 = Y1./max(Y1);

Y2 = 2*abs(fft(y2, nfft));
Y2(1) = .5*Y2(1);
Y2 = Y2./max(Y2);

有关信号的时域和频域图,请参见下图

在此处输入图像描述

我有两个问题:

  1. 怎样才能直观地看出指数的线性组合(1)将产生一个具有绝对幅度的总和|X(fk)|在基频处有最大值F = 20 Hz? 或者这是绝对幅度谱的优化问题|X(fk)|写给fk? 或者这只是DFT?

  2. 非矩形窗口的信号逐渐变细如何解释(1)最大值的更好定位在F = 20 Hz(即基频附近的泄漏最小化) ? 我很了解窗户,泄漏的理论。我更多的是寻找基于方程式的解释(1). 从目视检查y1y2在上图中,如何从线性组合中得出结论(1), 形状的差异 |Y1(fk)||Y2(fk)|?

0个回答
没有发现任何回复~