Costas环在无线扩频定位系统中的应用技术研究(1)
- 格式:pdf
- 大小:295.92 KB
- 文档页数:5
第14卷第2期 中国惯性技术学报 2006年4月 Vol.14
No.2 Journal of Chinese Inertial Technology Apr. 2006
文章编号:1005-6734(2006)02-0056-04
Costas环在无线扩频定位系统中的应用技术研究
宋建材,赵小明,杨兴文,赵素文,吴俊杰,张 鹏
(天津航海仪器研究所,天津300130)
摘要:分析了扩频定位系统中载波恢复和数据解调数学模型,设计了环路中应用的Butterworth
低通滤波器和产生正交两路信号的Hilbert变换滤波器,在Simulink中建立了Costas环仿真
模型。结果表明:该Costas环能够准确地实现扩频定位系统的载波恢复和发送端数据解调输
出,理论仿真和实际工程实验基本一致。
关 键 词:科斯特斯环;载波恢复;扩频;无线电定位
中图分类号:U666.1 文献标识码:A
Application of Costas loop in wireless spread spectrum positioning system
SONG Jian-cai, ZHAO Xiao-ming, YANG Xing-wen, ZHAO Su-Wen, WU Jun-Jie, ZHANG Peng
(Tianjin Navigation Instrument Research Institute, TianJin 300130, China)
Abstract: The mathematical models of Costas Loop for carrier recovery and data demodulation are
analyzed. Then the Butterworth low-pass filter and Hilbert transform filter for generating two channel
quadratic signal are designed, and a Costas simulation model in Simulink is developed. The
simulation results show that the carrier recovery and data demodulation can be accurately realized,
and are consisted with practical engineering experiments.
Key words: Costas loop; carrier recovery; spread spectrum; wireless position
0 引 言
无线扩频定位系统采用直接序列扩频技术,抑制载波平衡调制为BPSK方式,设载波为
0cosAtω
,
伪码信号为()
ct
,平衡调制信号为:
()()()()()
()
0000cos1cosct
ftAtcttAttωπϕωϕ
=⎡+×+⎤=−⎡+⎤
⎣⎦⎣⎦ (1)
式中,
0ω
为载频,A为载波振幅,
0()tϕ
为载波初相,()
ct
为伪随机码序列。为简便起见,设1A=
,
0()0tϕ
=
,则有:
()
()0
0cos=0
()
cos=1tct
ft
tctω
ω⎧
⎪
=
⎨
−
⎪
⎩ (2)
我们采用平衡Gold码序列对中频载波信号进行扩频,相当于进行BPSK调制,扩频信号的频谱如图
基金项目:国防科技应用基础研究项目
收稿日期:2006-02-21
作者简介:宋建材(1979—)男,天津航海仪器研究所硕士研究生,从事导航、制导与控制(扩频导航定位技术研究)。
电子邮箱:songjiancai_tj@163.com 第2期 宋建材等:Costas环在无线扩频定位系统中的应用技术研究 57
1所示。由图1可知,它不含有载波频率分量,用普通的锁相环路无法提取,需设计特殊的锁相环路—抑
制载波环对接收信号进行非线性处理后,可以重新
产生载波分量,实现载波的跟踪与同步。
Riter[1]
已经证明,跟踪低信噪比的抑制载波信号
的最佳装置是科斯特斯环(Costas Loop)及平方环
(Square Loop),J.K.Homes[2]
曾证明这两种环路是等
效的,但其实现方法完全不同。由于平方环可能有“理
想平方器”的问题,而且在本系统的中频频率上,Costas环在电路实现上比较方便,故而选用Costas环。
1 Costas环数学模型
Costas环又称为“IQ−”环,是用来解调抑制载波信号,也是二相或四相移相键控信号解调的专用
环路,其工作频率是载波频率。Costas环是将输入信号平方,形成离散的谱线,然后通过这些离散的谱
线跟踪输入信号。图2是“线性”形式的科斯特斯环。
压控振荡器(VCO)也用来产生载波参考信号,它与输入信号同相相乘及相移90°再相乘,乘法器再
经低通滤波器输出,其中的关键在于增加了第三个乘法器,而两路低通滤波器的输出都加到第三个乘法
器上,再用第三个乘法器的输出经环路滤波后去控制环路的压控振荡器(VCO)。
在考虑抑制载波信号加噪声的跟踪问题时,我们假设Costas环接收到的信号为:
()()()
ytstnt=+
(3)
其中,有用信号:
()
st
为
0()()cos()
sstActtωϕ
=+ (4)
加性窄带噪声()
nt
为:
00()()cos()()sin()
csssntnttnttωϕωϕ
=+++
(5)
在()
yt
的作用下,乘法器I
的输出为
IZ
经低通滤波后为
IU;乘法器Q的输出为
QZ,经低通滤波
后输出的
QU。
IU
与
QU在第三个乘法器相乘输出为()
det
。乘法器I的作用相当于鉴相器,其输出为:
[]
I0()()()cos()
rZtstnttωϕ
=++
(6)
经低通滤波器后成为: []
I11
()()()cos()()sin()
22csrssrUtActntntϕϕϕϕ
=+−+−
(7)
同理 []
Q0()()()sin()
rZtstnttωϕ
=++
[]
Q11
()()()sin()()cos()
22csrssrUtActntntϕϕϕϕ
=−+−+−
(9)
第三个乘法器输出为:
[][]2
2111
()()()sin2()()()()cos2()()sin2()
848dcsrscsrssretActntntActntntϕϕϕϕϕϕ
=−+−++−+−
(10) 式(10)右边第一项为环路误差控制电压,其它项为相位抖动,这正是Costas环的随机微分方程(SDE)。
2 Costas环的锁定检测器
Costas环的锁定检测器不止一种,图3是比较通用的相干锁定检测器,又称为对称锁定检测器。 图1 BPSK调制扩频信号的频谱图 600
400
200
0800
幅度
0 0.5 1 1.5 2 2.5 3 3.5
频率
×104
(8)低通滤波
低通滤波VCO
90o()
ytI解调输出数据
相干载波()
Fs
Q
图2 “线性”形式Costas环基本结构图58 中国惯性技术学报 2006年4月
假定相差:
0
srφϕϕ
=−=
(11)
同时Costas环路接收到的信号
如式(3)所示。将式(11)代入公
式(7)和(9)可得经过滤波器后
的信号:
()()()
I1
2cUtActnt=⎡+⎤
⎣⎦ (12)
()()
Q1
2sUtnt=
(13)
因此锁定检测器产生的结果为:
()()()
222222
IQ1
2
4scsztUtUtAcAcnnn⎡⎤
=−=++−
⎣⎦ (14)
3 计算机仿真及结果分析
由以上分析可知,在Costas环中每个乘法器输出都包含有高频分量,必须加入低通滤波器将其滤除。
本系统采用二阶Butterworth模拟低通滤波器的设计方法,通过采用双线性Z变换得到IIR数字低通滤波
器。设计指标是:通带衰减1 dB,阻带衰减40 dB,截止频率是采样频率的0.1倍,图4是二阶Butterworth
数字低通滤波器的幅频特性曲线和相频特性曲线。由图4可以看出幅频特性曲线满足技术要求,有很好
的带通和带阻特性。
Costas环中同相—正交两路信号
采用Hilbert变换滤波器来实现,理论
上Hilbert变换在不改变输入信号幅度
的前提下可以将其相位移动π/2−,
Hilbert变换的传递函数为:
()() 0
sgn
0Hj
Hj
jω
ωω
ω−>
⎧
=−=
⎨
<
⎩
(15)
具有传递函数()
sgnjω
−
的滤波器称为希尔伯特滤波器。希尔波特变换有多种实现方法,这里采用
一个30阶的滤波器来实现,滤波器的抽头系数由MATLAB函数remez产生[5]
:
% 得到将输入信号中频率为()
0.1-0.9
sf×的分量移项
% 90度所用的Hilbert滤波器的系数
Hilb = remez(30, [.1 .9], [1 1], 'Hilbert');
% 得到输入信号延迟15个单
%元所用离散数字滤波器系数
Delay = [zeros(1, N/2), 1];
Hilbert滤波器的结构及其仿真
结果如图5所示。
由图5可知,输入信号经过离
低通滤波低通滤波VCO
90o()
ytI
Q恢复载波()
Fs
低通滤波门限检测
−+
图3 带相干锁定检测器的Costas环框图
图4 Butterworth低通滤波器幅频特性曲线和相频特性曲线 幅频响应(dB)
相位响应
幅度/dB
相位/degree
归一化频率/归一化频率/
Delay(Z)
1
Hilb(Z)
1 Delay N/2
-90 dgg phase shift121
Delayed input
Hilbert TF of inputinput
时间/s 幅度
图5 Hilbert变换滤波器结构及仿真输出