第7章 傅里叶变换与滤波器形状

合集下载

第7章图像处理 课后答案

第7章图像处理 课后答案

7.1.1 图像金字塔
一系列以金字塔形状排列的分辨率逐步降低的图集合。 金字塔的底部是待处理图像的高分辨率表示,而顶部 是低分辨率的近似。基级J的尺寸是2J×2J或N×N (J=log2N), 中间级j的尺寸是2j×2j ,其中0<= j <=J。
图7.2b表示,各级的近似值和预测残差金字塔都是以 一种迭代的方式进行计算的。第一次迭代和传递时, j = J ,并且2J×2J的原始图像作为J级的输入图像,从 而产生J-1级近似值和J级预测残差,而J-1级近似值又 作为下一次迭代的输入,得到J-2级近似值和J-1级预 测残差。 迭代算法:
1, 0 x 0.5 ψ( x) 1, 0.5 x 1 0,在,有了尺度函数和小波函数,可以正式定义小 波变换了,它包括:一般小波序列展开、离散小波 变换和连续小波变换。
7.3.1 小波序列展开
首先根据小波函数ψ( x)和尺度函数 ( x)为函数f(x)定 义小波序列展开:
高斯近似值和预测残差金字塔
基级,第9级
第8级 第7级 第6级
图像重建
7.1.2 子带编码
另一种与多分辨率分析相关的重要图像技术是子带 编码。在子带编码中,一幅图像被分解成为一系列 限带分量的集合,称为子带,它们可以重组在一起 无失真地重建原始图像。最初是为语音(一维信号) 和图像压缩而研制的,每个子带通过对输入进行带 通滤波而得到(相当于分解一个频段为若干个子频 段)。因为得到的子带的带宽要比原始图像的带宽 小,子带可以无信息损失的抽样。 原始图像的重建可以通过内插、滤波和叠加单个子 带来完成。
k
V Spk an{k ( x)}
7.2.2 尺度函数
现在来考虑由整数平移和实数二值尺度、平方可积 函数 ( x) 组成的展开函数集合,即集合{ j ,k ( x)} : j/2 j j ,k ( x) 2 (2 x k ) 式7.2.10 k决定了 j ,k ( x)在x轴的位置(平移k个单位),j决定 了 j ,k ( x) 的宽度,即沿x轴的宽或窄的程度,而2j/2 控制其高度或幅度。由于 j ,k ( x)的形状随j发生变化, ( x) 被称为尺度函数。 如果为赋予一个定值,即j = j0,展开集合 { j0 ,k ( x)} 将是 { j ,k ( x)}的一个子集,一个子空间:

有限冲激响应滤波器的设计

有限冲激响应滤波器的设计
Ⅱ型线性相位滤波器的幅度函数和相位函数的特点: 幅度函数的特点: (1)当 时, =0,也就是说 在 处必然有一个零点; (2) 对 呈奇对称,对 呈偶对称。 相位函数的特点:同I型线性相位滤波器。
3.Ⅲ型线性相位滤波器
由于Ⅲ型线性相位滤波器关于 奇对称,且 为整数,所以,其频率响应可以表示为 (7.16) 其中
幅度函数为 (7.17) 相位函数为 (7.18)
2
则 就具有线性相位, 。
3
并考虑| |=1,等效地指定
4
(7.48)
5
根据DFT的性质可知,为保证 是实序列,应满足下列对称关系
6
(7.49)
7
由于
(7.50)
当N为偶数时, ;当N为奇数时, 。这样当N为偶数时,若按式(7.48)对 赋值,就不能满足式(7.49)的对应关系。由此,按如下原则对 赋值。
选取一个满意的窗函数,令
(7.41)
则 即为要设计的滤波器的单位抽样响应。
按上述方法设计的滤波器,由于满足了 的对称关系,因此都具有线性相位。
由DFT定义,得
(7.42)

(7.43)
频率抽样法是从频域出发,把给定的理想频率响应 加以等间隔抽样,即
Ⅳ型线性相位滤波器
型线性相位滤波器关于 奇对称,且N为偶数,所以 为非整数。其频率响应可以表示为 (7.19) 其中
1
2
幅度函数为 (7.20) 相位函数为 (7.21)
Ⅳ型线性相位滤波器的幅度函数和相位函数的特点: 幅度函数的特点: (1) 在 处必为零,也就是说 在 处为零点 ; (2) 在 处呈奇对称,在 处呈偶对称 相位函数的特点:同Ⅲ型线性相位滤波器。

