第
36
卷
第
11
期 系统工程与电子技术
Vol.36
No.1 1
2014
年
11
月
S
y
stems En
g
ineerin
g
and Electronics November 2014
文章编号:
1001-506X
(
2014
)
11-2149-07
网址:
www.s
y
s-ele.com
收稿日期:
2013 11 27
;修回日期:
2014 04 24
;网络优先出版日期:
2014 06 05
。
网络优先出版地址:
htt
p
:
∥
www.cnki.net
/
kcms
/
detail
/
11 .2422 .TN.20140605 .1743 .003 .html
基金项目:国家 自然 科学 基金(
61 271322
)资助课题
衍生于全相位
FFT
的双子段相位估计法
黄翔东,余
佳,孟天伟,王兆华
(天津大学电子信息工程学院,天津
300072
)
摘
要:为提高给定样本数量下 的 相位 估 计精 度,提出 双 子 段测 相 法。该 方法 是 从 全相 位 快 速傅 里 叶 变 换
(
fast Fourier transform
,
FFT
)测相法衍生的,只需取出全相位
FFT
包含的
N
路子数据段中的首、尾两个子段,再
综合这两个子段的
FFT
谱提取峰值谱相位信息即可完成高精度测 相。 理论 分析 了 在噪 声环 境中,所提方 法测 相
精度高于全相位
FFT
测相法的原因。仿真 实验 表 明,在 信 噪比 不 太 低的 场 合,双子 段 测 相 法 的 相位 估 计 方差 处
于全相位
FFT
估计法和两参数克拉美罗限之间,因而具有广泛应用前景。
关键词:双子段测相;全相位傅里叶变换;峰值谱;克拉美罗限
中图分类号:
TM 933 .3
;
TN 911 .72
文献标志码:
A
DOI
:
10.3969
/
j
.issn.1001-506X.2014.11 .07
Double-subse
g
ment measurin
g p
hase method derived from all-
p
hase FFT
HUANG Xian
g
-don
g
,
YU Jia
,
MENG Tian-wei
,
WANG Zhao-hua
(
School o
f
Electronic In
f
ormation En
g
ineerin
g
,
Tian
j
in Universit
y
,
Tian
j
in 300072
,
China
)
Abstract
:
To enhance the
p
hase measurement accurac
y
with a fixed amount of sam
p
les
,
a double-subse
g
-
ment measurin
g p
hase method is
p
ro
p
osed
,
which is derived from all-
p
hase fast Fourier transform
(
FFT
)
.In the
p
ro
p
osed method
,
onl
y
the first subse
g
ment and the last subse
g
ment amon
g
the
N
subse
g
ments included in all-
p
hase FFT are taken out
,
and accurate
p
hase measurement can be finished throu
g
h s
y
nthesizin
g
the
p
hase infor-
mation of the
p
eak s
p
ectra from these two se
g
ments FFT s
p
ectrums.Theoretic anal
y
sis is also conducted to an-
al
y
ze the reason wh
y
the
p
ro
p
osed method has hi
g
her
p
hase measurement accurac
y
than that of all-
p
hase FFT in
nois
y
circumstance.Moreover
,
simulation also shows that
,
when the si
g
nal-to-noise rate is not too low
,
the
p
hase estimation variance of the
p
ro
p
osed method is between that of all-
p
hase FFT and Cramer-Rao lower bound
CRLB of two
p
arameters
,
which endows the
p
ro
p
osed method with wide a
pp
lication
p
ros
p
ects.
Ke
y
words
:
double-subse
g
ment measurin
g p
hase
;
all-
p
hase fast Fourier transform
(
FFT
);
p
eak s
p
ectra
;
Cramer-Rao lower bound
(
CRLB
)
0
引
言
相位估计是信号处理的经典研究热 题
[
1 3
]
,广 泛应 用于
振动分析
[
4
]
、频率估计
[
5 6
]
、通信
[
7
]
、光 学 工 程
[
8
]
等 领域。 快
速傅里叶变换(
fast Fourier transform
,
FFT
)含有 丰富相位信
息,常作为测相手段。而克拉美罗限(
Cramer-Rao lower bound
,
CRLB
)则是衡量测相方法抗噪性能好坏的尺度
[
9 10
]
。
相位估计方法的测 相精 度与 信号 参数 估计 模型 紧 密相
关
[
9 1 1
]
。对振动信号最基本 成分———经典复指数信号
a
·
ex
p
[
j
(
ω
0
n
+
θ
0
)]测相时,由于涉及估计
3
个参数[
a
,
ω
0
,
θ
0
]
T
,因 而
相位估计方差由文献[
9
]所推导的
3
参数
CRLB
(
CRLB
3
)所
限定。现有大 多 数 测 相 方 法,如 希 尔 伯 特 变 换 法
[
12
]
、比 值
法
[
1 3
]
、相位差法
[
14 1 5
]
、余弦 窗 校 正 法
[
1 6
]
都 是 在 先 估 计 频 率
的前提下,再依据频率估计结果去估计相位,这就引 起误 差
扩散问题。正如文献 [
16
]所 指 出,频 率 估 计 的 偏 差 会 带 入
到相位估计中,从 而 增 加 了 相 位 估 计 的 不 确 定 性。 为 弥 补
相位估计依赖于频 率 估 计 的 缺 陷,文 献 [
17 18
]提 出 全 相
位
FFT
(
all-
p
hase FFT
,
a
p
FFT
)测 相 法,该 方 法 把 “
ω
0
n
+
θ
0
”作为整体研 究 对 象,因 而 相 位 估 计 不 依 赖 于 频 率 估 计,
文献[
19
]还推导出
a
p
FFT
测相 方 差理 论表 达 式,文 献[
11
]
证明了
a
p
FFT
测相 符 合 两 参 数 (而 不 是 传 统 的
3
参 数)估
计模型,其相 位 估 计 方 差 突 破 了 传 统 的
CRLB
3
,因 而 测 相
精度高于现有的符合
3
参数模型的测相法
[
12 1 6
]
。
然而作为新的符合 两 参 数 模 型 的 测 相 法,
a
p
FFT
测 相
精度仍有待提高。文 献[
11
]指 出,
a
p
FFT
测 相 方 差 离 两 参
数
CRLB
(
CRLB
2
)还 有 一 定 距 离。 为 此,在 深 入 剖 析