CN106814601B - A kind of Noise Filter method of signal specific - Google Patents

A kind of Noise Filter method of signal specific Download PDF

Info

Publication number
CN106814601B
CN106814601B CN201710182251.XA CN201710182251A CN106814601B CN 106814601 B CN106814601 B CN 106814601B CN 201710182251 A CN201710182251 A CN 201710182251A CN 106814601 B CN106814601 B CN 106814601B
Authority
CN
China
Prior art keywords
observation
kind
filter method
value
signal specific
Prior art date
Application number
CN201710182251.XA
Other languages
Chinese (zh)
Other versions
CN106814601A (en
Inventor
傅德锦
刘胜
方军涛
李俊
王毅
Original Assignee
厦门九华通信设备厂
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 厦门九华通信设备厂 filed Critical 厦门九华通信设备厂
Priority to CN201710182251.XA priority Critical patent/CN106814601B/en
Publication of CN106814601A publication Critical patent/CN106814601A/en
Application granted granted Critical
Publication of CN106814601B publication Critical patent/CN106814601B/en

Links

Classifications

    • GPHYSICS
    • G04HOROLOGY
    • G04RRADIO-CONTROLLED TIME-PIECES
    • G04R20/00Setting the time according to the time information carried or implied by the radio signal
    • G04R20/02Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS

Abstract

A kind of Noise Filter method of signal specific, according to observation Am, it is assumed that observation AmIt is exactly true value Zm, the guess value T with certain precision can be released by mapping relationship f (m-n)m‑n=Am+ f (m-n), because of observed quantity AmIt has been superimposed random noise, has often led to guess value Tm‑nThere is biggish deviations, or because speculate the increasing deviation of the accumulation that points increase.The basic principle of filtering algorithm of the invention assumes observation AmEqual to Zm, then according to observation AmObtained a series of guess value and observed difference score value should be equal to random noise rm, amplitude 80% falls between [- R R], whether exceeds [- R R] to T according to difference valuemIncrement or decreasing amount adjustment are carried out, increases T with adjustment numbermIt can become closer to Zm

Description

A kind of Noise Filter method of signal specific

Technical field

The present invention relates to field of signal processing, the Noise Filter method of especially a kind of signal specific.

Background technique

In time service field, either Beidou time service, GPS time service, GLONASS time service, PTP time service, long-wave time service or short Wave time service carries out temporal frequency transmission using wireless or wire channel, because receiver clock is unstable, system clock deviation, passes The time service information of the factors cause receivers such as the delay jitter of defeated channel reduction often has certain randomized jitter.Passive fixed Position field, the arrival time of the time pulse in time domain received often because signal strength or weakness variation and there are biggish waves It is dynamic.

Traditional linear filter method can effectively filter out random noise, but the useful component of some gradual changes simultaneously It is passivated together, becomes insensitive, the gradual change of the time service information as caused by the clock deviation of time dissemination system, such as the Doppler in passive location The gradual change of pulse arrival time caused by effect.Traditional linear filter method such as Kalman filtering has ratio just for linear signal Preferable effect.

Summary of the invention

It is a primary object of the present invention to overcome drawbacks described above in the prior art, propose a kind of according to a series of observation Value obtains the Noise Filter method of the signal specific of more accurate guess value, can effectively weaken random noise influence.

The present invention adopts the following technical scheme:

A kind of Noise Filter method of signal specific, it is characterised in that: preset TmIndicate Discrete Periodic Signals M-th of supposition amount, AmFor m-th of observation with certain error, K is filtering depth, and △ is adjustment stepping, and R is discrete week The jitter amplitude of phase signal, filter step are as follows:

1) when initial, T is enabledm=Am, n=1;

2) according to current observation AmWith the relationship of mapping function f (m-n), guess value T is obtainedm-n, and by Tm-nWith M-n observation Am-nIt makes comparisons, if Tm-n+ R < Am-n, then according to adjustment stepping △ to TmCarry out increasing amount adjustment;If Tm-n+ R > Am-n, then to TmCarry out decreasing amount adjustment;N=n+1;

3) n < K is judged, if so, returning to 2);If it is not, then exporting Tm

Preferably, Z is presetmFor the true value of m-th of not no error, mapping relations: Z are obeyedm-n=Zm+f(m- n);In step 2), it is assumed that current observation is equal to true value, then the current observation and mapping function f (m-n) Relationship is Am-n=Am+ f (m-n), the guess value obtained are Tm-n=Tm+f(m-n)。

Preferably, it is described according to adjustment stepping △ to TmCarry out increasing amount adjustment, including the two add, is subtracted, is multiplied or Division operation.

Preferably, it is described according to adjustment stepping △ to TmCarry out decreasing amount adjustment, including the two add, is subtracted, is multiplied or Division operation.

Preferably, the R is related to the statistical properties of random noise.

Preferably, the R takes the 75%-85% of the random noise of the jitter amplitude of Discrete Periodic Signals.

By the above-mentioned description of this invention it is found that compared with prior art, the invention has the following beneficial effects:

