第3章 离散傅里叶变换(DFT)西电

合集下载

离散傅里叶变换(DFT)

离散傅里叶变换(DFT)

尾补L-M个零后,再形成第一行的循环倒相序列。
(2) 第1行以后的各行均是前一行向右循环移1位 形成的。 (3) 矩阵的各主对角线上的序列值均相等。
x( L 1) x( L 2) y (0)c x(0) y (1) x(1) x(0) x( L 1) c y (2)c = x(2) x(1) x(0) y ( L 1)c x( L 1) x( L 2) x( L 3) x(1) h(0) x(2) h(1) x(3) h(2) x (0) h( L 1)
主值序列 x(n)
DFT变换对
x(n)的长度为M点,N≥M
N点DFT 变换对
DFT [ x(n)] X (k ) x(n)WNkn
n 0 N 1
WN e
j
2 N
k 0,1,..., N 1 n 0,1,..., N 1
1 N 1 IDFT [ X (k )] x(n) X (k )WN kn N k 0
1 IDFT[ X (k )]N N
N 1
[ x(m)WNmk ]WN kn
k 0 m 0
N 1 N 1
1 x ( m) N m 0
1 N
WNk ( m n )
k 0
N 1
W
k 0
N 1
k ( mn ) N
1 N
e
k 0
N 1 j 2 k ( m n ) N
x(n)
L称为循环卷积区间长度,L≥max[N,M]。
用矩阵计算循环卷积的公式
L 1 yc (n) h(m) x((n m)) L RL (n) m0

第3章 离散傅里叶变换(DFT)C

第3章 离散傅里叶变换(DFT)C

(3.4.9)
def 1 ' 1 ' X (k ) X a f k k X a kF f = T T NT T
p
k 0,1, 2,, N 1
由此可得: ' kF =TX (k ) T DFT[ x(n)] X a N
k 0,1, 2,, N 1
解:
1 1 Tp 0.1 s F 10
因此Tp min=0.1 s。因为要求Fs≥2fc,所以
Tmax
N min
1 1 0.2 103 s 2 f c 2 2500 2 f c 2 2500 500 F 10
第3章 离散傅里叶变换(DFT)
为使用DFT的快速算法FFT,希望N符合2的整数幂,为此 选用N =512点。 为使频率分辨率提高1倍,即F=5 Hz,要求:
说明了X(k)与Xa(jΩ)的关系. 为了符合一般的频谱描述习惯,以频率f为自变量
第3章 离散傅里叶变换(DFT)
令:
X a' ( f ) X a j X a j2πf 2 πf ' 2πf Xa ( f ) X X a a 2 πf
第3章 离散傅里叶变换(DFT)
x ( n) 如果 ~ 的周期预先不知道,可先截取M点进行DFT,即
(n) RM (n) xM (n) x X M (k ) DFT[ xM (n)]
再将截取长度扩大1倍,截取
0 k M 1
(3.4.18)
x (n)的频谱结构,只是在k=im 由此可见,XM(k)也能表示 ~ (i) ,表示 ~ x (n) 的i次谐波谱线,其幅度扩 时,X (im) mX

第3章 离散傅利叶变换(DFT)

第3章 离散傅利叶变换(DFT)
X1[k]和X2[k]分别为x1(n)和x2(n)的N点DFT 说明: (1)如果x1(n)和x2(n)皆为N点,则aX1[k]+bX2[k]也 是N点序列。 (2)若x1(n)和x2(n)的点数不等,则ax1(n)+bx2(n)应 为N=max[N1,N2]点,故DFT必须按N计算。例如,若 N1<N2,则取N=N2,那么需要将x1(n)补上(N2-N1)个 零点后变为N2点的序列,然后作N2点的DFT。
离散性、谐波性
2π τ
X(k)

2π τ
-Ω0 0 Ω0 3Ω0
4π τ
Ω
第三章 离散傅利叶变换(DFT)
2. 连续时间、连续频率——傅利叶变换(FT) 非周期、连续时间信号通过连续付里叶变换(FT) 得到非周期、连续频谱密度函数。
X ( jΩ ) x(t ) 1 2π



x ( t )e
jΩ t
dt
jΩ t



X ( j Ω )e

时域连续函数造成频 域是非周期的谱, 而时域的非周期性造 成频域是连续的谱。
第三章 离散傅利叶变换(DFT)
3. 离散时间、连续频率——序列的傅利叶变换(DTFT)
X (e
j
)
n



