CN104333525A - GMSK (Gaussian minimum shift keying) modulating system synchronization method - Google Patents

GMSK (Gaussian minimum shift keying) modulating system synchronization method Download PDF

Info

Publication number
CN104333525A
CN104333525A CN201410572532.2A CN201410572532A CN104333525A CN 104333525 A CN104333525 A CN 104333525A CN 201410572532 A CN201410572532 A CN 201410572532A CN 104333525 A CN104333525 A CN 104333525A
Authority
CN
China
Prior art keywords
gmsk
data
reference signal
local reference
vector space
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
Application number
CN201410572532.2A
Other languages
Chinese (zh)
Other versions
CN104333525B (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.)
Guangzhou Haige Communication Group Inc Co
Original Assignee
Guangzhou Haige Communication Group Inc Co
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 Guangzhou Haige Communication Group Inc Co filed Critical Guangzhou Haige Communication Group Inc Co
Priority to CN201410572532.2A priority Critical patent/CN104333525B/en
Publication of CN104333525A publication Critical patent/CN104333525A/en
Application granted granted Critical
Publication of CN104333525B publication Critical patent/CN104333525B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

The invention discloses a GMSK (Gaussian minimum shift keying) modulating system synchronization method which includes the steps: performing differential operation for GMSK modulating signals generated by a training sequence of preset numerous burst types to obtain a local reference signal; receiving GMSK data, marking current GMSK data receiving points as sampling points when the length of a received GMSK data segment is as same as that of the local reference signal, and performing differential operation for the current GMSK data segment to obtain a GMSK differential signal; gaining vector space cosine similarity according to the local reference signal and the GMSK differential signal in the same burst type; judging whether the vector space cosine similarity exceeds a threshold or not, continuing to receive the GMSK data if not, and then selecting one synchronization position from alternative synchronization positions to perform system synchronization within preset GMSK data receiving time. If the vector space cosine similarity exceeds the threshold, the burst type is the burst type of the currently received GMSK data, and the current sampling points are the alternative synchronization positions. Reliable synchronization can be realized, and the requirement of engineering application is more effectively met.

Description