Filtering method of the invention is according to a series of observation AmObtained a series of guess value Tm, in conjunction with discrete periodic Guess value is compared by the jitter amplitude R of signal with observation, according to comparison result to TmIncrement or decreasing amount adjustment are carried out, with Adjustment number increase, TmIt can become closer to Zm, can effectively weaken random noise rmIt influences.In addition, filtering side of the invention Method has relatively good filter effect to linear signal and nonlinear properties (such as periodic signal).

Detailed description of the invention

Fig. 1 is the flow chart of filtering method of the invention;

Fig. 2 is linear signal filter effect figure of the invention;

Fig. 3 is cyclical signal filter effect figure of the invention.

Specific embodiment

Below by way of specific embodiment, the invention will be further described.

For signal specific, it is assumed that it is superimposed with the random noise that 75%-85% (preferably 80%) amplitude is [- R, R].It should The mathematical model of signal specific is Am=Zm+rm,rmFor 80% between section [- R ,+R] random signal m-th of value, AmIt is The m observations with certain error, ZmFor the true value of m-th of not error, ZmObey known mapping relations Zm-n= Zm+f(m-n).For example, indicating f (m-n)=C (m-n) when linear signal, C is constant;Indicate f (m-n)=0 when periodic signal, As m-n=NT, N is integer, T is the period.

Referring to Fig.1, the Noise Filter method of a kind of signal specific of the invention presets K as filtering depth, △ To adjust stepping, R is the jitter amplitude of Discrete Periodic Signals.Can basis signal sample fluctuation characteristic, choose suitable K, △, R, filter step are as follows:

1) when initial, T is enabledm=Am, n=1.

2) assume that current observation is equal to true value, then the current observation and the relationship of mapping function f (m-n) are Am-n=Am+ f (m-n), the guess value obtained are Tm-n=Tm+ f (m-n) (n=NT, N are integer when periodic signal, and T is the period). By Tm-nWith the m-n observation Am-nIt makes comparisons, if Tm-n+ R < Am-n, then according to adjustment stepping △ to TmIncreasing amount adjustment is carried out, That is Tm=Tm+△;If Tm-n+ R > Am-n, then to TmCarry out decreasing amount adjustment, i.e. Tm=Tm-△;N=n+1.

3) n < K is judged, if so, returning to 2);If it is not, then exporting Tm.With increasing for adjustment number, the T of outputmIt will be more Come closer to Zm, therefore observation A can be replacedmIt exports as a result.

Method of the invention, above-mentioned increasing amount adjustment and decreasing amount adjustment also can be used other operation modes such as multiplication, division real It is existing.

Applicating example

Such as there is the 1PPS shake model of navigation time dissemination system similar in the signal specific model with the application.Its Am=Zm+ rm, Zm-n=Zm+ (m-n) * 1*f, f are the precise frequencies of local clock, rmIt is random between section [- 100ns ,+100ns] Noise.AmFor the observation of m-th of arrival time of 1PPS, ZmFor the true value of m-th of arrival time of 1PPS, rmFor time service of navigating The random clock deviation of the 1PPS of system.It can obtain a series of being substantially better than 1PPS arrival time and observing by the filtering algorithm of the application Measure AmSupposition arrival time Tm, and then restore the 1PPS that precision is better than ± 100ns.

It is linear signal filter effect figure, the sample of signal are as follows: f (m-n)=0, r referring to Fig. 2mFor fluctuation in [- 10 10] Random noise be superimposed one [- 50 50] infrequent noise.Abscissa is signal serial number, ordinate signal value.

It is cyclical signal filter effect figure, sample of signal referring to Fig. 3 are as follows: the period is the cosine wave of 50 signaling points, i.e. f (m-n)=0, working as n=NT, N is integer, when T=50, rmRandom noise to fluctuate in [- 10 10] is superimposed one [- 50 50] Infrequent noise.Abscissa is signal serial number, ordinate signal value.

The above is only a specific embodiment of the present invention, but the design concept of the present invention is not limited to this, all to utilize this Design makes a non-material change to the present invention, and should all belong to behavior that violates the scope of protection of the present invention.

Claims (6)