x ( n )e
j n
正变换
x(n)
第三章 离散傅利叶变换(DFT)
第三章 离散傅利叶变换(DFT)
引言
• 有限长序列在数字信号处理中占有很重要的地位。 由于计算机容量的限制,只能对过程进行逐段分析。 由于有限长序列,引入DFT(离散付里叶变换)。 • DFT变换除了作为有限长序列的一种付里叶表示,在 理论上重要之外,而且由于存在着计算机DFT的有效 快速算法--FFT,因而使离散付里叶变换(DFT)得以实 现,它使DFT在各种数字信号处理的算法中起着核心 的作用。

第3章 离散傅里叶变换(DFT)

第3章 离散傅里叶变换(DFT)
26
【例3.2.1】计算下面给出的两个长度为4的序列h(n)与
x(n)的4点和8
h(n) h(0), h(1), h(2), h(3) 1, 2,3, 4
解 h(n)与x(n)的4点循环卷积矩阵形式为 x(n) x(0), x(1), x(2), x(3) 1,1,1,1

yc
(1)

2
1
0
0
0
0
4
3 1

3


yc yc yc
(2)
(3)

(4)

3 4 0
2 3 4
1 2 3
0 1 2
0 0 1
0 0 0
0 0 0
4 1 6
0 0
1

0

10 9
yc (0) 1 4 3 2 1 10

yc
(1)


2
1
4
3 1 10

yc yc
(2)
(3)

3 4
2 3
1 2
4 1
1

1
10 10
h(n)与x(n)的8点循环卷积矩阵形式为
yc (0) 1 0 0 0 0 4 3 2 1 1

yc (5) 0 0 4 3 2 1 0 0 0 7

yc
(6)

0
0
0
4
3
2
1
0 0
4

yc (7) 0 0 0 0 4 3 2 1 0 0
h(n)和x(n)及其4点和8点循

第3章 离散傅里叶变换(DFT)

第3章 离散傅里叶变换(DFT)

X(k)与x(n)均为有限长序列,但由于WknN 的周期性,X(k)隐含周 期性, 且周期均为N。 对任意整数m, 总有
k ( WN WNk mN ) , k, m, N
N 1 n 0
均为整数
( 所以,X(k)满足 X (k mN ) x(n )WNk mN ) n kn x(n )WN X (k ) n 0 N 1
k 1 X 1 x n e
n 0
2 1n 4
x n e
n 0
3

x n ( j ) n 2 2 j
n 0
3
k 2
X 2 x n e
n 0 3
3
j
2 2n 4 2 3n 4
x n e j n x n (1) n 2

DFT后的X(k)具周期性,周期为N
x(n)满足
x(n+mN)=x(n)
IDFT后的x(n)具周期性,周期为N
主值区间和主值序列
任何周期为N的周期序列 ~(n) 可以看作长度为N的有限 x
x 长序列x(n)的周期延拓序列, x(n)是 ~(n) 的一个周期。 ~(n) 中n=0到N-1的第一个周期为 ~(n) 的主值区间。 x x x 主值区间上的序列为 ~(n)的主值序列;
x((n))N表示x(n)以N为周期的周期延拓序列,
((n))N表示n对N求余,
如果 则 n=MN+n1, 0≤n1≤N-1, M为整数, ((n))N=n1
--此运算符表示n被N除,商为M,余数为n1。
(n1) 是((n))N 的解,或称作取余数,或称作n对N取模值, 或 简称为取模值,n模N。

第3章 离散傅里叶变换(DFT)及其快速算法(FFT)

第3章 离散傅里叶变换(DFT)及其快速算法(FFT)

1 x(n) IDFT [ X (k )]N N

k 0
N 1
X (k )WN k n , n 0, 1, , N 1
也可以表示为矩阵形式: x DN1 X
DN1
称为N点IDFT矩阵,定义为:
1 1 1 1 W 1 WN 2 N 1 1 WN 2 WN 4 N ( N 1) WN 2( N 1) 1 WN 1 WN ( N 1) 2( N 1) WN WN ( N 1)( N 1)
3.1.3 DFT的矩阵表示
X (k ) DFT [ x(n)]N

n 0
N 1
k x(n)WN n , k 0, 1, , N 1
也可以表示成矩阵形式: X DN x 式中,X是N点DFT频域序列向量:
X [ X (0) X (1) X ( N 2) X ( N 1)]T
2
N 1

