CN109286588B - Method for estimating initial phase of binary phase shift keying signal - Google Patents

Method for estimating initial phase of binary phase shift keying signal Download PDF

Info

Publication number
CN109286588B
CN109286588B CN201710598645.3A CN201710598645A CN109286588B CN 109286588 B CN109286588 B CN 109286588B CN 201710598645 A CN201710598645 A CN 201710598645A CN 109286588 B CN109286588 B CN 109286588B
Authority
CN
China
Prior art keywords
shift keying
module
initial phase
sampling
phase shift
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.)
Active
Application number
CN201710598645.3A
Other languages
Chinese (zh)
Other versions
CN109286588A (en
Inventor
段云鹏
王平
朱智明
何靖
王百涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Remote Sensing Equipment
Original Assignee
Beijing Institute of Remote Sensing Equipment
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Remote Sensing Equipment filed Critical Beijing Institute of Remote Sensing Equipment
Priority to CN201710598645.3A priority Critical patent/CN109286588B/en
Publication of CN109286588A publication Critical patent/CN109286588A/en
Application granted granted Critical
Publication of CN109286588B publication Critical patent/CN109286588B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0014Carrier regulation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0014Carrier regulation
    • H04L2027/0024Carrier regulation at the receiver end
    • H04L2027/0026Correction of carrier offset
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/18Phase-modulated carrier systems, i.e. using phase-shift keying
    • H04L27/22Demodulator circuits; Receiver circuits

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

The invention discloses a binary phase shift keying signal initial phase estimation method, which firstly constructs a binary phase shift keying signal initial phase estimation platform, and the constructed system comprises: the system comprises a signal sampling module, an FFT module, a carrier frequency estimation module and an initial phase estimation module, wherein the signal sampling module is used for sampling binary phase shift keying signals to form a sampling sequence, the FFT module is used for carrying out FFT processing on the sampling sequence, the carrier frequency estimation module is used for estimating a carrier frequency according to the FFT processing result, and the initial phase estimation module is used for estimating the initial phase of the binary phase shift keying signals. The method can estimate the initial phase of the binary phase shift keying signal, has high accuracy of the estimated initial phase and can meet the use requirement.

Description

Method for estimating initial phase of binary phase shift keying signal
Technical Field
The invention relates to a signal initial phase estimation method, in particular to a binary phase shift keying signal initial phase estimation method.
Background
After the spread spectrum receiver works, it needs to perform initial phase parameter estimation on the received binary phase shift keying signal, so as to perform subsequent signal processing work. The binary phase shift keying signal belongs to phase modulation, the phase of the modulated signal is directly controlled by the digital baseband signal, and the amplitude and frequency of the modulated signal are kept unchanged. The commonly used method for estimating the initial phase of the signal is as follows: sampling the binary phase shift keying signal to form a sampling sequence; carrying out carrier frequency estimation and initial position estimation on the sampling sequence; performing digital quadrature frequency mixing processing on the sampling sequence by using a quadrature digital local oscillator to form in-phase data and quadrature data; initial phase estimation is performed using the in-phase and quadrature data based on the estimated initial position. Since the initial position estimation accuracy of the binary phase shift keying signal by the method is not very high, the initial phase estimation performance is reduced.
Fast fourier transform, abbreviated FFT, can convert time domain data to frequency domain data.
Disclosure of Invention
The invention aims to provide a binary phase shift keying signal initial phase estimation method, which solves the problem of poor estimation performance of a common signal initial phase estimation method.
A method for estimating an initial phase of a binary phase shift keying signal comprises the following specific steps:
first step, an initial phase estimation platform is established
An initial phase estimation stage comprising: the device comprises a signal sampling module, an FFT module, a carrier frequency estimation module and an initial phase estimation module. The following steps:
the signal sampling module has the functions of: sampling the binary phase shift keying signal to form a sampling sequence;
the function of the FFT module is as follows: performing FFT processing on the sampling sequence;
the carrier frequency estimation module has the functions of: estimating carrier frequency according to the FFT processing result;
the initial phase estimation module functions as: an initial phase of the binary phase shift keyed signal is estimated.
The second step signal sampling module samples the binary phase shift keying signal to form a sampling sequence
The signal sampling module samples a binary phase shift keying signal x (t) to form a sampling sequence x (n) with a sampling frequency fsWith a sampling time interval of Ts,Ts=1/fs(ii) a t is a time parameter, N is a time domain point index value, N is 0,1, …, N-1, and N is the length of the sampling sequence x (N).
The third step is that the FFT module carries out FFT processing on the sampling sequence
The FFT module performs FFT processing on the sample sequence x (n), transforms the sample sequence x (n) from the time domain to the frequency domain, and obtains an FFT processing result x (k):
Figure BDA0001356527270000021
wherein k is a frequency domain point index value, k is 0,1, …, N-1; j is an imaginary unit, representing
Figure BDA0001356527270000022
e‐j2πkn/NExpressed as complex numbers, cos (-2 π kn/N) + jsin (-2 π kn/N).
The fourth step carrier frequency estimation module estimates the carrier frequency according to the FFT processing result
The carrier frequency estimation module screens out the maximum value in the FFT processing result X (k), thereby determining the frequency domain point index value k corresponding to the maximum value0According to the formula: f. of0=k0/(NTs) Determining k0Frequency f corresponding to a point0
The fifth step is that the initial phase estimation module estimates the initial phase of the binary phase shift keying signal
The initial phase estimation module estimates the initial phase of the binary phase shift keying signal
Figure BDA0001356527270000023
Thus, estimation of the initial phase of the binary phase shift keying signal is achieved.
Preferably, the length N of the sampling sequence in the second step is 1024.
Preferably, the initial phase of the binary phase shift keyed signal in the fifth step
Figure BDA0001356527270000024
Obtained by the following formula:
Figure BDA0001356527270000025
wherein angle {. cndot } represents the angle of a complex number,
Figure BDA0001356527270000026
in complex representation form, represent: cos (-4 pi f)0nTs)+jsin(‐4πf0nTs)。
The method solves the problem of poor estimation performance of the common signal initial phase estimation method. The initial phase of the binary phase shift keying signal can be estimated, the estimated initial phase is high in precision, and the use requirement can be met.
Detailed Description
A method for estimating an initial phase of a binary phase shift keying signal comprises the following specific steps:
first step, an initial phase estimation platform is established
An initial phase estimation stage comprising: the device comprises a signal sampling module, an FFT module, a carrier frequency estimation module and an initial phase estimation module. The following steps:
the signal sampling module has the functions of: sampling the binary phase shift keying signal to form a sampling sequence;
the function of the FFT module is as follows: performing FFT processing on the sampling sequence;
the carrier frequency estimation module has the functions of: estimating carrier frequency according to the FFT processing result;
the initial phase estimation module functions as: an initial phase of the binary phase shift keyed signal is estimated.
The second step signal sampling module samples the binary phase shift keying signal to form a sampling sequence
The signal sampling module samples a binary phase shift keying signal x (t) to form a sampling sequence x (n) with a sampling frequency fsWith a sampling time interval of Ts,Ts=1/fs(ii) a t is a time parameter, N is a time domain point index value, N is 0,1, …, N-1, N is the length of the sample sequence x (N), and N is 1024.
The third step is that the FFT module carries out FFT processing on the sampling sequence
The FFT module performs FFT processing on the sample sequence x (n), transforms the sample sequence x (n) from the time domain to the frequency domain, and obtains an FFT processing result x (k):
Figure BDA0001356527270000031
wherein k is a frequency domain point index value, k is 0,1, …, N-1; j is an imaginary unit, representing
Figure BDA0001356527270000032
e‐j2πkn/NExpressed as complex numbers, cos (-2 π kn/N) + jsin (-2 π kn/N).
The fourth step carrier frequency estimation module estimates the carrier frequency according to the FFT processing result
The carrier frequency estimation module screens out the maximum value in the FFT processing result X (k), thereby determining the frequency domain point index value k corresponding to the maximum value0According to the formula: f. of0=k0/(NTs) Determining k0Frequency f corresponding to a point0
The fifth step is that the initial phase estimation module estimates the initial phase of the binary phase shift keying signal
The initial phase estimation module uses the formula:
Figure BDA0001356527270000033
estimating an initial phase of a binary phase shift keyed signal
Figure BDA0001356527270000034
Wherein angle {. cndot } represents the angle of a complex number,
Figure BDA0001356527270000035
in complex representation form, represent: cos (-4 pi f)0nTs)+jsin(‐4πf0nTs)。
Thus, estimation of the initial phase of the binary phase shift keying signal is achieved.

