CN108254730A - Auto-correlation function processing method at radar zero time delay based on damping fitting - Google Patents
Auto-correlation function processing method at radar zero time delay based on damping fitting Download PDFInfo
- Publication number
- CN108254730A CN108254730A CN201810319869.0A CN201810319869A CN108254730A CN 108254730 A CN108254730 A CN 108254730A CN 201810319869 A CN201810319869 A CN 201810319869A CN 108254730 A CN108254730 A CN 108254730A
- Authority
- CN
- China
- Prior art keywords
- time delay
- parameter
- value
- fitting
- incremented
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
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/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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses auto-correlation function processing methods at a kind of radar zero time delay based on damping fitting, mainly solve existing incoherent scattering radar in ionization detection layer, it is impossible to accurate the problem of obtaining information at zero time delay.Its implementation is:1) combination actual measurement autocorrelation value L is calculated using algorithm;2) establish damping function model, and in setting model each parameter value range and step-size in search;3) combination and fitting autocorrelation value S expression formulas are obtained according to damping function model, the multiple Selecting All Parameters in the range of parameter value calculate its corresponding residual error quadratic sum by L and S, compare to obtain error sum of squares more than Minimum Residual;4) final damping function model is determined by the corresponding parameter value of error sum of squares more than Minimum Residual, so as to which zero time delay fitting autocorrelation value be calculated.Processing procedure of the present invention is simply easily realized, and remains the data information at zero time delay, compensates for deficiency of the existing incoherent scattering radar in terms of zero time delay auto-correlation function processing.
Description
Technical field
The invention belongs to technical field of information processing, more particularly to incoherent scattering radar estimate of autocorrelation algorithm and
Approximating method is damped, it is specially a kind of based on auto-correlation function processing method at the radar zero time delay for damping fitting, available for counting
Calculate the autocorrelation value at incoherent scattering radar zero time delay.
Background technology
Ionosphere as component part important near earth space environment, directly to weather monitoring, broadcast, radar fix,
The activities such as radionavigation have tremendous influence, therefore seem most important to the detection in ionosphere.Common ionospheric probing
Means have vertical sounding, oblique incidence souding, Coherent Scattering Radar detection and incoherent scattering radar detection etc..Wherein, it is incoherent to dissipate
Radar is penetrated with outstanding advantages of measurement parameter is more, covering spatial dimension is wide, spatial and temporal resolution is high, becomes current ground observation ionization
The most powerful means of layer.But since it is built, difficulty is big, operating cost is high, and only there are about ten incoherent scattering thunders in the whole world at present
It reaches.Head set incoherent scattering radars in China are tentatively built up at the beginning of 2012 in Qujing of Yunnan, the ionization to lower latitude in China
Sheaf space weather monitoring is of great significance with research.
Incoherent scattering radar is for a wide range of continuously distributed ionosphere target, signal processing method and conventional radar
There is very big difference.The echo-signal in ionosphere is mainly the incoherent scattering signal as caused by the scatterings such as electronics, ion, is
A kind of typical random signal, within the short time of a few minutes, Ionospheric Echo signal has stationarity, can be by calculating from phase
Function and power spectral density are closed to characterize the statistical property of signal;Ionosphere is typical soft object, and the radar of different height returns
The mutual aliasing of wave signal needs to eliminate range ambiguity by the design of effective Signal coding and special signal processing algorithm.
The purpose of incoherent scattering radar signal processing is the autocorrelation value or power spectrum obtained at different delay, is subsequent inverting work
It prepares.
In incoherent scattering signal processing, since the high blur at zero time delay is big, resolution ratio is low, lead to auto-correlation letter
Number harmonies at different delay are poor, influence estimation performance compose to ionospheric scattering, it is therefore desirable to oneself at zero time delay
Correlation function carries out extra process.Article " the The use of multipulse zero that Oulu, Finland university delivers at it
It is disclosed in lag data to improve incoherent scatter radar power profile accuracy "
A kind of method for obtaining zero time delay information, this method to data at the zero time delay that is generated in incoherent scattering signal processing into
Row discard processing, and the pulse either Barker code of different frequency is retransmitted to obtain the data information at zero time delay.It is aobvious
And be clear to, this method requirement incoherent scattering radar possesses the coding mode of pulse or Barker code.However, in Yunnan Province of China
The first set incoherent scattering radar that Qujing is built up now only exists long pulse coding and two kinds of coding modes of two-phase alternate coded, is obtaining
Message context at zero time delay is taken to have larger difficulty, the method that Oulun Yliopisto is proposed is not appropriate for the incoherent scattering thunder in China
It reaches.
Invention content
At in for the above-mentioned prior art to incoherent scattering radar zero time delay auto-correlation function
The deficiency of reason method proposes auto-correlation function processing method at a kind of radar zero time delay based on damping fitting, current to solve
Incoherent scattering radar is in ionization detection layer, it is impossible to accurate the problem of obtaining information at zero time delay.
Realizing the concrete thought of the object of the invention is:1) calculation delay one is obtained to the actual measurement auto-correlation function value of time delay eight
Autocorrelation value is surveyed to a combination thereof;2) damping function model, and the value range of each parameter and search step in setting model are established
It is long;3) according to damping function model, time delay one is obtained to the fitting auto-correlation function expression formula of time delay eight, in parameter value range
Different multiple remnants are calculated from fitting auto-correlation function expression formula by combination actual measurement autocorrelation value in interior multiple Selecting All Parameters
Error sum of squares compares to obtain error sum of squares more than Minimum Residual;4) according to the corresponding parameter value of error sum of squares more than Minimum Residual, really
Fixed final damping function model;So as to solve the auto-correlation function value at zero time delay.It realizes and the accurate of information at zero time delay is obtained
It takes.
The present invention realizes that above-mentioned purpose is as follows:
(1) combination actual measurement autocorrelation value is calculated:
(1.1) the ionospheric scattering raw radar data D in incoherent scattering radar receiver is obtained;
(1.2) raw radar data D is filtered using Gaussian filter, obtains filtering data D1;
(1.3) the second filtering data D2 is obtained after filtering data D1 being translated i symbol, is obtained according to the following formula at time delay i
Actual measurement autocorrelation value L (i):
L (i)=D1 × D2,
Wherein i=0,1,2 ..., m-1, m be transmitting signal element number and m >=9;
(1.4) time delay i is taken as 1 to 8, according to step (1.3) respectively calculation delay 1 to time delay 8 actual measurement autocorrelation value,
Obtain combination actual measurement autocorrelation value L:
L=[L (1), L (2) ..., L (8)];
(2) damping function model is established:
If ionospheric scattering spectrum is C (v):
Wherein v represents frequency, A1To scatter spectral amplitude ratio, ε is halfwidth of line, v0To compose centre frequency, Δ v is spectral peak frequency
The difference of rate and spectrum centre frequency;
Ionosphere auto-correlation function I (t), i.e. damping function model are obtained as inverse Fourier transform to C (v):
Wherein, IDFT represents inverse Fourier transform, and t represents time, A1To scatter spectral amplitude ratio, δ represents damped coefficient, ωr=
2 π/T represents damping vibration angular frequency, and T represents the damping vibration period;
According to A in combination actual measurement autocorrelation value setting damping function model1Parameter value range;According to ionospheric scattering
Spectrum is main to concentrate frequency separation setting ωrWith the parameter value range of δ;
(3) error sum of squares more than acquisition Minimum Residual:
(3.1) t=i is substituted into formula<2>, obtain fitting autocorrelation value q (i) expression formulas at time delay i:
Q (i)=I (i)=2A1e-δ|i|cos(ωri) <3>
(3.2) time delay i is taken as 1 to 8, obtains time delay 1 respectively according to step (3.1) to the fitting autocorrelation value of time delay 8
Expression formula obtains combination and fitting autocorrelation value S expression formulas:
S=[q (1), q (2) ..., q (8)]<4>
(3.3) according to the parameter A set in step (2)1、ωrWith δ value ranges, its step-size in search is set separately;
(3.4) in preset parameter A1、ωrIn δ value ranges, the minimum value of each parameter is chosen respectively as just
Beginning parameter value, according to formula<5>Acquire the first residual error quadratic sum;
Solve residual error quadratic sum σ2Formula it is as follows:
σ2=[S-L] [S-L]T <5>
Wherein, []TRepresent transposition operation;
To A1、ωrWith the initial parameter value of δ, it is incremented by respectively according to the step-size in search of setting;
(3.5) judge the parameter size after being incremented by:
If the parameter A after being incremented by1、ωrIts parameter value range maximum value is respectively less than with δ, then is entered step (3.6);
If the parameter A after being incremented by1、ωrIt is greater than or equal to its parameter value range maximum value with any one in δ, then enters
Step (3.7);
(3.6) according to formula<5>The residual error quadratic sum corresponding to the parameter after being incremented by is acquired, and to each after being incremented by
Parameter is incremented by again according to its step-size in search, is then back to step (3.5) and the parameter after being incremented by again is judged;
(3.7) compare all residual error quadratic sum sizes acquired, obtain error sum of squares σ more than Minimum Residualmin 2;
(4) zero time delay fitting autocorrelation value is calculated:
If error sum of squares σ more than Minimum Residualmin 2Corresponding parameter value is final scattering spectral amplitude ratio A'1, final damping vibration
Angular frequency 'rWith final damped coefficient δ ', formula is substituted into<2>Obtain final damping function expression formula:
I'(t)=2A'1e-δ|t|cos(ω'rT),
Zero time delay is calculated by above formula and is fitted autocorrelation value q (0):
Q (0)=I'(0)=2A'1。
The present invention has the following advantages that compared with prior art:
First, since the present invention does not need to carry out additional pulse coding or Barker code to incoherent scattering radar, from
And reduce the cost of Radar Design, while reduce radar access times, reduce operation and maintenance expense;
Second, since the processing procedure that the present invention uses is not related to abandoning data, remain the data at zero time delay
Information so that, ionosphere information more perfect to the signal processing of incoherent scattering radar obtains more fully;
Third is introduced into the signal processing of incoherent scattering radar since the present invention will damp to be fitted, passes through resistance
The mode of Buddhist nun's fitting is fitted time delay auto-correlation function, solves the computational problem of autocorrelation value at zero time delay, and existing
The processing method of auto-correlation function is compared, and processing procedure is more succinct and easy to implement.
Description of the drawings
Fig. 1 is the realization flow chart of the present invention;
Fig. 2 is fitting auto-correlation analogous diagram and actual measurement auto-correlation analogous diagram;
Fig. 3 is the auto-correlation analogous diagram handled with incoherent scattering signal processing method echo-signal;
Fig. 4 is the auto-correlation analogous diagram handled with the present invention echo-signal;
Fig. 5 is that the more height powers handled with incoherent scattering signal processing method echo-signal compose analogous diagram;
Fig. 6 is that the more height powers handled with the present invention echo-signal compose analogous diagram.
Specific embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, below in conjunction with the accompanying drawings and specific implementation
Example is described further the present invention.
With reference to Fig. 1, specific implementation step of the invention is as follows:
Step 1:Calculate combination actual measurement autocorrelation value:
The ionospheric scattering echo raw data acquired from incoherent scattering radar receiver is stored with document form,
Then the data in selecting file carry out incoherent scattering radar signal processing, it is intended to calculation delay 1 to time delay 8 to data
Combination actual measurement autocorrelation value, in later step damper model parameter determine, specific calculating process is as follows:.
(1.1) the ionospheric scattering raw radar data D in incoherent scattering radar receiver is obtained;
(1.2) raw radar data D is filtered using Gaussian filter, obtains filtering data D1;Wherein Gauss filters
The width of the shock response of wave device is equal to transmitting signal element width;
(1.3) the second filtering data D2 is obtained after filtering data D1 being translated i symbol, is obtained according to the following formula at time delay i
Actual measurement autocorrelation value L (i):
L (i)=D1 × D2,
Wherein i=0,1,2 ..., m-1, m be transmitting signal element number and m >=9;
(1.4) time delay i is taken as 1 to 8, according to step (1.3) respectively calculation delay 1 to time delay 8 actual measurement autocorrelation value,
Obtain combination actual measurement autocorrelation value L:
L=[L (1), L (2) ..., L (8)];
The processing method of incoherent scattering radar signal realizes that step can be divided into this step:Data are read, are selected box-like
Wave filter is filtered signal, to eliminate noise jamming, while determines auto-correlation time delay spacing;To filtered number
According to progress related operation, calculation delay section matrix;Scattered signal difference is calculated using time delay section matrix according to coding mode
Auto-correlation function in altitude range, and the auto-correlation function in different cycles is accumulated.Incoherent scattering radar is basic
Parameter is according to following setting:Tranmitting frequency 500MHz, frequency acquisition 6.25MHz, 16 two-phase alternate codeds, pulse width
480us, time delay spacing 30us, time delay number 16, fitting using time delay 1 to time delay 8, totally 8 time delay points.
Step 2:Establish damping function model:
If ionospheric scattering spectrum is C (v):
Wherein v represents frequency, A1To scatter spectral amplitude ratio, ε is halfwidth of line, v0To compose centre frequency, Δ v is spectral peak frequency
The difference of rate and spectrum centre frequency;
Ionosphere auto-correlation function I (t), i.e. damping function model are obtained as inverse Fourier transform to C (v):
Wherein, IDFT represents inverse Fourier transform, and t represents time, A1To scatter spectral amplitude ratio, δ represents damped coefficient, ωr=
2 π/T represents damping vibration angular frequency, and T represents the damping vibration period;
According to A in combination actual measurement autocorrelation value setting damping function model1Parameter value range:A1By combination actual measurement certainly
Correlation L=[L (1), L (2) ..., L (8)] is determined, only considers L (1) here, sets A1Value range be [0.3,2] L
(1);
Main concentration frequency separation setting ω is composed according to ionospheric scatteringrWith the parameter value range of δ:Parameter ωrValue
Ranging from [0,1.6];The value range of parameter δ is [0,0.79];
The exact value of all parameters is unknown in model, by the parameter area of setting, in fit procedure is damped below
It is determined.
Step 3:Error sum of squares more than acquisition Minimum Residual:
(3.1) t=i is substituted into formula<2>, obtain fitting autocorrelation value q (i) expression formulas at time delay i:
Q (i)=I (i)=2A1e-δ|i|cos(ωri) <3>
(3.2) time delay i is taken as 1 to 8, obtains time delay 1 respectively according to step (3.1) to the fitting autocorrelation value of time delay 8
Expression formula obtains combination and fitting autocorrelation value S expression formulas:
S=[q (1), q (2) ..., q (8)]<4>
(3.3) according to the parameter A set in step (2)1、ωrWith δ value ranges, A is set separately1Step-size in search be
0.05L (1), ωrStep-size in search is that 0.03, δ step-size in search is 0.03;
(3.4) in preset parameter A1、ωrIn δ value ranges, the minimum value of each parameter is chosen respectively as just
Beginning parameter, i.e.,:Take A1For 0.3L (1)=2.0809*109、ωrFor 0, δ 0, then the corresponding combination and fitting autocorrelation value of the parameter
S and combination actual measurement autocorrelation value L are respectively:
Solve residual error quadratic sum σ2Formula it is as follows:
σ2=[S-L] [S-L]T <5>
Wherein, []TRepresent transposition operation;
According to formula<5>The first residual error quadratic sum is acquired as 1.5*1020;
To A1、ωrWith the initial parameter value of δ, it is incremented by respectively according to the step-size in search of setting, obtains being incremented by scattering spectral amplitude ratio
A1b=2.4277*109, it is incremented by damping vibration angular frequencyrb=0.03, it is incremented by damped coefficient δb=0.03;
(3.5) judge the parameter size after being incremented by:
Parameter A after being incremented by1b< 1.3873*1010, ωrb< 1.6, δb< 0.79, i.e., each parameter are respectively less than its parameter and take
It is worth range maximum value, therefore, enters step (3.6);
(3.6) be incremented by after parameter corresponding to combination and fitting autocorrelation value S and combine survey autocorrelation value L be respectively:
According to formula<5>It is 1.4563*10 to acquire the residual error quadratic sum corresponding to the parameter after being incremented by20;
Each parameter after being incremented by according to its step-size in search is incremented by again, obtains being incremented by scattering spectral amplitude ratio A1b=2.7745*
109, it is incremented by damping vibration angular frequencyrb=0.06, it is incremented by damped coefficient δb=0.06, step (3.5) is then back to again
Parameter after being incremented by is judged that each parameter is still respectively less than its parameter value range maximum value, therefore, continues to solve current ginseng
The corresponding residual error quadratic sum of number;The process is recycled after multiple value, the parameter A after being incremented by1、ωrWith it is any one in δ
It is a to be greater than or equal to its parameter value range maximum value, then it enters step (3.7);
There are 27 by the residual error quadratic sum repeatedly solved;
(3.7) compare the 27 residual error quadratic sum sizes acquired, obtain error sum of squares σ more than Minimum Residualmin 2For:
σmin 2=3.8605*1018;
Step 4:Calculate zero time delay fitting autocorrelation value:
Error sum of squares σ more than Minimum Residualmin 2Corresponding parameter value is:Final scattering spectral amplitude ratio A1b=9.0171*109、
Final damping vibration angular frequency 'r=0.60 and final damped coefficient δ '=0.60, substituted into formula<2>Finally hindered
Buddhist nun's function expression:
I'(t)=1.8034*1010*e-0.6|t|Cos (0.6t),
Zero time delay is calculated by above formula and is fitted autocorrelation value q (0):
Q (0)=I'(0)=1.8034*1010。
The application effect of the present invention is further described with reference to following emulation:
First, simulated conditions:Under 7 environment of Windows, emulation experiment is carried out using software MATLAB.
2nd, emulation content and result:
Emulation 1 for the echo data in ionosphere being highly 402km, uses incoherent scattering radar signal processing side
Method calculation delay one and emulates it to the autocorrelation value (actual measurement autocorrelation value) of time delay eight using MATLAB softwares, imitates
True result is as shown in " circle " mark curve in Fig. 2.To the autocorrelation value of the time delay zero using damping function fitting to time delay eight
MATLAB emulation is carried out, simulation result is as shown in " asterisk " mark curve in Fig. 2.
From Figure 2 it can be seen that damping function curve trend and the auto-correlation function under the calculating of incoherent scattering signal processing method
Curve moves towards roughly the same.Although observing presence of the figure due to range ambiguity, the two is not completely superposed, its error is permitting
Perhaps in the range of.The simulation result shows completely correct using the thought of damping fitting auto-correlation function.In addition to this, abscissa is
The value of ordinate is the autocorrelation value at zero time delay when zero, it is seen that its calculating is extremely simple.
Emulation 2 for the echo data in ionosphere being highly 402km, uses incoherent scattering radar signal processing side
Method calculates the autocorrelation value at all time delays, and it is emulated using MATLAB softwares, simulation result such as Fig. 3.
As seen from Figure 3, auto-correlation function harmony at different delay is poor, this is because range ambiguity at zero time delay
Greatly, distance resolution is relatively low.The imbalance of auto-correlation function has seriously affected the detection in ionosphere, it is therefore necessary to zero time delay
The auto-correlation function at place carries out extra process.
Emulation 3 for the echo data in ionosphere being highly 402km, is corrected non-using method described in the invention
The zero time delay autocorrelation value that is obtained under Coherent Scattering Radar signal processing method, then with MATLAB softwares at all time delays from
Correlation is emulated, and simulation result is as shown in Figure 4.
From fig. 4, it can be seen that compared with Fig. 3, it is harmonious at different delay by the corrected auto-correlation function of damping function
It makes moderate progress, is conducive to the detection in ionosphere.
Emulation 4, for the echo data of height all in ionosphere, uses incoherent scattering radar signal processing method meter
The power spectrum chart of full-height is calculated, and it is emulated using MATLAB softwares, simulation result such as Fig. 5.
As seen from Figure 5, it with reference to Fig. 3, since autocorrelation value is far longer than the autocorrelation value at other time delays at zero time delay, leads
The dispersion of incoherent scattering spectrum energy is caused, secondary lobe is high, seriously affects the estimation of each physical parameter to ionosphere.
For the echo data of height all in ionosphere, method described in the invention is used in each height for emulation 6
The zero time delay autocorrelation value obtained under incoherent scattering radar signal processing method is corrected, makees Fourier transformation to it, then use
MATLAB softwares emulate transformation results, and simulation result is as shown in Figure 6.
As seen from Figure 6, compared with Fig. 5, by the corrected power spectrum of damping function, energy comparison is concentrated, and is conducive to electricity
The estimation of each physical parameter of absciss layer.
In summary all simulation results, using the present invention, when can effectively and easily calculate incoherent scattering radar zero
The autocorrelation value of Yan Chu, and the harmony of auto-correlation function is substantially improved, the estimation performance of incoherent scattering spectrum is improved,
Be conducive to the detection in ionosphere, fully demonstrate the feasibility of the present invention.
Unspecified part of the present invention belongs to common sense well known to those skilled in the art.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, it is clear that for this field
Professional for, all may be in the situation without departing substantially from the principle of the invention, structure after the content of present invention and principle has been understood
Under, various modifications and variations in form and details are carried out, but these modifications and variations based on inventive concept are still at this
Within the claims of invention.
Claims (4)
1. auto-correlation function processing method at a kind of radar zero time delay based on damping fitting, which is characterized in that including walking as follows
Suddenly:
(1) combination actual measurement autocorrelation value is calculated:
(1.1) the ionospheric scattering raw radar data D in incoherent scattering radar receiver is obtained;
(1.2) raw radar data D is filtered using Gaussian filter, obtains filtering data D1;
(1.3) the second filtering data D2 is obtained after filtering data D1 being translated i symbol, obtains the reality at time delay i according to the following formula
Survey autocorrelation value L (i):
L (i)=D1 × D2,
Wherein i=0,1,2 ..., m-1, m be transmitting signal element number and m >=9;
(1.4) time delay i is taken as 1 to 8, calculation delay 1 to the actual measurement autocorrelation value of time delay 8, obtains respectively according to step (1.3)
Combination actual measurement autocorrelation value L:
L=[L (1), L (2) ..., L (8)];
(2) damping function model is established:
If ionospheric scattering spectrum is C (v):
Wherein v represents frequency, A1To scatter spectral amplitude ratio, ε is halfwidth of line, v0For compose centre frequency, Δ v for spectrum peak frequency with
Compose the difference of centre frequency;
Ionosphere auto-correlation function I (t), i.e. damping function model are obtained as inverse Fourier transform to C (v):
Wherein, IDFT represents inverse Fourier transform, and t represents time, A1To scatter spectral amplitude ratio, δ represents damped coefficient, ωr=2 π/T
Represent damping vibration angular frequency, T represents the damping vibration period;
According to A in combination actual measurement autocorrelation value setting damping function model1Parameter value range;It is composed and led according to ionospheric scattering
Frequency separation is concentrated to set ωrWith the parameter value range of δ;
(3) error sum of squares more than acquisition Minimum Residual:
(3.1) t=i is substituted into formula<2>, obtain fitting autocorrelation value q (i) expression formulas at time delay i:
Q (i)=I (i)=2A1e-δ|i|cos(ωri) <3>
(3.2) time delay i is taken as 1 to 8, the fitting autocorrelation value for obtaining time delay 1 to time delay 8 respectively according to step (3.1) is expressed
Formula obtains combination and fitting autocorrelation value S expression formulas:
S=[q (1), q (2) ..., q (8)]<4>
(3.3) according to the parameter A set in step (2)1、ωrWith δ value ranges, its step-size in search is set separately;
(3.4) in preset parameter A1、ωrIn δ value ranges, the minimum value of each parameter is chosen respectively as initial ginseng
Numerical value, according to formula<5>Acquire the first residual error quadratic sum;
Solve residual error quadratic sum σ2Formula it is as follows:
σ2=[S-L] [S-L]T <5>
Wherein, []TRepresent transposition operation;
To A1、ωrWith the initial parameter value of δ, it is incremented by respectively according to the step-size in search of setting;
(3.5) judge the parameter size after being incremented by:
If the parameter A after being incremented by1、ωrIts parameter value range maximum value is respectively less than with δ, then is entered step (3.6);
If the parameter A after being incremented by1、ωrIt is greater than or equal to its parameter value range maximum value with any one in δ, then enters step
(3.7);
(3.6) according to formula<5>The residual error quadratic sum corresponding to the parameter after being incremented by is acquired, and to each parameter after being incremented by
It is incremented by again according to its step-size in search, is then back to step (3.5) and the parameter after being incremented by again is judged;
(3.7) compare all residual error quadratic sum sizes acquired, obtain error sum of squares σ more than Minimum Residualmin 2;
(4) zero time delay fitting autocorrelation value is calculated:
If error sum of squares σ more than Minimum Residualmin 2Corresponding parameter value is final scattering spectral amplitude ratio A'1, final damping vibration angular frequency
Rate ω 'rWith final damped coefficient δ ', formula is substituted into<2>Obtain final damping function expression formula:
I'(t)=2A'1e-δ|t|cos(ω'rT),
Zero time delay is calculated by above formula and is fitted autocorrelation value q (0):
Q (0)=I'(0)=2A'1。
2. according to the method described in claim 1, it is characterized in that:Incoherent scattering radar basic parameter in step (1) is set
It puts as follows:
Tranmitting frequency 500MHz, frequency acquisition 6.25MHz, 16 two-phase alternate codeds, pulse width 480us, time delay spacing
30us, time delay number 16, fitting using time delay 1 to time delay 8, totally 8 time delay points.
3. according to the method described in claim 1, it is characterized in that:Gaussian filter in step (1.2), shock response
Width is equal to transmitting signal element width.
4. according to the method described in claim 1, it is characterized in that:Damping function model parameter value range is set in step (3)
It is fixed as follows:
Parameter A1Value range for [0.3,2] L (1), step-size in search is 0.05L (1);Parameter ωrValue range for [0,
1.6], step-size in search 0.03;Parameter δ value ranges be [0,0.79], step-size in search 0.03.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810319869.0A CN108254730B (en) | 2018-04-11 | 2018-04-11 | Radar zero-time-delay autocorrelation function processing method based on damping fitting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810319869.0A CN108254730B (en) | 2018-04-11 | 2018-04-11 | Radar zero-time-delay autocorrelation function processing method based on damping fitting |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108254730A true CN108254730A (en) | 2018-07-06 |
CN108254730B CN108254730B (en) | 2021-07-20 |
Family
ID=62748146
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810319869.0A Active CN108254730B (en) | 2018-04-11 | 2018-04-11 | Radar zero-time-delay autocorrelation function processing method based on damping fitting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108254730B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090135045A1 (en) * | 2007-11-28 | 2009-05-28 | Camero-Tech Ltd. | Through-the-obstacle radar system and method of operation |
CN102253372A (en) * | 2011-04-14 | 2011-11-23 | 西安电子科技大学 | ISR (incoherent scatter radar) signal processing system based on alternate code modulation |
CN103913723A (en) * | 2014-03-18 | 2014-07-09 | 南昌大学 | Waveform design method for processing scattered signals of incoherent scatter radar |
CN104181509A (en) * | 2014-08-31 | 2014-12-03 | 西安电子科技大学 | Incoherent scattering radar signal processing method based on frequency hopping and polyphase alternating codes |
CN104185353A (en) * | 2014-09-05 | 2014-12-03 | 中国人民解放军陆军军官学院 | Fusion reactor plasma density and temperature diagnosing method based on Thomson scattering weak coherent technique |
CN105182304A (en) * | 2015-06-04 | 2015-12-23 | 南昌大学 | Ionized layer incoherent scattering radar incoherent scattering spectrum modeling method |
CN105242274A (en) * | 2015-10-26 | 2016-01-13 | 南昌大学 | Ionosphere non-coherent scatter radar differential phase detection method |
CN107632323A (en) * | 2017-08-08 | 2018-01-26 | 西安电子科技大学 | Radar echo signal processing method for Ionospheric F detection |
-
2018
- 2018-04-11 CN CN201810319869.0A patent/CN108254730B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090135045A1 (en) * | 2007-11-28 | 2009-05-28 | Camero-Tech Ltd. | Through-the-obstacle radar system and method of operation |
CN102253372A (en) * | 2011-04-14 | 2011-11-23 | 西安电子科技大学 | ISR (incoherent scatter radar) signal processing system based on alternate code modulation |
CN103913723A (en) * | 2014-03-18 | 2014-07-09 | 南昌大学 | Waveform design method for processing scattered signals of incoherent scatter radar |
CN104181509A (en) * | 2014-08-31 | 2014-12-03 | 西安电子科技大学 | Incoherent scattering radar signal processing method based on frequency hopping and polyphase alternating codes |
CN104185353A (en) * | 2014-09-05 | 2014-12-03 | 中国人民解放军陆军军官学院 | Fusion reactor plasma density and temperature diagnosing method based on Thomson scattering weak coherent technique |
CN105182304A (en) * | 2015-06-04 | 2015-12-23 | 南昌大学 | Ionized layer incoherent scattering radar incoherent scattering spectrum modeling method |
CN105242274A (en) * | 2015-10-26 | 2016-01-13 | 南昌大学 | Ionosphere non-coherent scatter radar differential phase detection method |
CN107632323A (en) * | 2017-08-08 | 2018-01-26 | 西安电子科技大学 | Radar echo signal processing method for Ionospheric F detection |
Non-Patent Citations (6)
Title |
---|
LIN LI 等: "A modified parameter model-based spectrum estimation method for incoherent scatter radar", 《2017 2ND INTERNATIONAL CONFERENCE ON IMAGE, VISION AND COMPUTING(ICIVC)》 * |
ZIQING LU 等: "An effective method for incoherent scattering radar"s detecting ability evaluation", 《20161212 RADIO SCIENCE》 * |
吴永宏 等: "非相干散射信号数值模拟", 《电波科学学报》 * |
吴永宏 等: "非相干散射的一种非线性反演方法", 《电波科学学报》 * |
梅俊: "非相干散射雷达信号处理***外场测试实验与算法改进", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
王玮: "非相干散射雷达信号处理算法研究与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN108254730B (en) | 2021-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Carey et al. | Ocean ambient noise: measurement and theory | |
US7336222B2 (en) | System and method for measuring characteristics of a continuous medium and/or localized targets using multiple sensors | |
AU2008307362B2 (en) | Ground clutter mitigation using a parametric time domain method | |
CN106872969B (en) | Radar target angle estimation method based on MTD pulse accumulation and sliding processing | |
Cao et al. | Scale transition from geostrophic motions to internal waves in the northern South China Sea | |
US7239983B2 (en) | System and method for measuring characteristics of continuous media and/or localized targets using at least one multi-frequency sensor | |
CN105259537B (en) | Doppler spectral center frequency estimation method based on frequency displacement iteration | |
Bitner-Gregersen | Uncertainties in data for the offshore environment | |
Zhao et al. | Cramer-Rao lower bounds for the joint delay-Doppler estimation of an extended target | |
CN103412287A (en) | Linear frequency modulation signal parameter evaluation method based on LVD (Lv's distribution) | |
Park et al. | Simulation and analysis of polarimetric radar signatures of human gaits | |
CN106772302A (en) | A kind of knowledge assistance STAP detection methods under complex Gaussian background | |
CN108254730A (en) | Auto-correlation function processing method at radar zero time delay based on damping fitting | |
Yu et al. | Measurements of midfrequency acoustic backscattering from a sandy bottom in the South Yellow Sea of China | |
Rosenberg et al. | Comparison of bi-modal coherent sea clutter models | |
Goncharenko et al. | Adaptive moving target indication in a windblown clutter environment | |
Altes | Some theoretical concepts for echolocation | |
Oyedokun | Sea clutter simulation | |
Zabotin et al. | Studies of wave activity in the thermosphere-ionosphere system using Dynasonde techniques | |
Wang et al. | High-resolution spectral estimation of HF radar data for current measurement applications | |
RU2708072C1 (en) | Method of determining transverse dimensions of a target based on data of two spaced positions in multi-position radar stations based on their polarization-scattering properties | |
Choudhury et al. | Computer simulation of weather radar signals | |
Leonidov | Acoustic sounding in the problem of finding and registration of thermodynamic perturbations in the surface layer of atmosphere | |
Wang et al. | Measurement of surface currents from a shipborne HFSWR | |
VASSILEVA et al. | SIMULATION TOOL FOR ROC ANALYSIS IN THE CONTEXT OF WEATHER RADAR CALIBRATION |
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 |