k
DFT与DTFT变换
DFT所表示的不是序列的频谱,而是对序列频谱的一个采样! 采样间隔为2/N;N越大,X(k)越能反映X()的形状。
(2)序列的N点DFT是序列的Z变换在单位圆上的N点等间隔采样, 频率采样间隔为2/N。
X (k ) X ( z )
z e
j 2 k N
M 1

n 0
比较前面三式,得到:X (k ) X (e j )
结论:
2 k N
, k 0, 1, 2,, N 1
(1)序列的N点DFT是序列的傅里叶变换(DTFT)在频率区间 [0,2]上的N点等间隔采样,采样间隔为2/N。
X (e j )
X (k )

第三章离散傅里叶变换(DFT)

第三章离散傅里叶变换(DFT)
西北大学信息科学与技术学院 2007年
3.1.1 有限长序列的离散频域表示
我们已学过三种傅里叶分析工具,它们 分别应用于不同性质的信号。
1. 应用于连续周期信号——傅里叶级数展开
j2 kt
xa t Cke T
k
Ck
1 T
T2 -T2

(t
-j2
)e T
kt
dt
其中,T是信号 xa t的周期,Ck 表示了xa (t)的
离散傅里叶变换定义为
X (k)
N 1
x nWNkn
n0
0
0 k N 1 其他
西北大学信息科学与技术学院 2007年
反变换公式为
x(n)
N 1
X
k
W kn N
0 n N 1
k 0
0
其他
DFT是借用了DFS,这样就假定了序 列的周期性,但定义式本身对区间作了强制 约束,以符合有限长特点,这种约束不改变 周期性的实质,或者说,DFT隐含了周期 性。
fc n xn yn
M 1 m0
x
m
y
n m
l
RL
n
M 1 m0
x
m
r
y
n
m
rL
RL
n
r
M 1 m0
x
m
y
n
rL
m
RL
n
r
f
n
rL
Rl
n
西北大学信息科学与技术学院 2007年
圆周卷积fc (n) 等于一个周期序列的主值 序列,该周期序列是线性卷积f (n)以L为周期 进行周期延拓的结果,因此,当L ≥ L1满足 时, fc (n)必然等于f (n),但是,如果L < L1 , 则fc (n)不等于f (n) 。

第3章离散傅里叶变换(DFT)

第3章离散傅里叶变换(DFT)

|X(ejω)| 4
3 2 1 0
0
|X(k)| 5
4 3 2 1 0
0
|X(k)| 5
4 3 2 1 0
0
|X(k)| 5
4 3 2 1 0
0
第3章 离散傅里叶变换(DFT)
(a)R4(n)的幅频特性图
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
(b)4点DFT的幅频特性图
0.5
1
1.5

j 3 k
e8
sin( 2
sin(
k) ,
k)
8
k 0,1, , 7
第3章 离散傅里叶变换(DFT)
设变换区间N=16, 则
X (k)
15
x(n)W1k6n
3
j 2 kn
e 16
n0
n0

j 3 k
e 16
sin( k)
4

,
sin( k)
16
k 0,1, ,15
由此可见,x(n)的离散傅里叶变换结果与变换区间 长度N的取值有关。
第3章 离散傅里叶变换(DFT)
3.1.2 DFT与傅里叶变换(DTFT)和Z变换的关系 设序列x(n)长度为M,其Z变换和N(N≥M)点DFT分别为:
M 1
X (z) ZT[x(n)] x(n)zn
第3章 离散傅里叶变换(DFT)
【例3.1.2】 设x(n)=R4(n),X(ejω)=FT[x(n)]。分别计算 X(ejω)在频率区间[0,2π]上的16点和32点等间隔采样, 并绘制X(ejω)采样的幅频特性图和相频特性图。 解:

第3章--离散傅里叶变换(DFT)

第3章--离散傅里叶变换(DFT)

