第二讲 Part3 离散傅里叶变换_难点
信息与通信第3章离散傅里叶变换PPT课件

n
~x(n)
1
(b)
0123456 7
n
~x(n)
x(n rN ) x(n%N) x((n))N
r
x(n) ~x (n)RN (n) x((n))N RN (n)
x(n)为周期序列的主值序列
第4页/共46页
| X~(k)|
(c)
01 2 3 4 5 6 7
k
| X(k)|
(d)
01 2 3 4 5 6 7
四、用MATLAB计算序列的DFT
• Xk = fft(xn,N) • xn = ifft(Xk,N)
第12页/共46页
【例3.1.2】
(a)16点 DFT的 幅 频 特 性 图 4
3
(b)16点 DFT的 相 频 特 性 图 2
相位
幅度
2
0
1
0
0
0.5
1
1.5
2
/
(c)32点 DFT的 幅 频 特 性 图 4
W ( N 1)1 N
W 02 N
W 12 N
W ( N 1)2 N
W 0( N 1) N
W 1( N 1) N
x(0) x(1)
WN(
N
1)(
N
1)
x(N 1)
x(0) x(1)
1 N
W 00 N
W 10 N
x(N 1)
WN(N 1)0
W 01 N
W 11 N
1
n0
DFT的物理意义2:X(k)为x(n)的傅里叶变换在区间[0, 2π] 上的N点等间隔采样。
第9页/共46页
三、DFT的隐含周期性
N 1
离散傅里叶变换ppt

