matlab的fft关于频谱泄露和加窗.pdfVIP

精通MATLAB语言，有MATLAB编程问题的朋友，请直接联系我！ 联系方式： QQ：993878382 Email：993878382@ MATLAB博客：d.html 更多学习资料和电子书，请直接访问我的博客！ 如果我QQ不在线，请将问题发到我的邮箱，或者直接给我QQ留言，第 一时间答复你！简单问题请直接留言，复杂问题可以提供解决思路！ 同探讨，共进步！ 还可以为公司和科研单位设计各类算法，欢迎前来咨询！ 郑重声明：本资源来源于网络，仅限交流使用，请勿用于商业用途！ 最近看到不少论坛上的朋友都总结了一些自己的所思所想，我最近也在做FFT方面的东西，有 很多感触：一方面对自己工作的总结；一方面和大家一起分享一下希望对大家会有一点帮助。 我们分析的信号，如果只含整数次谐波的话，用FFT 分析信号的频谱和相位是非常准确的，如果信号含有确 定的间谐波，比如信号含有60HZ和65HZ的频率，那我 们也可以准确的分析出信号的频谱和相位，我们只要用 矩形窗截取 10 个周波的信号就可以分析出 50HZ/10=5HZ 以及 5HZ 的整数倍的信号的频谱和相位 了，分析的相位和频谱都是非常准确的！如果我们需要 分析的信号含有不确定的间谐波，比如我们根本不知道 信号含有什么样的间谐波，那么此时用FFT分析，必然 会有频谱泄露！那么怎么样才能减少频谱泄露呢（注意： 这种情况下我们只能减少频谱泄露，而基本上不可能完 全消除频谱泄露）？ 我们可以有两种方法：方法一：增加采样的长度，方 法二：加窗函数。 增加采样长度可以分析出更多频率的信号，可以减 少频谱泄露，不过增加采样长度必然会对数据处理的实 时性造成影响！ 理想的窗函数是主瓣很窄，旁瓣衰减很快，矩形窗 的主瓣很窄，但是旁瓣衰减却很慢，hanning窗、hamming 窗、blackman 窗等的旁瓣衰减有了明显的改进，但是主 瓣却宽了很多，大概是矩形窗主瓣的二倍，blackman 窗 的主瓣还要宽，这就造成了信号频谱的频率识别率很低！ 什么叫频率识别率呢？简单举个例子解释一下：比 如说我取十个周波（请注意：这里我用十个周波举例子 是为了便于计算，实际中我们是不取10 个周波的，因为 基2 的FFT运算要求取样点数是2 的整数倍整数次幂， 我们实际操作是取4 、8 或者16 个周波等2 的整数次幂）， 我用一个含有 25HZ 间谐波的信号做分析，那么我加 hamming窗和hanning窗以后，在20HZ和30HZ处的频率 点上都将有幅值，而且有趣的是：我们将20HZ 、25HZ 和 30HZ频率点处的幅值相加以后基本上接近与信号在 25HZ处的真实幅值（这一点我还没搞清楚，是否隐含什 么关系），并且在频谱上看到的25HZ处的幅值要小于实 际信号在 25HZ处的幅值，所以我们要对加窗后的FFT 变换的幅值要乘以一个恢复系数，不同窗的恢复系数也 是不同，矩形窗的是1，hanning窗的是2 ，hamming窗的 好像是1.84 左右吧（不好意思，记不清了，大家可以在 网上查到）！这样我们在分析25HZ频率点处的幅值时， 对于20HZ和30HZ频率点处的幅值都是不可信的，所以 我们至少要求20HZ和30HZ附近是不能有信号的，这样 频率分辨率相对来说是不是就降低了？因为 20HZ和 30HZ是不能用的！（哈哈，语言组织能力差了点，有什 么不严谨的地方大家包含啊）！ 大家可想而知，如果我的原始信号在20HZ和30HZ 处本来就都有幅值，那么20HZ和30HZ处的频率也将反 过来影响 25HZ处的信号，这样测出来的幅值误差将会 很大，所以我们在加窗时由一个要求：那就是谐波和间 谐波，以及间谐波与间谐波之间要相隔较远，我觉得， 至少要相隔 2 条谱线以上我们才能获得较高精度的幅 值！blackman 窗要求相隔的谱线还要多！这是因为 blackman窗的主瓣还要宽！ 所以对于加窗函数我们要权衡利弊，根据具体需要 来选择！