A kind of GMSK modulation system synchronization method
Technical field
The present invention relates to digital communication research field, particularly a kind of GMSK modulation system synchronization method.
Background technology
GMSK (Gaussian filtered MSK) signal is the premodulated that have passed through gauss low frequency filter before generating msk signal, and it and MSK belong to permanent envelope digital modulation.Permanent envelope digital modulation can use the nonlinear amplifier that power efficiency is high, avoids non-constant enveloped modulation to use the secondary lobe regeneration of the filtering caused during the high nonlinear amplifier of power efficiency.
Msk signal can regard the continuous phase fsk signal that modulation index is 0.5 as, although it has the advantage such as bandwidth of permanent envelope, Phase Continuation, relative narrower, does not still meet the demands for its out-of-band radiation some communication system.For effective inhibition zone external radiation, therefore before generating msk signal, add the gauss low frequency filter of upper level premodulated.
The waveform of GMSK modulation signal is typically expressed as:
In above formula, f cit is carrier frequency; it is carrier phase; E is a code-element period T binterior signal energy, for random initial phase, remain unchanged during whole burst modulation, wherein:
In formula, a k=± 1,
q ( t ) = 1 2 T b ∫ - ∞ t g ( τ - T b 2 ) dτ ;
In above formula, g (t) is the response of Gaussian filter to rectangular pulse.
For coherent demodulation, because Crystal Oscillator Errors, Doppler frequency-shift etc. can cause estimating carrier frequencies error and the unknown initial phase exists, directly carry out related calculation cannot obtain accurately relevant peaks using as synchronization metric.Meanwhile, for a certain specific applied environment, suppose that the energy intensity of noise is relatively stable, the signal energy intensity received for different communication distances is indefinite, can cause the significantly change of correlation, and then dependent threshold also cannot provide definite value, unfavorable Project Realization.The particularly fast-changing application scenario of communication distance.
Such as, a kind of existing adaptive demodulation method is: carry out upset process to the I received, Q signal, relevant method of sliding is adopted to carry out channel estimating to the signal after upset, Delay Synchronization is carried out according to channel estimation value, determine Timing Advance, and according to the energy of channel parameter and the type determining channel, when the channel parameter on maximal correlation position energy and maximum energy and ratio E ratiowhen exceeding certain threshold value, adopt zero forcing equalization method to carry out demodulation, export soft-decision-value simultaneously, work as E ratiowhen being less than certain threshold value, Equilized demodulation method is adopted to carry out demodulation.The method does not eliminate the interference that frequency deviation and initial phase cause.
Again such as, a kind of GMSK signal acceptance method is in addition: carried out " precoding " raw information before transmitting terminal carries out GMSK modulation, the two Chirp signals sent for transmitting terminal carry out precoding, frequency difference and experiment is estimated at receiving terminal, then on the basis correcting large frequency-difference, adopt least square method to carry out channel impulse response estimation, finally according to the remaining frequency difference in channel impulse response correction signal and skew, realize the semi-coherent demodulation of GMSK signal.The method can less frequency difference and skew in erasure signal, but the degree of calculating is complicated, and need arrange multiple threshold value, to energy intensity sensitive, and the not actual use of utilizing works.
To sum up, the shortcoming existing for method of existing related synchronization: (1) does not eliminate the interference that frequency deviation and initial phase cause.(2) synchronization metric is to energy intensity sensitive, is unfavorable for Project Realization.
Summary of the invention
Main purpose of the present invention is that the shortcoming overcoming prior art is with not enough, a kind of GMSK modulation system synchronization method is provided, it may be used for the receiving terminal of GMSK modulation system, realize synchronous reliably, particularly be applied in the middle of coherent demodulation, can the interference that causes of frequency offset resistant and initial phase, and to the GMSK signal energy Strength Changes received without specific requirement, make decision threshold comparatively stable, more meet the requirement of engineer applied.
Object of the present invention is realized by following technical scheme: a kind of GMSK modulation system synchronization method, comprises the following steps:
(1) calculus of differences is done to the GMSK modulation signal that the training sequence presetting several outburst types generates, obtain the local reference signal corresponding with outburst type;
(2) GMSK data are received, when the length of the GMSK data segment received is identical with the length of described local reference signal, marking current GMSK data receiver point is sampled point, and calculus of differences is carried out to the GMSK data segment of current reception, obtain the GMSK differential signal corresponding with outburst type;
(3) according to the local reference signal under same outburst type and GMSK differential signal, the vector space cosine similarity that this outburst type is corresponding is asked for;
(4) judge whether described vector space cosine similarity exceedes thresholding corresponding to this outburst type, if exceed described thresholding, this outburst type is the outburst type of the GMSK data of current reception, current sampled point is alternative sync bit, enters step (5); If be no more than described pre-determined threshold, perform step (2);
(5) in the GMSK data receipt time preset, from alternative sync bit, one of them sync bit is selected to carry out system synchronization.
Preferably, in described step (5), the GMSK data receipt time preset is a symbol period, selects the sync bit that wherein vector space cosine similarity is maximum to carry out system synchronization.Symbol period modulates the symbol time interval, alternatively is-symbol length or code-element period exactly.
Concrete, in described step (3), the calculation procedure of vector space cosine similarity is as follows:
If a certain GMSK differential signal and local reference signal are respectively X and Y, data length is n, if i-th data is X (i)=a (i)+jb (i) in GMSK differential signal, in local reference signal, i-th data is Y (i)=c (i)+jd (i), then the vector space cosine similarity computing formula that this outburst type is corresponding is:
ρ = Σ i = 1 n a ( i ) · c ( i ) + b ( i ) · d ( i ) | X ( i ) | · | Y ( i ) | .
Preferably, in described step (4), according to the transmission quality of current outburst type and GMSK data receiver channel, a selected factor alpha, 0.7 < α < 1, threshold value using α n as vector space cosine similarity, wherein n is the length of local reference signal.
Preferably, described step (2), GMSK differential signal adopts the mode of stream process to carry out synchronous process and obtains, step is: establish different outburst types to there being m kind local reference signal, then design m shifting cache district, shifting cache section length is consistent with the local reference signal length of each outburst type, to the calculus of differences that the GMSK data segment of each shifting cache district storage walks abreast; After reception GMSK data segment carries out calculus of differences, last bit data of each shifting cache district is shifted out and is shifted to remaining data, then by first position of new GMSK data stored in each buffer area.
Preferably, described step (1), is solidificated in the receiving terminal of signal by the local reference signal that different outburst type generates.Because once the standard that system adopts is determined, so midamble code is also just determined, the local reference signal of generation is also determined, so directly signal is solidificated in receiving terminal, need not calculate, can shorten the time of signal transacting thereupon at every turn.
Compared with prior art, tool has the following advantages and beneficial effect in the present invention:
1, the GMSK data of the present invention to the GMSK modulation signal generated and reception do calculus of differences, as can be seen from the generation block diagram of GMSK signal, what do that the signals revivification after calculus of differences goes out is by the signal after gauss low frequency filter, before integration, with this signal as the interference asking the local reference signal of synchronization metric can suppress first phase and frequency deviation to cause.
2, the present invention adopts vector space cosine similarity to judge the outburst type of GMSK data, the cosine value of what vector space cosine similarity was asked the is space angle of two complex signals, carry out measuring from the angle of two complex signals, hardly by the impact of energy intensity size, be only subject to the impact of 2 π Δ fT and signal to noise ratio.Thus the present invention can make synchronous thresholding be set to definite value, greatly improves the robustness of thresholding, solves Received signal strength energy intensity and change the problem that the synchronous thresholding caused cannot be set to definite value, be beneficial to Project Realization.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of the inventive method.
Fig. 2 is that in the present embodiment, GMSK signal generates schematic diagram.
Embodiment
Below in conjunction with embodiment and accompanying drawing, the present invention is described in further detail, but embodiments of the present invention are not limited thereto.
Embodiment 1
As shown in Figure 1, a kind of GMSK modulation system synchronization method of the present embodiment, its implementation process specifically comprises the following steps:
(1), according to the midamble code of the different outburst types of the different systems defined of use GMSK modulation, generate corresponding GMSK modulation signal respectively, as shown in Figure 2, its modulated process is classical GMSK orthogonal demodulation signal generating principle to its modulated process, therefore is not repeated herein.Then do calculus of differences to the GMSK modulation signal generated, as can be seen from the generation block diagram of GMSK signal, what do that the signals revivification after calculus of differences goes out is by the signal after gauss low frequency filter, before integration.
In the baseband signal that local carrier carries out after frequency spectrum shift to received signal, there is difference frequency partially and first phase.If do not consider noise, the phase place of baseband signal can simply be expressed as:
Δ f represents frequency deviation, represent phase noise.After difference processing, the phase place of baseband signal then becomes:
Wherein &theta; k - i = &pi; 2 &Integral; kT - T kT g ( &tau; - iT ) d&tau; .
T is a code-element period, normally a very little value, therefore when frequency deviation is not very large, 2 π Δ fT are very little on the impact of overall phase place.Therefore, with this signal as the interference asking the local reference signal of synchronization metric can suppress first phase and frequency deviation to cause.Once the standard that system adopts is determined, so midamble code is also just determined, the local reference signal of generation is also determined thereupon.Therefore, the local reference signal that different outburst type generates can be solidificated in the receiving terminal of signal.
(2), GMSK data are received, when the length of the GMSK data segment received is identical with the length of its corresponding local reference signal, marking current GMSK data receiver point is sampled point, and calculus of differences is carried out to current GMSK data segment, obtain the GMSK differential signal corresponding with outburst type.The mode of this example employing parallel processing realizes the calculus of differences to current GMSK data segment, concrete steps are: establish different outburst types to there being m kind local reference signal, then can design m shifting cache district, shifting cache section length is consistent with the local reference signal length of each outburst type.To the calculus of differences that the GMSK data segment of each shifting cache district storage walks abreast.After reception GMSK data segment carries out calculus of differences, last bit data of each shifting cache district is shifted out and is shifted to remaining data, then by first position of new GMSK data stored in each buffer area.
(3), according to the local reference signal under same outburst type and GMSK differential signal, the vector space cosine similarity that this outburst type is corresponding is asked for.
If a certain GMSK differential signal and local reference signal are respectively X and Y, data length is n, if i-th data is that in X (i)=a (i)+jb (i), Y, i-th data is Y (i)=c (i)+jd (i) in X.The vector space cosine similarity computing formula that then this outburst type is corresponding is:
&rho; = &Sigma; i = 1 n a ( i ) &CenterDot; c ( i ) + b ( i ) &CenterDot; d ( i ) | X ( i ) | &CenterDot; | Y ( i ) | .
As can be seen from formula, the cosine value of what vector space cosine similarity was asked the is space angle of two complex signals carries out measuring from the angle of two complex signals.If do not consider the impact of noise, if the complex signal direction that I, Q two paths of signals is formed is consistent, then represent that these two complex signal similarities are very large, difference is the size of energy intensity.Now, the vector space cosine similarity of trying to achieve is cos (0)=1.Therefore, vector space cosine similarity value is larger, represents that the similarity of two signals is larger.
(4), judge whether vector space cosine similarity exceedes thresholding corresponding to this outburst type, if cross thresholding, illustrate that this outburst type is the outburst type of the GMSK data of current reception, current sampled point is alternative sync bit, performs step (5); If but thresholding, perform step (2).
Choosing of threshold value will decide according to the length of local reference signal.Do not considering under noise situations, if two complex data directions are consistent, namely angle is 0, and its vector space cosine value reaches maximum 1.Therefore be the burst of n for local reference signal length, the maximum of its ρ is n.Consider the impact of the factor such as void inspection rate, loss index of noise, system requirements, comprehensively a factor alpha can be selected when reality uses, 0.7 < α < 1, the thresholding using α n as vector space cosine similarity.
(5), wait for a symbol period, the sync bit selecting vector space cosine similarity maximum from several the alternative sync bits obtained carries out system synchronization.
Above-described embodiment is the present invention's preferably execution mode; but embodiments of the present invention are not restricted to the described embodiments; change, the modification done under other any does not deviate from Spirit Essence of the present invention and principle, substitute, combine, simplify; all should be the substitute mode of equivalence, be included within protection scope of the present invention.

