窗函数的特性分析(DOC)

合集下载

数字信号处理窗函数特性

数字信号处理窗函数特性

信息科学与技术学院实验报告课程名称: 数字信号处理实验项目: 窗函数的特性实验地点:博西105 指导教师: 日期: 2013年5月7日实验类型:综合性实验(验证性实验综合性实验设计性实验)专业: 电子信息班级:姓名: 学号:一、实验目的及要求1.分析各种窗函数时域和频域特性。

2.为灵活运用窗函数分析信号频谱和设计FIR数字滤波器奠定基础。

二、实验仪器、设备或软件计算机 MATLAB软件三、实验步骤(或过程)(一):实验程序:N=51;beta=2;w1=boxcar(N);W1=fft(w1,256);w2=hanning(N);W2=fft(w2,256);w3=hamming(N);W3=fft(w3,256);w4=blackman(N);W4=fft(w4,256);w5=bartlett(N);W5=fft(w5,256);w6=Kaiser(N);W6=fft(w6,256);figure(1);subplot(2,1,1);stem([0:N-1],w1);title('Rectangle 窗');Y1=abs(fftshift(W1));subplot(2,1,2);plot([-128:127],Y1);title('幅度谱');figure(2);subplot(2,1,1);stem([0:N-1],w2);title('Hanning窗');Y2=abs(fftshift(W2));subplot(2,1,2);plot([-128:127],Y2);title('幅度谱');figure(3);subplot(2,1,1);stem([0:N-1],w3);title('Hamming窗');Y3=abs(fftshift(W3));subplot(2,1,2);plot([-128:127],Y3);title('幅度谱');figure(4);subplot(2,1,1);stem([0:N-1],w4);title('Blacekman窗');Y4=abs(fftshift(W4));subplot(2,1,2);plot([-128:127],Y4);title('幅度谱');figure(5);subplot(2,1,1);stem([0:N-1],w5);title('Bartlett窗');Y5=abs(fftshift(W5));subplot(2,1,2);plot([-128:127],Y5);title('幅度谱');figure(6);subplot(2,1,1);stem([0:N-1],w6);title('Kaiser 窗');Y6=abs(fftshift(W6));subplot(2,1,2);plot([-128:127],Y6);title('幅度谱');(二)程序设计实验:1.利用fft 函数分析常用窗函数的频谱特性,并从主瓣宽度和旁瓣相对幅度两个角度分析比较。

实验三 窗函数的特性分析报告

实验三 窗函数的特性分析报告

数字信号处理及实验实验报告实验题目窗函数的特性分析MYT 组别班级学号【实验目的】分析各种窗函数的时域和频率特性,灵活运用窗函数分析信号频谱和设计FIR数字滤波器。

【实验原理】在确定信号谱分析、随机信号功率谱估计以及FIR数字滤波器设计中,窗函数的选择对频谱分析和滤波器设计都起着重要的作用。

在确定信号谱分析和随机信号功率谱估计中,截短无穷长的序列会造成频率泄漏,影响频谱分析的精度和质量。

合理选取窗函数的类型,可以改善泄漏现象。

在FIR数字滤波器设计中,截短无穷长的系统单位脉冲序列会造成FIR滤波器的幅度特性产生波动,且出现过渡带。

【实验结果与数据处理】1、分析并绘出常用窗函数的时域特性波形。

程序如下:clc,clear,close allN=50figure(1)W1=boxcar(N);stem([0:N-1],W1);figure(2)W2=hanning(N);stem([0:N-1],W2);figure(3)W3=hamming(N);stem([0:N-1],W3);figure(4)W4=blackman(N);stem([0:N-1],W4);figure(5)W5=bartlett(N);stem([0:N-1],W5);figure(6)W6=kaiser(N,2*N);stem([0:N-1],W6);时域波形图如下:图 1 矩形窗图 2 汉宁窗图 3 汉明窗图 4 布莱克曼窗图 5 Bartlett窗图 6 凯泽窗2、研究凯泽窗(Kaiser)的参数选择对其时域和频域的影响。

(1)固定beta=4,分别取N=20,60,110。

clc,clear,close allN1=20;N2=60;N3=110;beat=4;figure(1)subplot(3,2,[1,2])W=kaiser(N1,beat);stem([0:N1-1],W);subplot(3,2,[3,4]);Ww=kaiser(N2,beat);stem([0:N2-1],Ww);subplot(3,2,[5,6]);WW=kaiser(N3,beat);stem([0:N3-1],WW);figure(2)subplot(3,2,[1,2])W1=fft(W,N1)plot([0:N1-1],abs(fftshift(W1)))subplot(3,2,[3,4]);W2=fft(Ww,N2)plot([0:N2-1],abs(fftshift(W2)))subplot(3,2,[5,6]);W3=fft(WW,N3)plot([0:N3-1],abs(fftshift(W3)))图7 凯泽窗频域图图8 凯泽窗时域图(2)固定N=60,分别取beta=1,5,11。

实验三窗函数的特性分析

实验三窗函数的特性分析

实验三窗函数的特性分析窗函数是在时间域上对信号进行加权的一种方法。

它在信号处理领域中应用广泛,用于去除频谱泄露和减少频谱波动。

窗函数可以改变信号的频谱特性,有助于减小频谱波动,提高频谱分析的准确性。

本实验将分析三种不同类型的窗函数:矩形窗、汉明窗和布莱克曼窗。

1.矩形窗:矩形窗是一种简单的窗函数,它将输入的信号乘以常数1、它在时间域上呈现出矩形的形状,频域上表现为sinc函数。

矩形窗的特点是具有较宽的主瓣,但是有很高的边瓣衰减,对于频谱泄露较为敏感。

它适用于信号频谱比较窄的情况,可以提供较好的分辨率。

2.汉明窗:汉明窗是一种平滑且对称的窗函数,它在时间域上具有一对对称的凸边,频域上表现为sinc-squared函数。

汉明窗的特点是在频域上拥有较窄的主瓣和较小的边瓣泄露。

这使得它在频谱分析中具有较好的分辨率和较低的波动。

它适用于信号频谱分析的大多数情况。

3.布莱克曼窗:布莱克曼窗是一种设计用于音频处理的窗函数,它在时间域和频域上都具有较好的性能。

它的形状和汉明窗类似,但有更宽的底部。

布莱克曼窗的特点是具有更强的边瓣抑制能力,相对于汉明窗能够更好地抑制频谱波动和频谱泄露。

它适用于对频谱准确性要求较高的信号处理任务。