Claims (1)

1. A binary phase shift keying signal initial phase estimation method is characterized by comprising the following specific steps:
firstly, constructing an initial phase estimation platform of binary phase shift keying signals
A binary phase shift keying signal initial phase estimation stage comprising: the device comprises a signal sampling module, an FFT module, a carrier frequency estimation module and an initial phase estimation module; the following steps:
the signal sampling module has the functions of: sampling the binary phase shift keying signal to form a sampling sequence;
the function of the FFT module is as follows: performing FFT processing on the sampling sequence;
the carrier frequency estimation module has the functions of: estimating carrier frequency according to the FFT processing result;
the initial phase estimation module functions as: estimating an initial phase of the binary phase shift keying signal;
the second step signal sampling module samples the binary phase shift keying signal to form a sampling sequence
The signal sampling module samples a binary phase shift keying signal x (t) to form a sampling sequence x (n) with a sampling frequency fsWith a sampling time interval of Ts,Ts=1/fs(ii) a t is a time parameter, N is a time domain point index value, N is 0,1, …, N-1, and N is the length of the sampling sequence x (N);
the third step is that the FFT module carries out FFT processing on the sampling sequence
The FFT module performs FFT processing on the sample sequence x (n), transforms the sample sequence x (n) from the time domain to the frequency domain, and obtains an FFT processing result x (k):
Figure FDA0002756039750000011
wherein k is a frequency domain point index value, k is 0,1, …, N-1; j is an imaginary unit, representing
Figure FDA0002756039750000012
e-j2πkn/NExpressed in complex number representation, representing cos (-2 π kn/N) + jsin (-2 π kn/N);
the fourth step carrier frequency estimation module estimates the carrier frequency according to the FFT processing result
The carrier frequency estimation module screens out the maximum value in the FFT processing result X (k), thereby determining the frequency domain point index value k corresponding to the maximum value0According to the formula: f. of0=k0/(NTs) Determining k0Frequency f corresponding to a point0
The fifth step is that the initial phase estimation module estimates the initial phase of the binary phase shift keying signal
The initial phase estimation module estimates the initial phase of the binary phase shift keying signal
Figure FDA0002756039750000013
Obtained by the following formula:
Figure FDA0002756039750000014
wherein angle {. cndot } represents the angle of a complex number,
Figure FDA0002756039750000015
in complex representation form, represent: cos (-4 pi f)0nTs)+jsin(-4πf0nTs)
Thus, estimation of the initial phase of the binary phase shift keying signal is achieved.
CN201710598645.3A 2017-07-21 2017-07-21 Method for estimating initial phase of binary phase shift keying signal Active CN109286588B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710598645.3A CN109286588B (en) 2017-07-21 2017-07-21 Method for estimating initial phase of binary phase shift keying signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710598645.3A CN109286588B (en) 2017-07-21 2017-07-21 Method for estimating initial phase of binary phase shift keying signal

