第三章离散傅立叶变换.
数字信号第三章 离散傅里叶变换
第三章离散傅里叶变换DFT: Discrete Fourier Transform第三章学习目标z理解傅里叶变换的几种形式z掌握离散傅里叶变换(DFT)及性质,圆周移位、共轭对称性,掌握圆周卷积、线性卷积及两者之间的关系z掌握频域抽样理论z掌握DFT的应用引言DFT要解决两个问题:一是频谱的离散化;二是算法的快速计算(FFT)。
这两个问题都是为了使计算机能够实时处理信号。
Fourier变换的几种可能形式时间函数频率函数连续时间、连续频率—傅里叶变换连续时间、离散频率—傅里叶级数离散时间、连续频率—序列的傅里叶变换可以得出一般的规律:一个域的离散对应另一个域的周期延拓;一个域的连续必定对应另一个域的非周期。
−jwndw e jwn 时域离散、非周期频域连续、周期z 时域周期化→频域离散化z 时域离散化→频域周期化离散连续周期性非周期性引言Fourier变换的几种可能形式时间函数频率函数连续时间、连续频率—傅里叶变换连续时间、离散频率—傅里叶级数离散时间、连续频率—序列的傅里叶变换离散时间、离散频率—周期序列的傅里叶级数由DTFT到DFS离散时间、离散频率的傅立叶级数(DFS)由上述分析可知,对DTFT,要想在频域上离散化,那么在时域上必须作周期延拓。
对长度为M的有限长序列x(n),以N为周期延拓(N≥M)。
注意:周期序列的离散傅里叶级数(DFS)只对有限长序列作周期延拓或周期序列成立。
……四种傅里叶变换形式的归纳时间函数频率函数连续和非周期非周期和连续连续和周期(T0)非周期和离散(Ω=2π/T)离散(T)和非周期周期(Ωs=2π/T)和连续离散(T)和周期(T0)周期(Ωs=2π/T)和离散(Ω=2π/T)在进行DFS 分析时,时域、频域序列都是无限长的周期序列周期序列实际上只有有限个序列值有意义长度为N 的有限长序列可以看成周期为N 的周期序列的一个周期(主值序列)借助DFS 变换对,取时域、频域的主值序列可以得到一个新的变换—DFT ,即有限长序列的离散傅里叶变换3.1 离散傅里叶变换(DFT )的定义及物理意义——有限长序列的离散频域表示x(n)的N 点DFT 是¾x(n)的z 变换在单位圆上的N 点等间隔抽样;¾x(n)的DTFT 在区间[0,2π)上的N 点等间隔抽样。
第三章 离散傅立叶变换
第三章 离散傅立叶变换一、离散傅立叶级数计算题:1.如果)(~n x 是一个周期为N 的周期序列,那么它也是周期为2N 的周期序列。
把)(~n x 看作周期为N 的周期序列有)(~)(~1k X n x ↔(周期为N );把)(~n x 看作周期为2N 的周期序列有)(~)(~2k X n x ↔(周期为2N );试用)(k X 1~表示)(k X 2~。
二、离散傅立叶变换定义填空题2.某DFT 的表达式是∑-==10)()(N k kl M Wk x l X ,则变换后数字频域上相邻两个频率样点之间的间隔是( )。
3.某序列DFT 的表达式是∑-==10)()(N k kl M W k x l X ,由此可看出,该序列的时域长度是( ),变换后数字频域上相邻两个频率样点之间隔是( )。
4.如果希望某信号序列的离散谱是实偶的,那么该时域序列应满足条件( )。
5.采样频率为Hz F s 的数字系统中,系统函数表达式中1-z 代表的物理意义是 ),其中时域数字序列)(n x 的序号n 代表的样值实际位置是( );)(n x 的N 点DFT )k X (中,序号k 代表的样值实际位置又是( )。
6.用8kHz 的抽样率对模拟语音信号抽样,为进行频谱分析,计算了512点的DFT 。
则频域抽样点之间的频率间隔f ∆为_______,数字角频率间隔w ∆为 _______和模拟角频率间隔∆Ω ______。
判断说明题:7.一个信号序列,如果能做序列傅氏变换对它进行分析,也就能做DFT 对它进行分析。
( )计算题8.令)(k X 表示N 点的序列)(n x 的N 点离散傅里叶变换,)(k X 本身也是一个N 点的序列。
如果计算)(k X 的离散傅里叶变换得到一序列)(1n x ,试用)(n x 求)(1n x 。
9.序列}{0,0,1,1)(=n x ,其4点DFT )(k x 如下图所示。
现将)(n x 按下列(1),(2),(3)的方法扩展成8点,求它们8点的DFT ?(尽量利用DFT 的特性)(1)⎩⎨⎧-=)4()()(1n x n x n y 7~43~0==n n(2)⎩⎨⎧=0)()(2n x n y 7~43~0==n n(3)⎪⎩⎪⎨⎧=0)2()(3n x n y 奇数偶数==n n 10.设)(n x 是一个2N 点的序列,具有如下性质:)()(n x N n x =+另设)()()(1n R n x n x N =,它的N 点DFT 为)(1k X ,求)(n x 的2N 点DFT )(k X 和)(1k X 的关系。
离散傅里叶变换(DFT)
第3章 离散傅里叶变换(DFT)
3.2 离散傅里叶变换的基本性质
一. 基本概念
1. 序列的圆周移位 序列x(n),长度为N,则x(n)的圆周移位定义为:
y(n) x((n m))N RN (n) circshift(a,[0,-1])
循环移位过程:
x(n) 周期延拓 x(n) x((n))N 左移m位
x(n), n 0,1, , N 1
N 1
X (z) ZT [x(n)] x(n)zn
n0 N 1
, X (k) DFT[x(n)] x(n)WNkn
n0
比较上面二式可得关系式:
0 k N-1
X (k ) X ( z) , j2 k ze N
0 k N-1
(3.1.3)
序列x(n)的N点DFT是 x(n)的Z变换在单位圆上的N点等间隔采样
(4) 周期为N 的离散周期信号
DFS
N 1
j 2 nk
X (k) x(n)e N
n0
x(n)
1
N 1
j 2 nk
X (k)e N
N k0
k ~ n ~
时域离散周期频域周期离散。频谱特点:周期为N的离散谱
第3章 离散傅里叶变换(DFT)
四种傅立叶变换:
1. 连续非周期 2. 连续周期 3. 离散非周期 4. 离散周期
m0
yc (n) y(n qN )RN (n)
q
序列的N点圆周卷积是序列线性卷积(以N为周期)周
期延拓序列的主值序列。故,当N≥[N1+N2-1]时,线性 卷积与圆周卷积相同。
圆周卷积 是针对DFT引出的一种表示方法
两序列长度必须等,不等时按要求补零
第三章离散傅里叶变换及其快速计算方法(DFT、FFT)
X (e jw )
(2)Z 变换 -- 提供任意序列的 z 域表示。
n
x( n)e jnw
X (z)
n
x ( n) z n
这两种变换有两个共同特征:
(1)变换适合于无限长序列 (2)它们是连续变量 ω 或 z 的函数
华北电力大学自动化系
3
3.1 问题的提出:可计算性
X (z)
而对于
n
x ( n) z n
n
x ( n) z n
找不到衰减因子使它绝对可和(收敛)。为此,定义新函 数,其 Z 变换:
华北电力大学自动化系
15
DFS 定义:正变换
X ( z)
n
x ( n) z n ~ ( n ) z n x
华北电力大学自动化系
6
3.1 问题的提出:傅里叶变换的四种形式 (3)
2. 周期连续时间信号:傅里叶级数 FS
~ (t ) x X (n 0 )
t T
时域周期频域离散
0
2 T
x(t)
~
n -
X(n 0 )e jn0t
时域连续函数造成频域是非周期的谱。 频域的离散对应时域是周期函数。
X (e jT )
T T
X (e jT )e jnT d
取样定理
n
x(nT )e jnT
1 X ( 0 ) T n
时域的离散化造成频域的周期延拓 时域的非周期对应于频域的连续
华北电力大学自动化系
8
《离散傅里叶变换-第三章》
n0 0 = kn 8 7
3
3
2π − j kn 8
3 − j kπ 8
(2) 设变换区间N=16, 则
X(k) = ∑ x(n)W
n= 0
3π k −j 16
π
N= 0 = n0 0
2 = ∑ e, k = 0,1, ⋅ ⋅ ⋅, 7 π N =0 sin( k ) 8
2. 时域循环移位定理 设x(n)是长度为N的有限长序列,y(n)为x(n)的循环移位,即: y(n)=x((n+m))NRN(n) 则: Y(k)=DFT[y(n)]=W-kmNX(k) 其中:X(k)=DFT[x(n)], 0≤k≤N-1
kn 证明: Y ( k ) = DFT [ y (n )] = x (( n + m )) N RN (n )WN ∑ N− 令n+m=n′,则有1 n =0 N −1
~
~ ∞
x (n ) =
m =−∞
∑
x ( n + mN )
(3.1.5)
(3.1.6) ••
~
x (n ) ••
0
••
N-1
•
n
x (n ) = x ( n ) ⋅ RN (n )
~
~
••
••
~(n ) x
•• •
0
••
•
••
•• •
~
••
N-1
•
n
一般定义周期序列 x(n) 中从n=0到N-1的第一个周期为 x(n)的主 n) x(n) (3.1.7) x( 值区间,而主值区间上的序列称为x(n) 的主值序列。(3.1.7) x(n)
第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 1.1) (3. k 0 N 1
所以, 在变换区间上满足下式:
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, 则
X (k ) x (n )W8kn e
n 0 N 0
Y (k ) DFT [ y ( n )]
kn x (( n m)) N RN ( n )WN n 0 kn x (( n m)) N WN n 0 N 1 N 1
令n+m=n’, 则有
Y (k )
N 1 m n m
k x((n)) NWN ( nm ) N 1 m n m
N 1
0 k N-1
比较上面二式可得关系式
X (k ) X ( z )
z e
j
2 k N
, ,
Hale Waihona Puke 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的隐含周期性
第三章 离散傅立叶变换.
第三章 离散傅立叶变换一、离散傅立叶级数计算题:1.如果)(~n x 是一个周期为N 的周期序列,那么它也是周期为2N 的周期序列。
把)(~n x 看作周期为N 的周期序列有)(~)(~1k X n x ↔(周期为N );把)(~n x 看作周期为2N 的周期序列有)(~)(~2k X n x ↔(周期为2N );试用)(k X 1~表示)(k X 2~。
二、离散傅立叶变换定义填空题2.某DFT 的表达式是∑-==10)()(N k kl M Wk x l X ,则变换后数字频域上相邻两个频率样点之间的间隔是( )。
3.某序列DFT 的表达式是∑-==10)()(N k kl M W k x l X ,由此可看出,该序列的时域长度是( ),变换后数字频域上相邻两个频率样点之间隔是( )。
4.如果希望某信号序列的离散谱是实偶的,那么该时域序列应满足条件( )。
5.采样频率为Hz F s 的数字系统中,系统函数表达式中1-z 代表的物理意义是 ),其中时域数字序列)(n x 的序号n 代表的样值实际位置是( );)(n x 的N 点DFT )k X (中,序号k 代表的样值实际位置又是( )。
6.用8kHz 的抽样率对模拟语音信号抽样,为进行频谱分析,计算了512点的DFT 。
则频域抽样点之间的频率间隔f ∆为_______,数字角频率间隔w ∆为 _______和模拟角频率间隔∆Ω ______。
判断说明题:7.一个信号序列,如果能做序列傅氏变换对它进行分析,也就能做DFT 对它进行分析。
( )计算题8.令)(k X 表示N 点的序列)(n x 的N 点离散傅里叶变换,)(k X 本身也是一个N 点的序列。
如果计算)(k X 的离散傅里叶变换得到一序列)(1n x ,试用)(n x 求)(1n x 。
9.序列}{0,0,1,1)(=n x ,其4点DFT )(k x 如下图所示。
现将)(n x 按下列(1),(2),(3)的方法扩展成8点,求它们8点的DFT ?(尽量利用DFT 的特性)(1)⎩⎨⎧-=)4()()(1n x n x n y 7~43~0==n n(2)⎩⎨⎧=0)()(2n x n y 7~43~0==n n(3)⎪⎩⎪⎨⎧=0)2()(3n x n y 奇数偶数==n n 10.设)(n x 是一个2N 点的序列,具有如下性质:)()(n x N n x =+另设)()()(1n R n x n x N =,它的N 点DFT 为)(1k X ,求)(n x 的2N 点DFT )(k X 和)(1k X 的关系。
第3章离散傅里叶变换(DFT)09-10-1
§3.2 离散傅里叶变换的基本性质
一. 线性性质
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
➢再 反 转 形 成 x2((-m))N , 取 主 值 序 列 则 得 到 x2((m))NRN(m),通常称之为x2(m)的圆周反转; ➢对x2(m)的圆周反转序列圆周右移n,形成
x2((n-m))NRN(m); ➢当n=0,1,2,…,N-1时,分别将x1(m)与x2((n-m))NRN(m)相 乘,并在m=0到N-1区间内求和,便得到其循环卷积y(n)。
y(n) x((n m))N RN (n)
则循环移位后的DFT为
Y (k) DFT [ y(n)] DFT [x((n m))N RN (n)] WNmk X (k)
证:利用周期序列的移位性质加以证明
DFS [x((n m)) N ] DFS [~x (n m)] WNmk X~(k)
x1(n)
0
N-1
~x2 (n)
0
N-1
n n
~x2 (m)
x2 0 mN RN (m)
0
m
x2 1 mN RN (m)
0
x2
2
mN
RN
(m)
m
0
m
x2 3 mN RN (m)
0
m
y(n) x1(n) N x2 (n) ➢两个长度
第三章 离散傅里叶变换(DFT)
− N
n
)*
W
n N
=
W
n N
+iN
3. 可约性 4. 正交性
W i⋅n N
= WNn / i
∑ ∑ 1
N
N −1
W
nk N
(WNmk
)
*
k =0
=
1 N
N −1
W (n−m)k N
k =0
=
⎧1, ⎩⎨0,
n − m = iN n − m ≠ iN
3.3 周期序列的离散傅里叶级数
z 可以看出,当0≤k≤N-1 时,X~(k) 是对X(z)在Z平面单 位圆上的N点等间隔采样,在此区间之外随着k的变 化,X~ (k ) 的值呈周期变化。
了。所以这种无穷长序列实际上只有N个序列值的信息是 有用的,因此周期序列与有限长序列有着本质的联系。
3.3 周期序列的离散傅里叶级数
z X~(k) ↔ ~x (n) 是一个周期序列的离散傅里叶 级数(DFS)变换对,这种对称关系可表示为:
∑ X
(k )
=
D F S [ x (n)]
=
N −1
x
10
X (k) =
|X(ejω)|
X (e jω ) ω= 2π k 10
=
− j 4π k
e 10
sin(π k / 2) sin(π k /10)
5
…
o
π
…
2π
3π
4π
ω
3.3 周期序列的离散傅里叶级数
例2 已知周期序列x (n),求X (k )。并讨论 X~ (k)与 X (e jω ) 的关系
将n和k互换,有 ∑ Nx (-k ) = N-1 X (n)WNkn n=0
第三章 离散傅里叶变换(DFT)
~ 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
第三章离散傅里叶变换
不变,F减小N增加,又因增加 因此,和N可按下面两式选择 例1 有一频谱分析用FFT处理器,抽样点数为2的幂,假定没有采用 任何 特殊的数据处理,已给条件为 ①频率分辨率 ②信号的最高频率 求:①最小记录长度 ②抽样点的最大间隔T ③在一个记录中最小点数N 解: ① ② ③ 取 (2)频域泄露(截短产生误差)
●任何有限长序列都可以表示成共轭对称分量和共轭反对称分量 之和,即 ………… ……….(3-2) 对(3-2)式n换成N-n,并取复共轭得 (3-3) 联立(3-2),(3-3)可得:
●任何序列也可以表示实部和虚部 (3-4) 其中 (3-5) (3-6) (3)DFT的共轭对称性 ●对(3-4)进行DFT得: (3-7) ① 对(3-5)进行DFT得: .(3-8) ② 对(3-6)进行DFT得 (3-9) 结论:由(3-7),(3-8),(3-9)可得 其中 ● 任何序列可以表示为共轭对称和共轭反对称分量: (3-10) (3-11) (3-12) ① 对(3-10)进行DFT得 ② 对(3-11)进行DFT得 ③ 对(3-12)进行DFT得 结论: 其中 ●是长度为N的实序列,且,则 ① 共轭对称,即
2 (a) n,m 3 1 0
(b) 1 2 3 n,m
-2 6 5
2 1 -3 N=4 (c) m
m 3 2 n=0 (d)
ቤተ መጻሕፍቲ ባይዱ
m 3 0 n=1 (e)
m 1 0 n=2 (f)
2 m 1 n=3 (g)
2 3 2 m 1 (h) 1
图4
4、复共轭序列的DFT
设是的复共轭序列,长度为N,则 (3-1) 且。 证明:根据DFT的唯一性,只要证明(3-1)式右边等于左边即可。 又由的隐含周期性有 。 同理可证 。
离散傅里叶变换(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
第3章 离散傅里叶变换(DFT)
时域循环移位定理表明:有限长序列的循环移位,在离散 频域中相当于引入一个和频率成正比的线性相移WN-mk 频域循环移位定理表明:时域序列的调制(相移)等效于频域 的循环移位
(3.1.7)
注:若x(n)实际长度为M,延拓周期为N,则当N<M时,(3.1.5) 式仍表示以N为周期的周期序列,但(3.1.6)和 (3.1.7)式仅对 N≥M时成立。
第3章 离散傅里叶变换(DFT)
图3.1.2(a)中x(n)实际长度M=6,
x (n) 如图 当延拓周期N=8时,~
3.1.2(b)所示。
DTFT:X(e )= x( n)e
M 1 n0
N (n) RN (n) xN ( n) x
(k ) x N (n)WNkn DFS : X
DFT与ZT关系:
k
z e
j k N
X (k ) X ( z )
k ,, ,..., N k ,, ,..., N
第3章 离散傅里叶变换(DFT)
(2)时/频域循)] X (k )
k 0,1,..., N 1
则
且
mk DFT [ x(( n m)) N RN (n)] WN X (k )
nl IDFT [ X (( k l )) N RN (k )] WN x ( n)
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
mk kn [ x ( m ) W ] W N N k 0 m 0 k ( mn ) W N k 0 N 1
数字信号处理第三章离散傅里叶变换(DFT)及其快速算法(FFT)
周期
2
s、fs N
分辨率
2
N
fs N
返回
回到本节
DFT和DFS之间的关系:
周期延拓
取主值
有限长序列
周期序列
主值区序列
有限长序列 x(n) n 0,1, 2, M 1
周期序列 xN (n) x(n mN ) x((n))N m 0 n0 N 1 n mN n0 ((n))N n0
四种傅立叶变换
离散傅立叶变换(DFT)实现了信号首次在频域 表示的离散化,使得频域也能够用计算机进行处理。 并且这种DFT变换可以有多种实用的快速算法。使信 号处理在时、频域的处理和转换均可离散化和快速 化。因而具有重要的理论意义和应用价值,是本课程 学习的一大重点。
本节主要介绍
3.1.1 DFT定义 3.1.2 DFT与ZT、FT、DFS的关系 3.1.3 DFT的矩阵表示
• X(k)为x(n)的傅立叶变换 X (e j ) 在区间 [0, 2 ]上的N
点等间隔采样。这就是DFT的物理意义。
j ImZ
2பைடு நூலகம்3
4
5 6
1 2
N
k=0 ReZ
7 (N-1)
DFT与z变换
X(ejω)
X(k)
0
o
2
0
N 1 k
DFT与DTFT变换
回到本节
变量
、f k
之间的某种变换关系.
• 所以“时间”或“频率”取连续还是离 散值,就形成各种不同形式的傅里叶变换 对。
3.1 离散傅里叶变换的定义及物理意义
时间域
t:连续
模拟域
3离散傅里叶变换解析
k
F (k )e
0
jk0t
4.离散、周期时域信号 f p (n) ←映射→周期、离散频域信号 F p (k ) ,它由离散傅里叶级数变 换构成映射关系,即
F p (k )
n 0
N 1
nk f p (n)WN
1 f p ( n) N
F (k )W
p n 0
N 1
Re[ x(n)] Re[ x(n)]
Im[ x(n)] Im[ x(n)]
则称为x(n)共轭反对称序列(conjugate antisymmetric sequence), 通常表示为: x (n) x * (n)
0 o
任何序列x(n)都可以表示为共轭对称序列与共轭反对称序列之和:
x ( n) x e ( n) x o ( n)
其中
xe (n) 1 [ x(n) x * (n)] 2
xo (n)
1 [ x(n) x * (n)] 2
4
若
F x(n) X (e j )
则
F x* (n) X * (e j )
上式说明共轭序列的傅里叶变换等于原序列傅里叶变换的共轭函数的 反函数。
f
n 0
N 1
p
(n)e jn0 r NFr
8
f 以上分析表明,系数 F 可以严格地由
r
N 1 n 0
p (n)e
jn0 r
NFr
式求出,也就是说
f p (n) Fk e jn0k
k 0
N 1
式表述的关系是存在的。
将
f p (n) Fk e
k 0
以上二式说明复指数 e
第3章离散傅里叶变换
第3章 离散傅里叶变换
二.序列的圆周移位 1.定义 一个有限长序列 x( n )的圆周移位定义为
xm (n) xn mN RN n
这里包括三层意思: ~ 先将 x( n )进行周期延拓 x (n) xn N 再进行移位 ~ x (n m) x n m N 最后取主值序列:
第3章 离散傅里叶变换
3.共轭对称特性之一
如果X ( k ) DFT [ x( n )],则 DFT [ x* ( n )] X * (( k ))N RN ( k )
证明:
X * (( N k ))N RN ( k )
nk DFT [ x ( n )] x* ( n )WN RN ( k ) * n 0 nk * nk * [ x( n )WN ] RN ( k ) [ x( n )WNNnWN ] RN ( k ) n 0 n 0 ( N k )n * [ x( n )WN ] RN ( k ) X * (( N k ))N RN ( k ) n 0 N 1 N 1 N 1 N 1
*复数序列实部的DFT 该序列DFT的圆周共轭对称分量。
5.共轭对称 特性之三
第3章 离散傅里叶变换
6.共轭对称 如果 X ( k ) DFT [ x( n )],则 DFT{ j Im [x( n )]} 特性之四 1 [ X (( k ))N X * (( N k ))N ] RN ( k ) X op ( k ) 2 1 * j Im [ x ( n )] [ x ( n ) x ( n )] 证明: 2 1 DFT{ j Im [x( n )]} { DFT [ x( n )] DFT [ x* ( n )]} 2 1 [ X ( k ) X * (( N k ))N RN ( k )] 2 1 [ X (( k ))N X * (( N k ))N ] RN ( k ) X op ( k ) 2
第3章 离散傅里叶变换(DFT)
M为整数 M为整数
x (n ) =
m = −∞
∑
∞
x ( n + mN )
(3.1.5) (3.1.6)
x (n ) = x (n ) ⋅ RN (n )
~
~
x(n)=x((n))N,
% X (k ) =
m =− ∞
∑ X (k + mN )
∞
% X (k ) = X (k ) RN (k )
回到本节
N k=0
k =0 N
为DFT变换 长度N≥M, , N 为DFT变换 长度N≥M, WN = e DFT 有限长 离散序列 有限长 离散序列
−j
2π N
第三章 离散傅里叶变换DFT
例1
解:
已知 x(n) = R4 (n),分别求N = 8和N =16 时的X (k)。
N = 8时
N−1 n=0 nk N
第三章 离散傅里叶变换DFT
式中x((n))N表示x(n)以N为周期的周期延拓序列, ((n))N 表示n对N求余, 即如果 n=MN+n1, 0≤n1≤N-1, 则 ((n))N=n1 例如 N = 5, x N (n) = x((n))5 则有
~
M为整数,
x (5) = x ((5))5 = x (0) x (6) = x ((6))5 = x (1)
∑e
n=0
k =0 8, = 0, k = 1, 2, 3, 4, 5, 6, 7
x(n)的16点DFT为
k 1 − W168 1 − e k X (k ) = W16 n = = k 2π −j k 1 − W16 n=0 1 − e 16 π 7π sin k −j k 2 = e 16 , k = 0,1, 2,L ,15 π sin k 16
第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