Claims (5)

1. a GMSK modulation system synchronization method, is characterized in that, comprises the following steps:
(1) calculus of differences is done to the GMSK modulation signal that the training sequence presetting several outburst types generates, obtain the local reference signal corresponding with outburst type;
(2) GMSK data are received, when the length of the GMSK data segment received is identical with the length of described local reference signal, marking current GMSK data receiver point is sampled point, and calculus of differences is carried out to the GMSK data segment of current reception, obtain the GMSK differential signal corresponding with outburst type;
(3) according to the local reference signal under same outburst type and GMSK differential signal, the vector space cosine similarity that this outburst type is corresponding is asked for;
(4) judge whether described vector space cosine similarity exceedes thresholding corresponding to this outburst type, if exceed described thresholding, this outburst type is the outburst type of the GMSK data of current reception, current sampled point is alternative sync bit, enters step (5); If be no more than described pre-determined threshold, perform step (2);
(5) in the GMSK data receipt time preset, from alternative sync bit, one of them sync bit is selected to carry out system synchronization.
2. GMSK modulation system synchronization method according to claim 1, it is characterized in that, in described step (5), the GMSK data receipt time preset is a symbol period, selects the sync bit that wherein vector space cosine similarity is maximum to carry out system synchronization.
3. GMSK modulation system synchronization method according to claim 1, is characterized in that, in described step (3), the calculation procedure of vector space cosine similarity is as follows:
If a certain GMSK differential signal and local reference signal are respectively X and Y, data length is n, if i-th data is X (i)=a (i)+jb (i) in GMSK differential signal, in local reference signal, i-th data is Y (i)=c (i)+jd (i), then the vector space cosine similarity computing formula that this outburst type is corresponding is:
&rho; = &Sigma; i = 1 n a ( i ) &CenterDot; c ( i ) + b ( i ) &CenterDot; d ( i ) | X ( i ) | &CenterDot; | Y ( i ) | .
4. GMSK modulation system synchronization method according to claim 1, it is characterized in that, in described step (4), according to the transmission quality of current outburst type and GMSK data receiver channel, a selected factor alpha, 0.7 < α < 1, the threshold value using α n as vector space cosine similarity, wherein n is the length of local reference signal.
5. GMSK modulation system synchronization method according to claim 1, it is characterized in that, described step (2), GMSK differential signal adopts the mode of stream process to carry out synchronous process and obtains, step is: establish different outburst types to there being m kind local reference signal, then design m shifting cache district, shifting cache section length is consistent with the local reference signal length of each outburst type, to the calculus of differences that the GMSK data segment of each shifting cache district storage walks abreast; After reception GMSK data segment carries out calculus of differences, last bit data of each shifting cache district is shifted out and is shifted to remaining data, then by first position of new GMSK data stored in each buffer area.
CN201410572532.2A 2014-10-23 2014-10-23 A kind of GMSK modulation system synchronization method Active CN104333525B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410572532.2A CN104333525B (en) 2014-10-23 2014-10-23 A kind of GMSK modulation system synchronization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410572532.2A CN104333525B (en) 2014-10-23 2014-10-23 A kind of GMSK modulation system synchronization method