频域信号 周期的 离散的
*时域是周期为Tp函数,频域的离散间隔为0
2
Tp
;
时域的离散间隔为T ,频域的周期为s
2
T
.
§ 3-1 周期序列的DFS
一.周期序列DFS的引入 导出周期序列DFS的传统方法是从连
续的周期信号的复数傅氏级数开始的:
~x (t) X~ ( jk0 )e jk0t k
对上式进行抽样,得:
n0
x(n)
IDFT X (k)
1 N
N 1
X (k )WNn,k
k 0
0nN-1
或者: X (k) X~(k)RN (k) x(n) ~x (n)RN (n)
练习题
参考答案
TP 1/ f 0.1(s) T 1/ 2 fh 1/ 8kHz 0.125(ms) N 2 fh / f 800
证明:
DFS[WNmn~x (n)]
N
1
WNmn
~x (n)WNkn
n0
N 1 ~x (n)WN(km)n n0
X~(k m)
WNmn
j 2 mn
eN
j 2 nm
eN
(e
j
2
N
n
)
m
时域乘以虚指数(
j 2
eN
n
)的m次幂,频域搬移m,调制特性。
四.周期卷积和
1.如果 Y~(k) X~1(k)X~2(k)
所以
DFS[~x (n
m)]
N 1m~x (i)WNik
W mk N
im
W mk N
N
1
~x (i)WNik
W mk N
~x (k
离散时间傅里叶变换

离散时间傅⾥叶变换1. 离散时间傅⾥叶变换的导出针对离散时间⾮周期序列,为了建⽴它的傅⾥叶变换表⽰,我们将采⽤与连续情况下完全类似的步骤进⾏。
考虑某⼀序列x[n],它具有有限持续期;也就是说,对于某个整数N1和N2,在 −N1⩽以外,x[n]=0。
下图给出了这种类型的⼀个信号。
由这个⾮周期信号可以构成⼀个周期序列\tilde x[n],使x[n]就是\tilde x[n]的⼀个周期。
随着N的增⼤,x[n]就在⼀个更长的时间间隔内与\tilde x[n]相⼀致。
⽽当N\to \infty,对任意有限时间值n⽽⾔,有\tilde x[n]=x[n]。
现在我们来考虑⼀下\tilde x[n]的傅⾥叶级数表⽰式\tag{1}\tilde x[n] = \sum_{k=(N)}a_ke^{jk{(2\pi/N)}n}\tag{2}a_k = \frac{1}{N} \sum_{n=(N)} \tilde x[n]e^{-jk{(2\pi/N)}n}因为在-N_1 \leqslant N \leqslant N_2区间的⼀个周期上\tilde x[n]=x[n],因此我们将上式的求和区间就选在这个周期上\tag{3}a_k = \frac{1}{N} \sum_{n=-N_1}^{N_2} x[n]e^{-jk{(2\pi/N)}n} = \frac{1}{N} \sum_{n=-\infty}^{+\infty} x[n]e^{-jk{(2\pi/N)}n}现定义函数\tag{4}X(e^{j\omega})=\sum_{n=-\infty}^{+\infty}x[n]e^{-j\omega n}可见这些系数a_k正⽐于X(e^{j\omega})的各样本值,即\tag{5}a_k = \frac{1}{N}X(e^{jk\omega_0})式中,\omega_0=2\pi/N⽤来记作在频域中的样本间隔。
离散傅里叶变换(DFT)63244ppt课件

令k=m
n0
精选ppt
ak
1
N1
j2kn
x(n)e N
Nn0
7
第3章 离散傅里叶变换(DFT)
ak
1
N1
j2kn
x(n)e N
Nn0
令X(k) Nak
N 1
j 2 nk
X (k) x(n)e N
n0
X
(k
N
)
N
1
x(n)e
j
2 N
n(k
N
)
N 1
j 2 nk
x(n)e N
X (k)
n0
(2)频域周期序列看作是有限长序列X(k)的周期延拓
(3)把周期序列DFS的定义式(时域、频域)各取主值 区间,就得到关于有限长序列时频域的对应变换对。
(前面已证:时域上周期序列的离散傅里叶级数在频域上仍是同 周期序列)
精选ppt
10
第3章 离散傅里叶变换(DFT) (1)周期序列的主值区间与主值序列
FT
X ( j) x(t)e jtdt
x(t) 1 X ( j)e jtd
2
时域非周期频域连续。频谱特点:连续非周期谱
精选ppt
3
第3章 离散傅里叶变换(DFT)
(3)离散非周期信号
X (e j ) x(e j )e jnd
2
1. 连续非周期 2. 连续周期 3. 离散非周期 4. 离散周期
连续非周期() FT
离散非周期 () FS
连续周期( ) DTFT
离散周期
DFS
切实理解四种FT之间的对应关系
精选ppt
5
第3章 离散傅里叶变换(DFT)
离散时间傅立叶变换(DTFT)

| X (e j ) | sin(N / 2) sin( / 2)
arg[ X (e j )] (N 1) arg[sin(N / 2)]
2
sin( / 2)
当N=4时,序列x(n)及其幅度谱与相位谱如下图示。
程序清单
clc; clear; y=[1 1 1 1]; x=0; n=[0:3]; w=0:0.01:2*pi; subplot(311); stem(n,y); xlabel('n'); ylabel('x(n)'); for n=0:3
xe (n) xe (n)
xo (n) xo(n)
xe (n)
1 2
[x(n)
x(n)]
xo (n)
1 2
[x(n)
x(n)]
(4)对序列x(n)旳X(ejω)
X(ejω)=Xe(ejω)+Xo(ejω)
Xe(ejω)=X*e(e-jω) Xo(ejω)=-X*o(e-jω)
X e (e j
)
对比上面两公式, 左边相等, 所以得到 xer(n)=xer(-n) xei(n)=-xei(-n)
(2)共轭反对称序列: 若满足下式: xO(n)=-x*O(-n) 则称xO(n)为共轭反对称序列。
共轭反对称序列旳性质:实部是奇函数, 虚部是偶函数。
例:共轭对称序列 共轭反对称序列
5-j -5+j
d
5、时域卷积定理
设
y(n)=x(n)*h(n),
则 Y(ejω)=X(ejω)·H(ejω)
时域卷积, 频域乘法
证明:
令k=n-m
y(n) x(m)h(n m)
m
Y (e j ) FT[ y(n)]
第3章 离散傅里叶变换 (2)

中国矿业大学信息与电气工程学院
二、DFT和Z变换的关系
进行 对比
X ( z ) ZT[ x( n)] x( n) z n
n 0
N 1
X (k ) DFT[ x( n)] x( n)WN
n 0
N 1
nk
经过抽样、截断和延拓后,信号时域和频域都是离散、周期的。
中国矿业大学信息与电气工程学院 3
学 习 方 法
从工程需要出发,理解信号频谱分析的实际问题。即
在实践中领悟处理原理的意义
从解决问题出发,理解各种信号处理方法的目的。即
在矛盾中思考工程实现的背景
在解决问题的过程中感受知识的力量、体会学习的快乐
N 1 k 0
j( ) nk ~ N X (k )e
2
~ 式中,乘以系数1/N是为了下面计算的方便。X (k ) 为k次 谐波的系数。 2 j 将上式两边同乘 e N rn ,并从n=0到N-1求和,得到:
N 1 n 0
~(n)e x
j(
2 ) rn N
1 N
N 1 N 1
X (k ) X R (k ) jX I (k )
或 X (k ) X (k ) e j ( k )
例3.1 求有限长序列的DFT,其中a=0.8,N=8。
a n , 0 n N 1 x(n) 其它 0,
解
X ( k ) x( n)W8nk a n e
中国矿业大学信息与电气工程学院
4
本章学习要点
理解傅里叶变换的四种形式的意义 了解离散傅里叶级数(DFS)的定义、基本性质 掌握离散傅里叶变换(DFT:Discrete Fourier Transform) 的定义、基本性质以及与Z变换和DTFT的关系,理解隐含 周期性的意义,掌握圆周卷积的计算 掌握频域采样理论的意义、分析过程和结论 掌握DFT在计算线性卷积、线性相关和谱分析等方面的应用
chapter3 离散时间傅里叶变换----DTFT.

a 1
n jn n jn a e a e n0 1
n
n
其中
n n jn n jn a e a e n 1 1
ae j 1 ae j
所以
2 j 1 a ae 1 X (e j ) j j 1 2a cos a 2 1 ae 1 ae
二、离散时间傅里叶变换的举例
1、单边指数序列 于是
X (e )
j n
x ( n)
a0
x(n) a u(n)
n
a 1
0
n jn
x ( n)e
DTFT
jn
a e
n0
1 1 ae j
0
1 2 3 4 5
n
x ( n)
a0
1 2
即
1 a u (n) 1 ae j
3、矩形窗序列 x(n) RN (n) u(n) u(n N )
X (e j )
jn jn x ( n ) e e n0 N 1
x(n) R4 (n)
1
n
0
1 2 3 4 5
n
1 e X (e ) 1 e j
j
jN
当x(n)是实序列,即 则
x(n) x* (n)
X (e j ) X * (e j )
X (e j ) e j( ) X R () jX I () X R () jX I () X (e j ) e j( )
即实序列的离散时间傅里叶变换,实部是偶对称的,虚部是奇 对称的,模是偶对称的,相位是奇对称的。 当x(n)是实偶序列,即
《离散傅里叶变换》课件

离散傅里叶级数
探索离散傅里叶级数的定义、性 质和计算方法以及在数字信号处 理中的应用。
离散傅里叶变换
仔细研究离散傅里叶变换的离散 性质和变换公式,揭示其在信号 分析中的独特优势。
离散傅里叶变换的性质
探索离散傅里叶变换的对称性、 线性性以及快速计算算法,解开 其工程应用的奥秘。
离散傅里叶变换实践1海明窗函数图像处理
探索离散傅里叶变换在图像滤波、增强和压缩中的重要作用。
视频编码
揭示离散傅里叶变换在视频编码和压缩领域的关键应用和优化策略。
总结
离散傅里叶变换的优点与缺点
离散傅里叶变换未来的发展趋势
2
深入了解海明窗函数的定义和特性,以
及在信号处理中的应用场景。
3
快速傅里叶变换算法
介绍快速傅里叶变换算法的基本原理和 实现方法,让你轻松掌握高效算法的使 用。
离散傅里叶变换与信号处理实例
通过实际案例演示离散傅里叶变换在语 音信号和图像信号处理中的应用与效果。
离散傅里叶变换应用
语音信号处理
深入研究离散傅里叶变换在语音信号分析、压缩和合成中的广泛应用。
《离散傅里叶变换》PPT 课件
本课件介绍离散傅里叶变换(Discrete Fourier Transform),让你轻松理解 该概念及其应用。从基本理论到实践应用,一网打尽。
简介
什么是离散傅里叶变换
深入探索离散傅里叶变换的定义、原理和作用,为你打开全新的数学世界。
应用领域
探索离散傅里叶变换在信号处理、图像处理、视频编码等领域的广泛应用。
傅里叶理论基础
1 傅里叶级数
揭秘傅里叶级数的概念和 原理,了解它在周期信号 分析中的作用。
2 傅里叶变换
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第三讲 Part3 DFT 的理论难点
1、抽样定理
连接离散信号与连续信号的桥梁。
()(){
()()j t a a j j n
s
n X j x t e dt
X e x nT e
ω
ω∞
-Ω-∞
∞
-=-∞
Ω==
⎰∑
根据频域卷积定理推导 ()
()()()
{1()()()()()2j j j j j y n x n h n Y e X e H e X e H e d πωωωθωθπ
θ
π--==*=⎰ 得到:1
()()j a
s k s
X e X
j jk T ω
∞
=-∞
=
Ω-Ω∑
2、FT 中的待研究的理论难点与关键之处
2.1 DFT 与DTFT 的关系
两种论述方法:
方法1:书P119-P120的论述;请同学看书后,上黑板叙述推演相关的过程。
方法2:书P121,连续频谱的抽样也必然使原来的时域信号变成周期的。
2.2 DFT 的()X k 是“()x n 的傅里叶变换”的某种程度上的近似。
用DFT 对连续信号和离散信号进行谱分析的基本原理和方法
2.2.1 怎样理解DFT 对FT 的近似?
由于用DFT 对连续信号做频谱分析的过程中隐含了频域和时域的两个周期延拓,又由于信号时宽和带宽的制约关系,因此,做DFT 得到的()N X k ,及由()N X k 做IDFT 得到的
()N x n 都是对原()a X j Ω及()a x t 的某种近似。
如果s T 选得足够小,则式1
()|()s j a T a
s l s
X e X
j jl T ω
ω∞
=Ω=-∞
=
Ω-Ω∑ 中将避免或大大减轻
频域的混叠。
如果N 选得足够大,一方面可以减轻式()()*()j j j a X e X e D e ω
ω
ω
=的窗口效应,另一方面也会减轻式()(),0,1, (1)
l x n x
n lN n N ∞
=-∞
=
+=-∑的时域混叠。
结论:在这两个条件均满足的情况下,上述的近似误差将减小到可接受的程度,从而
使()N x n 和()N X k 都是()a x t 和()a X j 的极好近似。
如何理解上述的陈述与结论?
DFT 对连续信号进行谱分析必然是近似的,其近似程度与信号带宽、采样频率和截取长度有关?
课本2之P123 “DFT 的图形解释”。
2.2.2另外一种近似的解释(西交大):
设连续信号()a x t 持续时间为p T ,最高频率为c f 。
(符合爆振信号的特点) 信号()a x t 的频谱分析是2()[()]()j ft a a a X jf FT x t x t e dt π∞
--∞
==
⎰
,
对()a x t 以采样间隔为s T 采样得()()a s x n x nT =,设有N 个采样点,s t nT = 并对()a X jf 作零阶近似,得到:1
20
()()s
N j fnT s
a
s
n X jf T x nT e
π--==∑
对f 的连续周期函数()a X jf 在区间[0,]s f 上等间隔采样N 点,采样间隔为F 。
参数,,,s p f T N F 满足关系式1
s s
f F N NT =
=
,另有s p NT T = 所以有1p
F T =
而且f kF =,代入得到:21
()(),01N j
kn N
s
a
s
n X jkF T x nT e
k N π--==≤≤-∑
令()(),()()a a s X k X jkF x n x nT ==,则有21
()()[()]N j kn N
a s s n X k T x n e
T DFT x n π--===∑
反过来,逆变换有1
()[()]a s
x n IDFT X k T =
结论1:()a X k 看不到()a X jf 的全部频谱特性,只能看到N 个离散采样点的谱特性,这就是所谓的栅栏效应。
结论2: 如果()a x t 持续时间无限长,进行截断处理,则会产生频率混叠和泄漏现象,使谱分析产生误差。
2.2.3
02/2/()()()()()()..................................................()()()*()s
s s
a a a N t nT j j j a a a T NT x t x n x n d n x n x n FT DTFT DTFT DFT DFS
x j X e X e D e ωωωππ=Ω=Ω=−−−−→−−−→−−−→−−→←−−−−Ω−−−−→−−−→s 周期延拓
抽样截短取一个周期周期延拓卷积抽样()()N X k X k −−−−→−−−−−→←−−−−周期延拓取一个周期
抽样定理:()s x nT 与()x t 的关系?()j X e ω
和()X j Ω的关系。
此乃数字信号处理中的基本问题。
()()j t a a X j x t e dt ∞
-Ω-∞
Ω=⎰,这里2f πΩ=为角频率,()a X j Ω为频谱密度。
()()j j n
s
n X e x nT e
ω
ω∞
-=-∞
=
∑
/1
()()|()s j s T a
s k s
X e X j X
j jk T ω
ω∞
Ω==-∞
=Ω=
Ω-Ω∑,本式的推导见P116-P117,即周期延拓
相对频率Ω,周期为2/2s s s T f ππΩ==; 相对圆频率ω,周期为2π。
变成周期的方法是将()a X j Ω在频率轴上以s Ω为周期移位后再叠加,并除以s T ,这种现象
称为频谱的周期延拓。
频域的“混叠”现象:一个周期中的()s X j Ω不等于()a X j Ω。
所谓混叠是指()j D e ω的主瓣宽度4/B N π=过大,无法区分欲区分的()j X e ω谱线? 混叠与泄漏是一回事吗?
2.3 DFT 泄漏
泄漏的定义?
对确定频率的正弦信号进行频谱分析,按理其频谱图上只有确定频率处有一根谱线,但实际上由于截取有限长度的信号,形成在其频谱图上除主要频率分量外,还出现了许多附加频率分量,造成能量不是集中于确定频率,部分能量泄漏到其他频率上。
从而给傅里叶变换带来误差,这种误差称为泄漏误差。
(1)对周期信号的整周期截取。
(2)对周期信号的非整周期截取。
(3)对任意信号的有限截取。
谱间干扰的定义?
例子:用DFT 计算理想低通滤波器频响曲线。
2.4 频率分辨率
信号()T x t 的长度为T 秒,()T X j Ω得频率分辨率1/f T ∆=
/s M T T =
()M x n 是无穷长离散信号()x n 和一宽度为M 的矩形窗相乘的结果,频率分辨率限制
为:/s f f M ∆=
DFT 的两根谱线间的距离为/s f f N ∆=
2.5 加窗
为了减少由截断所造成的泄漏误差,选取某些比矩形窗函数泄漏小的其他形状的窗函数,称为加窗处理。
减少泄漏与提高分辨率是矛盾的;
(1)要求精确读出峰值的频率而不考虑幅值的精度,则可直接截断,即矩形窗。
(2)对存在强干扰的窄带信号,选用旁瓣幅度小的窗函数。
2.6 补零
补零的作用?
数据过短时泄漏将严重影响对原频谱的辨认,为什么补零可在一定程度上克服这一现象?
实际例子描述。
Part2:DFT 的实际工程应用展示
(1) 通过KNOCK 模式识别应用案例,展示理论到实际的应用过程
领域专业背景知识—〉核心问题的理论实质—〉具体算法设计—〉完整的项目实现
(2) 工程应用中揭示的理论难点分析,可以改进的余地。
(3) DFT 应用背景
背景1:以卷积和相关运算的实现为依据;
背景2:以DFT 作为连续傅里叶变换的近似,对连续信号或序列进行谱分析。
Part3 :FFT 文档解释
3.1 概述
DFT 和卷积是信号处理中两个最基本也是最常用的运算。
21
()(),0,1,...,1,N j
nk
N
N
N n X k x n
W k N W e
π
--===
-=∑。