null - 程序员宅基地

FIR---Matlab 如何选择窗函数和 FFT 的长度_matlab 窗函数-程序员宅基地

技术标签： matlab 算法信息与通信 Matlab数字信号处理信号处理

FIR---Matlab 如何选择窗函数和 FFT 的长度

各种窗函数的介绍，见：

数据截断、频谱泄漏与窗函数的选择_布莱克曼窗_yishuihan-Oliver的博客-程序员宅基地

1、常用的四种窗函数

对于实际信号序列，如何选取窗函数呢？一般来说，选择第一旁瓣衰减大，旁瓣峰值衰减快的窗函数，有利于緩解截断过程中产生的频谱泄漏问题。但具有这两个特性的窗函数，其主瓣寛度较大，相应会帯来一些副作用，应用中需根据具体情况折中地选择。

设信号中包含 $f_{a}$ 和 $f_{b}$ 两个频率分量，窗函数的选择与两个频率分量的间距以及两个频率分量的幅度比例密切相关。窗函数选择的一般准则如表2-2-2所列。

主瓣窄的窗函数一般旁瓣泄漏大，频谱泄漏主要集中在旁瓣范围内。旁瓣衰减大的窗函数，一般主瓣较宽，泄漏主要集中在主瓣范围内。

当选择加窗DFT时，已知采祥长度 N 的选择与最小頻域分辨率有关，而从表2-2-1中看到采长度 N 还与窗函数的主瓣寛度相关，所以 N 的选择应为： $N\geqslant \frac{f_{s}}{\Delta f_{min}}*K$ (式 2-2-26)

式中：K 为窗函数的主瓣寛度与矩形窗的主瓣寛度之比。

根据窗函数对数据处理的影响，可参照下列原則选取理想的窗函数：

①窗函数频谱的主瓣尽可能地窄，以提高谱估计时的頻域分辨率和减小泄漏。
②尽量减小窗函数频谱的最大旁瓣的相对幅度，以使旁瓣高度随頻率尽快衰减。

若以上两条不能同时得到満足，则往往是増加主瓣寛度以換取旁瓣的抑制。

总之，在应用窗函数时，除要考窗函数本身的特性外，还应充分考虑被分析信号的特点及具体处理要求。

2、实例

例：已知信号 $x(t)$ 中 50Hz 频率分量的幅值为 311，46Hz 频率分量的幅值为 1.55，采样频率 $f_{s}=8000Hz$ 。要求 46Hz 信号的幅度分析进度不低于 5%，试问：

选择何种类型的窗函数比较合适？
采样长度 N 应为多少？
分析信号的实际频谱。

分析：由于 50Hz 频率分量幅度远大于 46Hz 频率分量，需要防止 50Hz 频率分量对 46Hz 频率分量的泄漏，同时考虑 46Hz 分量分析精度的要求，允许 50Hz 频率分量的最大泄漏为 $20*\lg \left | \frac{1.55}{311} *0.05\right |=-80dB$

以上介绍的窗函数中没有一个窗函数的第一旁瓣衰减能达到 -80dB；但可以选择旁瓣衰减大、高频衰减速度快的窗函数，以满足实际衰减要求。根据表2-2-1，选择布莱克曼窗比较合适。布莱克曼窗第3个旁瓣衰减为（-58-18*3）= -112dB。

此时，采样长度的选择不但要考虑窗函数的主瓣宽度，还要考虑旁瓣位置：

$N\geqslant \frac{f_{s}}{\Delta f_{min}}*(K+M)=\frac{8000}{50-46}*(3+3)=12000$

式中：K为窗函数的主瓣宽度与矩形窗的主瓣宽度之比，M为旁瓣位置。按表2-2-1中信息，矩形窗的主瓣宽度为 $\frac{4\pi }{N}$ ，而布莱克曼窗的主瓣宽度为 $\frac{12\pi }{N}$ ，所以 K=3；又取布莱克曼窗第3个旁瓣，M=3。

Matlab 程序如下：

% Matlab 如何选择窗函数和 FFT 的长度
% 
% pr2_2_11  
clear; clc; close all;

f1 = 50; a1 = 311; % 设置第1个分量的频率与幅值
f2 = 46; a2 = 1.57; % 设置第2个分量的频率与幅值
N = 12000; % 设置数据长度N
fs = 8000; % 设置采样频率fs
t = (0:N-1)/fs; % 设置时间刻度
x = a1*cos(2*pi*f1*t)+a2*cos(2*pi*f2*t); % 设置信号
freq = (0:N/2)*fs/N; % 设置频率刻度
wind = blackman(N)'; % 给出布莱克曼窗函数
X = fft(x.*wind); % FFT

% 作图
plot(freq,20*log10(abs(X(1:N/2+1))),'k');
grid; xlim([0 100])
xlabel('频率/Hz');
ylabel('幅值/dB');
title('信号谱图');
set(gcf,'color','w');

本文链接：https://blog.csdn.net/yishuihanq/article/details/132760292

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

心脏滴血漏洞HeartBleed CVE-2014-0160深入代码层面的分析_heartbleed代码分析-程序员宅基地

文章浏览阅读1.7k次。心脏滴血漏洞HeartBleed CVE-2014-0160 是由heartbeat功能引入的，本文从深入码层面的分析该漏洞产生的原因_heartbleed代码分析

java读取ofd文档内容_ofd电子文档内容分析工具（分析文档、签章和证书）-程序员宅基地