综上所述,不同的窗函数在频域上具有不同的特性。

矩形窗适用于频谱较窄的信号,提供较好的分辨率;汉明窗适用于大多数频谱分析的情况,具有较低的波动;布莱克曼窗能够更好地抑制频谱波动和泄露,适用于对准确性要求较高的任务。

在实际应用中,选择窗函数需要根据具体的信号特性和分析需求来进行。

需要折衷考虑分析的准确性和频谱泄露问题,并选择合适的窗函数来优化频谱分析的结果。

简述窗函数法的特点

简述窗函数法的特点

简述窗函数法的特点大凡讲高中数学的人都会听说过“窗函数法”,但这一名词从未进入高中教材,其含义也并不为多数老师所了解。

在初等函数里,我们经常可以用到许多函数方程,如二次函数、指数函数等,这些方程的实际意义往往比较抽象,有的甚至难以理解。

在实际应用中,我们可以借助某种简便手段把这些抽象函数化成具体的形式,使得实际问题更易处理。

这就是“转化”,即将研究对象由初等函数转化成其相应的次级数学表达式或微分方程。

在实际运用中,人们发现当一个函数值域很大时,将这个函数表示为有限个基本函数的乘积的形式,要比用有限个次级表达式去近似它更为简洁和准确。

因此人们便引入了“窗”( window)函数的概念,并逐渐将研究对象由基本函数转化成窗函数,再利用窗函数求出有限个基本函数的乘积的近似值,这样既简化了问题又节省了计算时间。

首先,“窗”函数有两个要素:一是函数的变化范围,一般是区间端点;二是函数表达式。

其次,所谓“窗函数法”就是在“转化”思想指导下建立起来的,它有两个层次的含义: 1.将函数转化成适当的次级数学表达式, 2.求函数f(x)在闭区间[a,b]上的近似值。

这种转化思想,就是“转化”。

“窗函数法”的特点是:①省去了繁琐的计算过程,保留了计算结果;②通过利用窗函数探索研究区间上连续函数的性质。

“窗函数法”中对应的函数方程是,在函数区间端点附近存在一个内接正实数x,使得。

而在开区间(a,b),则定义为。

也可以将窗函数定义为,当时,而且则。

当时,使得。

“窗函数法”的特点是,它只需作少量的初等变换,就能把数学模型的微分方程由原形式推广为非齐次、甚至可能不存在齐次解的形式。

同时,只要作较小的变换,使区间的范围缩小到函数表达式的允许误差范围之内,或改变解的形式,便可直接得到原微分方程的解。

此外,由于这种方法是依据初等函数的性质,运用解析法的一些基本技巧而构造出来的,因此,比之解析法,能得到更多、更好的近似解。

这些近似解和原形式之间的偏差是很小的。

FIR窗函数_信息与通信_工程科技_专业资料.doc

FIR窗函数_信息与通信_工程科技_专业资料.doc

窗函数性能的三个频域指标:Aw=—(1) 3db 带宽B ,单位为 N (最大可能的频率分辨力) (2) 最大旁瓣峰值A(dB) , A 越小,由旁瓣引起的谱失真越小 (3) 旁瓣谱峰渐进衰减速度D ( dB/oct )一个理想的窗函数,应该有最小的3dB 带宽、最小的旁瓣峰值及最大的旁 瓣衰减。

(1)矩形窗w(ri) = rlo, else其振幅响应为sinWr(p ))-——-sin在对过渡带和阻带衰减进行精确分析的过程中,对窗振幅响应进行连续积分 幅响应)是^必要的,即矩形窗函数w(n)以及它的振幅响应、累积振幅响应如下图所示。

振幅响应右W= W1Ns搭(2)必=M»1s 灯处具有第一个零点,因而主瓣的宽度为2厶3,所以过渡带宽也近似为2A 3。

其3dB 帯宽为B=0・886A3。

大约在w = 3^/N 出现第一个旁瓣(即主旁瓣),其幅度为:将它与主瓣振幅N 比较,则最大静峰值A (dB )为A=-13. 46 db.累积振幅相应的第一个旁瓣为2】dB,这个21dB 的阻蒂衰滅与窗长度N 无关。

根揭最小阻带衰滅,可以精确地计算出过渡带宽为:1.&T _ w” = ---旁瓣谱峰渐进衰滅速度D=-6db/octo(2)三角窗(或巴特利特Bartlett 窗)由于矩形窗从0到1(或1到0)有一个突变的过渡带,造成了吉布斯现象o Bartlett 提出了一种三角窗形式,它是两个矩形窗的卷积。

2刃 … N --- , YI — 0,1,...,— N 2一防,n =—2B=1.28A w,A=-27db, D=-12db/oct ,近似过渡带宽 8 /N,精确过渡带宽 6.R/N,最 小阻带衰减25dBo•1frequency in pi uiN-\AccumJW^IWs-aqo4°1......................I IB c i b e ls N J-1 0 1 frequency in pi unisAmplitude Response indB 13tloii!Amplitude^esponse : N=45 :巴dthF27«£«•triang (n)Forn odd:iv[k] = *2kn + 1' 2(n-k + 1)72 + 1iv [k 十 1]2k2 2佻) 一 n -Yfrequmcy BI piuiteTrimgultf Window :N=452246¥Air|Auh Re^ome in dB在Matlab 屮,函数bartlett (n)和triang(n)用来计算三角窗w 二 bartlett (n)• Forn oddForn even:♦ For n even〔2伙)2(n - k-1) n -1bartlett 函数返回的序列两端总是0,因此,对于奇数n,语句bartlett (n+2) 的中间部分等于triang(n);对于偶数n, bartlett 仍然是两个矩形序列的卷积。

常用窗函数的特点

常用窗函数的特点

常⽤窗函数的特点 转载⾃:1.矩形窗矩形窗相当使信号突然截断所乘的窗函数,它的旁瓣较⼤,且衰减较慢,旁瓣的死⼀个负峰值为主瓣的21%,第⼀个正峰值为主瓣的12.6%,第⼆个负负峰值为主瓣的9%,故巨星唱效果不适很好,泄漏较⼤。

2.汉宁窗汉宁窗的频谱时间上是由三个矩形窗经相互平移叠加⼆乘,汉宁窗的第⼀旁瓣幅值是主瓣的0.027%,这样旁瓣可以最⼤限度地互相抵消,从⽽达到加强主瓣的作⽤,使泄漏得到较为有效的抑制。

采⽤汉宁窗可以是主瓣加宽,倍频程衰减为18dB/otc,虽然平率分辨率⽐矩形窗稍有下降,但频谱幅值精度⼤为提⾼,因此,对要求显⽰不同频段上各频率成分的不同贡献⽽不关⼼频率分辨率的问题式时,建议使⽤汉宁窗。

3.海明窗海明窗与汉宁窗同属于余弦窗函数,它⽐汉宁窗在减⼩旁瓣幅值⽅⾯效果较好,但主瓣⽐汉宁窗也稍微宽⼀些。

海明窗的最⼤旁瓣⽐汉宁窗低,约为汉宁窗的1/5,其主瓣衰减率可达40dB/otc,这是海明窗⽐汉宁窗的优越之处。

但是海明窗的旁瓣衰减不及汉宁窗迅速,这是海明窗的缺点。

4.布莱克曼窗布莱克曼窗和汉宁窗及海明窗⼀样同属于⼴义余弦窗函数。

在与汉宁窗及海明窗相同长度的条件下,布莱克曼窗的主瓣稍宽,旁瓣⾼度稍低。

5.三⾓窗三⾓窗旁瓣较⼩,且⽆负值,衰减较快,但主瓣宽度加⼤,且使信号产⽣畸变。

6.余弦坡度窗余弦坡度窗是振动信号处理中常⽤的⼀种窗函数,是由矩形窗加汉宁窗组合⽽成。

它的窗函数曲线⼤部分持续时间⾥很平,如同矩形窗那样,之后加⼀段汉宁窗,平滑衰减到阶段处。

余弦坡度窗的有点介于矩形窗和汉宁窗之间。

因为矩形窗的频率主瓣窄,谱值衰减⼩,⽽汉宁窗的旁瓣⼩,主瓣宽。

因此,把两者结合起来取长补短,达到既有较窄频率主瓣,⼜有较好抑制谱泄漏效果。

7.帕曾窗帕曾窗是⼀种⾼次幂窗,但主瓣⽐汉宁窗窄,主瓣幅值⾼⼀些。

8.指数窗指数窗常⽤与结构冲击实验的数据处理。

当系统收到瞬态激励时,往往要做⾃由衰减运动,如果过结构的阻尼很⼩,幅值衰减的时间就越长,在进⾏有限点采样时因时域阶段⽽产⽣的能量泄漏就越⼤,频谱产⽣的畸变也就越严重。

(实验三窗函数的特性分析)

(实验三窗函数的特性分析)

实验报告实验课程:数字信号处理实验开课时间:2020—2021 学年秋季学期实验名称:窗函数的特性分析实验时间:2020年9月16日星期三学院:物理与电子信息学院年级:大三班级:182 学号:1843202000234 姓名:武建璋一、实验预习(2)固定N=60,分别取beta=1,5,11。

clc,clear,close allbeat1=1;beat2=5;beat3=11;N=60;figure(1)subplot(3,2,[1,2])W=kaiser(N,beat1);stem([0:N-1],W);subplot(3,2,[3,4]);Ww=kaiser(N,beat2);stem([0:N-1],Ww);subplot(3,2,[5,6]);WW=kaiser(N,beat3);stem([0:N-1],WW);figure(2)subplot(3,2,[1,2])W1=fft(W,N)plot([0:N-1],abs(fftshift(W1))) subplot(3,2,[3,4]);W2=fft(Ww,N)plot([0:N-1],abs(fftshift(W2))) subplot(3,2,[5,6]);W3=fft(WW,N)plot([0:N-1],abs(fftshift(W3)))4、某序列为x[k] = (11πk/20) + cos(9πk/20),使用fft函数分析其频谱。

(1) 利用不同宽度N的矩形窗截短该序列,N分别为20,40,160,观察不同长度N 的窗对谱分析结果的影响。

clc,clear,close allN1=20;N2=40;N3=160;k1=0:N1;k2=0:N2;k3=0:N3;X1=0.5.*cos((11*pi*k1)/20)+cos((9*pi*k1)/20)X2=0.5.*cos((11*pi*k2)/20)+cos((9*pi*k2)/20)X3=0.5.*cos((11*pi*k3)/20)+cos((9*pi*k3)/20)figure(1)subplot(3,2,[1,2])W1=fft(X1,N1)plot([0:N1-1],abs(fftshift(W1)))subplot(3,2,[3,4]);W2=fft(X2,N2)plot([0:N2-1],abs(fftshift(W2)))subplot(3,2,[5,6]);W3=fft(X3,N3)plot([0:N3-1],abs(fftshift(W3)))figure(2)subplot(3,2,[1,2])W=abs(fftshift(W1))stem([0:N1-1],W);subplot(3,2,[3,4]);Ww=abs(fftshift(W2))stem([0:N2-1],Ww);subplot(3,2,[5,6]);WW=abs(fftshift(W3))stem([0:N3-1],WW);(2) 利用汉明窗重做(1)。

如何选择窗函数窗函数的分析比较

如何选择窗函数窗函数的分析比较

如何选择窗函数窗函数的分析比较窗函数在信号处理和频谱分析中起着重要的作用,用于改善信号的频谱性质,以便更好地分析信号。

选择适合的窗函数可以提高信号的频域分辨率和抑制频谱泄漏。

首先,需要了解窗函数的基本概念和特性,以便更好地进行选择和分析。

1.窗函数的定义:窗函数是定义在有限时间和频率范围内的函数,用于将信号在时间和频域上进行截断。

常见的窗函数包括矩形窗、汉宁窗、汉明窗、布莱克曼窗等。

2.窗函数的性质:不同的窗函数具有不同的性质,如频域主瓣宽度、旁瓣衰减、频域泄漏等。

选择窗函数时需要考虑这些性质,以满足实际需求。

在选择窗函数时,需要考虑以下几个方面:1. 频域主瓣宽度:频域主瓣宽度反映了窗函数的频域分辨能力,即能否准确地分辨出信号的频率。

主瓣越窄,频率分辨能力越高。

因此,在需要高频率分辨率的应用中,应选择主瓣宽度较窄的窗函数,如Kaiser 窗、Slepian窗等。

2. 旁瓣衰减:窗函数的旁瓣衰减反映了窗函数对于频域旁瓣的抑制能力。

旁瓣越低,表示频域泄漏越小,能更好地抑制邻近频率的干扰。

因此,在需要高频域抑制能力的应用中,应选择旁瓣衰减较大的窗函数,如Blackman窗、Nuttall窗等。

3.时域响应:窗函数的时域响应直接影响波形的平滑程度和能否准确地表示信号的时域特征。

时域响应平滑的窗函数可以减小信号的突变,但也会造成时间分辨率的损失。

因此,在需要准确表示信号时域特征的应用中,应选择合适的时域响应窗函数,如Gaussian窗、Dolph-Chebyshev 窗等。

4.计算效率:窗函数的计算效率也是选择的重要因素。

复杂的窗函数可能需要更多的计算资源和消耗更多的时间。

因此,在需要实时处理和高效率计算的应用中,应选择计算效率较高的窗函数,如矩形窗和汉宁窗。

综合考虑以上因素,可以根据不同应用需求选择合适的窗函数。

在实际应用中,也可以通过试验和比较不同窗函数的效果,选择最符合要求的窗函数。

需要注意的是,窗函数的选择并没有绝对的标准,要根据具体的应用需求来进行选择,并对选择的窗函数进行分析和评估。

实验三窗函数特性分析

实验三窗函数特性分析

实验三窗函数特性分析窗函数特性分析是信号处理领域中一个重要的研究方向,通过对窗函数的分析可以有效地应用于噪声抑制、频谱分析等方面。

下面我们来详细分析几个常见的窗函数特性。

1.矩形窗矩形窗函数也被称为哈曼窗,其表达式为:w(n)={1(n∈[0,N-1])0otherwise(1)其中,N表示窗口长度。

矩形窗函数在频域上等效为一个 sinc 函数,其主瓣宽度与窗口长度成反比。

由于矩形窗函数在主瓣两侧具有较深的零点,因此具有较高的频率分辨率。

然而,由于其旁瓣较大,矩形窗函数容易产生假响应和泄露现象。

2.汉宁窗汉宁窗函数是一种改进的矩形窗函数,通过在矩形窗函数的基础上增加两个旁瓣,以减小旁瓣电平并抑制假响应。

汉宁窗函数的表达式为:w(n)=0.5−0.5cos⁡(2πnN−1)(2)其中,N表示窗口长度。

与矩形窗函数相比,汉宁窗函数的主瓣宽度增加了,旁瓣电平也较低。

在保持较高频率分辨率的同时,减小了假响应的可能性。

3.哈曼窗哈曼窗函数是一种基于最小旁瓣电平为目标的窗函数,通过调整汉宁窗函数的系数,使得旁瓣电平最小。

哈曼窗函数的表达式为:w(n)=0.4935N+0.4834cos⁡(2πnN−1)+0.0133cos⁡(4πnN−1)(3)其中,N表示窗口长度。

哈曼窗函数在主瓣两侧具有较深的零点,同时旁瓣电平较低,具有较高的频率分辨率和较小的假响应。

4.高斯窗高斯窗函数是一种基于高斯函数的窗函数,具有平滑的旁瓣衰减和较小的旁瓣电平。

高斯窗函数的表达式为:w(n)=e−n2/(2σ2)(4)其中,σ表示高斯函数的方差,N表示窗口长度。

高斯窗函数的主瓣宽度与窗口长度成反比,旁瓣电平随着远离主瓣而逐渐增大。

由于其旁瓣衰减较慢,高斯窗函数容易产生交叉干扰现象。

通过对以上常见窗函数的特性分析可知,不同的窗函数具有不同的频率响应特性。

在应用中需要根据具体需求选择合适的窗函数。

例如,当需要高频率分辨率时,可以选择矩形窗函数;当需要抑制假响应时,可以选择汉宁窗函数或哈曼窗函数;当需要平滑的旁瓣衰减时,可以选择高斯窗函数。

窗函数

窗函数

窗函数(window function)窗函数是频谱分析中一个重要的部分,CoCo包含了所有通用的窗函数以及冲击测试中的受迫/指数(force/exponential)窗。

窗函数修正了由于信号的非周期性并减小了频谱中由于泄露而带来的测量不准确性。

快速傅里叶变换假定了时间信号是周期无限的。

但在分析时,我们往往只截取其中的一部分,因此需要加窗以减小泄露。

窗函数可以加在时域,也可以加在频域上,但在时域上加窗更为普遍。

截断效应带来了泄漏,窗函数是为了减小这个截断效应,其设计成一组加权系数。

例如,一个窗函数可以定义为:w(t)=g(t) -T/2<t<T/2w(t)=0 其他g(t)是窗函数,T是窗函数的时间待分析的数据x(t)则表示为:x(t)=w(t)*x(t)'x(t)'表示原始信号x(t)表示待分析信号。

加窗在时域上表现的是点乘,因此在频域上则表现为卷积。

卷积可以被看成是一个平滑的过程。

这个平滑过程可以被看出是由一组具有特定函数形状的滤波器,因此,原始信号中在某一频率点上的能量会结合滤波器的形状表现出来,从而减小泄漏。

基于这个原理,人们通常在时域上直接加窗。

大多数的信号分析仪一般使用矩形窗(rectangular),汉宁(hann),flattop和其他的一些窗函数。

矩形窗函数:w(k)=1汉宁窗:w(k)=0.5*(1-cos(2*pi*k/(N-1))) 0<=k<=N-1由于加窗计算中衰减了原始信号的部分能量,因此对于最后的结果还需要加上修正系数。

在线性谱分析中,一般使用幅度系数(amplitude correction),在功率谱中,一般使用能量系数(energy correction)。

具体请看下以章节。

泄露效应对于简单的信号,比如一个单频率的正弦波,泄露就表现为不在其频率点上仍然会有能量的出现。

离其本身的频率越近的频率,泄露的情况越严重,而离的越远,则情况则会好一些。

窗函数的特性分析

窗函数的特性分析

窗函数的特性分析
窗函数技术是滤波器设计的重要部分。

它主要用来控制信号滤波器的
频率响应特性。

窗函数包括矩形窗,三角窗,汉宁窗,汉明窗,Hamming 窗,Kaiser窗等。

本文通过分析各种窗函数的特性,从而指导滤波器设
计的实现。

一、矩形窗函数的特性
矩形窗函数的特性是信号量和宽度恒定,即信号量不随时间变化,宽
度也不变,如下形式所示:
w[n]=1(0≤n≤N-1)
矩形窗的经典应用是定义时间信号的加权数,即叠加N个信号之和,
是滤波器设计的最基本的窗函数,但其窗函数的频率响应特性比较差。

二、三角窗函数的特性
三角窗函数是矩形窗函数的改进,其特性是信号量和宽度随时间变化,即信号量随时间变化,宽度也随时间变化,如下形式所示:
w[n]={1-,n-(N-1)/2,/(N-1)/2}(0≤n≤N-1)
三角窗函数的频率响应特性比矩形窗函数略好,同时在设计滤波器时
可以使用它,如果在误差允许的范围内的话。

三、汉宁窗函数的特性
汉宁窗函数是三角窗函数的一种变形函数,其特性是信号量和宽度随
时间变化,但信号量只允许有限的值,如下形式所示:
w[n]=1-{1-,2n/N-1,}^2(0≤n≤N-1)
汉宁窗函数的频率响应特性比三角窗函数略好。

【最新推荐】tukey窗函数-精选word文档 (9页)

【最新推荐】tukey窗函数-精选word文档 (9页)

本文部分内容来自网络整理,本司不为其真实性负责,如有异议或侵权请及时联系,本司将立即删除!== 本文为word格式,下载后可方便编辑和修改! ==tukey窗函数篇一:窗函数的实现及分析1窗函数1.1基本概念在实际进行数字信号处理时,往往需要把信号的观察时间限制在一定的时间间隔内,只需要选择一段时间信号对其进行分析。

这样,取用有限个数据,即将信号数据截断的过程,就等于将信号进行加窗函数操作。

而这样操作以后,常常会发生频谱分量从其正常频谱扩展开来的现象,即所谓的“频谱泄漏”。

当进行离散傅立叶变换时,时域中的截断是必需的,因此泄漏效应也是离散傅立叶变换所固有的,必须进行抑制。

而要对频谱泄漏进行抑制,可以通过窗函数加权抑制DFT的等效滤波器的振幅特性的副瓣,或用窗函数加权使有限长度的输入信号周期延拓后在边界上尽量减少不连续程度的方法实现。

而在后面的FIR滤波器的设计中,为获得有限长单位取样响应,需要用窗函数截断无限长单位取样响应序列。

另外,在功率谱估计中也要遇到窗函数加权问题。

窗函数的基本概念。

设x(n)是一个长序列,w(n)是长度为N的窗函数,用w(n)截断x(n),得到N点序列xn(n),即xn(n) = x(n) w(n)在频域上则有XNej????1πj?j?????Xe?Wed???π2π????由此可见,窗函数w(n)不仅仅会影响原信号x(n)在时域上的波形,而且也会影响到频域内的形状。

1.2设计原理窗函数设计法的基本原理是用有限长单位脉冲响应序列h?n?逼近hd?n?。

由于hd?n?往往是无限长序列,而且是非因果的,所以用窗函数??n?将hd?n?截断,并进行加权处理,得到:h?n??hd?n???n?h?n?就作为实际设计的FIR数字滤波器的单位脉冲响应序列,其频率响应函数Hej???为He????h?n?ej?n?0N?1j?n式中,N为所选窗函数??n?的长度。

用窗函数法设计的滤波器性能取决于窗函数??n?的1类型及窗口长度N的取值。

实验三窗函数的特性分析

实验三窗函数的特性分析

实验三窗函数的特性分析
一.窗函数的概念
窗函数是一种算法,它是一种带有其中一种形状的函数,通过对信号
进行处理,可以增强信号的一些特征,从而改善信号的可检测性和抑制噪声。

窗函数的定义:它在一些时间段上取特定的值,而在此之外的时间段上,则取零。

在细分时间段上,都按照固定的函数变换来求取取值,以保
证窗函数满足频率应答的要求。

二.常用窗函数
1)矩形窗函数:即矩形窗,也称为方形窗,最简单的窗函数形式,
是通过将脉冲在时间上延伸,而延伸后的脉冲形态则形成了“矩形”这样
一种特殊形状,从而被称为矩形窗。

2)凯廷窗:也称为汉明窗,是在矩形窗的基础上,进一步改进的一
种窗函数形式,是最常用的窗函数之一,它采用对称的函数形式,使得其
在频率响应上比矩形窗更加接近极低通滤波器的频率响应,从而有效地提
高了信号抑制噪声的能力,同时也保持了信号的清晰度。

3)高斯窗:又称为高斯滤波器,是一种基于高斯分布特性的滤波器,它的函数形状完全符合高斯分布的概率分布,在低噪声、低失真的环境中,效果最佳,是非常常用的窗函数。

4)黎曼窗:又叫黎曼汉明窗,它的特点是连续非均匀。

常用窗函数的特性与选用---Erwin

常用窗函数的特性与选用---Erwin

常用窗函数的特性与选用首先应该感谢Erwin站长,是他发起了这个好帖,我来整理下,供大家方便阅读学习和讨论!还是先列个提纲,慢慢补充内容。

1 什么是窗?2 为什么要窗?3 常用窗函数的时频特性与适用范围4 窗函数的综合比较与选用1 什么是“窗”?这个要从傅里叶分析说起。

从傅里叶分析本身定义看,它是对连续函数进行的,此时是没有窗的概念的。

但是傅里叶分析在理论上具有无限的完美性,但实用时却遇到很大的困难,因为它是一个积分表达式。

1963年,两位牛人提出了FFT的基本思想可以看做是傅里叶分析实用领域的一大突破(当然FFT计算量很大,直到计算机得到高速发展后FFT才有了广泛的应用)。

FFT有一个基本概念就是block,也就是一个数据块,FFT是对一个数据块的数据按照蝶型算法进行的。

那么,如何从一个连续的信号得到一个block以便进行FFT呢?这就需要一个窗从连续信号上截取一个block下来。

“窗”就是这样一个工具,用来从连续时间信号中提出一段有限的数据。

2 为什么要“窗”?答案很简单,加窗的目的有两个: 1)减小泄露; 2)改善栅栏效应;名词解释:泄露(leakage)在从一个连续信号中抽取一个block的过程,如果不加窗,实际上就是默认加了一个矩形窗,如下图示。

这样数据抽取的结果,就是使得原来连续信号中集中在f0处的能量被分散到两个较宽的频带中去了,这种现象称之为频谱能量泄漏(Leakage)。

栅栏效应(Picket Fence Effect)对信号做FFT时,得到的是一系列离散的谱线,如果信号中的频率成份位于谱线之间而不是正好落在谱线上,此时就会造成幅值和相位上的偏差。

离散的一条条谱线就象一个栅栏,因此称这种现象为栅栏效应。

栅栏效应可以形象地做一个比喻。

正如在栅栏外走过一个美女,你目不转睛去看,但因为栅栏效应总有一些关键部位被挡住,使得美女在一定程度上有失真。

