FFT分析问题

yunyi009

回复 8# BeiYangMan 问题就在这啊  不知道为什么会差这么多。。。晚点把序列发上来

quantus

回复 6# yunyi009

推导过程要看DSP的书。你说的凸起跟我说的不是一回事，你的凸起太大了。我的凸起只有一点点。你这种凸起肯定是混入了直流分量，你看看你是不是在信号没有稳定的时候就做了fft。如果是那么你可以把仿真时间加长，等稳定之后再取输出点做fft。

quantus

还有你0.003fs和0.004fs处有两个tone，你要查查为什么出现这两个tone

quantus

回复 8# BeiYangMan

加窗底噪不一样才正常

BeiYangMan

回复 14# quantus


    是的，我的理解有问题。我只是感觉差别不应该那么大。

BeiYangMan

回复 6# yunyi009

频谱就是这样的，你的结果没错。
我也尝试不加窗看频谱，跟你的一样。应该跟序列突然中止导致的谱的泄漏，由于低频处的谱本来就很低，所以低频处就平了。

quantus

回复 16# BeiYangMan
你对加窗的理解可能有点偏差。加窗是让信号慢慢被截断到0.

BeiYangMan

   

回复  BeiYangMan
你对加窗的理解可能有点偏差。加窗是让信号慢慢被截断到0.
quantus 发表于 2014-11-7 13:09

                               
登录/注册后可看大图

你说的没错，我在我的文档中没有指出这点。
严格的说加窗的谱分析 是让序列的一端为0，另外一段慢慢衰减趋于0，但不会等于0。这样整个序列做想象中的延拓，才没有“突变”。
这是加窗的目的。
所以matlab中hann这个函数，有两个option，一个是 对称，一个是 周期。
做窗函数滤波器，当然用对称
而做序列谱分析，就用周期。
Schreier，他自己定义了一个hann，其实是那个周期的窗。是为了用着方便，因为matlab默认是对称的窗。

加窗的效果 就是降低了频率的分辨率，泄漏还是会有的。比如hann，频谱还是会泄漏的，，对于一个单频信号，其能量只是只往左右一个bin的地方泄漏而已。这是我个人的理解。

另外，他的谱上的非线性跟窗没关系，不加窗时，那两个谐波几乎被泄漏的低频噪声淹没。加窗之后，显现出来了。非线性通常是dac的mismatch或者电路的本身的非线性，与窗无关。尽管加窗是非线性操作。

BeiYangMan

   

回复  BeiYangMan
你对加窗的理解可能有点偏差。加窗是让信号慢慢被截断到0.
quantus 发表于 2014-11-7 13:09

                               
登录/注册后可看大图

其实我是在解释不加窗时，噪声低为什么有一段是平直的。
不加窗就是加矩形窗，序列做延拓时，两端对接的地方有突变，没有平滑过渡，导致突变。
频谱上就出现所谓泄漏。

yunyi009

所谓噪底的大小究竟跟什么因素有关系呢？我也觉得加窗与不加窗噪底的变化有些过大了。而且我的DFT输入序列正好具有整数m个周期，即输入信号的频率正好在k=m频率单元的中心，这样一来抽样点的位置正好是幅值为0的位置，泄漏现象被掩盖，如图所示。那么加不加窗的结果就不应该差这么多了吧？