Publications (2)

Publication Number Publication Date
CN109286588A CN109286588A (en) 2019-01-29
CN109286588B true CN109286588B (en) 2021-02-26

Family

ID=65185143

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710598645.3A Active CN109286588B (en) 2017-07-21 2017-07-21 Method for estimating initial phase of binary phase shift keying signal

Country Status (1)

Country Link
CN (1) CN109286588B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1463084A (en) * 2003-06-18 2003-12-24 中国人民解放军理工大学通信工程学院 Method and device of iterative demodulation and decode for BPSK modulating system by Turbo encoding
CN1543088A (en) * 2003-11-07 2004-11-03 中国人民解放军理工大学通信工程学院 Iterative demodulating-decoding method and apparatus for code modulation system having output external information
CN1812295A (en) * 2005-11-28 2006-08-02 北京星河亮点通信软件有限责任公司 Method for testing 3G terminal radio frequency consistency
CN101751375A (en) * 2008-12-12 2010-06-23 普天信息技术研究院有限公司 Fast calculation method and device of discrete Fourier transformation (DFT)/inverse discrete Fourier transform (IDFT)
CN102571668A (en) * 2012-01-06 2012-07-11 合肥东芯通信股份有限公司 Phase compensation method and system in LTE (long term evolution) system
CN105099978A (en) * 2015-07-14 2015-11-25 华中科技大学 Method for removing phase modulation information

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7899126B2 (en) * 2007-07-27 2011-03-01 Mediatek Inc. Signal receiver and method for estimating residual doppler frequencies of signals thereof

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1463084A (en) * 2003-06-18 2003-12-24 中国人民解放军理工大学通信工程学院 Method and device of iterative demodulation and decode for BPSK modulating system by Turbo encoding
CN1543088A (en) * 2003-11-07 2004-11-03 中国人民解放军理工大学通信工程学院 Iterative demodulating-decoding method and apparatus for code modulation system having output external information
CN1812295A (en) * 2005-11-28 2006-08-02 北京星河亮点通信软件有限责任公司 Method for testing 3G terminal radio frequency consistency
CN101751375A (en) * 2008-12-12 2010-06-23 普天信息技术研究院有限公司 Fast calculation method and device of discrete Fourier transformation (DFT)/inverse discrete Fourier transform (IDFT)
CN102571668A (en) * 2012-01-06 2012-07-11 合肥东芯通信股份有限公司 Phase compensation method and system in LTE (long term evolution) system
CN105099978A (en) * 2015-07-14 2015-11-25 华中科技大学 Method for removing phase modulation information