设x(n)是一种长度为M旳有限长序列, 则定义x(n)旳N点
离散傅里叶正变换为
N 1
j 2 nk
X (k ) DFT[x(n)] x(n)e N
N 1
x(n)WNnk
n0
n0
离散傅里叶逆变换为
离散傅里叶变换对
x(n)
IDFT[ X (k )]
1 N
N 1
j 2 nk
X (k )e N
3.2 离散傅里叶变换旳基本性质
1 线性性质 假如x1(n)和x2(n)是两个有限长序列,长度分别为N1和N2。 y(n)=ax1(n)+bx2(n) 式中a、 b为常数, 即N=max[N1, N2],
则y(n)旳N点DFT为 Y(k)=DFT[y(n)]=aX1(k)+bX2[k], 0≤k≤N-1(3.2.1) 其中X1(k)和X2(k)分别为x1(n)和x2(n)旳N点DFT。 若N1<N2,则N=N2,那么需将x1(n)补上N2-N1个零值点后变
k 2 k f f s k
N
N
以上所讨论旳三种频率变量之间旳关系,在对模 拟信号进行数字处理以及利用模拟滤波器设计数 字滤波器乃至整个数字信号处理中十分主要,望 同学们高度注重。
第三章 离散傅里叶变换DFT
3.1.2 DFT旳隐含周期性------ DFT与 DFS旳关系
DFT变换对中,x(n)与X(k)均为有限长序列,但因为WknN旳周
第三章 离散傅里叶变换DFT
例2 : x(n) R8 (n),分别计算x(n)旳8点、16点DFT。
解: x(n)旳8点DFT为
X (k)
7 n0
R8 (n)W8k n
7 j2k n

数字信号处理:第3章 离散傅里叶变换(DFT)

数字信号处理:第3章  离散傅里叶变换(DFT)

N 1
X (k mN ) x(n)WN(kmN )n
n0
N 1
x(n)WNkn X (k)
n0
同理可证明(3.1.2)式中
x(n+mN)=x(n)
~
实际上, 任何周期为N的周期序列 x 都可以看
作长度为N的有限长序列x(n)的周期延拓序列, 而x(n)
则是
~
x
的一个周期, 即
~
x(n) x(n mN )
式和(3.2.10)式代入得到
x*(N-n)=x*ep(N-n)+x*op(N-n) =xep(n)-xop(n)
xep(n)=1/2[x(n)+x*(N-n)] xop(n)=1/2[x(n)-x*(N-n)]
(3.2.12) (3.2.13) (3.2.14)
2. DFT的共轭对称性 (1) 如果x(n)=xr(n)+jxi(n) 其中 xr=Re[x(n)]=1/2[x(n)+x*(n)] jxi(n)=jIm[x(n)]=1/2[x(n)-x*(n)] 由(3.2.7)式和(3.2.13)式可得 DFT[xr(n)]=1/2DFT[x(n)+x*(n)]
n0
X(k)的离散傅里叶逆变换为
X (k) DFT[x(n)]
1 N
N 1 n0
X (n)WNkn ,
k=0, 1, &, N-1
(3.1.2)
式中,
j 2
eN
, N称为DFT变换区间长度N≥M,
通常称(3.1.1)式和(3.1.2)式为离散傅里叶变换对。 下面
证明IDFT[X(k)]的唯一性。
则 ((n))N=n1
~
例如, N 5, x(n) x(n)5,

第3章离散傅里叶变换(DFT)

第3章离散傅里叶变换(DFT)

(3.1.10)
X (k ) X ((k )) N
第3章 离散傅里叶变换(DFT)
3.2 离散傅里叶变换的基本性质
• 线性性质 • 循环移位性质(序列的圆周移位) • 循环卷积定理 • DFT的共轭对称性
第3章 离散傅里叶变换(DFT)
3.2.1 线性性质
y(n) ax1 (n) bx2 (n)
第3章 离散傅里叶变换(DFT)
x(2)
x(1) x(0)
x(3) N=6 x(4) x(5)
左移 顺时
第3章 离散傅里叶变换(DFT)
2、时域循环移位定理
如果
X (k ) DFT [ x(n)] ,
0 k N 1
y(n) x((n m)) N RN (n)

Y (k ) DFT[ y(n)] WN km X (k ),
论上有重要意义,而且在各种信号的处理中亦起着
核心作用。
第3章 离散傅里叶变换(DFT)
(一)教学内容 1、离散傅立叶变换的定义 2、离散傅立叶变换的基本性质
3、频率域采样
4、DFT的应用举例 (二)基本要求 1、了解离散傅立叶变换概念及物理意义。 2、掌握离散傅立叶变换性质、频率域采样、DFT的应用。
y(n)的N点DFT为
N max{N1 , N2}
0≤k≤ N-1
Y (k ) DFT[ y(n)] aX1 (k ) bX 2 (k )
其中,X1 (k ), X 2 (k ) 分别是 x1 (n), x2 (n) 的N点的DFT 特别注意 注意:如果N1和N2不相等,则以N为DFT变换长度时,其中相对
第3章 离散傅里叶变换(DFT)
3.2.3

第三章 离散傅里叶变换(DFT)

