CN101436877B - Method for capturing multi-path interference resistant PN code self-adapting threshold - Google Patents

Method for capturing multi-path interference resistant PN code self-adapting threshold Download PDF

Info

Publication number
CN101436877B
CN101436877B CN2008102333579A CN200810233357A CN101436877B CN 101436877 B CN101436877 B CN 101436877B CN 2008102333579 A CN2008102333579 A CN 2008102333579A CN 200810233357 A CN200810233357 A CN 200810233357A CN 101436877 B CN101436877 B CN 101436877B
Authority
CN
China
Prior art keywords
filter
sign indicating
indicating number
threshold
code
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.)
Expired - Fee Related
Application number
CN2008102333579A
Other languages
Chinese (zh)
Other versions
CN101436877A (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN2008102333579A priority Critical patent/CN101436877B/en
Publication of CN101436877A publication Critical patent/CN101436877A/en
Application granted granted Critical
Publication of CN101436877B publication Critical patent/CN101436877B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a method for capturing a self-adapting threshold of an anti-multichannel interference PN code in the technical field of communication, which comprises the following steps: firstly constructing a two-way parallel digital matched filter, achieving correlative operation of a spread-spectrum signal and a local PN code, then summing two-way correlative values after being squared, and obtaining synchronously verified decision statistic of the PN code after taking a square root. The method is characterized in that the statistic is added to a sequence converting circuit and an input end of a median filter simultaneously, the statistic output by the sequence converting circuit is processed by a two-way parallel form filter and then output values are summed, a mean value is taken to obtain a self-adapting judgment threshold, a judging device compares the judgment statistic with an input self-adapting judgment threshold till the judgment amount is larger than a threshold value, and capture of the PN code is decided to be successful at this moment. The method has the advantages that the method can adapt to a multichannel environment better, has excellent performance of anti-multichannel interference, has strong practicability, and can be applied to various mobile communication systems using direct sequence spread-spectrum technique.

Description

The PN code self-adapting threshold catching method that a kind of anti-multipath disturbs
Technical field:
The present invention relates to communication technical field, relate in particular under the multipath channel environment, realize that in direct sequence spread spectrum communication system (DSSS) the local code sequence is slightly synchronous with the sign indicating number sequence that receives, a kind of PN code capture method that anti-multipath disturbs.
Background technology:
The PN sign indicating number be the most basic most important problem in the direct sequence spread spectrum communication system (DSSS) synchronously.Synchronizing process generally comprises two parts: catch and follow the tracks of.Catching is coarse adjustment, can make the local code sequence slightly synchronous with the sign indicating number sequence that receives, and general phase error is in 1/2 chip lengths.Tracking is fine tuning, and sending and receiving end symbol phases error is further dwindled, and reaches synchronous accurately and in whole transmission course, keeps this synchronous.
Forefathers have carried out a large amount of research for catching of PN sign indicating number, but mainly concentrate on the catching method based on fixed threshold.In the GSM of reality; Because user's fast moving adds the influence of channel fading and dynamic disturbance, the power that causes receiving signal can change in great dynamic range; Traditional fixed threshold can not obtain good acquisition performance, can't adapt to actual applied environment.For this reason, people have proposed multiple adaptive threshold catching method.Document [1] [Kim, C.J., Lee; H.J., and Lee, H.S.; Adaptive acquisition of PN sequences for DSSS communications, IEEE Trans.Commun., 1998; 46 (8): 993-996] invariable false alerting in the Radar Signal Processing (CFAR) inspection technology is incorporated in the PN sign indicating number acquisition procedure, has proposed average statistics (MLAP) and two kinds of methods of sequencing statistical (OSAP), utilize a plurality of forward directions and back Background Noise Power to be estimated to observed quantity.Document [2] [Xue Wei, Xiang Jingcheng, Zhou Zhizhong; The threshold adaptive method of estimation that a kind of PN sign indicating number is caught; Electronic letters, vol, 2003,31 (12): 1870-1873] and MLAP, OSAP method in the document [1] all be to utilize the maximum likelihood (ML) of the sample equivalent Gaussian noise power of judgment variables to estimate.Improve estimated accuracy and must increase and handle window length, but under multi-path environment, the relevant peaks of multi-path signal-component possibly fall into handles window, cause the noise power estimated value higherly, detection probability significantly descends.Document [3] [Choi Kwonhue; Cheun Kyungwhoon; And Jung Taejin, Adaptive PN code acquisition using instantaneous power-scaled detection threshold under rayleigh fading and pulsed gaussian noise jamming, IEEE Trans.on Commun.; 2002; 50 (8): 1232-1235] proposed a kind of adaptive threshold scheme, exist under the situation of narrow band interference, thereby estimated value also can higherly cause decreased performance based on real-time Energy Estimation.
Summary of the invention:
The objective of the invention is to: provide a kind of and can adapt to multi-path environment, background noise is estimated in real time, the PN code self-adapting threshold catching method that anti-multipath disturbs.
The objective of the invention is to realize through the enforcement following technical proposals:
1, a kind of PN code self-adapting threshold catching method of anti-multipath interference, realize by following step:
A, at first make up the parallel digital matched filter F of two-way I, F Q, F is set according to the value of corresponding position in the local PN sign indicating number IAnd F QIn the weight coefficient of multiplier: the corresponding position of PN sign indicating number is 0, then weight coefficient get+1; The corresponding position of PN sign indicating number is 1, and then weight coefficient gets-1.
B, with the in-phase component r of the spread-spectrum signal that receives IWith quadrature component r QBe input to digital matched filter F respectively IAnd F QIn, the related operation between realization and the local PN sign indicating number.
C, to digital matched filter F I, F QHomophase and the quadrature two-way correlation of output carries out square asking mould, then the root of making even after two modulus additions obtained PN sign indicating number verification of synchronization decision statistic amount Z n
It is characterized in that: d, the verification of synchronization decision statistic amount Z that is used to adjudicate nBe added to the input of sequence change-over circuit and median filter simultaneously, the sequence change-over circuit is exported X according to median filter nValue, to the Z of each input nCarry out like the conversion that acts, with list entries Z nConversion is output as a new sequence Y nIf: Z nValue greater than 5 times of X n, then the sequence change-over circuit is exported Y n=X nOtherwise Y n=Z n
Median filter adopts shift register structure, and filter window length is 2N+1, and N is a positive integer here, and being located at i, to import the sampling point of sample sequence in filter window constantly be x I-N... X i..., x I+N, this moment, filter was exported y so i=med (x I-N... X i..., x I+N), sample values all in the window of med () expression are here got the computing that it comes a middle value output by after the sequence arrangement from small to large.
E, the parallel morphological filter of structure two-way, the one tunnel is form open-close filter F Oc, another road is that form is closed-Kai filter F Co, with sequence Y nBe added to F OcAnd F CoInput, it is carried out further smothing filtering.
The related definition of morphological filter is following:
If being respectively, f (n) and g (m) be defined in F={0,1 ...., N-1} and G={0,1 ...., the discrete function that M-1} is last, and N>>M.If f (n) is a list entries, g (m) is a structural element.F (n) is defined as respectively about corrosion and the expansion of g (m):
Figure GSB00000657068500021
(n=0,1,......,N-M)
( f ⊕ g ) ( n ) = max m = 0,1 , . . . , M - 1 { f ( n - m ) + g ( m ) } - - - ( 2 )
(n=0,1,......,N+M-2)
Where the symbols
Figure GSB00000657068500023
and
Figure GSB00000657068500024
denote erosion and dilation.
F (n) opens to close respectively with form about the form of g (m) and is defined as:
Figure GSB00000657068500025
Symbol о and represent that respectively form opens and the form closed operation in the formula.
Oc (f (n)), co (f (n)) represent that respectively form open-close filter and form close-output of Kai filter, and form open-close filter and form close-the Kai filter adopts the cascade combining form of modal of morphological opening and closing operation, definition as follows:
oc(f(n))=(fоg·g)(n)(5)
co(f(n))=(f·gоg)(n)(6)
F, with morphological filter F OcAnd F CoThe addition of two-way output valve, get its mean value then and obtain self-adaptive decision threshold V T
Because the anti-dilatancy of opening operation and the dilatancy of closed operation, there is the statistical bias phenomenon in the output of form OC and CO filter, uses separately somely can't obtain good effect.Therefore the present invention to open-close with close-output valve of Kai filter made average treatment.
G, the PN sign indicating number verification of synchronization decision statistic amount Z of decision device to importing nSelf-adaptive decision threshold V with input TRelatively big or small, as verification of synchronization decision statistic amount Z nBe lower than threshold value V TThe time, show local PN code phase and receive signal not synchronously that the decision statistic amount and the thresholding that then produce corresponding next phase place are made comparisons again, up to the corresponding verification of synchronization decision statistic amount Z of one of them phase place nSurpass thresholding V T, be judged to be acquisition success this moment, decision device output PN sign indicating number acquisition success index signal.
The invention has the advantages that: owing to used the estimator that is made up of sequence change-over circuit, median filter and morphological filter that background noise has been carried out real-time estimation, so the present invention can carry out the self adaptation adjustment according to the noise size by the detection threshold that is used for prize judgment that estimator obtains when catching carrying out the PN sign indicating number.Simultaneously, owing to carried out the sequence conversion to received signal, rejected the exceptional value in the sample; Eliminated the influence that possibly fall into the relevant peaks of handling the multi-path signal-component in the window, reduced the influence that multipath disturbs greatly, accomplished preliminary treatment data; Therefore; Adaptive threshold catching method proposed by the invention can be good at adapting to the multipath channel environment, and the anti-multipath jamming performance is excellent, can be applied to use the various GSMs of direct sequence spread spectrum skill.
Description of drawings
Fig. 1 catches the implementation method sketch map for PN sign indicating number of the present invention
Fig. 2 realizes receiving spread frequency signal and local PN sign indicating number for digital matched filter and carries out the sketch map of related operation
Fig. 3 is the false alarm probability P of system FaBe respectively 0.01,0.001 and at 0.0001 o'clock, the graph of a relation of PN sign indicating number Acquisition Detection probability of the present invention and signal to noise ratio
Mark among the figure: r IBe the in-phase component of the spread-spectrum signal that receives, r QBe the quadrature component of the spread-spectrum signal that receives, F IAnd F QBe digital matched filter, () 2Be square operation, symbol
Figure GSB00000657068500031
Be add operation, symbol
Figure GSB00000657068500032
Be extracting operation, F Oc, F CoRepresent the form open-close respectively and close-the Kai filter, symbol 1/2 expression to input value divided by 2 computing, Expression is to Z nAnd V TCompare the computing of size, symbol The computing to a plurality of input variable summations is represented in expression multiplying, ∑.
Embodiment
In an embodiment of the present invention; Receiving signal is the Direct swquence spread spectrum signal; The modulation system of spread spectrum system adopts difference binary phase shift keying (DBPSK; Differential Binary Phase Shift Keying), the code length of spread spectrum PN sign indicating number is L=1023 (a m sequence), and spreading gain is about 30dB.Data rate is R b=9.6kbit/s confirms that thus bit rate is R c=9.6L=9820.8k chip/s, the about 20MH of required channel width ZAdopt data bit and the synchronous mode of PN sign indicating number during spread spectrum, promptly the PN sign indicating number aligns with the transfer point of data bit, and each data bit width equals a PN sign indicating number cycle.
In an embodiment, digital matched filter F I, F QThe delay tap number be 1023, equal the code length of the PN sign indicating number that adopts, promptly digital matched filter will carry out related operation complete cycle to receiving the PN sign indicating number.Median filter sliding window length gets 9, the form open-close with close-the used structural element of Kai filter get g (m)=0,1,2,1,0}.So far, accomplished parameter setting to PN sign indicating number capture systems self-adaptive decision threshold.
The spread-spectrum signal that receives is divided into homophase (R through after the non-coherent demodulation I) and quadrature (r Q) two paths of signals, be input to digital matched filter F respectively IAnd F QIn, the related operation between realization and the local PN sign indicating number.
In order to overcome the influence of the phase deviation that the incoherent reception mode brings, to digital matched filter F I, F QHomophase and the quadrature two-way correlation of output carries out square asking mould, then the root of making even after two modulus additions obtained PN sign indicating number verification of synchronization decision statistic amount Z n
The verification of synchronization decision statistic amount Z that is used to adjudicate nBe added to the input of sequence change-over circuit and median filter simultaneously, the sequence change-over circuit is exported X according to median filter nValue, to the Z of each input nCarry out like the conversion that acts, with list entries Z nConversion is output as a new sequence Y nIf: Z nValue greater than 5 times of X n, then the sequence change-over circuit is exported Y n=X nOtherwise Y n=Z n
With sequence Y nBe added to parallel form open-close filter F OcClose with form-Kai filter F CoInput, it is carried out further smothing filtering.With morphological filter F OcAnd F CoThe addition of two-way output valve, get its mean value then and obtain decision device test statistics is adjudicated required self-adaptive decision threshold V T
Decision device is with the PN sign indicating number verification of synchronization decision statistic amount Z of input nSelf-adaptive decision threshold V with input TCompare, as verification of synchronization decision statistic amount Z nBe lower than threshold value V TThe time, show local PN code phase and receive signal not synchronously that the decision statistic amount and the thresholding that then produce corresponding next phase place are made comparisons again, up to the corresponding verification of synchronization decision statistic amount Z of one of them phase place nSurpass thresholding V T, be judged to be acquisition success this moment, decision device output PN sign indicating number acquisition success index signal.
In an embodiment, the PN sign indicating number is caught all and is realized by hardware.In practical application, the designer can select corresponding software realization mode as the case may be.
Fig. 3 has provided the false alarm probability P of system FaBe respectively 0.01,0.001 and at 0.0001 o'clock, the PN sign indicating number Acquisition Detection probability of present embodiment and the relation of signal to noise ratio, along with signal to noise ratio increases, detection probability progressively increases, and has shown validity of the present invention.

