CN110346772A - A kind of higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method - Google Patents
A kind of higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method Download PDFInfo
- Publication number
- CN110346772A CN110346772A CN201910777757.4A CN201910777757A CN110346772A CN 110346772 A CN110346772 A CN 110346772A CN 201910777757 A CN201910777757 A CN 201910777757A CN 110346772 A CN110346772 A CN 110346772A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- time
- ionosphere
- significantly
- 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.)
- Pending
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/03—Details of HF subsystems specially adapted therefor, e.g. common to transmitter and receiver
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/418—Theoretical aspects
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention provides a kind of higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection methods, the contamination model of sea clutter is broadened using high frequency day ground wave radar as background constructing, and significantly ionosphere phase diameter Disturbance Rejection method is proposed based on Generalized parametering time frequency analysis (GPTF) algorithm, and further pass through the validity of the proposed method of simulating, verifying.Theory analysis and emulation experiment show that the processing method can effectively inhibit significantly ionosphere phase diameter disturbance, and have the advantages that high resolution, without cross term.
Description
Technical field
According to the present invention is a kind of higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method, and the technology is especially suitable
The target detection under sea clutter background is broadened by ionosphere phase diameter disturbing influence for high-frequency sky-wave over-the-horizon radar.
Background technique
High-frequency sky-wave over-the-horizon radar has many advantages, such as that sky wave propagation distance is remote, coverage area is big, Anti-antiradiation missile, by
To the extensive attention of domestic and international research institution.However the influence that sky-wave OTH radar will receive the disturbance of ionosphere phase diameter leads to sea
Foreign echo spectrum widening, especially significantly the disturbance of ionosphere phase diameter will lead to sea clutter spectrum widening seriously and complex shape, seriously
Influence the detection of Ship Target.Phase diameter disturbance in ionosphere can cause sea clutter Bragg video stretching, influence ship target at a slow speed and examine
It surveys.Most of ocean clutter cancellation algorithm obtains in the case where being all based on sea clutter instantaneous frequency as steady state value at present, electricity
Sea clutter Bragg frequency bandspread caused by absciss layer phase diameter disturbs will lead to the failure of ocean clutter cancellation algorithm or reduce clutter recognition effect
Fruit.Therefore the influence for needing to eliminate the disturbance of ionosphere phase diameter before carrying out ocean clutter cancellation first, sharpens the sea clutter of broadening
Spectrum, i.e. ionosphere phase diameter Disturbance Rejection are handled, and are then carried out ocean clutter cancellation processing again, are conducive to the detection of Ship Target.
Current ionosphere phase diameter Disturbance Rejection method is mentioned by the variation of estimation broadening echo spectrum instantaneous frequency
Frequency modulation(PFM) function is taken, to obtain ionospheric phase variation correction function, then echo is believed using the correction function obtained
Number correction, can so that broadening echo spectrum sharpen, improve the target detection performance of radar.Such method mainly has maximum entropy at present
Power estimation method, PGA method, minimum entropy search method, the solution pollution method based on feature decomposition and piecewise polynomial modeling etc..This
A little methods require the peak Bragg for extracting broadening first to do instantaneous Frequency Estimation.When phase perturbation amplitude is smaller, sea clutter
Two single order Bragg peak stretching degree are smaller in spectrum, can use band-pass filtering method at this time and filter out one of peak Bragg
Estimate disturbing function.However when significantly phase diameter disturbs in ionosphere, the Bragg summit of two broadenings is overlapping together,
It is difficult effectively to extract single Bragg signal using based on band-pass filtering method at this time, so as to cause the processing of phase diameter Disturbance Rejection
Poor effect.
The case where for the disturbance of significantly phase diameter, time frequency analysis is the relatively good method of target.It is currently used to be mainly
Imparametrization Time-Frequency Analysis Method, as Short Time Fourier Transform (STFT), Winger-Will are distributed (WVD), puppet WVD distribution
(PWVD) and smooth-switch method distribution (SPWVD) etc..But such method is more sensitive for signal-to-noise ratio and by cross term interference
With the influence of resolution ratio limitation.The Major Technology detected to folded Clutter in Skywave Radars ship target is ocean clutter cancellation, representative
Algorithm include sea clutter opposition method, subspace class method, based on singular value decomposition (SVD) ocean clutter cancellation algorithm and from
Adaptive filter method etc., but these algorithms are difficult to the broadening sea clutter to non-stationary and are effectively inhibited.When with imparametrization
Frequency analysis method is compared, and Parametric Time-frequency Analysis method selects suitable description non-stationary by the model of introducing apriori signals
The core of signal can effectively improve time frequency resolution, typical method has adaptively when kernel form meets with analyzed signal
Chirplet decomposition, Atomic Decomposition and multinomial Fourier transformation etc., but such method mainly uses polynomial kernel, it is uncomfortable
Faster strong time-varying non-stationary signal is changed over time in analysis.
Patent CN201510607745.9 discloses a kind of time-variable filtering parameter based on time frequency analysis and generates with realizing
System and method, this method are related to communication and field of signal processing, in order to preferably filter out multi-user interference and noise, lead to
Signal can be filtered according to the variation of the bandwidth of signal by crossing this method, preferably filter out multi-user interference and noise.
Patent CN201611139350.1 discloses a kind of start and stop vehicle fault signature extraction based on Parametric Time-frequency Analysis
And diagnostic method, this method decompose rotor start and stop vehicle vibration signal using Parametric Time-frequency Analysis and acquire each frequency component
Then accurate amplitude and phase combine with holography spectrum and efficiently extract rotor fault feature, it is right under steady-state speed to avoid
Mistaken diagnosis, omitting diagnosis when rotor progress fault diagnosis, breaching conventional photographic spectral technology can be only applied to stationary vibration signal
Limitation, expanded its application range, can be applied to fault diagnosis and status monitoring under large rotating machinery non-stationary operating condition.
Patent CN201611169629.4 discloses a kind of high-order PPS being fitted based on time frequency analysis and instantaneous frequency profile
Modulated parameter estimating method, this method can effectively inhibit the cross term interference of high-order PPS signal, estimate not under low signal-to-noise ratio
Know the PPS signal phase parameter of order, there is preferable parameter Estimation performance, overcome conventional method by time-frequency cross term interference
Influence, be of great significance to the subsequent processing and signature analysis of non-stationary signal.
Patent CN201510922006.9 discloses a kind of time frequency analysis side of improvement generalized S-transform based on parameter optimization
One is introduced into the window function of S-transformation by method, this method by the function of first order of independent variable of frequency, so that window function form spirit
Work is adjustable, and is optimized by parameter of the time-frequency concentration class to function of first order, can improve the energy accumulating of signal time-frequency distributions
Property, the estimated accuracy of signal transient parameter is improved, and calculation amount is small, suitable for communication, radar, earthquake and biomedicine signals
Analysis and processing.
In journal article " the instantaneous micro-doppler frequency extraction method of skirt target is bored in the precession based on Parametric Time-frequency Analysis ",
The method that author proposes is directed to the cone skirt target echo of multicomponent data processing composition, is interfered using coherent signal list range Doppler
(CSRDI) method estimates coning frequency, and then utilizes the micro-doppler curve of Parametric Time-frequency Analysis estimation scattering point, Zhi Houli
The scattering point echo-signal obtained with bandstop filter separation estimation, based on the effective of the mentioned method of Electromagnetic Simulation data verification
Property.
When journal article " the rotor full working scope dynamic balance method based on Parametric Time-frequency Analysis " is proposed based on parametrization
The rotor full working scope dynamic balance method of frequency analysis.This method can be convenient, quickly determine rotor unbalance amount and unbalance side
Position, is effectively reduced rotor-support-foundation system unbalance vibration, while reducing in equilibrium process and opening train number number.
Summary of the invention
Ship detection is difficult under broadening echo spectrum background to solve the problems, such as higher-frequency radar ionosphere contamination, and the present invention mentions
A kind of significantly ionosphere phase diameter Disturbance Rejection method includes following procedure out:
Step 1 carries out Generalized parametering time frequency analysis using formula to sea clutter signal and obtains GPTF time-frequency distributions GPTF1
[·]；
The GPTF that step 2 obtains step 11[] carries out feature decomposition, and each to the matrix reconstruction after feature decomposition
Signal xi(n), i=1,2 ..., wherein xiIt (n) is the positive and negative Bragg instantaneous frequency of sea clutter；
Each x after step 3 pair reconstructi(n), i=1,2 ... it is obtained by Generalized parametering time frequency analysis again respectively
Corresponding GPTF time-frequency distributions GPTF2[·]；
The each GPTF of step 42[] carries out optimal path algorithm detection respectively, extracts the instantaneous frequency of ionospheric phase variation
Rate；
Step 5 integrates the ionosphere modulating frequency of extraction to obtain phase correction function, correction pollution sea clutter echo letter
Number.
Optionally, step 5 further includes following procedure:
Indicate the peak the sea clutter Bragg instantaneous frequency extracted, fBIt is the peak First-order sea clutter Bragg frequency, then it is electric
Frequency modulation(PFM) caused by absciss layer are as follows:
Therefore, ionospheric phase variation function are as follows:
The compensated echo-signal result of phase perturbation are as follows:
In formula, A (t) is echo amplitude, and θ (t) is phase of echo when ionosphere is steady, and m (t) is ionospheric disturbance phase
Position.
Optionally, optimal path described in step 4 shows on time-frequency domain for the sea clutter by ionospheric phase variation
The instantaneous frequency distribution curve changed over time out, optimal path detection algorithm are adaptive in such a way that signal energy is accumulated
It searches for the time-frequency distributions curve of energy peak signal in time-frequency figure with answering, and extracts the instantaneous frequency profile of signal.
Optionally, the Generalized parametering time frequency analysis contaminated sea clutter signal carried out, and to each of after reconstruct
The Generalized parametering time frequency analysis that signal carries out respectively includes: by pull-in frequency rotation and translation operator, being rotated using frequency
Operator rotates the time-frequency characteristics of non-stationary signal, and signal is made to tend to be steady, and is then Fu in short-term to postrotational signal
In leaf transformation, signal time-frequency characteristics are finally moved into true Location of ridge axis using frequency translation operator.
Optionally, the Generalized parametering time frequency analysis in step 1 and step 3 is to carry out following operation to corresponding signal:
Wherein, sr(τ) indicates that instantaneous frequency is the complex signal of arbitrary function, κP(τ) indicates Generalized parametering time frequency analysis
Transformation kernel, P indicate transformation nuclear parameter,Indicate frequency rotation operator,Indicate frequency translation operator, gσ
(τ) indicates the window function of time-frequency conversion,Indicate gσThe conjugation of (τ), f indicate the instantaneous frequency of signal.
Optionally, Generalized parametering time frequency analysis further comprises: by transformation kernel parameters precision in conjunction with time-frequency concentration degree,
Transformation nuclear parameter is iterated and asks optimal, i.e., cyclic approximation refinement is carried out to signal time-frequency characteristics and obtains most suitable transformation kernel
Parameter determines frequency rotation operator and translation operator, using determining frequency rotation operator and frequency translation operator to non-stationary
Signal carries out rotation translation, and signal is made to tend to be steady.
Optionally, significantly ionosphere phase diameter Disturbance Rejection method further includes pollution sea clutter number to the higher-frequency radar
Input signal mould according to generation, when broadening sea clutter spectrum model after establishing ionosphere contamination is as ionosphere phase diameter Disturbance Rejection
Type；
Broadening sea clutter spectrum model after the ionosphere contamination is
In formula, N is beam position number, and K is distance unit number, and l is range gate, θiFor beam position,For sea
The amplitude of clutter,For first of range gate, θiThe First-order sea clutter of a beam position, k-th distance unit
Bragg frequency, φ (t) are ionosphere phase diameter disturbing function.
Optionally, significantly ionosphere phase diameter Disturbance Rejection method is significantly ionized the higher-frequency radar suitable for processing
False-alarm and false dismissal situation caused by the peak sea clutter Bragg caused by layer phase diameter disturbs significantly is broadened or divided.
Optionally, significantly ionosphere phase diameter Disturbance Rejection method is suitable for HF skywave radar to the higher-frequency radar.
Present invention implementation provides a kind of higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method.With high frequency day earthwave
Radar is background, establishes the contamination model of broadening sea clutter；Significantly ionosphere phase diameter disturbance is proposed based on GPTF algorithm
Suppressing method；Pass through the validity of the proposed method of simulating, verifying.
Technical solution provided by the invention, has at least the following technical effects or advantages:
(1) compared with the Time-frequency methods such as classical STFT and SPWVD, the significantly ionosphere phase diameter Disturbance Rejection based on GPTF
The time frequency resolution of method is high, not by cross term interference, can disturb and carry out to significantly ionosphere phase diameter under low miscellaneous noise ratio
Effectively inhibit.
(2) compared with the ionosphere phase diameter Disturbance Rejection method with traditional based on time frequency analysis, the when frequency division based on GPTF
Analysis method is suitable for analysis and changes over time faster strong time-varying non-stationary ionosphere contamination echo-signal.
Detailed description of the invention
Fig. 1 is based on GPTF significantly ionosphere phase diameter Disturbance Rejection method and step
Fig. 2 transformation kernel parameter evaluation method flow chart
Fig. 3 is not by the sea clutter spectrogram of ionosphere phase diameter disturbance pollution
Sea clutter time-frequency figure of the Fig. 4 by the disturbance pollution of ionosphere phase diameter
The pollution function that Fig. 5 is extracted based on GPTF Time-Frequency Analysis Method
The error precision for the pollution function that Fig. 6 is extracted based on GPTF Time-Frequency Analysis Method
Specific embodiment
The present invention provides a kind of higher-frequency radar based on Generalized parametering time frequency analysis, significantly phase diameter in ionosphere is disturbed
Suppressing method.Firstly, this method using high frequency day ground wave radar as background, establishes the contamination model of broadening sea clutter；Secondly, base
Significantly ionosphere phase diameter Disturbance Rejection method is proposed in GPTF algorithm；Finally, passing through having for mentioned method of simulating, verifying
Effect property.Theory analysis and emulation experiment show that this method can effectively inhibit significantly ionosphere phase diameter disturbance, and have resolution ratio
The advantages of high, without cross term.
With reference to the accompanying drawings and embodiments, the present invention will be described in detail.
(1) broadening sea clutter modeling
The present invention is using high frequency day ground wave radar as application background.High frequency day ground wave radar is a kind of New System early warning radar,
Wherein ionosphere is important propagation medium, and phase diameter disturbance in ionosphere influences mainly by Ionospheric Parameters the broadening that sea clutter is composed
Variation range and parameter perturbation amplitude determine that the phase perturbation function under different ionospheric disturbance states can use multinomial phase
Phase perturbation model describes, therefore the sea clutter spectrum model after ionosphere contamination can be described as:
In formula, N is beam position number, and K is distance unit number, and l is range gate, θiFor beam position,For sea
The amplitude of clutter,For first of range gate, θiA beam position, k-th distance unit First-order sea clutter Bragg
Frequency, φ (t) are ionosphere phase diameter disturbing function.
Broadening sea clutter model after ionosphere contamination is that the ionosphere phase diameter proposed by the present invention based on GPFT disturbs suppression
Algorithm processed provides basic input signal model.
(2) basic principle of GPTF
Time frequency analysis can indicate the relationship between the Time And Frequency of non-stationary signal, become and analyze having for such signal
Power tool.Short Time Fourier Transform (STFT) is one of fairly simple method.This method is the base in Fourier transformation
To signal windowing process on plinth, signal is approximate steady in default windows, it can portray signal in time domain and frequency domain simultaneously, thus instead
Signal spectrum has been reflected to change with time characteristic.Due to regular length window function cannot signal acquisition frequency in time variation,
Therefore the time-frequency representation concentration degree of short time discrete Fourier transform is poor, is unable to the time-frequency characteristics of such signal of accurate description.Generalized Parameters
Change Time-Frequency Analysis Method is rotated by pull-in frequency and translation operator, using frequency rotation operator that the time-frequency of non-stationary signal is special
Sign is rotated, and so that signal is tended to be steady, is then STFT to postrotational signal, finally utilizes frequency translation operator by signal
Time-frequency characteristics move to true Location of ridge axis.Since the analyzed signal of STFT is approximate stationary signal, so Generalized parametering time-frequency
Analysis can effectively improve time frequency resolution, and cross term interference is not present.
1. the basic principle of GPTF
Assuming that signal transient frequency is arbitrary function fm(t), then complex signal form at this time is
Signal, which does Generalized parametering time frequency analysis, to be expressed as form:
Wherein, sr(τ) indicates that instantaneous frequency is the complex signal of arbitrary function, κP(τ) indicates Generalized parametering time frequency analysis
Transformation kernel, P indicate transformation nuclear parameter,Indicate frequency rotation operator,Indicate frequency translation operator, gσ(τ) is indicated
The window function of time-frequency conversion,Indicate gσThe conjugation of (τ), f indicate the instantaneous frequency of signal.According to formula (1), Generalized Parameters
The basic process for changing time frequency analysis can be described as: the rotating signal first in time frequency plane, i.e., by signal transient frequency fm(τ) subtracts
Remove κP(τ), then with κP(t0) it is that frequency increment translates signal on time frequency plane, finally signal is done using window function
Short time discrete Fourier transform.Work as κP(τ)=fmWhen (τ), signal srShown in the signal form such as formula (2) of (τ) after rotation translation:
At this point, signalInstantaneous frequency κP(t0) τ does not change at any time, signalFor stationary signal.At this time to letter
NumberShort Time Fourier Transform is done, compared to non-stationary signal, time-frequency concentration degree is improved significantly.It can be with from (1), (2)
Find out, by choosing suitable transformation nuclear parameter, so that it is determined that frequency rotation operator and translation operator, are revolved using determining frequency
Turn operator and frequency translation operator carries out rotation translation to non-stationary signal, so that signal be made to tend to be steady, reaches raising time-frequency
The purpose of resolution ratio.As can be seen that the time frequency resolution and transformation kernel parameters precision of Generalized parametering time frequency analysis are closely related.
From be analyzed above it is found that the key of Generalized parametering time frequency analysis be estimate transformation kernel parameter.
2. converting the estimation method of nuclear parameter
The working mechanism of GPTF is actually based on the short time discrete Fourier transform of frequency rotation operator and translation operator.It is elected
When the transformation kernel selected more is matched with the time-frequency characteristics of analyzed signal, the signal after rotation translation is done just closer to steady to signal
Signal, the time frequency resolution after doing Short Time Fourier Transform is also higher, also more accurate to portraying for non-stationary signal time-frequency characteristics.
It can be seen that the precision of transformation nuclear parameter determines the performance of Generalized parametering Time-frequency method, so that it is non-flat to affect its analysis
The performance of steady signal.Time-frequency characteristics approximation theory points out that, when transformation kernel parameter Estimation is more accurate, time-frequency concentration degree is higher, time-frequency
Concentration degree is higher, and it is more accurate that time-frequency curve extracts, then the transformation nuclear parameter estimated is more accurate, higher so as to obtain
Time-frequency concentration degree.
Generalized parametering Time-Frequency Analysis Method is by conjunction with time-frequency concentration degree, realizing transformation kernel parameters precision to change
The iteration for changing nuclear parameter asks optimal.Above-mentioned principle is actually based on Generalized parametering time frequency analysis and carries out to signal time-frequency characteristics
Cyclic approximation refinement, to obtain most suitable transformation nuclear parameter.As shown in Fig. 2, converting nuclear parameter according to time-frequency approximation theory
The detailed step of estimation is as follows:
Step H1 determines transformation kernel form according to signal form priori；
Step H2, the number of iterations i=0, initialization transformation nuclear parameter P0；
Step H3 is Parametric Time-frequency Analysis TF (t to signal according to formula (1)0,f；Pi)；
Step H4 detects the peak ridge in time-frequency figure
Step H5 selects suitable approximating method, fitting crestal line estimation transformation nuclear parameter according to crestal line form
Step H6, the number of iterations i increase by 1, update transformation nuclear parameter
Step H7 calculates termination conditionWherein,Indicate i-th
The peak ridge of iteration,Indicate the peak ridge of (i-1)-th iteration, it should be pointed out that
Termination condition parameter ζ is arranged in step H8；
Step H9 compares the size of Λ and ζ, as Λ > ζ, return step H3；As Λ≤ζ, step H10 is carried out；
Step H10, output transform nuclear parameter obtain time-frequency curve.
(3) the significantly ionosphere phase diameter Disturbance Rejection algorithm based on GPTF
Firstly, the sea clutter signal after acquisition or generation ionosphere contamination, then utilizes the time frequency analysis side based on GPTF
Method analyzes pollution echo-signal, recycles the pollution function extracted in time-frequency figure based on optimal path detection algorithm, most
The pollution function compensation pollution echo-signal extracted is utilized afterwards.
As shown in Figure 1, a kind of higher-frequency radar of the present invention significantly disturb by ionosphere phase diameter Disturbance Rejection method, the phase diameter
Dynamic suppressing method is suitable for the peak sea clutter Bragg caused by processing is disturbed by significantly ionosphere phase diameter and significantly broadens or divide
Caused false-alarm and false dismissal situation include following procedure:
Step 1 carries out Generalized parametering time frequency analysis to contaminated sea clutter signal and obtains GPTF time-frequency distributions GPTF1
[·]；
The GPTF that step 2 obtains step 11[] carries out feature decomposition, and each to the matrix reconstruction after feature decomposition
Signal xi(n), i=1,2 ..., wherein xiIt (n) is the positive and negative Bragg instantaneous frequency of sea clutter；
Each signal x after step 3 pair reconstructi(n), i=1,2 ... it is obtained again by Generalized parametering time frequency analysis respectively
To its corresponding GPTF time-frequency distributions GPTF2[·]；
The each GPTF of step 42[] carries out optimal path algorithm detection respectively, extracts the instantaneous frequency of ionospheric phase variation
Rate；The optimal path be by ionospheric phase variation changing over time of being shown on time-frequency domain of sea clutter it is instantaneous
Curve of frequency distribution, optimal path detection algorithm are adaptively to search for energy in time-frequency figure in such a way that signal energy is accumulated
The time-frequency distributions curve of peak signal, and extract the instantaneous frequency profile of signal；
Step 5 integrates the ionosphere modulating frequency of extraction to obtain phase correction functionCorrection pollution sea clutter
Echo-signal.Assuming thatThe peak the sea clutter Bragg instantaneous frequency extracted is indicated, then frequency modulation(PFM) caused by ionosphere
Are as follows:
In formula, fBIt is the peak First-order sea clutter Bragg frequency.Therefore, ionospheric phase variation function are as follows:
Therefore, the compensated echo-signal result of phase perturbation are as follows:
In formula, A (t) is echo amplitude, and θ (t) is phase of echo when ionosphere is steady, and m (t) is ionospheric disturbance phase
Position.
Simulation analysis is carried out according to above step.
Fig. 3 show original unpolluted sea clutter spectrogram, and Fig. 4 is shown after by the disturbance pollution of ionosphere phase diameter
Sea clutter time-frequency figure (positive Bragg frequency).Wherein, ionosphere contamination function is φ (t)=exp (jm2 π sin (2 π f0·
t))·exp(j2πfBT), wherein fBFor sea clutter Bragg frequency, m and f0It is the parameter of ionosphere contamination function.Fig. 5 is shown
The ionosphere contamination function extracted based on context of methods.Fig. 6 show the ionosphere contamination function extracted based on context of methods
Error curve.As seen from Figure 5, the pollution function curve that algorithm proposed by the invention extracts is kissed substantially with ideal curve
It closes, and Fig. 6 can be seen that the root mean square between the algorithm proposed by the invention pollution function curve extracted and ideal curve
Error is small than the root-mean-square error of traditional Short Time Fourier Transform (STFT) method.
It is discussed in detail although the contents of the present invention have passed through above preferred embodiment, but it should be appreciated that above-mentioned
Description is not considered as limitation of the present invention.After those skilled in the art have read above content, for of the invention
A variety of modifications and substitutions all will be apparent.Therefore, protection scope of the present invention should be limited to the appended claims.
Claims (9)
1. a kind of higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method, which is characterized in that including following procedure:
Step 1 carries out Generalized parametering time frequency analysis to contaminated sea clutter signal and obtains GPTF time-frequency distributions GPTF1[·]；
The GPTF that step 2 obtains step 11[] carries out feature decomposition, and to each signal x of matrix reconstruction after feature decompositioni
(n), i=1,2 ..., wherein xiIt (n) is the positive and negative Bragg instantaneous frequency of sea clutter；
Each signal x after step 3 pair reconstructi(n), i=1,2 ... it is obtained by Generalized parametering time frequency analysis again respectively
Corresponding GPTF time-frequency distributions GPTF2[·]；
The each GPTF of step 42[] carries out optimal path algorithm detection respectively, extracts the instantaneous frequency of ionospheric phase variation；
Step 5 integrates the ionosphere modulating frequency of extraction to obtain phase correction function, correction pollution sea clutter echo-signal.
2. higher-frequency radar as described in claim 1 significantly ionosphere phase diameter Disturbance Rejection method, which is characterized in that step 5 into
One step includes following procedure:
Indicate the peak the sea clutter Bragg instantaneous frequency extracted, fBIt is the peak First-order sea clutter Bragg frequency, then ionosphere
Caused frequency modulation(PFM) are as follows:
Therefore, ionospheric phase variation function are as follows:
The compensated echo-signal result of phase perturbation are as follows:
In formula, A (t) is echo amplitude, and θ (t) is phase of echo when ionosphere is steady, and m (t) is ionospheric disturbance phase.
3. higher-frequency radar as described in claim 1 significantly ionosphere phase diameter Disturbance Rejection method, which is characterized in that in step 4
The optimal path is the instantaneous frequency changed over time shown on time-frequency domain by the sea clutter of ionospheric phase variation
Rate distribution curve, optimal path detection algorithm are adaptively to search in time-frequency figure energy most in such a way that signal energy is accumulated
The time-frequency distributions curve of big signal, and extract the instantaneous frequency profile of signal.
4. higher-frequency radar as described in claim 1 significantly ionosphere phase diameter Disturbance Rejection method, which is characterized in that contaminated
The Generalized parametering time frequency analysis that carries out of sea clutter signal, and to the Generalized parametering that each signal after reconstruct carries out respectively
Time frequency analysis includes: using frequency rotation operator that the time-frequency of non-stationary signal is special by pull-in frequency rotation and translation operator
Sign is rotated, and so that signal is tended to be steady, is then done Short Time Fourier Transform to postrotational signal, and frequency translation is finally utilized
Signal time-frequency characteristics are moved to true Location of ridge axis by operator.
5. higher-frequency radar as described in claim 1 significantly ionosphere phase diameter Disturbance Rejection method, which is characterized in that step 1 and
Generalized parametering time frequency analysis in step 3 is to carry out following operation to corresponding signal:
Wherein, sr(τ) indicates that instantaneous frequency is the complex signal of arbitrary function, κPThe transformation of (τ) expression Generalized parametering time frequency analysis
Core, P indicate transformation nuclear parameter,Indicate frequency rotation operator,Indicate frequency translation operator, gσ(τ) is indicated
The window function of time-frequency conversion,Indicate gσThe conjugation of (τ), f indicate the instantaneous frequency of signal.
6. higher-frequency radar as claimed in claim 4 significantly ionosphere phase diameter Disturbance Rejection method, which is characterized in that Generalized Parameters
Changing time frequency analysis further comprises: by transformation kernel parameters precision in conjunction with time-frequency concentration degree, being iterated and asks to transformation nuclear parameter
It is optimal, i.e., cyclic approximation refinement is carried out to signal time-frequency characteristics and obtain most suitable transformation nuclear parameter, determines frequency rotation operator
With translation operator, rotation translation is carried out to non-stationary signal using determining frequency rotation operator and frequency translation operator, makes letter
It number tends to be steady.
7. higher-frequency radar as described in claim 1 significantly ionosphere phase diameter Disturbance Rejection method, which is characterized in that further packet
The sea clutter data containing pollution generate, and the broadening sea clutter spectrum model after establishing ionosphere contamination is as ionosphere phase diameter Disturbance Rejection
When input signal model；
Broadening sea clutter spectrum model after the ionosphere contamination is
In formula, N is beam position number, and K is distance unit number, and l is range gate, θiFor beam position,For sea clutter
Amplitude,For first of range gate, θiA beam position, k-th distance unit First-order sea clutter Bragg frequency
Rate, φ (t) are ionosphere phase diameter disturbing function.
8. higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method, feature as described in any one of claim 1-7
It is, significantly ionosphere phase diameter Disturbance Rejection method is suitable for processing by the disturbance of significantly ionosphere phase diameter to the higher-frequency radar
False-alarm caused by the peak caused sea clutter Bragg is significantly broadened or divided and false dismissal situation.
9. higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method, feature as described in any one of claim 1-7
It is, significantly ionosphere phase diameter Disturbance Rejection method is suitable for HF skywave radar to the higher-frequency radar.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910777757.4A CN110346772A (en) | 2019-08-22 | 2019-08-22 | A kind of higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910777757.4A CN110346772A (en) | 2019-08-22 | 2019-08-22 | A kind of higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110346772A true CN110346772A (en) | 2019-10-18 |
Family
ID=68181059
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910777757.4A Pending CN110346772A (en) | 2019-08-22 | 2019-08-22 | A kind of higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110346772A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112114334A (en) * | 2020-09-01 | 2020-12-22 | 华东交通大学 | GNSS ionosphere abnormal disturbance parameter extraction method |
| CN112130138A (en) * | 2020-07-06 | 2020-12-25 | 电子科技大学 | Ionospheric phase pollution correction method for over-the-horizon radar |
| CN112231987A (en) * | 2020-11-04 | 2021-01-15 | 华东交通大学 | Ionosphere forecasting method based on VMD and Elman neural network |
| CN112834993A (en) * | 2020-12-31 | 2021-05-25 | 中国地质大学（武汉） | LMSCT time-frequency analysis method for IPIX radar signal target detection |
| CN113901741A (en) * | 2021-10-20 | 2022-01-07 | 哈尔滨工业大学（威海） | Ocean-ionosphere dissipation dynamics simulation method, simulation system, equipment and application |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102540162A (en) * | 2011-12-12 | 2012-07-04 | 中国船舶重工集团公司第七二四研究所 | Method for estimating low-altitude electromagnetic wave propagation characteristic on basis of sea clutter |
| CN102788969A (en) * | 2012-07-04 | 2012-11-21 | 中国人民解放军海军航空工程学院 | Sea surface micromotion target detection and feature extraction method based on short-time fractional Fourier transform |
| CN104155646A (en) * | 2014-08-22 | 2014-11-19 | 西安电子科技大学 | Sea surface low-speed motion object detection method based on Doppler spectrum signature |
| CN104215951A (en) * | 2014-08-20 | 2014-12-17 | 上海无线电设备研究所 | System and method for detecting low-speed small target under sea cluster background |
| CN104391279A (en) * | 2014-11-24 | 2015-03-04 | 哈尔滨工业大学 | Ionosphere propagation characteristic based phase diameter disturbance suppression method |
| CN106597408A (en) * | 2016-12-16 | 2017-04-26 | 重庆邮电大学 | Method for estimating high-order PPS signal parameter based on time-frequency analysis and instantaneous frequency curve-fitting |
| EP3239656A1 (en) * | 2016-04-29 | 2017-11-01 | Thales | Method for optimising the detection of sea targets and airborne radar implementing such a method |
| CN109063613A (en) * | 2018-07-20 | 2018-12-21 | 东北大学 | Based on the synchronous nonstationary random response method for extracting transformation of Generalized parametering |
| US20190159703A1 (en) * | 2017-11-29 | 2019-05-30 | Samsung Electronics Co., Ltd. | System and method for obtaining blood glucose concentration using temporal independent component analysis (ica) |
-
2019
- 2019-08-22 CN CN201910777757.4A patent/CN110346772A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102540162A (en) * | 2011-12-12 | 2012-07-04 | 中国船舶重工集团公司第七二四研究所 | Method for estimating low-altitude electromagnetic wave propagation characteristic on basis of sea clutter |
| CN102788969A (en) * | 2012-07-04 | 2012-11-21 | 中国人民解放军海军航空工程学院 | Sea surface micromotion target detection and feature extraction method based on short-time fractional Fourier transform |
| CN104215951A (en) * | 2014-08-20 | 2014-12-17 | 上海无线电设备研究所 | System and method for detecting low-speed small target under sea cluster background |
| CN104155646A (en) * | 2014-08-22 | 2014-11-19 | 西安电子科技大学 | Sea surface low-speed motion object detection method based on Doppler spectrum signature |
| CN104391279A (en) * | 2014-11-24 | 2015-03-04 | 哈尔滨工业大学 | Ionosphere propagation characteristic based phase diameter disturbance suppression method |
| EP3239656A1 (en) * | 2016-04-29 | 2017-11-01 | Thales | Method for optimising the detection of sea targets and airborne radar implementing such a method |
| CN106597408A (en) * | 2016-12-16 | 2017-04-26 | 重庆邮电大学 | Method for estimating high-order PPS signal parameter based on time-frequency analysis and instantaneous frequency curve-fitting |
| US20190159703A1 (en) * | 2017-11-29 | 2019-05-30 | Samsung Electronics Co., Ltd. | System and method for obtaining blood glucose concentration using temporal independent component analysis (ica) |
| CN109063613A (en) * | 2018-07-20 | 2018-12-21 | 东北大学 | Based on the synchronous nonstationary random response method for extracting transformation of Generalized parametering |
Non-Patent Citations (3)
| Title |
|---|
| Y. YANG; Z. K. PENG; X. J. DONG; W. M. ZHANG; G. MENG: "《General Parameterized Time-Frequency Transform》", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
| 李亚军: "《天发地收高频雷达展宽海杂波建模与抑制技术研究》", 《中国博士学位论文全文数据库 信息科技辑》 * |
| 肖金国: "《基于广义参数化时频分析的空间目标微多普勒提取与参数估计》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112130138A (en) * | 2020-07-06 | 2020-12-25 | 电子科技大学 | Ionospheric phase pollution correction method for over-the-horizon radar |
| CN112130138B (en) * | 2020-07-06 | 2022-10-14 | 电子科技大学 | Ionized layer phase pollution correction method for over-the-horizon radar |
| CN112114334A (en) * | 2020-09-01 | 2020-12-22 | 华东交通大学 | GNSS ionosphere abnormal disturbance parameter extraction method |
| CN112231987A (en) * | 2020-11-04 | 2021-01-15 | 华东交通大学 | Ionosphere forecasting method based on VMD and Elman neural network |
| CN112231987B (en) * | 2020-11-04 | 2022-05-17 | 华东交通大学 | Ionosphere forecasting method based on VMD and Elman neural network |
| CN112834993A (en) * | 2020-12-31 | 2021-05-25 | 中国地质大学（武汉） | LMSCT time-frequency analysis method for IPIX radar signal target detection |
| CN112834993B (en) * | 2020-12-31 | 2023-06-30 | 中国地质大学（武汉） | LMSCT time-frequency analysis method for IPIX radar signal target detection |
| CN113901741A (en) * | 2021-10-20 | 2022-01-07 | 哈尔滨工业大学（威海） | Ocean-ionosphere dissipation dynamics simulation method, simulation system, equipment and application |
| CN113901741B (en) * | 2021-10-20 | 2023-08-18 | 哈尔滨工业大学（威海） | Ocean-ionosphere dissipation dynamics simulation method, simulation system and equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110346772A (en) | A kind of higher-frequency radar significantly ionosphere phase diameter Disturbance Rejection method | |
| Zheng et al. | Radar high-speed target detection based on the scaled inverse Fourier transform | |
| Huang et al. | Ground maneuvering target imaging and high-order motion parameter estimation based on second-order keystone and generalized Hough-HAF transform | |
| Wang et al. | Integrated cubic phase function for linear FM signal analysis | |
| Li et al. | Separation of micro-Doppler signals based on time frequency filter and Viterbi algorithm | |
| CN106597408B (en) | High-order PPS signal parameter estimation method based on time-frequency analysis and instantaneous frequency curve fitting | |
| CN111142105A (en) | ISAR imaging method for complex moving target | |
| CN109324322A (en) | A kind of direction finding and target identification method based on passive phased array antenna | |
| Tong et al. | Parameter estimation of FH signals based on STFT and music algorithm | |
| CN108333568B (en) | Broadband echo Doppler and time delay estimation method based on Sigmoid transformation in impact noise environment | |
| Lu et al. | Ionospheric decontamination and sea clutter suppression for HF skywave radars | |
| Li et al. | Cascaded method for ionospheric decontamination and sea clutter suppression for high‐frequency hybrid sky‐surface wave radar | |
| Wang et al. | Approach for high‐resolution inverse synthetic aperture radar imaging of ship target with complex motion | |
| Asghari et al. | Doa estimation of noncircular signals under impulsive noise using a novel empirical characteristic function-based music | |
| CN107861115A (en) | A kind of OTHR maneuvering target method for parameter estimation based on instantaneous autocorrelation matrix Its Sparse Decomposition | |
| Lei et al. | Transient interference suppression method based on an improved TK energy operator and fuzzy reasoning | |
| Raju et al. | MST radar signal processing using iterative adaptive approach | |
| Gholami et al. | Two-stage estimator for frequency rate and initial frequency in LFM signal using linear prediction approach | |
| Xu et al. | Parameter estimation of underwater moving sources by using matched Wigner transform | |
| Orović et al. | A class of highly concentrated time-frequency distributions based on the ambiguity domain representation and complex-lag moment | |
| CN112014813B (en) | Method and system for correcting ionosphere pollution of sky wave radar | |
| CN112731292B (en) | Low-altitude flight target signal time delay estimation method with local IMF energy weighting | |
| CN112824927B (en) | Sky wave over-the-horizon radar ionospheric phase pollution correction method based on sparse filtering | |
| Li et al. | Computational implementation and asymptotic statistical performance analysis of range frequency autocorrelation function for radar high-speed target detection | |
| Bi et al. | Sea clutter suppression algorithm in low SCR based on improved fractional Fourier transform |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191018 |
|
| RJ01 | Rejection of invention patent application after publication |