第三章  离散傅里叶变换(DFT)
N 1
~ X ( k ) N k ( r pn)
k 0
N 1
~ NX ( r pN ) ~ NX ( r )
j 2 nr N
1 ~ 因此, X (r ) N
~ ( n )e x
n 0
N 1
将r换成k则有 1 ~ X (k ) N

n 0
则有
~ ~ ~ (n) b~ (n) aX (k ) bX (k ) DFSax1 x2 1 2
其中,a,b为任意常数。
二.序列的移位
~ ~(n) X (k ) 如果 DFSx
则有:
~ ~(n m) W mk X (k ) DFSx N e
2 j mk N
即:
N 1 n 0 j 2 kn N
~ ~( n )e X (k ) x ~( n ) 1 x N
N 1 k 0
~ X ( k )e
2 j kn N
~ X (k ) 的周期性 2 N 1 j ( k mN ) n ~ 周期性: ( k m N) ~( n )e N X x
) X (k )
0
0 20
N 0 N

k
四.离散时间、离散频率的傅氏变换--DFT
x(nT)=x(n)
1 2 T0 F0 0
T0 NT
0
x (e
j k 0T
T 2T
1 2
( N 1) ( N 1)
NT N
0
)
2 T s 1 T 2
x(k )
n 0 N 1 j 2 nk N
~ ( n )W nk x N
N 1 n 0

离散傅里叶变换(DFT)

离散傅里叶变换(DFT)