Claims (1)

1. the PN code self-adapting threshold catching method that disturbs of an anti-multipath, realize by following step:
A, at first make up the parallel digital matched filter F of two-way I, F Q, F is set according to the value of corresponding position in the local PN sign indicating number IAnd F QIn the weight coefficient of multiplier: the corresponding position of PN sign indicating number is 0, then weight coefficient get+1; The corresponding position of PN sign indicating number is 1, and then weight coefficient gets-1;
B, with the in-phase component r of the spread-spectrum signal that receives IWith quadrature component r QBe input to digital matched filter F respectively IAnd F QIn, the related operation between realization and the local PN sign indicating number;
C, to digital matched filter F I, F QHomophase and the quadrature two-way correlation of output carries out square asking mould, then the root of making even after two modulus additions obtained PN sign indicating number verification of synchronization decision statistic amount Z n
It is characterized in that: d, the verification of synchronization decision statistic amount Z that is used to adjudicate nBe added to the input of sequence change-over circuit and median filter simultaneously, the sequence change-over circuit is exported X according to median filter nValue, to the Z of each input nCarry out like the conversion that acts, with list entries Z nConversion is output as a new sequence Y nIf: Z nValue greater than 5 times of X n, then the sequence change-over circuit is exported Y n=X nOtherwise Y n=Z n
E, the parallel morphological filter of structure two-way, the one tunnel is form open-close filter F Oc, another road is that form is closed one and opened filter F Co, with sequence Y nBe added to F OcAnd F CoInput, it is carried out further smothing filtering;
F, to morphological filter F OcAnd F CoThe addition of two-way output valve, getting its mean value then is self-adaptive decision threshold V T
G, the PN sign indicating number verification of synchronization decision statistic amount Z of decision device to importing nSelf-adaptive decision threshold V with input TRelatively big or small, as verification of synchronization decision statistic amount Z nBe lower than threshold value V TThe time, show local PN code phase and receive signal not synchronously that the decision statistic amount and the thresholding that then produce corresponding next phase place are made comparisons again, up to the corresponding verification of synchronization decision statistic amount Z of one of them phase place nSurpass thresholding V T, be judged to be acquisition success this moment, decision device output PN sign indicating number acquisition success index signal.
CN2008102333579A 2008-12-16 2008-12-16 Method for capturing multi-path interference resistant PN code self-adapting threshold Expired - Fee Related CN101436877B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008102333579A CN101436877B (en) 2008-12-16 2008-12-16 Method for capturing multi-path interference resistant PN code self-adapting threshold

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008102333579A CN101436877B (en) 2008-12-16 2008-12-16 Method for capturing multi-path interference resistant PN code self-adapting threshold