1. a kind of Noise Filter method of signal specific, it is characterised in that: preset TmIndicate Discrete Periodic Signals m A supposition amount, AmFor m-th of observation with certain error, K is filtering depth, and △ is adjustment stepping, and R is discrete periodic letter Number jitter amplitude, filter step is as follows:
1) when initial, T is enabledm=Am, n=1;
2) according to current observation AmWith the relationship of mapping function f (m-n), guess value T is obtainedm-n, and by Tm-nWith m-n Observation Am-nIt makes comparisons, if Tm-n+ R < Am-n, then according to adjustment stepping △ to TmCarry out increasing amount adjustment;If Tm-n+ R > Am-n, then To TmCarry out decreasing amount adjustment;N=n+1;
3) n < K is judged, if so, returning to 2);If it is not, then exporting Tm
2. a kind of Noise Filter method of signal specific as described in claim 1, it is characterised in that: preset ZmFor The true value of m-th of not no error obeys mapping relations: Zm-n=Zm+f(m-n);In step 2), it is assumed that current observation Value is equal to true value, then the current observation and the relationship of mapping function f (m-n) are Am-n=Am+ f (m-n), obtains Guess value is Tm-n=Tm+f(m-n)。
3. a kind of Noise Filter method of signal specific as described in claim 1, it is characterised in that: in step 2), It is described according to adjustment stepping △ to TmIncreasing amount adjustment is carried out, including the two add, is subtracted, is multiplied or division operation.
4. a kind of Noise Filter method of signal specific as described in claim 1, it is characterised in that: in step 2), It is described according to adjustment stepping △ to TmDecreasing amount adjustment is carried out, including the two add, is subtracted, is multiplied or division operation.
5. a kind of Noise Filter method of signal specific as described in claim 1, it is characterised in that: the R and random The statistical properties of noise are related.
6. a kind of Noise Filter method of signal specific as described in claim 1, it is characterised in that: the R takes discrete The 75%-85% of the random noise of the jitter amplitude of periodic signal.
CN201710182251.XA 2017-03-24 2017-03-24 A kind of Noise Filter method of signal specific CN106814601B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710182251.XA CN106814601B (en) 2017-03-24 2017-03-24 A kind of Noise Filter method of signal specific

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710182251.XA CN106814601B (en) 2017-03-24 2017-03-24 A kind of Noise Filter method of signal specific

Publications (2)

Publication Number Publication Date
CN106814601A CN106814601A (en) 2017-06-09
CN106814601B true CN106814601B (en) 2019-06-04

Family

ID=59115500

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710182251.XA CN106814601B (en) 2017-03-24 2017-03-24 A kind of Noise Filter method of signal specific

Country Status (1)

Country Link
CN (1) CN106814601B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1266343A (en) * 1999-03-03 2000-09-13 株式会社日立制作所 Mobile terminal, base station and its synchronizing control method
CN103454905A (en) * 2013-09-13 2013-12-18 烟台东方英达康自动化技术有限公司 Method for synchronizing time
CN103873223A (en) * 2014-02-25 2014-06-18 上海微小卫星工程中心 Satellite internal equipment timing method and system
WO2018093046A1 (en) * 2016-11-15 2018-05-24 전자부품연구원 Utc time synchronization method of device using gps that improves abnormal operating situation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1266343A (en) * 1999-03-03 2000-09-13 株式会社日立制作所 Mobile terminal, base station and its synchronizing control method
CN103454905A (en) * 2013-09-13 2013-12-18 烟台东方英达康自动化技术有限公司 Method for synchronizing time
CN103873223A (en) * 2014-02-25 2014-06-18 上海微小卫星工程中心 Satellite internal equipment timing method and system
WO2018093046A1 (en) * 2016-11-15 2018-05-24 전자부품연구원 Utc time synchronization method of device using gps that improves abnormal operating situation

Also Published As

Publication number Publication date
CN106814601A (en) 2017-06-09

Similar Documents

Publication Publication Date Title
US6646587B2 (en) Doppler radar apparatus
EP0518442B1 (en) FMCW radar range calibration appts
EP2659239B1 (en) Radar level gauging using frequency modulated pulsed wave
EP0797105B1 (en) Method for measuring the time of flight of electric, electromagnetic or acoustic signals
CN1159592C (en) Digital measuring method of frequency and phase
EP1744176B1 (en) Distance measuring device, distance measuring method and distance measuring program
AU596497B2 (en) Doppler radar method and apparatus for measuring muzzle velocity
US6628229B1 (en) Stabilization of oscillators in a radar level transmitter
US4766436A (en) Radar device for measuring the distance of the device to a surface
FR2722005A1 (en) Apparatus and method for attenuating ambiguits in pulse doppler radars
GB1109710A (en) Self-calibrating, self-testing radio altimeter
Chereskin et al. Identifying and screening filter skew and noise bias in acoustic Doppler current profiler measurements
CN100434869C (en) Flow meter filter system and method
US20070210955A1 (en) Error corrector for radar timing systems
US2418538A (en) Measurement of distance by frequency-modulated carrier wave
US4061030A (en) Method and device for correcting the output signal from a digital transducer for measuring a physical magnitude or variable
EP0764277A1 (en) Process for determining the intermediate frequency deviation in frequency pulsed radar systems
RU2269206C2 (en) Method for estimation of doppler distribution/speed in mobile wireless communication devices, variants
EP2833548A3 (en) Digitally controlled edge interpolator (DCEI) for digital-to-time converters (DTC)
EP1798568B1 (en) Systems and methods for self-test of a radar altimeter
US8730093B2 (en) MFPW radar level gauging with distance approximation
US7752010B2 (en) Distance measuring apparatus capable of controlling range and resolution
EP0501566A1 (en) Linearizing a swept-frequency radar
US9026403B2 (en) Frequency measurement device and electronic device
Roscoe et al. Frequency and fundamental signal measurement algorithms for distributed control and protection applications

Legal Events

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