CN104267413B - Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search - Google Patents
Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search Download PDFInfo
- Publication number
- CN104267413B CN104267413B CN201410436400.7A CN201410436400A CN104267413B CN 104267413 B CN104267413 B CN 104267413B CN 201410436400 A CN201410436400 A CN 201410436400A CN 104267413 B CN104267413 B CN 104267413B
- Authority
- CN
- China
- Prior art keywords
- signal
- wavelet
- self adaptation
- max
- algorithm
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The present invention relates to a kind of Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search, its technical characteristics is: utilize differential coherence Cumulate algorithm to capture faint received global navigation satellite system signal;Faint received global navigation satellite system signal is analyzed by application Lifting Wavelet decomposition;Dual threshold in Lifting Wavelet denoising is selected to be optimized by application signal intensity self adaptation tabu search algorithm;Utilize optimized dual threshold that the faint received global navigation satellite system signal of capture is carried out post processing.The present invention is reasonable in design, Weak GNSS signal is captured by differential coherence Cumulate algorithm, the Lifting Wavelet dual threshold denoising method of application signal intensity self adaptation TABU search optimization carries out post processing to the Weak GNSS signal obtained, achieve relatively low calculating time cost and higher accuracy, improve the signal to noise ratio of signal output.
Description
Technical field
The invention belongs to Weak Signal Processing technical field, especially one search based on the taboo of signal intensity self adaptation
The Lifting Wavelet dual threshold Denoising Algorithm of rope.
Background technology
In recent years, in order to obtain more accurate navigation information, such as height, speed and position, worldwide navigation is defended
Star system (GNSS) technology is widely used in military field, and vehicle mounted guidance, personal hand-held terminal and
Multiple civil area all plays very important role.But, at ratio under relatively rugged environment, such as township
Village, forest, valley and indoor environment, owing to the reason of propagation loss and multipath fading makes the letter that receives
Number signal to noise ratio the lowest, compared with open environment, signal can produce the decay of 15dB-20dB.It addition,
The power of aeronautical satellite small-signal is far below the scope of application of general operation of receiver power.Therefore, tradition
GNSS signal receiver can not capture gps signal in the environment of high request, not to mention follow the tracks of and fixed
Position.In sum, need to increasingly focus on the capture of Weak GNSS signal, and receiver is at capture amplitude
Small-signal aspect that is relatively low and that be embedded in noise requires the sensitiveest, meanwhile, accumulated time, higher
The lifting of signal to noise ratio and sensitivity is all to need research and the direction paid close attention to.
The method of Weak absorption is different from the detection of general signal, and traditional signal detecting method is typically paid close attention to
Reduction and the raising of signal to noise ratio in noise.But when processing small-signal, these conventional methods often lost efficacy,
Because conventional method can not be effectively improved the signal to noise ratio of the small-signal corroded by noise, and at traditional letter
Number amplify during, small-signal easily torsional deformation.It addition, what some traditional small-signals obtained
Method (such as low frequency phase sensitive filter) needs some to have the specific information of OFF signal, such as phase place,
Frequency etc..Accordingly, it would be desirable to invent some new methods obtaining small-signal and use it for solving practical problems.
Summary of the invention
It is an object of the invention to overcome the deficiencies in the prior art, it is provided that one is reasonable in design, precision is high and the time
The Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search of low cost.
The present invention solves existing technical problem is that and takes techniques below scheme to realize:
A kind of Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search, including following
Step:
Step 1, utilize differential coherence Cumulate algorithm capture faint received global navigation satellite system signal;
Faint received global navigation satellite system signal is analyzed by step 2, application Lifting Wavelet decomposition;
Dual threshold in Lifting Wavelet denoising is selected by step 3, application signal intensity self adaptation tabu search algorithm
It is optimized;
Step 4, utilize optimized dual threshold to capture faint received global navigation satellite system signal carry out after locate
Reason.
And, the implementation method of described step 1 is: obtain small-signal by differential coherence accumulation, then
Accumulation results is inputted Wavelet Denoising Method wave filter to find its code phase and Doppler frequency shift, thus complete capture
Process.
And, the concrete step that processes of described step 2 is:
(1) decomposition step, by original signal sequence SjIt is decomposed into two subdivisions according to parity
evenj-1And oddj-1:
(oddj-1,evenj-1)=Split (Sj)
(2) prediction steps, is predicted by following prediction machine:
dj-1=oddj-1-P(evenj-1)
Wherein, P (evenj-1) it is to predict odd number value, specifying information d by even itemj-1Pass through oddj-1With in advance
Survey P (evenj-1Between), error represents;
(3) update step, update average as the following formula:
Sj-1=evenj-1+ U (dj-1)。
And, the concrete step that processes of described step 3 is:
(1) upper threshold and lower limit T ' are supposedminWith T 'max
(2) the initial interval of threshold optimization is set to [Tmin,Tmax], by upper limit TmaxIt is set to Tmax=3 σ,
Wherein σ is estimated as: σ ≈ Mid/0.6745, in formula, Mid be Decomposition order minimum time wavelet coefficient by
Median after numerical ordering, lower limit TminIt it is the absolute value of minimum coefficient than 0 in jth layer wavelet decomposition;
(3) at T 'minWith T 'maxBetween produce a candidate solution subset, then according to following fitness function
It is ranked up:
Wherein, α is the auto-adaptive parameter that the strength information according to faint received global navigation satellite system signal is chosen, fdFor
Reconstruction signal after denoising, f is reference signal;
Wherein, σωT () is the root-mean-square value of signal,For system average noise power;
(4) scan for by signal intensity self adaptation TS algorithm, be modified current solution obtaining feasible
Solving, if moving to its neighbouring solution to be not better than the most current optimal solution, the most adjacent solution will be accepted;Then draw
Enter to avoid list and avoid circulation;
(5) step (3) is returned until meeting stopping criterion;
(6) utilize equation below calculate optimal threshold:
T′max=y1j×100+y2j×10-1+y3j×10-2+y4j×10-3
T′m in=y5j×100+y6j×10-1+y7j×10-2+y8j×10-3。
And, the concrete processing method of described step 4 is: at optimal double threshold value [T 'min,T′max] be determined after,
Kth wavelet coefficient W in j layer wavelet decompositionj,kProcess as the following formula:
Advantages of the present invention and good effect be:
1, the present invention is during capture small-signal, uses differential coherence accumulation to obtain small-signal, energy
Enough reducing Squared Error Loss, reduce the amplification of noise, the raising to signal to noise ratio is bigger simultaneously.
2, the present invention is by small-signal capture and the combination of Wavelet Denoising Method, signal to noise ratio can be brought up to some
Specific degree, utilizes Wavelet noise-eliminating method to improve output signal-to-noise ratio.
3, the present invention uses tabu search algorithm, it is possible to the optimization feature of simulating human memory function, and passes through
Local neighbor seaching mechanism and corresponding taboo criterion avoid roundabout search, it is possible to discharge some outstanding taboos
Avoid the multiformity that state ensures to have efficient search, can suppress to fall into lacking of easy precocious in genetic algorithm
Fall into, reach final optimization.
4, the present invention is reasonable in design, captures Weak GNSS signal, application letter by differential coherence Cumulate algorithm
The Lifting Wavelet dual threshold denoising method that number Self-adaptive strength TABU search the optimizes Weak GNSS signal to obtaining
Carry out post processing, it is achieved that relatively low calculates time cost and higher accuracy, improve signal output
Signal to noise ratio.
Accompanying drawing explanation
Fig. 1 is the signal intensity self adaptation TABU search process schematic of the present invention;
Fig. 2 is relation comparison diagram between the different wavelet basis wavelet decomposition number of plies and output signal-to-noise ratio;
Fig. 3 is present invention denoising effect comparison diagram compared with primary signal;
Fig. 4 is the output signal-to-noise ratio comparison diagram under the present invention and other three kinds of algorithms before and after noise reduction.
Detailed description of the invention
Below in conjunction with accompanying drawing, the embodiment of the present invention is further described.
A kind of Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search, including following
Step:
Step 1, utilize differential coherence Cumulate algorithm capture faint received global navigation satellite system signal.
In the presence of a harsh environment, such as genuine and valley, prior art proposed different technological means and overcame
A difficult problem for faint GLONASS (GNSS) signal capture.In order to improve the sensitivity of receiver,
And obtaining Weak GNSS signal more accurately, traditional algorithm includes coherent accumulation (COH), non-coherent accumulation
(NCOH) and differential coherence accumulation (DFC).In terms of Weak GNSS signal capture, COH is a kind of straight
Connect and effective method, the signal to noise ratio of the signal detected can be improved.But, the navigation data of 50bps and
Frequency shift (FS) will double to limit the time of coherent accumulation.Use the performance of differential coherence accumulation (DFC) algorithm
The ratio performance mean height 1.5dB of non-coherent accumulation (NCOH), DFC algorithm can obtain the processing gain of maximum.
Therefore, this step utilizes differential coherence accumulation (DFC) to capture faint received global navigation satellite system signal.
During capture small-signal, first implement differential coherence accumulation and obtain small-signal, then will accumulation
Result input Wavelet Denoising Method wave filter, to find its code phase and Doppler frequency shift, has captured subsequently.Differential
Coherent accumulation can reduce Squared Error Loss, and it is conjugated multiplication in adjacent coherent integrator catch cropping.Therefore, this
Method is less to the amplification of noise, and the raising to signal to noise ratio is bigger.
The digital intermediate frequency signal of global navigation satellite system GNSS can be represented as:
Wherein, PiIt is signal amplitude, DiIt is navigation data, CiIt is the C/A when moment k with time delay τ
Spreading code.Having the dependency between the satellite-signal of same pseudo-random noise PRN and local C/A code can
To be expressed as:
Wherein, sIT () is in-phase component part, and sQT () is quadrature component part.
If correlation time is N number of cycle, then correlated results is:
Wherein, A represents that signal amplitude, R (τ) are auto-correlation function, noise n2,IWith n2,QGaussian distributedN () represents Gauss distribution.
When using the accumulation of real differential, the differential cumulative process of a signal is as follows:
If differential cumulative number is M step, then result is:
Wherein Δ τkWith Δ τk-1Represent Delay Estima-tion error during moment k and k-1, Δ ω respectivelyd,kWith Δ ωd,k-1Point
Do not represent the Doppler frequency estimation error of moment k and k-1, Δ ωdRepresent average to moment k from moment k-1
Doppler frequency estimation error.May certify that and have when code phase and Doppler frequency shift are perfectly aligned
Δτk=Δ τk-1=0 and Δ ωd,k=Δ ωd,k-1=0.
Wherein, noise n3,IAnd n3,QGaussian distributedN () represents Gauss distribution, and γ is big
In the constant coefficient of zero, its value depends on the time of integration.
Faint received global navigation satellite system signal is analyzed by step 2, application Lifting Wavelet decomposition.
This step application Lifting Wavelet is decomposed Weak GNSS signal is analyzed including following three step:
(1) decomposition step:
In decomposition, original signal sequence SjIt is broken down into two subdivisions even according to parityj-1With
oddj-1
(oddj-1,evenj-1)=Split (Sj)
The lazy wavelet transformation that this step is referred to as in lifting scheme.Signal decomposition is simply two by this step the most simply
Part, can not promotion signal.Next step promoted will identify that the two sequence is to reduce its dependency.
(2) prediction steps:
Prediction machine is defined as
dj-1=oddj-1-P(evenj-1)
Wherein, P (evenj-1) it is to predict odd number value by even item, and specifying information dj-1Pass through oddj-1With
Prediction P (evenj-1Between), error represents.The dual lifting that this step is referred to as in lifting scheme.Phase when signal
When closing property is the highest, it was predicted that be very effective.
(3) step is updated:
In order to make similarity signal Sj-1Maintain primary signal SjSome characteristics, update average the most as the following formula:
Sj-1=evenj-1+ U (dj-1)
The original lifting that this step is referred to as in lifting scheme.
If three steps of above-mentioned decomposition are all used for decomposed similarity signal Sj-1, then in a number of iteration time
Can obtain primary signal S after numberjMultilamellar decompose.It is found that use lifting scheme to carry out small echo
One of the sharpest edges decomposed are that wavelet transformation can be decomposed into some basic steps and carries out, and theirs is inverse
Conversion is easier to find.
For the setting identical with wavelet filter, lifting process can improve the speed of wavelet transformation doublely.
The computation complexity of wavelet transformation is defined as when ground floor decomposes the multiplication calculated required for a pair coefficient
Number and addition number of times.For lifting scheme, it is assumed that the length of predicted operation coefficient and the length of renewal coefficient of performance
Degree M and N respectively.Can obtain:
Wherein, Md=M/2-1, Nd=N/2-1.
Can draw from above formula, the computation complexity of lifting scheme is 2 (M+N), including (M+N) secondary taking advantage of
The addition that method is secondary with (M+N).It addition, the computation complexity of wavelet transform (DWT) is 4 (M+N)
+2。
Therefore, for longer wave filter, the amount of calculation of application lifting scheme is about the one of traditional algorithm
Half.
Dual threshold in Lifting Wavelet denoising is selected by step 3, application signal intensity self adaptation tabu search algorithm
It is optimized.
The present invention is small-signal capture and Wavelet Denoising Method to be combined, and wherein Wavelet Denoising Method extensively should
For a lot of fields.Signal to noise ratio can be brought up to some specific degree by Wavelet Denoising Method, therefore at satellite
Navigation application also becomes a study hotspot.When Wavelet Denoising Method proposes first for target be reduce
White noise, is then applied to ultrasonic signal and obtains and noise suppressed aspect.As for field of satellite navigation, carry
Go out a kind of method based on wavelet transform (DWT) and be applied to the capture of GNSS signal.Further
Ground, utilizes Wavelet noise-eliminating method to improve output signal-to-noise ratio and has been also applicable in weak GPS signals acquisition detection
Post processing aspect.
Tabu search algorithm is a kind of overall situation Neighborhood-region-search algorithm, and the method simulates the optimization of human memory's function
Feature, and avoid roundabout search by local neighbor seaching mechanism and corresponding taboo criterion.Meanwhile, logical
Crossing the standard of breaking a taboo, it can also discharge the multiformity that some outstanding taboo states ensure to have efficient search.
It addition, TABU search can suppress to fall into the defect of easy precocious in genetic algorithm, reach final
Optimization.
Wavelet Denoising Method mode includes following two mode: hard-threshold mode and soft-threshold mode.The former is by low
All coefficients in bottom threshold are set to zero, and other coefficient all keeps constant;The latter is to will be less than threshold value equally
All coefficients of lower limit are set to zero, but other coefficient shrinks to zero.Scheme, heuristic calculation are determined for threshold value
Method, after the iteration carrying out certain step number, can choose the threshold value of optimum in theory.
Weak GNSS signal comprises high-frequency components and low frequency component.Owing to the ratio of both compositions is not solid
Fixed, if therefore easy choice hard-threshold mode or soft-threshold mode are certain by inevitably losing
The information of frequency range, because they delete high-frequency components the most rambunctiously.Therefore, two kinds above-mentioned
Algorithm all can not obtain gratifying performance during Wavelet Denoising Method.
The present invention apply signal intensity self adaptation Tabu search algorithm carry out Wavelet Denoising Method.Due to tabu search algorithm
Being a kind of heuritic approach, it is prior to the heuritic approach combined by random algorithm and local investigation algorithm.
The signal intensity self adaptation Tabu search algorithm process of the present invention is as follows:
(1) upper threshold and lower limit T ' are supposedminWith T 'maxIt is defined as four position effective digitals, wherein three positions
After arithmetic point, such as [y5,y6,y7,y8], [y1,y2,y3,y4]。
(2) the initial interval of threshold optimization is set to [Tmin,Tmax], upper limit TmaxIt is set to:
Tmax=3 σ
In formula, σ is estimated as σ ≈ Mid/0.6745.
Mid be Decomposition order minimum time wavelet coefficient by the median after numerical ordering.
Lower limit TminIt it is the absolute value of minimum coefficient than 0 in jth layer wavelet decomposition.
To retrain T ' as followsminWith T 'maxInitialize:
Tmin≤T′min< T 'max≤Tmax
(3) at T 'minWith T 'maxBetween produce a candidate solution subset, then according to following fitness function
It is ranked up:
Wherein, α is the auto-adaptive parameter that the strength information according to faint received global navigation satellite system signal is chosen,
It can be according to the value of the actual information self-adaptative adjustment fitness function of the small-signal of capture so that adapt to
Degree function can be made for different small-signals and judging more accurately;fdFor the reconstruction signal after denoising, f
For reference signal, and the value of s is the biggest, represents that fitness is the highest;
Wherein, σωT () is the root-mean-square value of signal,For system average noise power.
(4) according to the process of the TABU search of signal intensity self adaptation shown in Fig. 1, current solution is carried out once letter
Single correction obtains feasible solution, and such process is referred to as a moved further.If moving to its neighbouring solution wait
Choosing solves T* and is not better than current optimal solution, and in order to avoid local optimum, candidate solution T* will be accepted regardless of it
Whether it is optimal solution.Then introduce taboo list and avoid circulation.Taboo list stores all can not be by
It is accepted as the mobile step currently solved.The mobile step meeting taboo rule will be stored in taboo list.Taboo
The use of list reduces the probability of circulation, because it prevent that return in the recent period in certain iterative steps
The solution just accessed.
(5) (3) are returned until meeting stopping criterion.
(6) utilize equation below calculate optimal threshold:
T′max=y1j×100+y2j×10-1+y3j×10-2+y4j×10-3
T′min=y5j×100+y6j×10-1+y7j×10-2+y8j×10-3。
Step 4, utilize optimized dual threshold to capture faint received global navigation satellite system signal carry out after locate
Reason.
Optimal double threshold value [T ' in step 3min,T′max] be determined after, as the following formula by kth in jth layer wavelet decomposition
Individual wavelet coefficient Wj,k:
Achieved by above four steps be applied to Weak GNSS signal capture based on signal intensity self adaptation
The Lifting Wavelet dual threshold Denoising Algorithm that TABU search optimizes, hence it is evident that reduce time complexity and improve
Make an uproar accuracy.
In order to illustrate effect of the present invention, the mode using Computer Simulation below is faint to being applied to
The present invention of GNSS signal capture is modeled, and achieves the simulation to real scene by assignment.Numeral
Intermediate-freuqncy signal is as the input signal of emulation, and detailed process point following four aspect is analyzed:
(1) wavelet basis and the selection of Decomposition order
As shown in Figure 2, when decomposition level rises to 3, signal to noise ratio keeps stable and even performs better than.
Haar wavelet transform base is put up the best performance in all candidate's wavelet basiss, and owing to signal is square wave, has minor fluctuations
Wavelet basis more applicable, therefore select Haar wavelet transform base as in inventive algorithm application wavelet basis.
It addition, also contains the discussion selected for Decomposition order in fig. 2.If Decomposition order is too much,
Owing to calculating time shortage can make efficiency reduce.And, the risk losing effective information can cause one relatively
The signal to noise ratio result of difference.But, when Decomposition order deficiency, optimum can not be obtained equally.Therefore,
The present invention selects three layers of decomposition.
(2) denoising effect
The line being located below in Fig. 3 represents primary signal, and another represents and uses based on signal intensity self adaptation
Lifting Wavelet dual threshold Denoising Algorithm (present invention) that TABU search optimizes carries out the signal after denoising.Permissible
It is evident that along with the raising of input signal power, not only the output signal-to-noise ratio of primary signal improves, denoising
The output signal-to-noise ratio of signal improves the most therewith.And, the difference between them is held essentially constant, this meaning
The system of wearing can process this input signal power small-signal between-176dBw to-168dBw, at letter
Make an uproar than the gain that above can provide 8dB.
(3) present invention and the contrast of other algorithm
Fig. 4 shows the comparison diagram of output signal-to-noise ratio in the case of four kinds.Line with round dot represents this
The denoising effect of the algorithm of bright middle proposition, the Lifting Wavelet i.e. optimized based on signal intensity self adaptation TABU search
Dual threshold Denoising Algorithm.It is obvious that fact proved, it is possible to the method obtaining optimal output signal-to-noise ratio is this
The algorithm of bright middle proposition, and be used alone Lifting Wavelet carry out the effect of denoising be better than be used alone tradition from
Dissipate wavelet transformation (DWT).Meanwhile, along with the increase of input signal power, it is used alone DWT and individually makes
The snr gain being obtained in that with Lifting Wavelet denoising is gradually lowered, on the contrary, adaptive based on signal intensity
The performance answering the Lifting Wavelet dual threshold Denoising Algorithm of TABU search optimization is stable and gratifying.
(4) application wavelet decomposition lifting scheme and the contrast on time complexity of the DWT scheme
In the following table, compared for the calculating time of two kinds of denoising methods, i.e. avoiding based on signal intensity self adaptation
The DWT small echo dual threshold Denoising Algorithm of chess game optimization and based on signal intensity self adaptation TABU search optimization carry
Rise small echo dual threshold Denoising Algorithm.Should be noted that, the specific calculating time can not represent real reality
Time, but this contrast can be regarded as computation complexity and the reflection of enforcement difficulty.Therefore, it can
Reaching a conclusion, DWT needs to consume the time far above application wavelet decomposition lifting scheme.
It is emphasized that embodiment of the present invention is illustrative rather than determinate, therefore originally
Invention includes the embodiment that is not limited to described in detailed description of the invention, every by those skilled in the art according to
Other embodiments that technical scheme draws, also belong to the scope of protection of the invention.
Claims (5)
1. a Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search, its feature
It is to comprise the following steps:
Step 1, utilize differential coherence Cumulate algorithm capture faint received global navigation satellite system signal;
Faint received global navigation satellite system signal is analyzed by step 2, application Lifting Wavelet decomposition;
Dual threshold in Lifting Wavelet denoising is selected by step 3, application signal intensity self adaptation tabu search algorithm
It is optimized;
Step 4, utilize optimized dual threshold to capture faint received global navigation satellite system signal carry out after locate
Reason.
Lifting Wavelet dual threshold based on signal intensity self adaptation TABU search the most according to claim 1
Denoising Algorithm, it is characterised in that: the implementation method of described step 1 is: obtain micro-by differential coherence accumulation
Weak signal, then inputs Wavelet Denoising Method wave filter to find its code phase and Doppler frequency shift by accumulation results,
Thus complete acquisition procedure.
Lifting Wavelet dual threshold based on signal intensity self adaptation TABU search the most according to claim 1
Denoising Algorithm, it is characterised in that: the concrete step that processes of described step 2 is:
(1) decomposition step, by original signal sequence SjIt is decomposed into two subdivisions according to parity
evenj-1And oddj-1:
(oddj-1,evenj-1)=Split (Sj)
(2) prediction steps, is predicted by following prediction machine:
dj-1=oddj-1-P(evenj-1)
Wherein, P (evenj-1) it is to predict odd number value, specifying information d by even itemj-1Pass through oddj-1With in advance
Survey P (evenj-1Between), error represents;
(3) update step, update average as the following formula:
Sj-1=evenj-1+U(dj-1)。
Lifting Wavelet dual threshold based on signal intensity self adaptation TABU search the most according to claim 1
Denoising Algorithm, it is characterised in that: the concrete step that processes of described step 3 is:
(1) upper threshold and lower limit T are supposedm'inAnd Tm'ax
(2) the initial interval of threshold optimization is set to [Tmin,Tmax], by upper limit TmaxIt is set to Tmax=3 σ,
Wherein σ is estimated as: σ ≈ Mid/0.6745, in formula, Mid be Decomposition order minimum time wavelet coefficient by
Median after numerical ordering, lower limit TminIt it is the absolute value of minimum coefficient than 0 in jth layer wavelet decomposition;
(3) at T 'minWith T 'maxBetween produce a candidate solution subset, then according to following fitness function
It is ranked up:
Wherein, α is the auto-adaptive parameter that the strength information according to faint received global navigation satellite system signal is chosen, fdFor
Reconstruction signal after denoising, f is reference signal;
Wherein, σωT () is the root-mean-square value of signal,For system average noise power;
(4) scan for by signal intensity self adaptation TS algorithm, be modified current solution obtaining feasible
Solving, if moving to its neighbouring solution to be not better than the most current optimal solution, the most adjacent solution will be accepted;Then draw
Enter to avoid list and avoid circulation;
(5) step (3) is returned until meeting stopping criterion;
(6) utilize equation below calculate optimal threshold:
T′max=y1j×100+y2j×10-1+y3j×10-2+y4j×10-3
T′min=y5j×100+y6j×10-1+y7j×10-2+y8j×10-3。
Lifting Wavelet dual threshold based on signal intensity self adaptation TABU search the most according to claim 1
Denoising Algorithm, it is characterised in that: the concrete processing method of described step 4 is: at optimal double threshold value [T 'min,T′max]
After being determined, kth wavelet coefficient W in jth layer wavelet decompositionj,kProcess as the following formula:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410436400.7A CN104267413B (en) | 2014-08-29 | 2014-08-29 | Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410436400.7A CN104267413B (en) | 2014-08-29 | 2014-08-29 | Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104267413A CN104267413A (en) | 2015-01-07 |
CN104267413B true CN104267413B (en) | 2016-10-19 |
Family
ID=52158953
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410436400.7A Expired - Fee Related CN104267413B (en) | 2014-08-29 | 2014-08-29 | Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104267413B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108703741B (en) * | 2018-03-30 | 2019-04-09 | 西安交通大学 | A kind of shifty NIRS Interference Detection and minimizing technology |
CN113391302B (en) * | 2020-03-13 | 2023-03-28 | 成都锦江电子系统工程有限公司 | Dynamic double-threshold echo detection method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101950024A (en) * | 2010-08-11 | 2011-01-19 | 中国电子科技集团公司第二十研究所 | Code carrier consistency detection method applied to local area augmentation system |
CN102594472A (en) * | 2012-03-22 | 2012-07-18 | 北京邮电大学 | Method and system for wireless channel measurement based on wavelet decomposition threshold de-nosing |
CN103778921A (en) * | 2013-12-28 | 2014-05-07 | 中国人民解放军兰州军区乌鲁木齐总医院 | Method for eliminating nonuniform noise in speech collected by radar |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3334226B1 (en) * | 2009-02-05 | 2020-03-25 | Apple Inc. | Method and system for user equipment location determination on a wireless transmission system |
-
2014
- 2014-08-29 CN CN201410436400.7A patent/CN104267413B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101950024A (en) * | 2010-08-11 | 2011-01-19 | 中国电子科技集团公司第二十研究所 | Code carrier consistency detection method applied to local area augmentation system |
CN102594472A (en) * | 2012-03-22 | 2012-07-18 | 北京邮电大学 | Method and system for wireless channel measurement based on wavelet decomposition threshold de-nosing |
CN103778921A (en) * | 2013-12-28 | 2014-05-07 | 中国人民解放军兰州军区乌鲁木齐总医院 | Method for eliminating nonuniform noise in speech collected by radar |
Also Published As
Publication number | Publication date |
---|---|
CN104267413A (en) | 2015-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104267379B (en) | A kind of active radar and passive radar based on Waveform Design works in coordination with anti-interference method | |
CN102638888B (en) | Indoor positioning fingerprint grouping method based on signal statistics | |
CN110348288A (en) | A kind of gesture identification method based on 77GHz MMW RADAR SIGNAL USING | |
CN104597435B (en) | Correction frequency domain compensation and fractional order Fourier transformation based multi-frame coherent TBD method | |
CN106125053A (en) | Pulse Doppler radar polarization anti jamming method | |
CN102879766A (en) | Method and apparatus for detecting and tracking faint target of high frequency ground wave radar | |
CN1831558A (en) | Single-channel synthetic aperture radar moving-target detection method based on multi-apparent subimage paire | |
CN104597434B (en) | Improve the multiframe coherent TBD methods of envelope shift compensation and Fourier Transform of Fractional Order | |
CN101893698B (en) | Noise source test and analysis method and device | |
CN107861123A (en) | A kind of through-wall radar is under complex environment to the method for multiple mobile object real-time tracking | |
CN107153208A (en) | A kind of method of gps carrier phase cycle-slip detection and repair | |
CN107607937A (en) | Radar target distance-finding method based on time reversal | |
CN102353952A (en) | Line spectrum detection method by coherent accumulation of frequency domains | |
CN102323601A (en) | Method for capturing GNSS-BOS (Global Navigation Satellite System-Binary Offset Carrier) modulation signal | |
CN104267413B (en) | Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search | |
CN102565825B (en) | Received signal strength fiduciary level decision maker, method and code phase error calculation method | |
CN113406671A (en) | GNSS forwarding type deception jamming detection method based on C/N0-MV | |
CN103926605B (en) | GPS weak signal catching method based on difference circulation coherent integration | |
Zhu et al. | Iterative time-reversal mirror method for imaging the buried object beneath rough ground surface | |
CN103809193B (en) | Improve the capture systems of Weak GNSS signal processing gain | |
Closas et al. | Evaluation of GNSS direct position estimation in realistic multipath channels | |
CN116930963A (en) | Through-wall imaging method based on wireless communication system | |
CN102903083A (en) | Context model and dual-tree complex wavelet transform-based denoising method for underwater sonar image | |
CN105652256A (en) | High-frequency surface wave radar TBD method based on polarization information | |
CN115311532A (en) | Ground penetrating radar underground cavity target automatic identification method based on ResNet network model |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161019 Termination date: 20170829 |