Publications (2)

Publication Number Publication Date
CN101436877A CN101436877A (en) 2009-05-20
CN101436877B true CN101436877B (en) 2012-05-23

Family

ID=40711125

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008102333579A Expired - Fee Related CN101436877B (en) 2008-12-16 2008-12-16 Method for capturing multi-path interference resistant PN code self-adapting threshold

Country Status (1)

Country Link
CN (1) CN101436877B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101231983B1 (en) * 2011-06-28 2013-02-08 엘에스산전 주식회사 A method of obtaining synchronization in a real time locating system
CN102801441B (en) * 2012-07-27 2014-08-27 北京邮电大学 Synchronous clock extraction method for receiving end in fully digital spread spectrum communication system
CN106301416B (en) * 2016-08-26 2019-05-28 成都九洲迪飞科技有限责任公司 Bpsk signal and msk signal hybrid spread spectrum system
CN106357300B (en) * 2016-08-26 2019-04-16 成都九洲迪飞科技有限责任公司 The anti-impulse disturbances method of BPSK spread spectrum system
CN109633704B (en) * 2018-12-28 2021-04-27 四川安迪科技实业有限公司 Maximum value-based satellite communication capturing method and system
CN110535620B (en) * 2019-08-22 2022-03-08 重庆大学 Signal detection and synchronization method based on decision feedback
CN114759949A (en) * 2022-03-22 2022-07-15 陕西烽火电子股份有限公司 Fast code capture method in BPSK modulation direct sequence spread spectrum communication

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5724384A (en) * 1994-07-14 1998-03-03 Samsung Electronics Co., Ltd. PN code sync device using an adaptive threshold
CN1168249C (en) * 2000-09-29 2004-09-22 华为技术有限公司 Self-adaptive threshold pseudo-random code catching and identifying method
CN1622482A (en) * 2003-11-27 2005-06-01 中国电子科技集团公司第三十研究所 A PN code capture method having decision threshold adaptive estimation function
CN101282139A (en) * 2008-03-11 2008-10-08 电子科技大学 Real time pseudo-code period synchronous technique for non-cooperation low signal-noise ratio direct-sequence spread-spectrum signal

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5724384A (en) * 1994-07-14 1998-03-03 Samsung Electronics Co., Ltd. PN code sync device using an adaptive threshold
CN1168249C (en) * 2000-09-29 2004-09-22 华为技术有限公司 Self-adaptive threshold pseudo-random code catching and identifying method
CN1622482A (en) * 2003-11-27 2005-06-01 中国电子科技集团公司第三十研究所 A PN code capture method having decision threshold adaptive estimation function
CN101282139A (en) * 2008-03-11 2008-10-08 电子科技大学 Real time pseudo-code period synchronous technique for non-cooperation low signal-noise ratio direct-sequence spread-spectrum signal

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
何昭君,唐友喜.一种自适应门限的PN码捕获方案.《2007中国西部青年通信学术会议论文集》.2007,574-579. *
张春海,王世练,张尔扬.基于噪声功率估计的抗窄带干扰接收机自适应PN码捕获技术.《信号处理》.2007,第23卷(第2期),165-168. *