Publications (2)

Publication Number Publication Date
CN104333525A true CN104333525A (en) 2015-02-04
CN104333525B CN104333525B (en) 2017-09-22

Family

ID=52408179

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410572532.2A Active CN104333525B (en) 2014-10-23 2014-10-23 A kind of GMSK modulation system synchronization method

Country Status (1)

Country Link
CN (1) CN104333525B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104852876A (en) * 2015-05-18 2015-08-19 中国航空无线电电子研究所 Wireless aviation burst communication system
WO2016165080A1 (en) * 2015-04-15 2016-10-20 华为技术有限公司 Reference signal sending and receiving method and device
CN109462567A (en) * 2018-12-25 2019-03-12 成都天奥信息科技有限公司 A method of improving the anti-time slot collision ability of airborne AIS
CN109936520A (en) * 2018-12-20 2019-06-25 杭州电子科技大学 The GMSK of low complex degree demodulates time synchronization method in a kind of burst communication
CN111212006A (en) * 2019-12-25 2020-05-29 奉加微电子(上海)有限公司 Signal demodulation system and demodulation method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6449304B1 (en) * 1999-05-14 2002-09-10 The Aerospace Corporation Method of determining the carrier reference phase of coherent demodulation of frequency hopped continuous phase modulated signals with synch words arbitrarily located within the hop
CN1710898A (en) * 2005-06-29 2005-12-21 西安电子科技大学 Modulation and demodulation method for continuous phase signals
CN104092527A (en) * 2014-06-26 2014-10-08 广州海格通信集团股份有限公司 GMSK coherent demodulation stream processing synchronization method combined with difference

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6449304B1 (en) * 1999-05-14 2002-09-10 The Aerospace Corporation Method of determining the carrier reference phase of coherent demodulation of frequency hopped continuous phase modulated signals with synch words arbitrarily located within the hop
CN1710898A (en) * 2005-06-29 2005-12-21 西安电子科技大学 Modulation and demodulation method for continuous phase signals
CN104092527A (en) * 2014-06-26 2014-10-08 广州海格通信集团股份有限公司 GMSK coherent demodulation stream processing synchronization method combined with difference

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016165080A1 (en) * 2015-04-15 2016-10-20 华为技术有限公司 Reference signal sending and receiving method and device
US10389559B2 (en) 2015-04-15 2019-08-20 Huawei Technologies Co., Ltd. Reference signal sending method, reference signal receiving method, and apparatus
CN104852876A (en) * 2015-05-18 2015-08-19 中国航空无线电电子研究所 Wireless aviation burst communication system
CN104852876B (en) * 2015-05-18 2018-01-16 中国航空无线电电子研究所 A kind of aviation wireless burst communication system
CN109936520A (en) * 2018-12-20 2019-06-25 杭州电子科技大学 The GMSK of low complex degree demodulates time synchronization method in a kind of burst communication
CN109936520B (en) * 2018-12-20 2021-08-17 杭州电子科技大学 Low-complexity GMSK demodulation timing synchronization method in burst communication
CN109462567A (en) * 2018-12-25 2019-03-12 成都天奥信息科技有限公司 A method of improving the anti-time slot collision ability of airborne AIS
CN109462567B (en) * 2018-12-25 2021-07-06 成都天奥信息科技有限公司 Method for improving anti-time slot conflict capability of airborne AIS (automatic identification System)
CN111212006A (en) * 2019-12-25 2020-05-29 奉加微电子(上海)有限公司 Signal demodulation system and demodulation method
CN111212006B (en) * 2019-12-25 2022-07-15 奉加微电子(上海)有限公司 Signal demodulation system and demodulation method