数字信号处理第三版第七章

数字信号处理第三版第七章

对称,是满足式(7.1.9)的一组解,
因为cos[ω(n-τ)]关于n=τ偶对称,所以要求τ和h(n)满
足如下条件:

()
,
N1

2
2
h(n)h(N1n), 0≤ n≤ N1
(7.1.10)
2. 线性相位FIR滤波器幅度特性Hg(ω)的特点 实质上,幅度特性的特点就是线性相位FIR滤波
因为cos[ω(n-τ)]关于ω=0, π, 2π三点偶对称,所以由 式(7.1.11)可以看出,Hg(ω)关于ω=0, π, 2π三点偶对称。 因此情况1可以实现各种(低通、高通、带通、带阻)滤 波器。
情况2: h(n)=h(N-n-1), N为偶数。
仿照情况1的推导方法得到:
H ( e j ) H g () e j = N 1 h ( n ) e j n e j M 2 h ( n )c o s (( n ) )
第7章 有限脉冲响应数字滤波器的设计
7.1 线性相位FIR数字滤波器的条件和特点 7.2 利用窗函数法设计FIR滤波器 7.3 利用频率采样法设计FIR滤波器 7.4 利用等波纹最佳逼近法设计FIR滤波器 7.5 IIR和FIR数字滤波器的比较 7.6 几种特殊类型滤波器简介 7.7 滤波器分析设计工具FDATool
用情况3的推导过程可以得到:
M
Hg() 2h(n)sin[(n)] n0
(7.1.13)
N是偶数,τ=(N-1)/2=N/2-1/2。所以,当ω=0, 2π时,
sin[ω(n-τ)]=0;当ω=π时,sin[ω(n-τ)]=(-1)n- N/2, 为峰值点。而且sin [ω(n-τ)]关于过零点ω=0和
如何减少吉布斯效应的影响,设计一个满足要求的FIR滤波器呢? 直观上,增加矩形窗口的宽度(即加大N)可以减少吉布斯效应 的影响。N 时, 在主瓣附近, WRg(ω)近似为:

傅里叶变换与滤波器形状

傅里叶变换与滤波器形状

增益
增益
增益
高 通
增益
低 通
低 通
高 通
左列为线性幅度响应 右列为对数幅度响应
对数刻度可以将一个变化范围 很大的增益表示在同一个图上。
增益
带 通
带 通
增益
增益
带 阻
增益
带 阻
对于数字频率,只考虑 0 ~ π 的范围
幅 度 响 应
幅度响应表明该滤波器为高通滤波器
相 位 响 应
表明此脉冲响应 所描述的滤波器 为带通滤波器
第7章 傅里叶变换与滤波器形状
内容纲要
7.1 傅里叶变换基础
定义离散时间傅里叶变换
7.2
频率响应及其它形式
定义数字滤波器的频率响应;建立频率响应、 传输函数、差分方程和脉冲响应之间的联系
7.3
频率响应和滤波器形状
用离散时间傅里叶变换计算滤波器在正弦输入时的输出; 确定数字滤波器的幅度响应和相位响应;建立模拟频率 和数字频率的关系;说明极零点位置如何决定滤波器形 状;研究一阶系统和二阶系统实例
极点离单位圆距离 较远,因此滤波器的 选择性不好(不是很 好的滤波器)。

0.5
实部
实部
二阶滤波器
a b c d
极点越接近单 位圆,滤波器的 选择性越好。
e f g h
a
b
c
d
e
f
g
h
内容总结
增 大
减 小
低通滤波器 由极零点位置得到的幅值与由 频率响应得到的幅值一致。 滤波器形状
虚部
实部
滤波器形状
滤波器 (a) 虚部 滤波器 (b) 滤波器 (c)
实部
滤波器形状
0.8π

第7章 傅里叶变换与滤波器形状

第7章  傅里叶变换与滤波器形状

第7章 傅里叶变换与滤波器形状 7.1离散时间傅里叶变换基础离散时间傅里叶变换(DTFT )是数字信号分析的一个重要工具。