Also Published As

Publication number Publication date
CN101436877A (en) 2009-05-20

Similar Documents

Publication Publication Date Title
CN101436877B (en) Method for capturing multi-path interference resistant PN code self-adapting threshold
Tian et al. A GLRT approach to data-aided timing acquisition in UWB radios-Part I: Algorithms
CN100586115C (en) Modulation and demodulation method for continuous phase signals
US7796686B2 (en) Adaptive ultrawideband receiver and method of use
CN109167746A (en) Continuous wave and the quick identification device of pulse signal
Rabbachin et al. ML time-of-arrival estimation based on low complexity UWB energy detection
CN105871765A (en) Wireless communication carrier wave tracking method based on FFT assistant S-PLL
CN109088838B (en) Pseudo code-Doppler fast capturing method of direct sequence spread spectrum DPSK signal under high dynamic condition
CN100550662C (en) A kind of signal acceptance method that is used for pulse ultra-broadband communication system
CN105790788B (en) A kind of pseudo-code of direct expansion msk signal-Doppler's joint acquisition method
CN101951356A (en) Synchronization method of orthogonal frequency division multiplexing-ultra wideband (OFDM-UWB) system based on peak detection
CN109039573A (en) A method of based on CPM signal multiple-symbol detection
CN110224720A (en) A kind of CSS system non-coherent demodulation method based on multiplication cross second differnce structure
CN105445767B (en) BOC signal parameter blind estimating methods based on average Ambiguity Function
CN107370705A (en) FFT optimization method in the capture of high dynamic weakly continuous phase modulated signal
CN112702080B (en) Direct sequence spread spectrum signal and pseudo code estimation method based on K-means algorithm
CN105656511B (en) Differential correlation acquisition method suitable for environment with frequency offset and low signal-to-noise ratio
CN100420163C (en) Frequency domain receiving method for super broad band radio communication system
CN103501187B (en) Interference cancellation-based short wave multi-path signal synchronization method
CN101710889B (en) Blind symbol rate estimating method for multi-path linear hybrid digital modulation signals
CN108718461B (en) Frequency offset resistant burst CPM signal frame synchronization method
CN109104215A (en) A kind of Frequency Hopping Signal chip rate blind estimating method based on wavelet transformation
CN108401582B (en) A kind of GEO satellite mobile communication system initial frequency synchronization method adapting to TDMA/FDMA systems
CN100433563C (en) Correlated receiver and signal receiving method for super-broad band radio communication system
CN109361632B (en) Parameter blind estimation method of multipath BOC signal based on spectrum correlation

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120523

Termination date: 20141216

EXPY Termination of patent right or utility model