文章浏览阅读1.4k次。前言ofd是国家文档标准，其对标的文档格式是pdf。ofd文档是容器格式文件，ofd其实就是压缩包。将ofd文件后缀改为.zip，解压后可看到文件包含的内容。ofd文件分析工具下载：点我下载。ofd文件解压后，可以看到如下内容：对于xml文件，可以用文本工具查看。但是对于印章文件(Seal.esl)、签名文件(SignedValue.dat)就无法查看其内容了。本人开发一款ofd内容查看器，..._signedvalue.dat

基于FPGA的数据采集系统（一）_基于fpga的信息采集-程序员宅基地

文章浏览阅读1.8w次，点赞29次，收藏313次。整体系统设计本设计主要是对ADC和DAC的使用，主要实现功能流程为：首先通过串口向FPGA发送控制信号，控制DAC芯片tlv5618进行DA装换，转换的数据存在ROM中，转换开始时读取ROM中数据进行读取转换。其次用按键控制adc128s052进行模数转换100次，模数转换数据存储到FIFO中，再从FIFO中读取数据通过串口输出显示在pc上。其整体系统框图如下：图1：FPGA数据采集系统框图从图中可以看出，该系统主要包括9个模块：串口接收模块、按键消抖模块、按键控制模块、ROM模块、D.._基于fpga的信息采集

微服务 spring cloud zuul com.netflix.zuul.exception.ZuulException GENERAL-程序员宅基地

文章浏览阅读2.5w次。1.背景错误信息：-- [http-nio-9904-exec-5] o.s.c.n.z.filters.post.SendErrorFilter : Error during filteringcom.netflix.zuul.exception.ZuulException: Forwarding error at org.springframework.cloud..._com.netflix.zuul.exception.zuulexception

邻接矩阵-建立图-程序员宅基地

文章浏览阅读358次。1.介绍图的相关概念　　图是由顶点的有穷非空集和一个描述顶点之间关系-边（或者弧）的集合组成。通常，图中的数据元素被称为顶点，顶点间的关系用边表示，图通常用字母G表示，图的顶点通常用字母V表示，所以图可以定义为:　　G=(V,E)其中，V(G)是图中顶点的有穷非空集合，E(G)是V(G)中顶点的边的有穷集合1.1 无向图：图中任意两个顶点构成的边是没有方向的1.2 有向图：图中..._给定一个邻接矩阵未必能够造出一个图

IAR调试程序闪退问题_iar在debug时停止工作-程序员宅基地

文章浏览阅读1.7k次。问题:IAR调试STM32程序，点击调试按钮后软件自动关闭，并弹出报错提示框解决：将调试的接口模式改为SWD模式即可。我的原先设置为JTAG模式。_iar在debug时停止工作

随便推点

python--xlrd/xlwt/xlutils_xlutils模块可以读xlsx吗-程序员宅基地

文章浏览阅读219次。python–xlrd/xlwt/xlutilsxlrd只能读取，不能改,支持 xlsx和xls 格式xlwt只能改，不能读xlwt只能保存为.xls格式xlutils能将xlrd.Book转为xlwt.Workbook，从而得以在现有xls的基础上修改数据，并创建一个新的xls，实现修改xlrd打开文件import xlrdexcel=xlrd.open_workbook('E:/test.xlsx') 返回值为xlrd.book.Book对象,不能修改获取sheett_xlutils模块可以读xlsx吗

关于新版本selenium定位元素报错：‘WebDriver‘ object has no attribute ‘find_element_by_id‘等问题_unresolved attribute reference 'find_element_by_id-程序员宅基地

文章浏览阅读8.2w次，点赞267次，收藏656次。运行Selenium出现'WebDriver' object has no attribute 'find_element_by_id'或AttributeError: 'WebDriver' object has no attribute 'find_element_by_xpath'等定位元素代码错误，是因为selenium更新到了新的版本，以前的一些语法经过改动。..............._unresolved attribute reference 'find_element_by_id' for class 'webdriver

FIR---Matlab 如何选择窗函数和 FFT 的长度_matlab 窗函数-程序员宅基地

FIR---Matlab 如何选择窗函数和 FFT 的长度

1、常用的四种窗函数

2、实例

智能推荐

心脏滴血漏洞HeartBleed CVE-2014-0160深入代码层面的分析_heartbleed代码分析-程序员宅基地

java读取ofd文档内容_ofd电子文档内容分析工具（分析文档、签章和证书）-程序员宅基地

基于FPGA的数据采集系统（一）_基于fpga的信息采集-程序员宅基地

微服务 spring cloud zuul com.netflix.zuul.exception.ZuulException GENERAL-程序员宅基地

邻接矩阵-建立图-程序员宅基地

IAR调试程序闪退问题_iar在debug时停止工作-程序员宅基地

随便推点

python--xlrd/xlwt/xlutils_xlutils模块可以读xlsx吗-程序员宅基地

关于新版本selenium定位元素报错：‘WebDriver‘ object has no attribute ‘find_element_by_id‘等问题_unresolved attribute reference 'find_element_by_id-程序员宅基地

DOM对象转换成jQuery对象转换与子页面获取父页面DOM对象-程序员宅基地

什么是算法？-程序员宅基地

【网络安全】网络安全的标准和规范_网络安全标准规范-程序员宅基地

Windows上的巧克力味Chocolatey详解_chocolate怎么卸载-程序员宅基地

推荐文章

热门文章

相关标签