DTFT 把信号或滤波器从时域变换到频域,主要是为了研究信号或滤波器的频率特性。

该变换主要用于分析信号和滤波器的频谱性质。

对于信号,DTFT 提供的信息称为信号的频谱。

对于滤波系统,DTFT 得到的信息称为滤波器的频率响应(frequency response )。

它由两部分组成:幅度响应(magnitude response )和相位响应(phase response )。

幅度响应给出了滤波器的形状,通过它我们可以深入了解滤波器的工作特性。

信号x[n]的离散时间傅里叶变换定义为:()[]jn n X x n e∞-Ω=-∞Ω=∑,这里Ω为数字频率,单位弧度。

记为(){[]}X x n Ω=F利用欧拉公式,DTFT 变换为()[][](cos()sin())jn n n X x n ex n n j n ∞∞-Ω=-∞=-∞Ω==Ω-Ω∑∑变换()X Ω在每个不同的数字频率上可有不同的值,当信号x[n]与正弦或余弦“共振”时,最大。

也就是说,当x[n]以接近频率Ω变化时,()X Ω较大。

离散时间傅里叶变换反应了信号的频率。

例7.1 求如图信号的离散时间傅里叶变换注意,一般情况,DTFT 是复值。

例7.2 求信号x[n]=4(u[n]-u[n-3])的DTFT 。

离散时间傅里叶变换有两个重要的特性,时延特性和周期性。

00[]()[]()Fjn Fx n X x n n eX -Ω−−→Ω-−−→Ω(2)()X X πΩ+=ΩDTFT 是周期性的,周期为2π。

即离散时间傅里叶变换对所有的数字频率Ω,每2π重复一次,不断重复。