Also Published As

Publication number Publication date
CN109286588A (en) 2019-01-29

Similar Documents

Publication Publication Date Title
US7639750B2 (en) Phase tracking method and device thereof
CN106209701B (en) MFSK signal code rate-estimation method and device under Alpha Stable distritation noise circumstance
CN106059973B (en) Frequency deviation estimating method and system
CN111740932B (en) Short wave communication synchronization head capturing method under extremely low signal-to-noise ratio
WO2015070820A1 (en) Spread spectrum signal generating method, generating apparatus, receiving method and receiving apparatus
CN109067680B (en) Carrier frequency offset estimation method and device of baseband signal
CN110300079B (en) MSK signal coherent demodulation method and system
CN106788734B (en) Optical OFDM system adopting data-free auxiliary frequency offset estimation algorithm
CN108566353B (en) Continuously corrected carrier synchronization device and method
CN105162533B (en) Transmitter amplitude imbalance and phase imbalance measuring method
CN109286588B (en) Method for estimating initial phase of binary phase shift keying signal
CN108562918B (en) BOC (n, n) ambiguity-free capturing method and device based on correlation shift
CN102137053A (en) Method for estimating signal to noise ratio of BPSK (Binary Phase Shift Keying) signal
WO2017152438A1 (en) Systems and methods for time synchronization between transmitters and receivers in a communication system
EP2922262B1 (en) Electronic circuit with synchronous demodulator of a phase-modulated signal
CN106569182B (en) Phase-coded signal carrier frequency estimation method based on minimum entropy
CN109660481B (en) Quadrature amplitude modulation and demodulation method and constellation diagram optimization method thereof
WO2018057320A1 (en) Wideband residual sideband calibration
CN103873225A (en) Timing estimation method for burst communication
CN105743830B (en) A kind of demodulation method of MFSK signals
CN105406921A (en) MFSK water sound communication method and MFSK water sound communication system
CN108120974B (en) Pulse compression network generation method for multi-time-width linear frequency modulation signal
KR102625211B1 (en) Apparatus And Method for Calculating Receiving Time Of Wireless Communication Signal
CN105704074A (en) Maximum likelihood sequence detection in the phase domain
CN101646232A (en) Method, device and communication equipment for estimating frequency deviations

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant