CN103595405A - 一种基于粒子滤波的锁相环实现方法 - Google Patents
一种基于粒子滤波的锁相环实现方法 Download PDFInfo
- Publication number
- CN103595405A CN103595405A CN201310597377.5A CN201310597377A CN103595405A CN 103595405 A CN103595405 A CN 103595405A CN 201310597377 A CN201310597377 A CN 201310597377A CN 103595405 A CN103595405 A CN 103595405A
- Authority
- CN
- China
- Prior art keywords
- particle
- phase
- locked loop
- sigma
- centerdot
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
本发明涉及一种基于粒子滤波的锁相环实现方法,属于通信信号处理技术领域。本发明方法,具有比传统锁相环应用范围更加广泛的优势,可应用于非线性,非高斯的场合。
Description
技术领域
本发明涉及一种基于粒子滤波的锁相环实现方法,属于通信信号处理技术领域。
背景技术
锁相环技术应用于对宇宙中飞行目标的跟踪和遥测。锁相环以锁定输入载波信号的相位和频率为目标的一种载波环实现形式,然而传统锁相环中,环路滤波器采用的是低通滤波器,属于线性器件,只能解决线性问题,在本发明中采用粒子滤波器替代低通滤波器,可以使得锁相环的应用范围更加广泛
粒子滤波可以解决非线性,非高斯***的问题,因此应用范围广,日渐成为学术界的研究热点。近年来,粒子滤波逐渐引起信号处理和通信领域的重视。而本发明将粒子滤波方法应用于锁相环中来提高锁相环的性能,扩大锁相环的应用范围。
发明内容
本发明的目的是为改善现有锁相环不能应用于非线性非高斯的情况下的缺陷,提出一种基于粒子滤波的锁相环实现方法,扩大锁相环能够使用的范围。
本发明的目的是通过以下技术方案实现的。
本发明的一种基于粒子滤波的锁相环实现方法,其实现步骤如下:
步骤1、将输入的数字信号ui输入进入鉴相器,鉴相器另一个端口输入的是上一时刻锁相环输出的信号,记为uo。经鉴相器输出的结果为ud。若为初始状态,则uo可以为锁相环前端的器件粗略估计的信号。
步骤2、产生粒子。
按照设定的概率分布π(π通常选取一个均值为0,方差很大的Gaussian分布,可以依据具体情况而定,这种情况下,一般视为均匀分布),采样得到k时刻粒子集i=1,2,...,N,其中上标i表示样本序号,i=1,2,…,N,下标k为符号周期序号,k=1,2,…。
步骤3、确定状态转移方程与观测方程。
粒子滤波的状态方程如式(1)所示:
xk+1=Φxk+Ψek+Γuk (1)
其中, 一般xk代表此时滤波器的状态; 其中a1=T(T为环路滤波器更新时间间隔); 其中a2=ω3 np·T·k(ωnp为环路等效噪声的1.2倍,k为根据实际情况设定的系数),a3=2ω2 np·T·k; 和表示的是在环路滤波过程中产生的噪声。
粒子滤波的观测方程如式(2)所示:
其中,a4=2ωnp·k;vk+1表示在观测过程中所产生的噪声。
步骤4、粒子滤波初始化
步骤4.1根据式(1)k=0时,可以设x0=0,将步骤2里产生的粒子集合里的所有粒子带入到方程(1)中,得到N个粒子;
步骤4.3,将步骤4.1输出的每个粒子重要性权值进行归一化
步骤6、将步骤5的输出值输入至压控振荡器,将输出结果作为锁相环的结果,同时反馈给鉴相器。
步骤7、进行步骤1。
步骤8、粒子滤波:
步骤8.1、重复步骤2,得到粒子集合,根据式(1)可以将粒子集合里所有粒子带入到方程(1)中,得到N个粒子;
步骤8.2、计算第k时刻里,每一个粒子的重要性权值;
步骤8.3、权值归一化,当各个粒子权值计算完成后进行权值归一化;
有益效果
本发明方法,具有比传统锁相环应用范围更加广泛的优势,可应用于非线性,非高斯的场合。
附图说明
图1是本发明的方法粒子滤波部分的冲击响应图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明。
实施例
一种基于粒子滤波的锁相环实现方法,其实现步骤如下:
步骤1、将输入的数字信号ui=sin(2π×20t)输入进入鉴相器,鉴相器另一个端口输入的是上一时刻锁相环输出的信号,记为uo,作为初始时刻,没有上一时刻的输出信号uo,这是可根据前端的收集到的信息作为uo,我们在这里设uo=cos(2π×19.5t)。经鉴相器输出的结果为
步骤2、产生粒子。
按照设定的概率分布π(π通常选取一个均值为0,方差很大的Gaussian分布,可以依据具体情况而定,我们在这里设定方差为1000),采样得到k时刻粒子集i=1,2,...,N,其中上标i表示样本序号,i=1,2,…,N,这里N我们取值为100,即取1024个粒子;下标k为符号周期序号,k=1,2,…。
步骤3、确定状态转移方程与观测方程。
粒子滤波的状态方程如式(1)所示:
xk+1=Φxk+Ψek+Γuk (1)
其中, 一般xk代表此时滤波器的状态; 其中a1=T(T为环路滤波器更新时间间隔在这里,我们设这个值为0.01ms); 其中a2=ω3 np·T·k(ωnp为环路等效噪声的1.2倍,在这里我们设环路等效噪声带宽为50Hz,则ωnp=60Hz,k为根据实际情况设定的系数,我们将k值设为1),a3=2ω2 np·T·k; 和表示的是在环路滤波过程中产生的噪声,我们设定和都服从高斯分布N(0,0.012)。
粒子滤波的观测方程如式(2)所示:
其中,a4=2ωnp·k;vk+1表示在观测过程中所产生的噪声,我们设定vk+1服从高斯分布N(0,0.012)。
步骤4、粒子滤波初始化
步骤4.1根据式(1)k=0时,可以设x0=0,将步骤2里产生的粒子集合里的所有粒子带入到方程(1)中,得到1024个粒子;
步骤4.3,将步骤4.1输出的每个粒子重要性权值进行归一化
步骤6、将步骤5的输出值输入至压控振荡器,将输出结果作为锁相环的结果,同时反馈给鉴相器输出结果。
步骤7、进行步骤1。
步骤8、粒子滤波。
步骤8.1、重复步骤2,得到粒子集合,根据式(1)可以将粒子集合里所有粒子带入到方程(1)中,得到1024个粒子;
步骤8.2、计算第k时刻里,每一个粒子的重要性权值;
步骤8.3、权值归一化,当各个粒子权值计算完成后进行权值归一化;
Claims (1)
1.一种基于粒子滤波的锁相环实现方法,其特征在于步骤如下:
步骤1、将输入的数字信号ui输入进入鉴相器,鉴相器另一个端口输入的是上一时刻锁相环输出的信号,记为uo;经鉴相器输出的结果为ud;若为初始状态,则uo可以为锁相环前端的器件粗略估计的信号;
步骤2、产生粒子;
按照设定的概率分布π(π通常选取一个均值为0,方差很大的Gaussian分布,可以依据具体情况而定,这种情况下,一般视为均匀分布),采样得到k时刻粒子集i=1,2,...,N,其中上标i表示样本序号,i=1,2,…,N,下标k为符号周期序号,k=1,2,…;
步骤3、确定状态转移方程与观测方程;
粒子滤波的状态方程如式(1)所示:
xk+1=Φxk+Ψek+Γuk (1)
其中, 一般xk代表此时滤波器的状态; 其中a1=T(T为环路滤波器更新时间间隔); 其中a2=ω3 np·T·k(ωnp为环路等效噪声的1.2倍,k为根据实际情况设定的系数),a3=2ω2 np·T·k; 和表示的是在环路滤波过程中产生的噪声;
粒子滤波的观测方程如式(2)所示:
其中,a4=2ωnp·k;vk+1表示在观测过程中所产生的噪声;
步骤4、粒子滤波初始化
步骤4.1根据式(1)k=0时,可以设x0=0,将步骤2里产生的粒子集合里的所有粒子带入到方程(1)中,得到N个粒子;
步骤4.3,将步骤4.1输出的每个粒子重要性权值进行归一化
步骤6、将步骤5的输出值输入至压控振荡器,将输出结果作为锁相环的结果,同时反馈给鉴相器;
步骤7、进行步骤1;
步骤8、粒子滤波:
步骤8.1、重复步骤2,得到粒子集合,根据式(1)可以将粒子集合里所有粒子带入到方程(1)中,得到N个粒子;
步骤8.2、计算第k时刻里,每一个粒子的重要性权值;
步骤8.3、权值归一化,当各个粒子权值计算完成后进行权值归一化;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310597377.5A CN103595405B (zh) | 2013-11-22 | 2013-11-22 | 一种基于粒子滤波的锁相环实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310597377.5A CN103595405B (zh) | 2013-11-22 | 2013-11-22 | 一种基于粒子滤波的锁相环实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103595405A true CN103595405A (zh) | 2014-02-19 |
CN103595405B CN103595405B (zh) | 2016-06-29 |
Family
ID=50085412
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310597377.5A Active CN103595405B (zh) | 2013-11-22 | 2013-11-22 | 一种基于粒子滤波的锁相环实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103595405B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108594277A (zh) * | 2018-04-27 | 2018-09-28 | 北京邮电大学 | 一种相位差确定方法、装置、电子设备及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050047496A1 (en) * | 2003-08-13 | 2005-03-03 | Mcintire William K. | Modem with pilot symbol synchronization |
CN1866763A (zh) * | 2006-05-11 | 2006-11-22 | 上海交通大学 | 基于粒子滤波的信道估计方法 |
CN101055563A (zh) * | 2007-05-21 | 2007-10-17 | 北京理工大学 | 基于多建议分布的粒子滤波方法 |
CN101662433A (zh) * | 2009-06-23 | 2010-03-03 | 中山大学 | 一种基于粒子滤波修正的信道预测方法 |
CN103036669A (zh) * | 2012-12-24 | 2013-04-10 | 北京理工大学 | 一种基于粒子滤波的符号同步方法 |
-
2013
- 2013-11-22 CN CN201310597377.5A patent/CN103595405B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050047496A1 (en) * | 2003-08-13 | 2005-03-03 | Mcintire William K. | Modem with pilot symbol synchronization |
CN1866763A (zh) * | 2006-05-11 | 2006-11-22 | 上海交通大学 | 基于粒子滤波的信道估计方法 |
CN101055563A (zh) * | 2007-05-21 | 2007-10-17 | 北京理工大学 | 基于多建议分布的粒子滤波方法 |
CN101662433A (zh) * | 2009-06-23 | 2010-03-03 | 中山大学 | 一种基于粒子滤波修正的信道预测方法 |
CN103036669A (zh) * | 2012-12-24 | 2013-04-10 | 北京理工大学 | 一种基于粒子滤波的符号同步方法 |
Non-Patent Citations (2)
Title |
---|
叶斌丽: "雷达目标跟踪的粒子滤波算法研究及实现", 《中国优秀硕士学位论文全文数据库信息科技辑》, no. 4, 15 April 2011 (2011-04-15), pages 3 - 4 * |
王法胜,赵清杰: "一种用于解决非线l性滤波问题的新型粒子滤波算法", 《计算机学报》, vol. 31, no. 2, 29 February 2008 (2008-02-29) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108594277A (zh) * | 2018-04-27 | 2018-09-28 | 北京邮电大学 | 一种相位差确定方法、装置、电子设备及存储介质 |
CN108594277B (zh) * | 2018-04-27 | 2020-06-12 | 北京邮电大学 | 一种相位差确定方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103595405B (zh) | 2016-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103399202B (zh) | 可变速抽水蓄能机组控制***的锁相方法 | |
CN103490772B (zh) | 一种基于无功补偿的滑动加权单相软锁相方法 | |
CN101582692A (zh) | 一种提高数字锁相环性能的方法 | |
CN103472302B (zh) | 用单相光伏并网逆变器检测电网电压相位的方法 | |
CN107623522B (zh) | 一种基于d-q变换的双二阶广义积分锁相环控制方法 | |
CN109818370B (zh) | 一种带有超前矫正的二阶广义积分锁频环控制方法 | |
CN109245103A (zh) | 基于改进滑动平均值滤波器的软件锁相环实现方法及装置 | |
CN108494714A (zh) | 一种快速克服多普勒频移的gmsk相干解调的方法 | |
CN103152012B (zh) | 一种降低通带群延迟误差的fir滤波器设计方法 | |
CN104730339A (zh) | 一种三相电压不对称情况下数字锁相方法 | |
CN106209701A (zh) | Alpha稳定分布噪声环境下MFSK信号符号速率估计方法及装置 | |
CN104698254A (zh) | 一种电网基波正序电压提取方法及锁相方法 | |
CN101807918A (zh) | 基于同步坐标系的单相锁相环及其实现方法 | |
CN103986458A (zh) | 一种基于重复控制的微电网单相并网锁相环控制方法 | |
CN104410407A (zh) | 一种自适应数字锁相环及锁相方法 | |
CN102868163B (zh) | 一种自适应并网变流器三相软锁相环 | |
CN103078633A (zh) | 一种单相锁相环 | |
CN103647550A (zh) | 一种用于动态电压无功补偿的锁相环方法 | |
CN103595405A (zh) | 一种基于粒子滤波的锁相环实现方法 | |
CN103592152A (zh) | 确定模型水轮机转轮叶片进口部位空化的声学方法 | |
CN109981510A (zh) | 一种适用深空通信的极低码速率遥控可靠接收方法 | |
CN106130505A (zh) | 一种可变分数时延数字滤波器的优化方法 | |
CN106249586A (zh) | 一种支持单相和三相电压跟踪的锁相环方法 | |
CN104506188A (zh) | 基于带通滤波器的抗谐波干扰型锁相环方法 | |
CN201742387U (zh) | 一种基于同步坐标系的单相锁相环 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200430 Address after: 065001 New Asia R & D building, No.106, No.1 Road, Langfang Economic and Technological Development Zone, Hebei Province Patentee after: LANGFANG YINHE XINTONG TECHNOLOGY Co.,Ltd. Address before: 100081 No. 5, Zhongguancun South Street, Haidian District, Beijing Patentee before: BEIJING INSTITUTE OF TECHNOLOGY |
|
TR01 | Transfer of patent right |