k=floor((-Nw/2+0.5):(Nw/2+0.5)); %建立关于纵轴对称的频率相量
for r=0:3;
K=3*r+1;
% 1,4,7,10
nx=0:(K*Nx-1); x=xn(mod(nx,Nx)+1);
%周期延拓后的时间向量 %周期延拓后的时间信号x
Xk=x*(exp(-j*dw*nx'*k))/K; %DFS
0
DFT的提出:
离散傅里叶变换不仅具有明确的物理意义,相对于DTFT, 它更便于用计算机处理。但是,直至上个世纪六十年代,由 于数字计算机的处理速度较低以及离散傅里叶变换的计算量 较大,离散傅里叶变换长期得不到真正的应用,快速离散傅 里叶变换算法的提出,才得以显现出离散傅里叶变换的强大 功能,并被广泛地应用于各种数字信号处理系统中。近年来, 计算机的处理速率有了惊人的发展,同时在数字信号处理领 域出现了许多新的方法,但在许多应用中始终无法替代离散 傅里叶变换及其快速算法。
X (e j ) x(n)e jn n
x(n) 1 X (e j )e jnd
2
其中ω为数字角频率,单位为弧度。 注意:非周期序列,包含了各种频率的信号。
局限性:离散时间傅里叶变换(DTFT)是特殊的Z变换,在数学和信号分 析中具有重要的理论意义。但在用计算机实现运算方面比较困难。这是因为, 在DTFT的变换对中,离散时间序列在时间n上是离散的,但其频谱在数字角
§1、傅里叶级数
周期为N的序列 ~x(n) ~x(n rN), (r为整数)
j( 2 )n
基频序列为 e1(n) e N
k次谐波序列为
ek (n)
j( 2 )nk
e N

第三章离散傅里叶变换DFT一

第三章离散傅里叶变换DFT一

(1)
(2)
x1(n) (n 3)
x2 (n)
1 2
(n
1)
(n)
1 2
(n
1)
X1(e j )
(n 3)e jn e j3
n
X 2 (e j )
x2 (n)e jn
n
1 e j 2
1
1 e j 2
1 cos
3.2离散时间序列的傅里叶变换
(3) x3(n) anu(n), 0 a 1 (4) x4 (n) u(n 3) u(n 4)
3.4离散时间周期序列的傅里叶级数 (DFS)
3.4离散时间周期序列的傅里叶级数 (DFS)
2. 时域频域各取一个周期,得到DFT
x(t) X (k0 )e jk0t (1) k
( 0
2
T
)
x(nTs )
x(t) t nTs
X (k0 )e jk0nTs
k
X
(k
0
)e
jk
2 T
nTs
主要内容
3.1连续时间信号的傅里叶变换 3.2离散时间序列的傅里叶变换
(DTFT) 3.3连续时间信号的抽样 3.4离散时间周期序列的傅里叶级数
(DFS)
3.1连续时间信号的傅里叶变换
周期连续信号傅里叶级数展开
周期信号f(t)=f(t+nT) ,满足狄氏条件(有限区间逐 段光滑)时,可展成:
x(n)e N
N n0
(3)
3.4离散时间周期序列的傅里叶级数 (DFS)
习惯上将以上的式(2),(3)中的定标因子移 到反变换中,得到离散傅里叶变换( DFT ):
X
(k)
N 1 n0

第3章 离散傅里叶变换(DFT)09-10-

第3章  离散傅里叶变换(DFT)09-10-

N 1
X (k) DFT[x(n)]
x(n)W
nk N
,0
k
N
1
n0
DFT
x(n)
IDFT [X (k)]
1 N
N 1
X
(k
)W
nk N
k 0
,0
n
N
1
例:x(n)=R4(n) ,求x(n)的8点和16点DFT 解:设变换区间N=8, 则
7
X1(k)
3
x(n)W8kn
j 2 kn
n0
n0
N1
k n N 1
kn kN
x n WN2 x n N WN2 WN 2
n0
n0
N 1
x
kn
n WN2
1 e jk
n0
2
X
k 2
,
k 偶数
0, k 奇数
0 k 2N -1
~
~
N 1 ~
j 2 kn
X (k) DFS[x(n)] x(n)e N

yn, 试x求nY(NkR)=2DNFnT[y(n)]与
X(k)之间的关系。
解:
2 N 1
2 N 1
Yk
ynW2kNn
xnN
R2
N
n
W kn 2N
n0
n0
N 1
2 N 1
xnW2kNn xnW2kNn
n0
nN
N 1
N 1
x
n
W kn 2N
x n N W2kNnN
序列的DFS级数系数的主值序列!
§3.2 离散傅里叶变换的基本性质
一. 线性性质

第3章 3.1-3.2离散傅里叶变换(DFT)

第3章  3.1-3.2离散傅里叶变换(DFT)

n0
WNkm X (k)
第3章 离散傅立叶变换(DFT)
对比记忆:
循环时移:
x((n
m))
N
RN
(n)
W mkm N
X(k
)
线性时移:
x(n n0 ) e jn0 X(e j )
29
时域移位,频域相移
2020/4/5
第3章 离散傅立叶变换(DFT)
3. 频域循环移位定理 如果: X (k) DFT[x(n)], 0 k N 1 则 : Y (k) X ((k l))N RN (k)
e8
n0
n0
j 3k
e8
sin(
2
sin(
k) k)
,k
0,1,, 7
8
17 2020/4/5
第3章 离散傅立叶变换(DFT)
提高谱密度
18
图3.1.1 R4(n)的FT和DFT的幅度特性关系
2020/4/5
第3章 离散傅立叶变换(DFT)
3.3.2 DFT和DTFT、ZT的关系
设序列x(n)的长度为N, 其ZT、DTFT和
对任意整数m, 总有:
WNk WN(kmN) , k, m, N均为整数
所以(3.3.6)式中, X(k)满足:
N 1
X (k mN ) x(n)WN(kmN )n
n0
N 1
x(n)WNkn X (k)
n0
同理可证明(3.3.7)式中:
14 2020/4/5
x(n mN) x(n)
1.
设序列h(n)和x(n)的长度分别为N和M。h(n)与x(n)的
L点循环卷积定义为:L1
kn
e4
n0

(n) 第3章离散傅里叶变换(DFT)

(n) 第3章离散傅里叶变换(DFT)

X (e j ) DTFT [ x(n)] x( n)e j n
比较上面三式可得关系式
X (k ) X ( z )
z e
j 2 k N
, ,
0 k N -1 0 k N -1
(3.1.3) (3.1.4)
X
X (k ) X (e j )
2 k N
X
第3章 离散傅里叶变换(DFT)
重要公式:
1 N
e
k 0
N 1 j 2 k ( m n ) N
1 0
m n MN , M为整数 m n MN , M为整数
X
第3章 离散傅里叶变换(DFT)
例 3.1.1 x(n)=R4(n) ,求x(n)的8点和16点DFT 。 解: (1)变换区间N=8, 则
kn n 0 N 1
(3.1.1)
X(k)的离散傅里叶逆变换为
1 N 1 kn x(n) IDFT [ X (k )] X (k )WN , n 0,1,...... N 1 (3.1.2) N k 0 2 j 式中,WN e N
N称为DFT变换区间长度,N≥M.
X
第3章 离散傅里叶变换(DFT)
3.2.2
循环移位性质
1. 序列的循环移位 设x(n)为有限长序列, 长度为N, 则x(n)的循环 移位定义为 y(n)=x((n+m))NRN(n) (3.2.2)
X
第3章 离散傅里叶变换(DFT)
x(n)
1
0.5
n
0 -2 0 2 4 0 1 2 3 4 5 6 7
X
第3章 离散傅里叶变换(DFT)
左图为4点矩形 信号的Z变换, 红色的曲线为4 点矩形信号的 DTFT。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

~
实际上, 任何周期为N的周期序列
则是 x
~
x
都可以看
作长度为N的有限长序列x(n)的周期延拓序列, 而x(n)
~
的一个周期, 即
x(n)
m


x(n mN )
(3.1.5) (3.1.6)
x(n ) x(n ) RN (n )
~
~
为了以后叙述方便, 将(3.1.5)式用如下形式表示:
x (n ) x1 (n ) x2 (n ) x1 (m) x2 ((n m)) N RN (n )
m 0 N 1
由于
X (k ) DFT [ x(n )] X 1 (k ) X 2 (k ) X 2 (k ) X 1 (k )
所以
x (n ) IDFT [ X (k )] x1 (n ) x2 (n ) x2 (n ) x1 (n )
X(k)的离散傅里叶逆变换为
1 X (k ) DFT [ x(n)] N

n 0
N 1
X ( n)WN kn , k=0, 1, &, N-1 (3.1.2)
式中, e , N称为DFT变换区间长度N≥M, 通常称(3.1.1)式和(3.1.2)式为离散傅里叶变换对。 下面 证明IDFT[X(k)]的唯一性。 把(3.1.1)式代入(3.1.2)式有
第3章 离散傅里叶变换(DFT)
3.1 离散傅里叶变换的定义 3.2 离散傅里叶变换的基本性质 3.3 频率域采样 3.4 DFT的应用举例
3.1 离散傅里叶变换的定义
3.1.1 DFT的定义 设x(n)是一个长度为M的有限长序列, 则定义x(n) 的N点离散傅里叶变换为
kn X (k ) DFT [ x (n )] x (n )WN , k=0, 1, &, N-1 (3.1.1) n 0 N 1
即循环卷积亦满足交换律。
作为习题请读者证明频域循环卷积定理:
如果 则 x(n)=x1(n)x2(n)
1 X (k ) DFT [ x ( n )] X 1 ( k ) X 2 ( k ) N 1 N 1 X 1 (l ) X 2 (( k l )) N RN (k ) N l 0 1 X (k ) X 2 (k ) X 1 (k ) N 1 N 1 X 2 (l ) X 1 (( k l )) N RN (k ) N l 0
n 0 N 0
j
e
3 j k 8
sin( k ) 2 , k 0,1, ,7 sin( k ) 8

设变换区间N=16, 则
X (k ) x(n )W8kn e
n 0 N 0
7
3
j
2 kn 8
e
3 j k 8
sin( k ) 4 , k 0,1, ,15 sin( k ) 16
m
0 1 2 3 4 5 6 7 x2 ((2 -m))NRN(m) 1
m
0 1 2 3 4 5 6 7 x(n) 4 3 2 1 0 1 2 3 4 5 6 7
m
n
图3.2.2 循环卷积过程示意图
3.2.5 DFT的共轭对称性
1. 有限长共轭对称序列和共轭反对称序列 为了区别于傅里叶变换中所定义的共轭对称(或共 轭反对称)序列, 下面用xep(n)和xop(n)分别表示有限长 共轭对称序列和共轭反对称序列, 则二者满足如下定 义式: xep(n)=x*ep(N-n), xop(n)=-x*op(N-m), 0≤n≤N-1 0≤n≤N-1 (3.2.9) (3.2.10)
j
2 N
1 IDFT [ X (k )] N
m 0 N 1

k 0
N 1
mk [ x(m)WN ]WN kn m 0
N 1
1 x ( m) N

k 0
N 1
k WN ( mn )
1 N
W
k 0
N 1
k ( m n ) N源自{01m n MN , M M为整数 m n MN , M M为整数
则有
x (5) x ((5))5 x(0) x (6) x ((6))5 x (1)
所得结果附合图2.1.2所示的周期延拓规律。
~
~
如果x(n)的长度为N, 且

~
~
x (n)的离散傅里叶级数表示为
N 1 n 0 ~ kn N
~
x
(n)=x((n))N, 则可写
X (k ) x(n)W
3.2.2 循环移位性质
1. 序列的循环移位 设x(n)为有限长序列, 长度为N, 则x(n)的循环移 位定义为 y(n)=x((n+m))NRN(N) (3.2.2)
图 3.2.1
循环移位过程示意图
2. 时域循环移位定理
设x(n) 是长度为N的有限长序列, y(n)为x(n)的循 环移位, 即 y(n)=x((n+m))NRN(n) 则
所以, 在变换区间上满足下式:
IDFT[X(k)]=x(n), 0≤n≤N-1
由此可见, (3.1.2)式定义的离散傅里叶变换是唯一的。 例 3.1.1 x(n)=R4(n) ,求x(n)的8点和16点DFT 设变换区间N=8, 则
7 3 2 kn 8
X (k ) x (n )W8kn e
Y (k ) W
km N

n
N 1
kn x ((n)) N WN
W
km N
n 0

N 1
kn x (n)WN
WN km X (k )
3. 频域循环移位定理如果
X(k)=DFT[x(n)], 0≤k≤N-1
Y(k)=X((k+l))NRN(k) 则 y(n)=IDFT[Y(k)]=WnlNx(n) (3.2.4)
比较上面二式可得关系式
X (k ) X ( z )
z e
j
2 k N
, ,
0 k N-1 0 k N-1
(3.1.3) (3.1.4)
X ( k ) X ( z j )
2 k N
图 3.1.1 X(k)与X(e jω)的关系
3.1.3 DFT的隐含周期性
前面定义的DFT变换对中, x(n)与X(k)均为有限长 序列, 但由于WknN的周期性, 使(3.1.1)式和(3.1.2)式 中的X(k)隐含周期性, 且周期均为N。 对任意整数m, 总有
DFT [ x (n )]
又由X(k)的隐含周期性有X(N)=X(0) 用同样的方法可以证明 DFT[x*(N-n)]=X*(k) (3.2.8)
1
0 1 2 3 4 5 6 7 x2 (n) 1
n,m
0 1 2 3 4 5 6 7 x2 ((- m))NRN(m) 1
n
0 1 2 3 4 5 6 7 x2 ((1 -m))NRN(m) 1
k ( WN WNk mN ) , k , m, N 均为整数 所以(3.1.1)式中, X(k)满足
( X (k mN ) x(n )WNk mN ) n n 0 kn x(n )WN X (k ) N 1 n 0 N 1
同理可证明(3.1.2)式中
x(n+mN)=x(n)
3.2.3 循环卷积定理
有限长序列x1(n)和x2(n), 长度分别为N1 和N2 , N=max[ N1, N2 ]。 x1(n)和x2(n)的N点DFT分别为: X1(k)=DFT[x1(n)] X2(k)=DFT[x2(b)]
如果
X(k)=X1(k)·X2(k) 则
x(n) IDFT [ X (k )] x1 (m)((n m)) N RN (n) x(n) IDFT [ X (k )] x2 (m)((n m)) N RN (n)
3.2 离散傅里叶变换的基本性质
3.2.1 线性性质
如果x1(n)和x2(n)是两个有限长序列, 长度分别为 N1和N2。
y(n)=ax1(n)+bx2(n)
式中a、 b为常数, 即N=max[N1, N2], 则y(n) 的N点DFT为 Y(k)=DFT[y(n)]=aX1(k)+bX2[k], 0≤k≤N-1(3.2.1) 其中X1(k)和X2(k)分别为x1(n)和x2(n)的N点DFT。
x(n) x(n) N
(3.1.7)
图 3.1.2 有限长序列及其周期延拓
式中x((n))N表示x(n)以N为周期的周期延拓序列,
((n))N表示n对N求余, 即如果 n=MN+n1, 0≤n1≤N-1, 则 ((n))N=n1 例如, N 5, x(n) x(n)5 ,
~
M为整数,
N 1
kn′ , N
以N为周期, 所以对
其在任一个周期上求和的结果不变。 因此
kn X (k ) x1 (m)WN m 0
X 1 (k ) X 2 (k ),
0 k N 1
循环卷积过程中, 要求对x2(m)循环反转, 循环移 位, 特别是两个N长的序理的循环卷积长度仍为N。 显 然与一般的线性卷积不同, 故称之为循环卷积, 记为

3.1.2 DFT和Z变换的关系
设序列x(n)的长度为N, 其Z变换和DFT分别为:
X ( z ) ZT [ x (n )] x (n ) z n
n 0 kn X (k ) DFT [ x (n )] x (n )WN n 0 N 1
N 1
0 k N-1
~
x((n)) N W
相关文档
最新文档