7.2.1 频率响应和差分方程 对差分方程逐项求DTFT∑∑==-=-Mk kN k kk n x b k n y a 0][][][]2[]1[][][]2[]1[][210210M n x b n x b n x b n x b N n y a n y a n y a n y a M N -++-+-+=-++-+-+0101()()()()()()j jN N j jM M a Y a e Y a e Y b X b eX b eX -Ω-Ω-Ω-ΩΩ+Ω++Ω=Ω+Ω++Ω0101()()()()j jN j jM N M a a e a e Y b b e b e X -Ω-Ω-Ω-Ω+++Ω=+++Ω010010()()()Mjk j jM k M k Nj jN jk N k k b eY b b e b eH X a a e a ea e-Ω-Ω-Ω=-Ω-Ω-Ω=Ω+++Ω===Ω+++∑∑例7.3 求差分方程频率响应y[n]=-0.85y[n-1]+0.5x[n].例7.4 求差分方程频率响应y[n]+0.1y[n-1]+0.85y[n-2]=x[n]-0.3x[n-1]7.2.2 频率响应和传输函数10101010()()()MkMkM k NNkN k k b zb b z b zY z H z X z a a z a za z---=---=+++===+++∑∑010010()()()Mjk j jM k M k Nj jN jk N k k b eY b b e b eH X a a e a e a e-Ω-Ω-Ω=-Ω-Ω-Ω=Ω+++Ω===Ω+++∑∑例 7.5 求滤波器的频率响应,它的传输函数21210.2()10.50.9z H z z z ----=++2210.2()10.50.9j j j e H e e-Ω-Ω-Ω-Ω=++频率响应是脉冲响应的DTFT 。

傅里叶变换和频率域滤波的介绍

傅里叶变换和频率域滤波的介绍

傅立叶变换和频率域滤波的介绍序言.傅立叶变换的作用和意义1、下图中,最后一个波是由前面四个波组合而成的,是用傅立叶变换,可以很容易将其区分出来。

2、对比物理上对光谱的理解(赤橙黄绿青蓝紫、初中对光的三棱镜分解、高中对燃烧的钠元素所发出光的光谱分析实验),可以将傅立叶变换理解成“数学上的三棱镜”,傅立叶变换使我们能够通过频率成分来分析一个函数。

一.一维傅立叶变换及其反变换单变量连续函数f(x)的傅立叶变换F(u)定义为等式:(1)其中j=。

相反,给定F(u),通过傅立叶反变换可以获得f(x):(2)这两个等式组成了傅立叶变换对。

很明显,()F u 是一个复函数,即 ()()()F u R u jI u =+ (3)R(u)和I(u)分别是F(u)的实部和虚部,其中 (4) 称为傅立叶变换的幅度或频率谱,同时 (5) 称为变换的相角或相位谱。

在研究图像增强时,我们主要关心频率谱的性质。

所以需要定义的另一个量是功率谱,它被定义为傅立叶变换的平方:(6) 术语“谱密度”也用来指功率谱。

这些等式很容易扩展到两个变量u 和v 的情况:(7)类似地,反变换为:(8) 例1、已知()(0,0)x f x e x β-∂=>∂>,而且()0(0)f x t =<=。

则其傅立叶变换为: 2(2)00(2)(2)00222222()2(2)2(2)(2)2(2)2(2)(2)x j ux x j ux j u x j u xF u e e dx e dx e dx e j uj u j u j u j u j u u u j u u ππππββββπββππππββππββπππ+∞+∞-∂--∂++∞-∂+-∂++∞==-==∂+-∂-==∂+∂+∂-∂-=∂+∂=-∂+∂+⎰⎰⎰我们的兴趣在于离散函数,所以将不停留在这些数学定义中。

然而,在某些情况下,利用这些等式比利用它们的离散形式更容易证明二维傅立叶变换的性质。

傅里叶变换与滤波器设计

傅里叶变换与滤波器设计

傅里叶变换与滤波器设计引言:傅里叶变换是一种重要的数学工具,它可以将一个函数在时域中的描述转换为频域中的描述。

滤波器设计是一种基于傅里叶变换的方法,用于对信号进行处理和改变其频率特性。

本文将详细介绍傅里叶变换的原理和滤波器设计的基本概念。

一、傅里叶变换的概述傅里叶变换是将信号从时域转换到频域的一种数学变换。

通过傅里叶变换,我们可以将一个信号分解为一系列的正弦和余弦分量,每个分量对应于不同的频率。

傅里叶变换可以分为离散傅里叶变换(DFT)和连续傅里叶变换(CTFT)两种形式。

二、离散傅里叶变换(DFT)与快速傅里叶变换(FFT)离散傅里叶变换是将连续信号转换为离散信号的一种方法。

在离散傅里叶变换中,信号被看作是一系列离散的采样点,通过计算每个采样点的幅值和相位,可以得到信号在频域中的表示。

而快速傅里叶变换是一种高效的计算离散傅里叶变换的方法,通过减少计算复杂度,可以快速地得到信号的频域表示。

三、滤波器的基本概念滤波器是一种电子设备或算法,用于改变信号的频率特性。

滤波器可以实现信号的滤波、增强某些频率成分或抑制其他频率成分的功能。

常见的滤波器有低通滤波器、高通滤波器、带通滤波器和带阻滤波器等。

四、滤波器设计方法滤波器设计是根据滤波器的要求和设计规范,选择适当的滤波器类型并确定其参数的过程。

常见的滤波器设计方法有时域方法和频域方法。

时域方法包括窗函数法和脉冲响应法,频域方法包括理想滤波器法和波导滤波器法。

根据不同的需求和应用场景,可以选择合适的设计方法和滤波器类型。

五、应用实例:语音信号处理中的滤波器设计滤波器设计在语音信号处理中具有重要的应用。

例如,在语音通信系统中,为了提高语音质量和降低噪声干扰,常常需要设计滤波器对语音信号进行降噪和增强。

此外,在语音识别和语音合成等领域中,滤波器设计也扮演着关键的角色。

六、总结傅里叶变换是一种重要的数学工具,可以将信号从时域转换到频域,为信号处理和滤波器设计提供了基础。

数字信号处理讲义第7章滤波器的设计方法

数字信号处理讲义第7章滤波器的设计方法

第7章滤波器的设计方法教学目的1.掌握由连续时间滤波器设计离散时间IIR滤波器的方法,包括冲激响应不变法,双线性变换法等;2.了解常用的窗函数,掌握低通IIR滤波器的频率变换法、用窗函数法设计FIR滤波器的方法;3.掌握FIR滤波器的逼近原理与设计方法。

教学重点与难点重点:本章是本课程的重中之重,滤波器的设计是核心内容之一。

1.连续时间滤波器设计离散时间IIR滤波器的方法,包括冲激响应不变法,双线性变换法等;2.常用的窗函数,掌握低通IIR滤波器的频率变换法、用窗函数法设计FIR滤波器的方法;3.掌握FIR滤波器的逼近原理与设计方法。

难点:1.冲激响应不变法,双线性变换法2.用窗函数法设计FIR滤波器FIR滤波器的逼近原理与设计方法基本概念7.0.1 选频滤波器的分类数字滤波器是数字信号处理的重要基础。

在对信号的过滤、检测与参数的估计等处理中, 数字滤波器是使用最广泛的线性系统。

数字滤波器是对数字信号实现滤波的线性时不变系统。

它将输入的数字序列通过特定运算转变为输出的数字序列。

因此,数字滤波器本质上是一台完成特定运算的数字计算机。

我们已经知道,一个输入序列x(n),通过一个单位脉冲响应为h(n)的线性时不变系统后,其输出响应y(n)为∑∞-)(y))()()(n(nn=m*=xmhnhx将上式两边经过傅里叶变换,可得式中,Y (e j ω)、X (e j ω)分别为输出序列和输入序列的频谱函数, H (ejω)是系统的频率响应函数。

可以看出,输入序列的频谱X (e j ω)经过滤波后,变为X (e j ω)H (e j ω)。

如果|H (e j ω)|的值在某些频率上是比较小的,则输入信号中的这些频率分量在输出信号中将被抑制掉。

因此,只要按照输入信号频谱的特点和处理信号的目的,适当选择H (ej ω),使得滤波后的X (e j ω)H (e j ω)符合人们的要求,这就是数字滤波器的滤波原理。

和模拟滤波器一样,线性数字滤波器按照频率响应的通带特性可划分为低通、高通、带通和带阻几种形式。

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

第7章傅里叶变换与滤波器形状7.1离散时间傅里叶变换基础离散时间傅里叶变换(DTFT)是数字信号分析的一个重要工具。

DTFT把信号或滤波器从时域变换到频域,主要是为了研究信号或滤波器的频率特性。

该变换主要用于分析信号和滤波器的频谱性质。

对于信号,DTFT提供的信息称为信号的频谱。

对于滤波系统,DTFT得到的信息称为滤波器的频率响应(frequency response)。

它由两部分组成:幅度响应(magnitude response)和相位响应(phase response)。

幅度响应给出了滤波器的形状,通过它我们可以深入了解滤波器的工作特性。

信号x[n]的离散时间傅里叶变换定义为:,这里为数字频率,单位弧度。

记为利用欧拉公式,DTFT变换为变换在每个不同的数字频率上可有不同的值,当信号x[n]与正弦或余弦“共振”时,最大。

也就是说,当x[n]以接近频率变化时,较大。

离散时间傅里叶变换反应了信号的频率。

例7.1 求如图信号的离散时间傅里叶变换注意,一般情况,DTFT是复值。

例7.2 求信号x[n]=4(u[n]-u[n-3])的DTFT。

离散时间傅里叶变换有两个重要的特性,时延特性和周期性。

DTFT是周期性的,周期为。

即离散时间傅里叶变换对所有的数字频率,每重复一次,不断重复。

7.2 频率响应及其他形式7.2.1 频率响应和差分方程对差分方程逐项求DTFT例7.3 求差分方程频率响应y[n]=-0.85y[n-1]+0.5x[n].例7.4 求差分方程频率响应y[n]+0.1y[n-1]+0.85y[n-2]=x[n]-0.3x[n-1]7.2.2 频率响应和传输函数例 7.5 求滤波器的频率响应,它的传输函数7.2.3 频率响应和脉冲响应频率响应是脉冲响应的DTFT。

例7.6 数字滤波器的脉冲响应写出其频率响应。

7.3 频率响应与滤波器形状7.3.1 滤波器对正弦输入的作用由于复杂信号可以由各种频率和相位的正弦波叠加而成,我们先考虑单一频率即正弦信号的输入。

时域与频域的输出关系。

频率响应是个复数,可表示成是增益,无量纲,但可用分贝dB表示,此时增益为。

是相位差,单位是度或弧度。

增益是对输入的放大量,相位差是对输入的相移。

对于给定的频率,输出的幅度是滤波器的增益与输入幅度的乘积,输出相位是滤波器相位差与输入相位的和。

7.3.2 幅度响应和相位响应幅度响应是增益与频率的关系图。

相位响应是相位与频率的关系图。

例7.9 系统频率响应,每间隔pi/4弧度,计算相应的频率响应,绘制幅度响应和相位响应图。

幅度响应和相位响应是周期的,每2pi弧度重复一次。

幅度响应和相位响应是连续函数,在每个频率上有值。

幅度响应是偶函数,相位响应是奇函数。

由于负频率没有实际意义,在0~pi间已经包含了所有重要的信息。

采用分贝的优点是,在增益变化范围非常大时,可以方便的绘制在一个图上。

对数刻度实际上是对原图进行比例缩小。

弧度和度对相位响应形状没有什么影响。

例7.10 滤波器幅度响应和相位响应如图,可以观察到不同频率对应的滤波器增益和相位差。

例7.11 把信号加到具有如图所示频率响应的滤波器上。

增益、分贝表示的幅度响应和弧度、度表示的相位响应。

各种滤波器形状。

例7.13 滤波器频率响应。

例7.14 滤波器脉冲响应。

例7.15 滤波器的差分方程y[n]=1.5y[n-1]-0.85y[n-2]+x[n]例7.16 递归滤波器传输函数这是个梳状滤波器。

具有模拟回响的作用。

梳状滤波器的另一个例子第8章将介绍周期信号频谱,对于周期信号,其频谱是在固定间隔的尖峰上(谐波),梳状滤波器可以用来衰减谐波尖峰间的分量(噪声)。

7.3.3 模拟频率与数字频率不同采样频率对数字滤波器的影响。

以数字频率表示的滤波器特性,只有在采样频率选定后,才能确定。

例如,中心频率为0.6pi弧度的窄带滤波器,当采样频率为1kHz时,可通过300Hz的信号,而采样频率为4kHz时可通过1.2kHz的信号。

例7.18 采样频率分别为12kHz和30kHz时,低通滤波器截止频率分别为1800Hz和4500Hz。

例7.19 对于带通滤波器,采样频率分别为4kHz和10kHz时,观察它对1kHz输入信号的影响。

P196 例7.207.3.4 由极零点确定滤波器形状幅度与到极点的距离成反比,与到零点的距离成正比。

极点与单位圆距离越近的位置幅度越大,越远的位置幅度越小。

零点与单位圆距离越近的位置幅度越小,越远的位置幅度越大。

例7.21 滤波器传输函数例7.22 滤波器差分方程y[n]=x[n-1]+x[n-3]例7.23 比较滤波器形状7.3.5 一阶滤波器形状一阶系统:例7.24 7.25观察两个滤波器y[n]+0.5y[n-1]=x[n]和y[n]-0.5y[n-1]=x[n]的形状。

7.3.6 二阶滤波器形状二阶系统:我们给出8个滤波器的形状P203,其极零点图在第6章 P159 所示。

这其中包括了低通、高通与带通滤波器,极点距离单位圆近的滤波器选择性较好。

第8章数字信号频谱8.1 数字信号频谱的意义前章研究了滤波器的频率特性,本章则研究信号的频率特性。

在频域每个数字信号都有其典型的特征。

比如正弦波含单一频率,白噪声(white noise)包含所有频率分量。

信号的平稳变化源于它的低频分量,陡峭边缘和急剧变化源于它的高频分量。

例如方波,即包含产生平稳变化的低频分量,又包含了形成陡峭边缘的高频分量。

信号的频谱(spectrum)详细描述了信号所包含的频率分量。

了解信号频率特性对滤波器设计极为重要。

例如,要用扬声器再现音乐,那么音乐的频率分量决定了扬声器的特性。

语音识别系统中的麦克风必须有足够宽的频率响应,才能采集到语音的所有重要频率。

信号频谱由两部分构成,幅度频谱和相位频谱。

分非周期信号和周期信号考察其频谱特性。

8.2 非周期数字信号非周期信号就是不按固定间隔重复的信号。

计算非周期数字信号频谱,使用离散时间傅里叶变换(DTFT)。

DTFT计算需要信号所有采样值,当信号有无限多个递减的非零采样值时,就可以在幅度降低到某个合适的门限后截断信号,计算近似的DTFT。

非周期数字信号频谱是连续和周期性的,周期为。

当采样频率确定时对应的模拟频率为。

例8.1 求矩形脉冲x[n]=u[n]-u[n-4]的幅度响应和相位频谱。

幅度频谱是所有矩形脉冲具有的特性,称为sinc函数。

相位频谱仅用计算的点很难准确画出。

频谱是连续周期的。

0~п弧度的信息足以推出所有的频谱信息。

幅度频谱是偶函数,相位频谱是奇函数。

幅度频谱有零点。

频率的主要分量集中在低频部分。

采样频率确定时,才能获知频谱对应的模拟频率。

例8.2 求信号的频谱。

采样频率为15kHz。

由于信号衰减很快,前3个采样点足以近似信号的DTFT。

信号在低频和高频间的幅度差小于2dB,说明信号包含的各频率分量幅度近似。

0~2kHz之间频谱幅度较大。

例8.3 元音“eee”发音的一段发音,采样频率8kHz。

信号很有规律,但不是理想的周期信号,主要由200Hz和400Hz信号构成。

8.3 周期数字信号周期信号是在整个时域按固定间隔重复的信号。

每个重复周期中出现的采样点数称为信号的数字周期。

使用傅里叶级数(DFS)或称离散傅里叶级数来计算其频谱。

数字周期为N的信号,可以表示为正弦或余弦和,可以由欧拉公式展开为复指数形式:n为采样点编号,标号k从0~N-1变化,事实上由于信号的周期性,可以是任意相继的N个值。

c k为傅里叶系数:计算傅里叶系数需要x[n]的任意N个相继的采样值,因此c k具有k个主值,随k值增大,这k个主值不断重复。

这里含有信号幅度信息,由它随k的变化表示信号的幅度频谱;携带相位信息,由它随k的变化表示信号的相位频谱。

此式显示了傅里叶级数所有分量的频率。

由k=0可以获得信号的直流分量|c0|/N,它是信号的平均值。

k>0的其它频率称为周期信号的谐波。

k=1的频率为信号的一次谐波,或称为基频。

大小为f s/N,此频率十分重要,其倒数为信号完成一个完整循环的时间,即NT s,T s为采样间隔。

所有谐波为基频的整数倍。

N个傅里叶级数的系数位于从0直到接近采样频率的频率上。

幅度频谱和相位频谱都具有周期为N的周期性。

幅度频谱是k的偶函数,相位频谱是k的奇函数。

DFS与DTFT有着重要的区别:DTFT产生连续频谱,频谱在所有频率处都有值。

非周期信号的频谱是光滑连续的谱线。

DFS产生的频谱仅有N点的频谱,且仅含有有限个频率。

周期信号的频谱是线谱。

周期N=60的周期信号。

周期信号每60个采样点重复一次。

基频为f s/60Hz。

其它竖线表示其它谐波,位于基频整数倍处。

N-1次谐波位于59f s/60Hz。

并不是所有周期信号都含有全部谐波,有些信号只含有奇次谐波,偶次谐波为0,有些信号仅在某些谐波处的值为0。

例8.4 求如图所示方波信号的幅度频谱和相位频谱。

N=8。

该信号含有直流分量和奇次谐波。

图形表明,它是非周期信号频谱的采样形式P215。

例8.5 求信号x[n]=sin(nп/5)的频谱,采样频率为1kHz。

例8.6 求信号x[n]=1+sin(nп/2)+cos(nп/4) 的频谱。

例8.8 周期为256点的周期方波的频谱。

只含奇次谐波,第3、5、7次谐波的幅度分别为一次谐波幅度的1/3、1/5和1/7,进一步证实了方波可由奇谐正弦波构造出来。

例8.9 数字系统以每秒20k采样点的速率对信号进行采样,得到数字信号为x[n]=cos(nп/6)+cos(nп3/7)+cos(nп2/3),三个子信号的周期分别为12、14、3,整个信号的周期为84。

根据公式,,计算出三个子信号分别为1667Hz、4286Hz和6667Hz。

频谱图上k=7、18、28三条谱线对应这三个子信号,可由式算出。

相关文档
最新文档