七,用 DFT对模拟信号作频谱分析信号的频谱分析:计算信号的傅里叶变换
0
0
s
h
T
f
T
F
N
f
时 域 采 样 间 隔时 域 采 样 频 率信 号 记 录 长 度
( 频 率 分 辨 率 ) 频 域 采 样 间 隔采 样 点 数信 号 最 高 频 率
0
0
sTfN
TF
1/sfT?2shff? 001/TF?
0sf NF? 0T NT?
对连续时间非周期信号的 DFT逼近
( ) ( ) jtX j x t e dt
1() 2 jtx t X j e d
( ) ( ) ( )j t j nT
n
X j x t e dt x nT e T
n
t nT dt T dt T
1)将 在 轴上等间隔( T)分段()xt t
2)将 截短成有限长序列()xn 00~t T N?,个时域抽样点
N- 1
0
( ) ( ) j nT
n
X j T x nT e
002 F
0k
21
0
()
N j nk
N
n
T x n e
3)频域抽样:一个周期分 N段,采样间隔,时域周期延拓,
周期为
0F
001/TF?
0
N - 1
0
0
( ) ( ) jk nT
n
X jk T x nT e
002 / 2 /sT F f N
[ ( )]T D F T x n
0
1( ) ( )
2
s j nTx nT X j e d
0
1
00
0
1 ()
2
N
jk nT
k
X jk e?
0
1
0
0
N
k
d
d
21
00
0
()
N j nk
N
k
F X jk e
21
0
0
1 ()N j nkN
s
k
f X jk eN
1N N
01 / [ ( ) ]T I D F T X jk
对连续时间非周期信号的 DFT逼近过程
1)时域抽样
2)时域截断
3)频域抽样
0( ) [ ( ) ]X j k T D F T x n
0
1( ) [ ( ) ]x n I D F T X jk
T
近似逼近:
对连续时间周期信号的 DFS逼近
0 0
0 0
0
1( ) ( )T jk tX jk x t e d t
T
00() jk t
k
x t X jk e
0
1
0
00
1( ) ( )N jk nT
n
X jk x nT e TT
0
1
0 0
NT
n
t nT dt T dt T
1)将 在 轴上等间隔( T)分段()xt t
1 [ ( ) ]DFS x n
N
0 2/TN
0T NT?
21
0
1 ()N j nkN
n
x n eN
2)频域截断:长度正好等于一个周期
0
1
0
0
( ) ( )
N
jk nT
k
x nT X jk e
21
0
0
()
N j nk
N
k
X jk e
21
0
0
1 ()N j nkN
k
N X jk eN
0[ ( ) ]N I D F S X jk
0
1( ) [ ( ) ]X jk DFS x n
N
0( ) [ ( ) ]x n N I D F S X j k
近似逼近:
频率响应的混叠失真及参数的选择
0
0
sTfN
TF
2shff?时域抽样:
001/FT?频域抽样:
同时提高信号最高频率和频率分辨率,需增加采样点数 N。
0
0
sTfN
TF
hsff要 增 加 信 号 最 高 频 率 则
0NF当 给定 必,即分辨率
00
0
1FT
F要提高频率分辨率,即 则
shN T f f当 给定 则 要不产生混叠,必信号最高频率与频率分辨率之间的矛盾
hf信号最高频率 的确定
0 /2htT?
0
11
2h hf Tt
FFT例:有一频谱分析用的 处理器,其抽样点数必须是2的 整数幂,假设没有采用任何的数据处理措施,已给条件为:
1 1 0
24
Hz
k Hz
)频率分辨率
)信号最高频率 0 1
2 T
N
T
试确定以下参量:
)最小记录长度
)抽样点间的最大时间间隔 (即最小抽样频率)
3 )在一个记录中最少点数
1解,)最小记录长度:
0
0
11 0,1
10TsF
2 2 1 /s h sf f f T)最大抽样间隔 ( )
3
11 0 1 2 5
2 2 4 1 0hT m sf,
3)最小记录点数
3
0
2 2 4 1 0 800
10
hfN
F
102 2 1 0 2 4 8 0 0mN取频谱泄漏改善方法:
对时域截短,使频谱变宽拖尾,称为泄漏
1)增加 x(n)长度 2)缓慢截短栅栏效应改善方法:
增加频域抽样点数 N(时域补零),使谱线更密
DFT只计算离散点(基频 F0的整数倍处)的频谱,
而不是连续函数频率分辨率提高频率分辨率方法:
增加信号实际记录长度补零并不能提高频率分辨率
001/FT?
0
0
s
h
T
f
T
F
N
f
时 域 采 样 间 隔时 域 采 样 频 率信 号 记 录 长 度
( 频 率 分 辨 率 ) 频 域 采 样 间 隔采 样 点 数信 号 最 高 频 率
0
0
sTfN
TF
1/sfT?2shff? 001/TF?
0sf NF? 0T NT?
对连续时间非周期信号的 DFT逼近
( ) ( ) jtX j x t e dt
1() 2 jtx t X j e d
( ) ( ) ( )j t j nT
n
X j x t e dt x nT e T
n
t nT dt T dt T
1)将 在 轴上等间隔( T)分段()xt t
2)将 截短成有限长序列()xn 00~t T N?,个时域抽样点
N- 1
0
( ) ( ) j nT
n
X j T x nT e
002 F
0k
21
0
()
N j nk
N
n
T x n e
3)频域抽样:一个周期分 N段,采样间隔,时域周期延拓,
周期为
0F
001/TF?
0
N - 1
0
0
( ) ( ) jk nT
n
X jk T x nT e
002 / 2 /sT F f N
[ ( )]T D F T x n
0
1( ) ( )
2
s j nTx nT X j e d
0
1
00
0
1 ()
2
N
jk nT
k
X jk e?
0
1
0
0
N
k
d
d
21
00
0
()
N j nk
N
k
F X jk e
21
0
0
1 ()N j nkN
s
k
f X jk eN
1N N
01 / [ ( ) ]T I D F T X jk
对连续时间非周期信号的 DFT逼近过程
1)时域抽样
2)时域截断
3)频域抽样
0( ) [ ( ) ]X j k T D F T x n
0
1( ) [ ( ) ]x n I D F T X jk
T
近似逼近:
对连续时间周期信号的 DFS逼近
0 0
0 0
0
1( ) ( )T jk tX jk x t e d t
T
00() jk t
k
x t X jk e
0
1
0
00
1( ) ( )N jk nT
n
X jk x nT e TT
0
1
0 0
NT
n
t nT dt T dt T
1)将 在 轴上等间隔( T)分段()xt t
1 [ ( ) ]DFS x n
N
0 2/TN
0T NT?
21
0
1 ()N j nkN
n
x n eN
2)频域截断:长度正好等于一个周期
0
1
0
0
( ) ( )
N
jk nT
k
x nT X jk e
21
0
0
()
N j nk
N
k
X jk e
21
0
0
1 ()N j nkN
k
N X jk eN
0[ ( ) ]N I D F S X jk
0
1( ) [ ( ) ]X jk DFS x n
N
0( ) [ ( ) ]x n N I D F S X j k
近似逼近:
频率响应的混叠失真及参数的选择
0
0
sTfN
TF
2shff?时域抽样:
001/FT?频域抽样:
同时提高信号最高频率和频率分辨率,需增加采样点数 N。
0
0
sTfN
TF
hsff要 增 加 信 号 最 高 频 率 则
0NF当 给定 必,即分辨率
00
0
1FT
F要提高频率分辨率,即 则
shN T f f当 给定 则 要不产生混叠,必信号最高频率与频率分辨率之间的矛盾
hf信号最高频率 的确定
0 /2htT?
0
11
2h hf Tt
FFT例:有一频谱分析用的 处理器,其抽样点数必须是2的 整数幂,假设没有采用任何的数据处理措施,已给条件为:
1 1 0
24
Hz
k Hz
)频率分辨率
)信号最高频率 0 1
2 T
N
T
试确定以下参量:
)最小记录长度
)抽样点间的最大时间间隔 (即最小抽样频率)
3 )在一个记录中最少点数
1解,)最小记录长度:
0
0
11 0,1
10TsF
2 2 1 /s h sf f f T)最大抽样间隔 ( )
3
11 0 1 2 5
2 2 4 1 0hT m sf,
3)最小记录点数
3
0
2 2 4 1 0 800
10
hfN
F
102 2 1 0 2 4 8 0 0mN取频谱泄漏改善方法:
对时域截短,使频谱变宽拖尾,称为泄漏
1)增加 x(n)长度 2)缓慢截短栅栏效应改善方法:
增加频域抽样点数 N(时域补零),使谱线更密
DFT只计算离散点(基频 F0的整数倍处)的频谱,
而不是连续函数频率分辨率提高频率分辨率方法:
增加信号实际记录长度补零并不能提高频率分辨率
001/FT?