这就是栅栏效应。

3 常用窗函数的时频特性与适用范围-- uniform widow(1)uniform窗,即又名矩形窗,具有较大的实时带宽,最大栅栏效应误差为3.92dB。

常用窗函数的特性与选用

常用窗函数的特性与选用

常用窗函数的特性与选用在数字信号处理领域,窗函数是一种非常重要的工具,用于改善信号的频谱特性。

它们在时域和频域中都有特定的作用,可以帮助我们更好地理解和分析信号。

本文将介绍几种常用的窗函数,并探讨它们各自的特性和选用方法。

一、矩形窗矩形窗是最简单的一种窗函数,其特性如下:1. 优点:计算简单,处理速度快。

2. 缺点:主瓣宽度较宽,旁瓣较大,导致频率分辨率较低,频谱泄露严重。

选用矩形窗的场景:当信号处理速度要求较高,且对频率分辨率要求不高时,可以选用矩形窗。

二、汉宁窗汉宁窗是一种常用的窗函数,其特性如下:1. 优点:主瓣宽度适中,旁瓣较小,频率分辨率较高,频谱泄露较少。

2. 缺点:计算相对复杂,处理速度较慢。

选用汉宁窗的场景:当信号处理速度要求适中,且对频率分辨率要求较高时,可以选用汉宁窗。

三、汉明窗汉明窗是汉宁窗的一种变体,其特性如下:1. 优点:主瓣宽度较窄,旁瓣较小,频率分辨率较高,频谱泄露较少。

2. 缺点:计算相对复杂,处理速度较慢。

选用汉明窗的场景:当信号处理速度要求适中,且对频率分辨率要求较高时,可以选用汉明窗。

四、布莱克曼窗1. 优点:主瓣宽度较窄,旁瓣较小,频率分辨率较高,频谱泄露较少。

2. 缺点:计算复杂,处理速度较慢。

选用布莱克曼窗的场景:当信号处理速度要求较低,且对频率分辨率要求较高时,可以选用布莱克曼窗。

五、凯泽窗凯泽窗是一种可调窗函数,其特性如下:1. 优点:通过调整参数,可以灵活控制主瓣宽度和旁瓣高度,以满足不同场景的需求。

2. 缺点:计算复杂,处理速度较慢。

选用凯泽窗的场景:当信号处理速度要求较低,且对频率分辨率和旁瓣高度有特殊要求时,可以选用凯泽窗。

根据信号处理速度、频率分辨率和旁瓣高度等需求,我们可以选择合适的窗函数。

在实际应用中,我们需要权衡各种窗函数的优缺点,以便在满足需求的前提下,提高信号处理的性能。

六、窗函数的选择与优化1. 了解信号特性:在选用窗函数之前,要了解信号的特性,包括频率成分、信号长度等。

常见的窗函数及基本参数

常见的窗函数及基本参数

常见的窗函数及基本参数一、概述在信号处理中,窗函数是一种用于减少频谱泄漏和增加频谱分辨率的技术。

它们通常用于傅里叶变换和相关算法中。

窗函数是一个非常重要的概念,因为它们可以帮助我们更好地理解信号处理中的许多问题。

在本文中,我们将介绍一些常见的窗函数及其基本参数。

二、矩形窗函数矩形窗函数是最简单的窗函数之一,也称为“盒形窗”。

它是一个由0和1组成的序列,其中1表示数据被保留在该位置上,0表示数据被舍弃。

它的数学表达式如下:w(n) = 1, 0 <= n <= N-1w(n) = 0, 其他情况其中N为序列长度。

三、汉明窗函数汉明窗函数是一种平滑过渡到零值的窗函数。

它可以减少频谱泄漏,并且具有较好的抑制旁瓣能力。

它的数学表达式如下:w(n) = a - (1-a) * cos(2*pi*n/(N-1)), 0 <= n <= N-1其中a为系数,通常取0.54。

四、汉宁窗函数汉宁窗函数是一种平滑过渡到零值的窗函数,与汉明窗函数类似。

它也可以减少频谱泄漏,并且具有较好的抑制旁瓣能力。

它的数学表达式如下:w(n) = 0.5 - 0.5 * cos(2*pi*n/(N-1)), 0 <= n <= N-1五、布莱克曼窗函数布莱克曼窗函数是一种平滑过渡到零值的窗函数,具有较好的抑制旁瓣能力。

它的数学表达式如下:w(n) = a0 - a1*cos(2*pi*n/(N-1)) + a2*cos(4*pi*n/(N-1)) -a3*cos(6*pi*n/(N-1)) + a4*cos(8*pi*n/(N-1)), 0 <= n <= N-1其中a0=0.42, a1=0.5, a2=0.08, a3=0.025, a4=0.01。

六、海明窗函数海明窗函数是一种平滑过渡到零值的窗函数,具有良好的旁瓣抑制能力。

它的数学表达式如下:w(n) = 0.54 - 0.46 * cos(2*pi*n/(N-1)), 0 <= n <= N-1七、升余弦窗函数升余弦窗函数是一种平滑过渡到零值的窗函数,具有较好的旁瓣抑制能力。

窗函数的实现及分析

窗函数的实现及分析

窗函数的实现及分析窗函数是指将理想的频谱截断成有限的频谱,并对信号进行加权的函数。

在信号处理中,窗函数被广泛应用于频谱分析、滤波器设计、波形合成和信号的时频分析等方面。

其作用是减小频谱泄漏、降低旁瓣干扰和改善频谱估计的准确性。

1. 直接实现法(Direct Approach):直接实现法是指通过直接计算窗函数的定义式来得到窗函数的采样值。

例如,常见的矩形窗函数可以通过以下公式计算得到:w(n)=1,0<=n<Nw(n)=0,其他情况其中,n为窗函数的采样序号,N为窗函数的长度。

类似地,其他窗函数如汉宁窗、汉明窗、布莱克曼窗等也可以通过相应的定义式计算得到。

直接实现法的优点是实现简单,计算速度快。

缺点是窗函数的采样点数需要提前确定,并且无法根据需要动态调整窗函数的长度。

此外,直接实现法在频率分辨率方面相对较差,易产生频谱泄漏现象。

2. 卷积实现法(Convolution Approach):卷积实现法是指利用卷积运算的性质,通过将序列信号和窗函数进行卷积来实现窗函数。

例如,矩形窗可以通过以下卷积运算实现:w(n)=RECT(n)=δ(n)*δ(n)其中,δ(n)为单位脉冲函数。

卷积实现法的优点是可以根据需要动态调整窗函数的长度和形状,适应不同的信号分析要求。

此外,卷积实现法拥有较好的频率分辨率和抗频谱泄漏能力。

对于窗函数的分析,可以从以下几个方面进行:1.主瓣宽度:主瓣宽度是指窗函数的主瓣在频谱中的宽度。

窗函数的主瓣宽度决定了频率分辨率的能力,主瓣宽度越窄,频率分辨率越高。

例如,矩形窗的主瓣宽度较宽,频谱分辨率相对较低;而汉宁窗、汉明窗等窗函数的主瓣宽度相对较窄,频谱分辨率较高。

2.旁瓣干扰:旁瓣干扰是指窗函数在频谱中产生的旁瓣能量。

窗函数的旁瓣干扰会引入频谱泄漏现象,降低频谱估计的准确性。

一般而言,窗函数的旁瓣干扰越低,频谱估计的准确性越高。

常见的窗函数如布莱克曼窗具有较低的旁瓣干扰能力。

窗函数类型和特点

窗函数类型和特点

窗函数有多种类型,它们各自具有不同的特性和用途。

以下是一些常见的窗函数类型及其特点:
1.矩形窗(Rectangular Window):矩形窗的幅度为常数1,对应于时域上的矩形函
数。

在频域上,矩形窗的频谱为sinc函数,具有宽主瓣和高旁瓣的特性。

矩形窗的优点是主瓣比较集中,缺点是旁瓣较高,并有负旁瓣,导致变换中带进了高频干扰和泄漏,甚至出现负谱现象。

2.汉明窗(Hamming Window):汉明窗是一种常用的窗函数,具有较低的旁瓣衰
减。

它在时域上的形式为0.54 - 0.46*cos(2πn/N),n为时域上的序号,N为窗口长度。

汉明窗在频域上具有较为平坦的幅度响应和较低的旁瓣,适用于需要减少旁瓣干扰的应用场景。

3.海宁窗(Hanning Window):海宁窗是汉明窗的一种特殊形式,也被称为汉宁
窗。

它在时域和频域上的特性与汉明窗类似,但是旁瓣衰减更快。

海宁窗在时域上表示为正弦函数和余弦函数的组合,其频谱具有较为尖锐的主瓣和较低的旁瓣。

适用于需要快速衰减旁瓣的应用场景。

4.凯泽窗(Kaiser Window):凯泽窗是一种具有可调参数的窗函数,其旁瓣衰减速
度和主瓣宽度可以根据参数进行调整。

凯泽窗在频域上具有较为尖锐的主瓣和快速衰减的旁瓣,适用于需要灵活控制窗函数特性的应用场景。

这些窗函数类型各有特点,应根据具体应用场景选择适合的窗函数类型。

窗函数

窗函数

窗函数(window function)窗函数是频谱分析中一个重要的部分,CoCo包含了所有通用的窗函数以及冲击测试中的受迫/指数(force/exponential)窗。

窗函数修正了由于信号的非周期性并减小了频谱中由于泄露而带来的测量不准确性。

快速傅里叶变换假定了时间信号是周期无限的。

但在分析时,我们往往只截取其中的一部分,因此需要加窗以减小泄露。

窗函数可以加在时域,也可以加在频域上,但在时域上加窗更为普遍。

截断效应带来了泄漏,窗函数是为了减小这个截断效应,其设计成一组加权系数。

例如,一个窗函数可以定义为:w(t)=g(t) -T/2<t<T/2w(t)=0 其他g(t)是窗函数,T是窗函数的时间待分析的数据x(t)则表示为:x(t)=w(t)*x(t)'x(t)'表示原始信号x(t)表示待分析信号。

加窗在时域上表现的是点乘,因此在频域上则表现为卷积。

卷积可以被看成是一个平滑的过程。

这个平滑过程可以被看出是由一组具有特定函数形状的滤波器,因此,原始信号中在某一频率点上的能量会结合滤波器的形状表现出来,从而减小泄漏。

基于这个原理,人们通常在时域上直接加窗。

大多数的信号分析仪一般使用矩形窗(rectangular),汉宁(hann),flattop和其他的一些窗函数。

矩形窗函数:w(k)=1汉宁窗:w(k)=0.5*(1-cos(2*pi*k/(N-1))) 0<=k<=N-1由于加窗计算中衰减了原始信号的部分能量,因此对于最后的结果还需要加上修正系数。

在线性谱分析中,一般使用幅度系数(amplitude correction),在功率谱中,一般使用能量系数(energy correction)。

具体请看下以章节。

泄露效应对于简单的信号,比如一个单频率的正弦波,泄露就表现为不在其频率点上仍然会有能量的出现。

离其本身的频率越近的频率,泄露的情况越严重,而离的越远,则情况则会好一些。

4窗函数及频谱分析

4窗函数及频谱分析

4窗函数及频谱分析窗函数是一种对信号进行截断的方法,通常用于频谱分析、滤波和信号重构等应用中。

在频谱分析中,窗函数可以降低谱泄漏效应,提高频谱分辨率。

常见的窗函数有矩形窗、汉明窗、汉宁窗和布莱克曼窗等。

矩形窗是最简单的窗函数,它将信号在截断窗口内的值保持不变,其频谱包含了原信号的所有频率成分。

然而,由于矩形窗的边界不连续,会引入频谱泄漏效应,导致在频域上原本存在的尖峰变得模糊不清。

汉明窗是一种典型的对称窗函数,它的特点是边界平滑,能够降低频谱泄漏效应。

汉明窗的形状类似于一个谐波振荡器,对信号进行截断时,它在边界处施加了平滑的衰减,从而减少了频谱泄漏。

汉宁窗是汉明窗的变种,它在衰减的斜率上更加陡峭,能够进一步减少频谱泄漏。

汉宁窗的形状类似于一个凸起的典型窗口,信号在截断窗口内的值按照窗函数的形状进行衰减。

布莱克曼窗是一种优化的窗函数,它在边界处的衰减斜率更加陡峭,能够最大程度地减少频谱泄漏效应。

布莱克曼窗在截断窗口内的信号衰减方式是非线性的,可以更好地适应信号的动态范围。

频谱分析是一种对信号进行频域分析的方法。

通过将信号转换到频率域,可以得到信号的频率成分和相应的幅度信息。

频谱分析可以用于信号的谱线显示、频谱检测和频谱估计等应用中。

频谱分析的基本原理是将信号转换到频域,通常使用快速傅里叶变换(FFT)算法完成。

FFT算法能够高效地计算信号的频谱,对于长度为N 的信号,FFT算法的时间复杂度为O(NlogN)。

在实际应用中,频谱分析常常使用窗函数对信号进行截断。

通过选择合适的窗函数,可以减少频谱泄漏效应,提高频谱分辨率。

频谱分析还可以通过窗函数长度的选择来调节频谱分辨率和频谱平滑度。

总结起来,窗函数是对信号进行截断的方法,常用于频谱分析中。

常见的窗函数有矩形窗、汉明窗、汉宁窗和布莱克曼窗等。

频谱分析是一种对信号进行频域分析的方法,通过将信号转换到频率域,可以得到信号的频率成分和相应的幅度信息。

使用合适的窗函数可以降低频谱泄漏效应,提高频谱分辨率。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

本科学生验证性实验报告
学号********* 姓名李开斌
学院物理与电子信息专业、班级11电子
实验课程名称窗函数的特性分析
教师及职称李宏宁
开课学期2013 至2014 学年下学期
填报时间2014 年03 月26 日云南师范大学教务处编印
1.实验现象与结果
1.分析并绘出常用窗函数的时域特性波形%矩形窗时域波形及频谱
N=51;
w=boxcar(N);
Y=fft(w,256);
subplot(2,1,1);
stem([0:N-1],w);
xlabel('w');
ylabel('y');
title('时域波形');
subplot(2,1,2);
Y0= abs(fftshift(Y));
plot([-128:127], Y0)
xlabel('W');
ylabel('Y0');
title('频谱图形');
N=51;
k=0:N-1;
w=0.54-0.46*cos(2*pi*k/(N-1)) Y=fft(w,256);
subplot(2,1,1);
stem([0:N-1],w);
xlabel('w');
ylabel('y');
title('时域波形');
subplot(2,1,2);
Y0= abs(fftshift(Y));
plot([-128:127], Y0)
xlabel('W');
ylabel('Y0');
title('频谱图形');
N=51;
k=0:N-1;
w=1/2*(1-cos(2*pi*k/(N-1))); Y=fft(w,256);
subplot(2,1,1);
stem([0:N-1],w);
xlabel('w');
ylabel('y');
title('时域波形');
subplot(2,1,2);
Y0= abs(fftshift(Y));
plot([-128:127], Y0)
xlabel('W');
ylabel('Y0');
title('频谱图形');
N=51;
w=bartlett(N);
Y=fft(w,256); subplot(2,1,1); stem([0:N-1],w); xlabel('w');
ylabel('y');
title('时域波形'); subplot(2,1,2);
Y0= abs(fftshift(Y)); plot([-128:127], Y0) xlabel('W');
ylabel('Y0');
title('频谱图形');
N=51;
beta=4;
w=Kaiser(N,beta); Y=fft(w,256); subplot(2,1,1); stem([0:N-1],w); xlabel('w');
ylabel('y');
title('时域波形'); subplot(2,1,2);
Y0= abs(fftshift(Y)); plot([-128:127], Y0) xlabel('W');
ylabel('Y0');
title('频谱图形');
plot([-128:127], Y0)
xlabel('W');
ylabel('Y0');
title('频谱图形beta=4,N=110');
%凯撒窗时域波形及频谱N=60,beta取不同值的波形比较N=60;
beta=1;
w=Kaiser(N,beta);
Y=fft(w,256);
subplot(3,2,1);
stem([0:N-1],w);
xlabel('w');
ylabel('y');
title('时域波形N=60,beta=1');
subplot(3,2,2);
Y0= abs(fftshift(Y));
plot([-128:127], Y0)
xlabel('W');
ylabel('Y0');
title('频谱图形N=60,beta=1');
N=60;
4.某序列为x[k]=0.5cos(k 2011π)+cos(k 20
9π),使用fft 函数分析其频谱。

%利用矩形窗分析序列
N=20;
k=0:N-1;
x=0.5*cos(11*pi*k/20)+cos(9*pi*k/20);
w=ones(1,N);
y=x.*w;
Y=fft(y,512);
subplot(3,2,1);
stem([0:N-1],y);
title('抽样信号');
xlabel('频率');
ylabel('幅值');
subplot(3,2,2);
Y0=abs(fftshift(Y));
plot([-256:255], Y0);
title('时域波形');
%利用汉明窗分析序列
N=20;
k=0:N-1;
x=0.5*cos(11*pi*k/20)+cos(9*pi*k/20); w=1/2*(1-cos(2*pi*k/(N-1)));
y=x.*w;
Y=fft(y,512);
subplot(3,2,1);
stem([0:N-1],y);
title('抽样信号');
xlabel('频率');
ylabel('幅值');
subplot(3,2,2);
Y0=abs(fftshift(Y));
plot([-256:255], Y0);
title('时域波形');
xlabel('频率');
ylabel('幅值');
N=40;
k=0:N-1;
x=0.5*cos(11*pi*k/20)+cos(9*pi*k/20); w=1/2*(1-cos(2*pi*k/(N-1)));
y=x.*w;
Y=fft(y,512);
subplot(3,2,3);
stem([0:N-1],y);
title('抽样信号');
xlabel('频率');
ylabel('幅值');
subplot(3,2,4);
Y0=abs(fftshift(Y));
plot([-256:255], Y0);
title('时域波形');
xlabel('频率');
ylabel('幅值');
N=160;
k=0:N-1;
x=0.5*cos(11*pi*k/20)+cos(9*pi*k/20); w=1/2*(1-cos(2*pi*k/(N-1)));
y=x.*w;
Y=fft(y,512);
subplot(3,2,5);
stem([0:N-1],y);
title('抽样信号');
xlabel('频率');
ylabel('幅值');
subplot(3,2,6);
Y0=abs(fftshift(Y));
plot([-256:255], Y0);
title('时域波形');
xlabel('频率');
ylabel('幅值');
subplot(3,2,5);
stem([0:N-1],y);
title('抽样信号');
xlabel('频率');
ylabel('幅值');
subplot(3,2,6);
Y0=abs(fftshift(Y));
plot([-256:255], Y0);
title('时域波形');
xlabel('频率');
ylabel('幅值');
2.实验总结
通过本次实验,我学会了分析常用窗函数的时域和频域特性,懂得了灵活运用窗函数分析信号频谱和设计FIR数字滤波器。

教师评语及评分:
签名:年月日。

相关文档
最新文档