Also Published As

Publication number Publication date
CN104333525B (en) 2017-09-22

Similar Documents

Publication Publication Date Title
CN104333525A (en) GMSK (Gaussian minimum shift keying) modulating system synchronization method
CN104852876B (en) A kind of aviation wireless burst communication system
CN108512791B (en) Satellite-borne AIS demodulation method based on timing frequency offset compensation
US8275077B1 (en) Coherent demodulation of ais-GMSK signals in co-channel
CA2995500C (en) A differential chaos shift keying (dcsk) based on hybrid chaotic system
CN102075472B (en) Method for spreading intermediate frequency of OQPSK and despreading and demodulating spread intermediate frequency
CN103023831B (en) A kind of carrier frequency bias estimation being applicable to burst waveform
CN101868949B (en) Wireless communication method, radio transmitter apparatus and radio receiver apparatus
US20170250792A1 (en) Signal transmission apparatus, carrier phase recovery apparatus and method
CN111935046B (en) Low-complexity frequency shift keying signal symbol rate estimation method
CN110300079B (en) MSK signal coherent demodulation method and system
CN105827551A (en) Method and system for estimating modulation index and frequency offset
US9847895B2 (en) Modulation index shift signaling
CN202906963U (en) A frequency deviation estimating system of a coherent demodulation frequency shift keying modulating signal
CN104363194A (en) PSK (phase shift keying) modulation recognition method based on wave form transformation
CN103905085A (en) Burst hybrid spread spectrum underwater sound covert communication method
US9722833B2 (en) Circuits and methods for frequency offset estimation in FSK communications
CN111147102B (en) Low signal-to-noise ratio code capturing method based on frequency modulation signal
CN106452652B (en) A kind of MPI suppression method based on chaos wireless communication system
CN106508104B (en) A kind of method of extension remote measurement coherent receiver frequency offset estimation range
CN104320243A (en) Analysis processing method of CD-FM-CDSK secret communication system
CN105407067A (en) Recovery method of coherent carrier wave and timing of burst mode GMSK
CN108123724A (en) A kind of communication system based on shortwave narrowband waveform
EP3035622B1 (en) Maximum likelihood sequence detection in the phase domain
CN103873225A (en) Timing estimation method for burst communication

Legal Events

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