CN101984613A - 一种低速率bpsk突发信号码速率估计方法 - Google Patents
一种低速率bpsk突发信号码速率估计方法 Download PDFInfo
- Publication number
- CN101984613A CN101984613A CN2010105650775A CN201010565077A CN101984613A CN 101984613 A CN101984613 A CN 101984613A CN 2010105650775 A CN2010105650775 A CN 2010105650775A CN 201010565077 A CN201010565077 A CN 201010565077A CN 101984613 A CN101984613 A CN 101984613A
- Authority
- CN
- China
- Prior art keywords
- signal
- bit rate
- low
- bpsk
- symbol period
- 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
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明公开了一种低速率BPSK突发信号码速率估计方法,先对接收到的采样频率为fs的BPSK信号的r(n)下变频得到基带信号,由迟滞比较器对基带信号波形进行整形,再对经过整形的信号进行差分运算并由此得到高、低电平持续期序列,然后由高、低电平持续期序列进行最小二乘直线拟合,得到拟合误差,由拟合误差的大小选出符号周期的估计值,对符号周期取倒数即得到码速率的估计值;其中:n=0,1,……,N-1;本发明原理简单,计算量小,易于实现;采用直线拟合,具有更高的稳健性,而且原理简单,不涉及复杂运算,适用于工程应用;通过对实际采集的信号分析表明,由该方法得到的估计值对常规的码速率估计算法进行辅助时,取得了优越的性能。
Description
技术领域
本发明涉及数字通信中一种BPSK(Binary Phase Shift Keying,二元数字相移键控)信号码速率估计方法,具体是低速率BPSK突发信号码速率的一种估计方法。
技术背景
BPSK调制因其抗噪性能良好而在卫星通信等数字通信系统中得到广泛应用。对数字通信系统来说,发送端所采用的码速率是接收端进行解调的一项必不可少的参数,而在非合作通信信号分析中,码速率只能由接收数据估计得到,因此,如何提高码速率估计的精度对非合作通信具有极其重要的意义。在突发数据通信中,当信号的码速率特别低时,接收数据中的符号个数往往非常少,只有十几个甚至几个,此时,传统的适用于符号个数比较多的场合下的估计方法往往会失效,需要由其他方法缩小码速率估计值的取值范围,以此辅助传统的估计方法,从而提高估计精度。
目前,低速率BPSK突发信号码速率估计中,用于约束码速率取值范围的方法主要有两种:一种是在频域由带宽辅助估计;另一种是在时域将基带信号高、低电平持续时间的最小值作为符号周期来辅助估计。对低码率信号,当符号个数比较少时,带宽估计自身就是一个难点,而基带信号高、低电平持续时间的最小值有可能会是符号周期的某个倍数,从而给出错误的辅助估计范围。
发明内容
本发明旨在提出一种适用于低码率BPSK突发信号的码速率估计方法,基于上述的时域辅助估计方法,采用直线拟合,比上述方法具有更高的稳健性,而且原理简单,不涉及复杂运算,适用于工程应用。
本发明的技术方案如下:
一种低速率BPSK突发信号码速率估计方法,其特征在于:先对接收到的采样频率为 的BPSK信号下变频得到基带信号,由迟滞比较器对基带信号波形进行整形,再对经过整形的信号进行差分运算并由此得到高、低电平持续期序列,然后由高、低电平持续期序列进行最小二乘直线拟合,得到拟合误差,由拟合误差的大小选出符号周期的估计值,此处,因为码速率和符号周期呈倒数关系,所以对符号周期取倒数即得到码速率的估计值;其中:n = 0,1,……,N-1。
对所述BPSK信号的载波频率进行估计后进行正交下变频处理,是对接收到的信号BPSK信号分别乘以一个信号和一个信号并进行低通滤波,然后低通滤波输出的两路信号分别称作I路(In-phase,同相)输出信号和Q路(Quadrature,正交)输出信号。其中,是载波频率的估计值,是载波相位的估计值。
在下变频处理时,对载波相位进行扫描指的就是载波相位的估计值获取的过程,也就是载波相位估计过程,具体实施时将等间隔划分为若干份,并进行扫描。因为当=且=时,信号将与发送端载波同频同相,此时的I路输出信号的能量达到最大值,所以,可以通过I路输出信号能量最大化作为准则对载波相位进行估计,实现载波相位的同步。
当发送端基带信号不进行成形处理时,在I路输出的信号是经过信道而产生失真并叠加了噪声的方波信号;其中:n = 0,1,……,N-1。
其中:是一个正的经验参数,是一个指示函数,当条件成立时,指示函数,;反之,当条件成立时,指示函数,。这里并没有定义情况下指示函数的取值,因为出现的概率为0,同时,将的情况并入到上述两种情况中的任意一个均可。
取值为2或者-2的位置是基带信号正负电平的交替处(因为定义指示函数的值等于+1或者-1,也就是经过整形的信号的取值只有+1或-1,而又是相邻两点的差值,所以取值只能是-2或0或者+2。等于0的地方代表没有符号交替,而等于-2或者+2的地方代表有符号交替),也就是符号发生更替的边界处。因此,由可以得到符号更替处的索引值序列 (其中:j = 0, 1, ……, L,这里记的元素个数为L+1),进而得到高、低电平的持续期序列:
其中:j = 0, 1, ……, L-1。
。
高、低电平持续期中的最小值,也就是的第一项的取值近似等于符号周期T(符号周期指的是数字通信中,用于代表一个发送符号的信号的持续时间)的1倍、2倍、3倍,或者4倍,而大于4倍的概率非常低。从实验中观察到,的取值通常都近似等于符号周期T的1倍,偶尔出现近似等于符号周期2倍的情况。
记:
对排序后的高、低电平持续期序列采用最小二乘直线拟合,则有:
对e关于T求导,并令导数等于0,则可以得到符号周期的估计值:
同时,得到拟合误差为:
这样,对i从1到n_max按如上步骤扫描,总结成伪代码如下:
for i = 1 to n_max
for j = 0 to L-1
end
end
则取对应拟合误差最小的符号周期估计值为符号周期的最终估计值,若拟合误差存在若干个最小值,则取符号周期长的,也就是上述i值小的符号周期估计值为码速率的最终估计值。
本发明的有益效果如下:
本发明原理简单,计算量小,易于实现;采用直线拟合,具有更高的稳健性,而且原理简单,不涉及复杂运算,适用于工程应用;通过对实际采集的信号分析表明,由该方法得到的估计值对常规的码速率估计算法进行辅助时,取得了优越的性能。
附图说明
图1为本发明的流程图
图2为本发明的迟滞比较器输入输出特性曲线示意图
图3 为本发明的基带信号及校正后的实例信号示意图。
具体实施方式
如图1所示,一种低速率BPSK突发信号码速率估计方法,先对接收到的采样频率为的BPSK信号下变频得到基带信号,由迟滞比较器对基带信号波形进行整形,再对经过整形的信号进行差分运算并由此得到高、低电平持续期序列,然后由高、低电平持续期序列进行最小二乘直线拟合,得到拟合误差,由拟合误差的大小选出符号周期的估计值,此处,因为码速率和符号周期呈倒数关系,所以对符号周期取倒数即得到码速率的估计值;其中:n = 0,1,……,N-1。
对所述BPSK信号的载波频率进行估计后进行正交下变频处理时,得到I路输出信号和Q路输出信号,同时对载波相位进行扫描,由I路输出信号能量最大化实现载波相位的同步;当发送端基带信号不进行成形处理时,在I路输出的信号是经过信道而产生失真并叠加了噪声的方波信号;其中:n = 0,1,……,N-1。
对所述失真并叠加噪声的方波信号,利用迟滞比较器的迟滞特性进行校正,得到校正后的信号;其中:n = 0,1,……,N-1。
其中:是一个正的经验参数,实验中得到,将设置为0.9时取得了较好的性能;是一个指示函数,当条件成立时,指示函数,;反之,当条件成立时,指示函数,。这里并没有定义情况下指示函数的取值,因为出现的概率为0,同时,将的情况并入到上述两种情况中的任意一个均可。
取值为2或者-2的位置是基带信号正负电平的交替处(因为定义指示函数的值等于+1或者-1,也就是经过整形的信号的取值只有+1或-1,而又是相邻两点的差值,所以取值只能是-2或0或者+2。等于0的地方代表没有符号交替,而等于-2或者+2的地方代表有符号交替),也就是符号发生更替的边界处。因此,由可以得到符号更替处的索引值序列 (其中:j = 0, 1, ……, L,这里记的元素个数为L+1),进而得到高、低电平的持续期序列:
其中:j = 0, 1, ……, L-1。
高、低电平持续期中的最小值,也就是的第一项的取值近似等于符号周期T(符号周期指的是数字通信中,用于代表一个发送符号的信号的持续时间)的1倍、2倍、3倍,或者4倍,而大于4倍的概率非常低。从实验中观察到,的取值通常都近似等于符号周期T的1倍,偶尔出现近似等于符号周期2倍的情况。
记:
对排序后的高、低电平持续期序列采用最小二乘直线拟合,则有:
对e关于T求导,并令导数等于0,则可以得到符号周期的估计值:
同时,得到拟合误差为:
这样,对i从1到n_max按如上步骤扫描,得到拟合误差最小值对应的符号周期估计值,即为符号周期的最终估计值,若拟合误差存在若干个最小值,则取符号周期长的,也就是上述i值小的符号周期估计值为码速率的最终估计值。
对i从1到n_max扫描时,采用最小二乘直线拟合得到符号周期估计值分别为:
对应的拟合误差分别为:
Claims (10)
1.一种低速率BPSK突发信号码速率估计方法,其特征在于:先对接收到的采样频率为 的BPSK信号下变频得到基带信号,由迟滞比较器对基带信号波形进行整形,再对经过整形的信号进行差分运算并由此得到高、低电平持续期序列,然后由高、低电平持续期序列进行最小二乘直线拟合,得到拟合误差,由拟合误差的大小选出符号周期的估计值,对符号周期取倒数即得到码速率的估计值;其中:n = 0,1,……,N-1。
10.根据权利要求9所述的一种低速率BPSK突发信号码速率估计方法,其特征在于:当得到的拟合误差存在若干个最小值,则比较符号周期的长度,取所述i值小的符号周期估计值为码速率的最终估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010565077 CN101984613B (zh) | 2010-11-30 | 2010-11-30 | 一种低速率bpsk突发信号码速率估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010565077 CN101984613B (zh) | 2010-11-30 | 2010-11-30 | 一种低速率bpsk突发信号码速率估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101984613A true CN101984613A (zh) | 2011-03-09 |
CN101984613B CN101984613B (zh) | 2013-06-05 |
Family
ID=43641784
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010565077 Expired - Fee Related CN101984613B (zh) | 2010-11-30 | 2010-11-30 | 一种低速率bpsk突发信号码速率估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101984613B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103441975A (zh) * | 2013-08-30 | 2013-12-11 | 天津理工大学 | 一种基于功率谱的二相编码信号参数估值方法 |
CN106209701A (zh) * | 2016-07-01 | 2016-12-07 | 中国人民解放军信息工程大学 | Alpha稳定分布噪声环境下MFSK信号符号速率估计方法及装置 |
CN106506414A (zh) * | 2016-11-17 | 2017-03-15 | 中国电子科技集团公司第四十研究所 | 一种基于峰值位置的相位调制码速率估计方法 |
CN112162152A (zh) * | 2020-08-31 | 2021-01-01 | 南京亿杰明信息技术有限公司 | 基于相位直线拟合的正弦波相参脉冲串信号频率估计方法 |
CN112162153A (zh) * | 2020-08-31 | 2021-01-01 | 南京亿杰明信息技术有限公司 | 基于相位直线拟合的正弦波信号频率估计方法 |
CN112260693A (zh) * | 2020-09-14 | 2021-01-22 | 中国工程物理研究院电子工程研究所 | 一种数字调制信号重构方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4887007A (en) * | 1987-02-18 | 1989-12-12 | U.S. Philips Corporation | DC-AC converter for supplying a gas and/or vapour discharge lamp |
EP0932259A1 (en) * | 1998-01-27 | 1999-07-28 | Lucent Technologies Inc. | Iterative decoding on demand |
CN1499755A (zh) * | 2002-11-08 | 2004-05-26 | 深圳市中兴通讯股份有限公司上海第二 | 优化的外环功率控制方法和装置 |
CN101162910A (zh) * | 2006-10-10 | 2008-04-16 | 北京六合万通微电子技术有限公司 | 本振泄漏自动消除装置 |
-
2010
- 2010-11-30 CN CN 201010565077 patent/CN101984613B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4887007A (en) * | 1987-02-18 | 1989-12-12 | U.S. Philips Corporation | DC-AC converter for supplying a gas and/or vapour discharge lamp |
EP0932259A1 (en) * | 1998-01-27 | 1999-07-28 | Lucent Technologies Inc. | Iterative decoding on demand |
CN1499755A (zh) * | 2002-11-08 | 2004-05-26 | 深圳市中兴通讯股份有限公司上海第二 | 优化的外环功率控制方法和装置 |
CN101162910A (zh) * | 2006-10-10 | 2008-04-16 | 北京六合万通微电子技术有限公司 | 本振泄漏自动消除装置 |
Non-Patent Citations (1)
Title |
---|
王甲峰: "SDPSK调制参数的盲估计", 《无线电工程》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103441975A (zh) * | 2013-08-30 | 2013-12-11 | 天津理工大学 | 一种基于功率谱的二相编码信号参数估值方法 |
CN103441975B (zh) * | 2013-08-30 | 2016-12-28 | 天津理工大学 | 一种基于功率谱的二相编码信号参数估值方法 |
CN106209701A (zh) * | 2016-07-01 | 2016-12-07 | 中国人民解放军信息工程大学 | Alpha稳定分布噪声环境下MFSK信号符号速率估计方法及装置 |
CN106209701B (zh) * | 2016-07-01 | 2019-08-02 | 中国人民解放军信息工程大学 | Alpha稳定分布噪声环境下MFSK信号符号速率估计方法及装置 |
CN106506414A (zh) * | 2016-11-17 | 2017-03-15 | 中国电子科技集团公司第四十研究所 | 一种基于峰值位置的相位调制码速率估计方法 |
CN106506414B (zh) * | 2016-11-17 | 2019-04-09 | 中国电子科技集团公司第四十一研究所 | 一种基于峰值位置的相位调制码速率估计方法 |
CN112162152A (zh) * | 2020-08-31 | 2021-01-01 | 南京亿杰明信息技术有限公司 | 基于相位直线拟合的正弦波相参脉冲串信号频率估计方法 |
CN112162153A (zh) * | 2020-08-31 | 2021-01-01 | 南京亿杰明信息技术有限公司 | 基于相位直线拟合的正弦波信号频率估计方法 |
CN112162152B (zh) * | 2020-08-31 | 2024-01-26 | 南京亿杰明信息技术有限公司 | 基于相位直线拟合的正弦波相参脉冲串信号频率估计方法 |
CN112260693A (zh) * | 2020-09-14 | 2021-01-22 | 中国工程物理研究院电子工程研究所 | 一种数字调制信号重构方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101984613B (zh) | 2013-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101984613B (zh) | 一种低速率bpsk突发信号码速率估计方法 | |
US9954701B2 (en) | Bluetooth signal receiving method and device using improved packet detection and symbol timing acquisition | |
CN109981506B (zh) | 一种基于开环结构的全数字解调方法 | |
KR102041342B1 (ko) | 다중 상관 검출을 이용하여 패킷 검출, 심볼 타이밍 포착 및 반송파 주파수 오프셋 추정을 동시에 병렬 수행하는 방법 및 그 블루투스 장치 | |
KR102269195B1 (ko) | 변조된 신호의 수신자의 초기 타이밍 동기화를 수행하는 시스템 및 방법 | |
CN107968757B (zh) | 一种频移键控调制信号的解调方法及系统 | |
CN103209143B (zh) | 一种数据辅助下的星载ais信号同步参数估计方法及系统 | |
CN101184071B (zh) | 基于伪误码率统计的盲信噪比估计方法 | |
CN103023831B (zh) | 一种适用于突发波形的载波频偏估计方法 | |
CN103117965A (zh) | 一种星载ais信号定时频偏联合估计方法及其实现系统 | |
CN107342960B (zh) | 一种适合幅度相移键控的非数据辅助频偏估计方法 | |
CN104852876A (zh) | 一种航空无线突发通信系统 | |
US7477707B2 (en) | Computationally efficient demodulation for differential phase shift keying | |
CN101087281A (zh) | 一种正交幅度调制的信噪比和幅噪比的测量方法及装置 | |
US20060209993A1 (en) | Demodulator and receiver for pre-coded partial response signals | |
CN111901269A (zh) | 可变调制指数的高斯频移键控调制方法、装置及系统 | |
CN113132285A (zh) | 一种数字解调系统及方法 | |
CN108833321B (zh) | 基于差分相位波形匹配的编码cpm信号码块同步方法 | |
CN102378354A (zh) | 一种用于短时突发通信的同步方法 | |
CN101404564A (zh) | 8psk格雷映射的一种软解调方法 | |
CN108718461B (zh) | 一种抗频偏的突发cpm信号帧同步方法 | |
US6874096B1 (en) | Apparatus and method for detecting packet arrival time | |
EP0789958A1 (en) | Method of, and apparatus for, symbol timing recovery | |
US9059835B1 (en) | Method and apparatus for demodulating a wireless signal | |
CN110535620B (zh) | 一种基于判决反馈的信号检测与同步方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130605 Termination date: 20131130 |