CN102445685A - Small spot radar signal decomposition method - Google Patents

Small spot radar signal decomposition method Download PDF

Info

Publication number
CN102445685A
CN102445685A CN2011102983263A CN201110298326A CN102445685A CN 102445685 A CN102445685 A CN 102445685A CN 2011102983263 A CN2011102983263 A CN 2011102983263A CN 201110298326 A CN201110298326 A CN 201110298326A CN 102445685 A CN102445685 A CN 102445685A
Authority
CN
China
Prior art keywords
parameter
waveform
radar signal
value
initial value
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
Application number
CN2011102983263A
Other languages
Chinese (zh)
Other versions
CN102445685B (en
Inventor
赖旭东
侯文广
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN 201110298326 priority Critical patent/CN102445685B/en
Publication of CN102445685A publication Critical patent/CN102445685A/en
Application granted granted Critical
Publication of CN102445685B publication Critical patent/CN102445685B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention relates to a small spot radar signal decomposition method. In the method, original radar signal waveform is taken as superposition of Gaussian waves, so as to estimate a Gaussian decomposition parameter initial value and optimize the value. The process of estimating the Gaussian decomposition initial value comprises the following steps of: extracting three parameters of one waveform with a maximum peak value in the current waveforms; subtracting corresponding component of the waveform from the original radar signal waveform; judging whether a termination condition effects, if so, exiting, otherwise, continuing to decompose the signal. The method can effectively decompose hybrid waveforms, particularly for the condition that the peaks are very closer; the method can accurately estimate the number of waveforms and the parameter initial value; in the method, pretreatments, such as filtering and the like, is not needed for the initial data; and the termination condition is easy to control. With the result of decomposition by the method, a more accurate and detailed digital ground model can be provided.

Description

Small light spot radar signal decomposition method
Technical field
The present invention relates to field of information processing, particularly a kind of decomposition method of small light spot radar signal.
Background technology
On the theoretical foundation that laser wave graphic data Gaussian function decomposes, Liu Feng etc. have proposed to utilize Generalized Gaussian pattern function match pulse waveform, and propose the method (Central-South Forestry University science and technology journal 30 is rolled up in August, 2,010 8 phases) of pulse waveform important parameter.Small light spot laser radar (LiDAR) data can provide more accurately detailed digital terrain model.To the small light spot laser radar data, horse flood is superfine then to have proposed a kind of position that improved EM pulse detection method obtains echo-pulse and higher waveform decomposition method (2009 1 phases of remote sensing journal) of dependable performance, precision of width of using.
These technology have all realized the decomposition of LiDAR Wave data preferably; Also obtained certain parameter, the location is inaccurate but these technology exist initial value, the also inaccurate problem of wavelet shape number estimation; So just cause decomposition result often not restrain, can produce error.
Summary of the invention
To the deficiency of prior art, the present invention proposes a kind of small light spot radar signal decomposition method.
Technical scheme of the present invention is a kind of small light spot radar signal decomposition method; Original radar signal waveform
Figure 728111DEST_PATH_IMAGE001
is considered as a plurality of high bass waves ,
Figure 588751DEST_PATH_IMAGE003
... The stack of
Figure 2011102983263100002DEST_PATH_IMAGE004
Figure 287059DEST_PATH_IMAGE005
， =1，2，… ；
Wherein,
Figure 2011102983263100002DEST_PATH_IMAGE008
expression waveform point reflection interval position; The truth of a matter of expression natural logarithm; The amplitude of parameter
Figure 2011102983263100002DEST_PATH_IMAGE010
expression high bass wave
Figure 690118DEST_PATH_IMAGE011
; The intermediate value of parameter
Figure 2011102983263100002DEST_PATH_IMAGE012
expression high bass wave
Figure 46144DEST_PATH_IMAGE011
, the variance of parameter
Figure 400902DEST_PATH_IMAGE013
expression high bass wave
Figure 444950DEST_PATH_IMAGE011
;
By iterating on Gaussian
Figure 139237DEST_PATH_IMAGE011
number
Figure 853115DEST_PATH_IMAGE007
, and each Gaussian parameter
Figure 275262DEST_PATH_IMAGE010
,
Figure 483521DEST_PATH_IMAGE012
and
Figure 94631DEST_PATH_IMAGE013
for estimation, parameter
Figure 494257DEST_PATH_IMAGE010
, and
Figure 212388DEST_PATH_IMAGE013
initial value, and then optimized to obtain initial parameters
Figure 678004DEST_PATH_IMAGE010
,
Figure 749996DEST_PATH_IMAGE012
and
Figure 566643DEST_PATH_IMAGE013
determine the value, according to a Gaussian
Figure 615239DEST_PATH_IMAGE011
number
Figure 13990DEST_PATH_IMAGE007
and parameters
Figure 506152DEST_PATH_IMAGE010
,
Figure 784993DEST_PATH_IMAGE012
and
Figure 388013DEST_PATH_IMAGE013
The determined value decomposition raw radar signal waveform
Figure 906850DEST_PATH_IMAGE001
,
Figure 2011102983263100002DEST_PATH_IMAGE014
Estimate that Gauss decomposes initial value and may further comprise the steps,
Step 1; If
Figure 428967DEST_PATH_IMAGE006
=1, initial current waveform are original radar signal waveform ;
Step 2; Extract the maximum point of peak value in the current waveform; The horizontal ordinate of this point
Figure 613534DEST_PATH_IMAGE015
is as the initial value of parameter
Figure 173829DEST_PATH_IMAGE012
; Peak value is as the initial value of parameter
Figure 555262DEST_PATH_IMAGE010
; Search peak is the horizontal ordinate of 0.5 * point respectively along this peak value left and right sides; Obtain
Figure 2011102983263100002DEST_PATH_IMAGE016
and ; Get the littler person of
Figure 193420DEST_PATH_IMAGE016
and
Figure 198285DEST_PATH_IMAGE017
middle distance and be designated as
Figure 2011102983263100002DEST_PATH_IMAGE018
, then the initial value of parameter is (
Figure 719287DEST_PATH_IMAGE018
-
Figure 567157DEST_PATH_IMAGE012
)/0.8326;
Step 3; Initial value according to step 2 gained parameter
Figure 554705DEST_PATH_IMAGE010
, and
Figure 417673DEST_PATH_IMAGE013
; From current waveform, deduct the corresponding component of high bass wave
Figure 764341DEST_PATH_IMAGE005
, obtain a new Wave data;
Step 4; Judge whether to satisfy end condition; Satisfy end condition and then make n= ; Withdraw from then; Otherwise make
Figure 752337DEST_PATH_IMAGE006
=
Figure 996236DEST_PATH_IMAGE006
+1; The Wave data that step 3 gained is new is as current waveform, repeating step 2 and step 3; The combination of one or more that said end condition is the following option,
(1) peak value of high bass wave
Figure 766003DEST_PATH_IMAGE011
is less than the peak value resolution of radar equipment;
(2) Gaussian
Figure 728142DEST_PATH_IMAGE011
of the area with raw radar signal waveform
Figure 294253DEST_PATH_IMAGE001
of the area is less than the preset value than the threshold area;
(3) Gaussian
Figure 143391DEST_PATH_IMAGE011
The peak Gaussian
Figure 831862DEST_PATH_IMAGE002
The peak of the peak ratio is less than a preset threshold value;
(4) the high bass wave number that has decomposited
Figure 218981DEST_PATH_IMAGE006
reaches the preset decomposition waveform number upper limit.
And, the initial value of parameter , and
Figure 650334DEST_PATH_IMAGE013
is optimized adopts memory-limited BFGS algorithm or expectation-maximization algorithm or non-linear least square algorithm.
The present invention can effectively decompose hybrid waveform, particularly corrugation pitch from very near situation; Can accurately estimate waveform number and parameter initial value; Need not carry out pre-service such as filtering to raw data; End condition is controlled easily.Decompose the result who obtains like this, more accurately detailed digital terrain model can be provided.
Description of drawings
Fig. 1 is the process flow diagram of the embodiment of the invention.
Fig. 2 is the oscillogram after the embodiment of the invention iteration first time.
Fig. 3 is the oscillogram after the embodiment of the invention iteration second time.
Fig. 4 is the embodiment of the invention oscillogram after the iteration for the third time.
Fig. 5 is the oscillogram after the 4th iteration of the embodiment of the invention.
Fig. 6 is the oscillogram after the 5th iteration of the embodiment of the invention.
Fig. 7 is the oscillogram after the embodiment of the invention parameter optimization.
Embodiment
Specify technical scheme of the present invention below in conjunction with accompanying drawing and embodiment.
Embodiment is considered as a plurality of high bass waves
Figure 245711DEST_PATH_IMAGE002
,
Figure 318710DEST_PATH_IMAGE003
with original radar signal waveform
Figure 72219DEST_PATH_IMAGE001
... The stack of
Figure 612898DEST_PATH_IMAGE004
，
Figure 30421DEST_PATH_IMAGE006
=1，2，…
Figure 161189DEST_PATH_IMAGE007
；
Wherein, expression waveform point reflection interval position; The truth of a matter of
Figure 945791DEST_PATH_IMAGE009
expression natural logarithm; The amplitude of parameter
Figure 507091DEST_PATH_IMAGE010
expression high bass wave
Figure 554681DEST_PATH_IMAGE011
; The intermediate value of parameter
Figure 864440DEST_PATH_IMAGE012
expression high bass wave
Figure 748214DEST_PATH_IMAGE011
, the variance of parameter
Figure 801620DEST_PATH_IMAGE013
expression high bass wave
Figure 703717DEST_PATH_IMAGE011
;
By iterating on Gaussian
Figure 498891DEST_PATH_IMAGE011
number
Figure 853649DEST_PATH_IMAGE007
, and each Gaussian parameter
Figure 586113DEST_PATH_IMAGE010
,
Figure 591984DEST_PATH_IMAGE012
and
Figure 243545DEST_PATH_IMAGE013
for estimation, parameter ,
Figure 559437DEST_PATH_IMAGE012
and
Figure 170547DEST_PATH_IMAGE013
initial value, and then optimized to obtain initial parameters
Figure 573103DEST_PATH_IMAGE010
,
Figure 840136DEST_PATH_IMAGE012
and determine the value, according to a Gaussian
Figure 317702DEST_PATH_IMAGE011
number
Figure 638962DEST_PATH_IMAGE007
and parameters
Figure 314663DEST_PATH_IMAGE010
,
Figure 379571DEST_PATH_IMAGE012
and The determined value decomposition raw radar signal waveform ,
Figure 498727DEST_PATH_IMAGE019
。
The Gauss who estimates is decomposed initial value to be optimized and can to adopt L-BFGS (memory-limited BFGS algorithm) or prior aries such as EM (expectation maximization) or LM (non-linear least square) algorithm.
Referring to Fig. 1, among the embodiment, estimate that Gauss decomposes initial value and may further comprise the steps,
Step 1; If
Figure 101747DEST_PATH_IMAGE006
=1, initial current waveform are original radar signal waveform
Figure 479638DEST_PATH_IMAGE001
;
Step 2; Extract the maximum point of peak value in the current waveform; The horizontal ordinate of this point
Figure 391968DEST_PATH_IMAGE015
is as the initial value of parameter
Figure 917628DEST_PATH_IMAGE012
; Peak value is as the initial value of parameter
Figure 75071DEST_PATH_IMAGE010
; Search peak is the horizontal ordinate of 0.5 *
Figure 938160DEST_PATH_IMAGE010
point respectively along this peak value left and right sides; Obtain
Figure 506545DEST_PATH_IMAGE016
and
Figure 722763DEST_PATH_IMAGE017
; Get the littler person of
Figure 418317DEST_PATH_IMAGE016
and
Figure 770801DEST_PATH_IMAGE017
middle distance
Figure 510087DEST_PATH_IMAGE015
and be designated as , then the initial value of parameter
Figure 945802DEST_PATH_IMAGE013
is (
Figure 215110DEST_PATH_IMAGE018
-
Figure 141609DEST_PATH_IMAGE012
)/0.8326;
Step 3; Initial value according to step 2 gained parameter
Figure 129156DEST_PATH_IMAGE010
,
Figure 353464DEST_PATH_IMAGE012
and
Figure 995054DEST_PATH_IMAGE013
; From current waveform, deduct the corresponding component of high bass wave
Figure 341722DEST_PATH_IMAGE005
, obtain a new Wave data;
Step 4; Judge whether to satisfy end condition; Satisfy end condition and then make n=
Figure 301719DEST_PATH_IMAGE006
; Withdraw from then; Otherwise make
Figure 329718DEST_PATH_IMAGE006
=
Figure 573617DEST_PATH_IMAGE006
+1; The Wave data that step 3 gained is new is as current waveform, repeating step 2 and step 3; The combination of one or more that said end condition is the following option,
(1) peak value of high bass wave
Figure 28870DEST_PATH_IMAGE011
; (being the initial value of step 2 gained parameter
Figure 240277DEST_PATH_IMAGE010
) is less than the peak value resolution of radar equipment;
(2) Gaussian
Figure 868704DEST_PATH_IMAGE011
of the area with raw radar signal waveform
Figure 904793DEST_PATH_IMAGE001
of the area is less than the preset value than the threshold area;
(3) Gaussian The peak Gaussian
Figure 842367DEST_PATH_IMAGE002
The peak of the peak ratio is less than a preset threshold value;
(4) the high bass wave number that has decomposited
Figure 274485DEST_PATH_IMAGE006
reaches the preset decomposition waveform number upper limit.
All high bass waves that decomposited before at high bass wave
Figure 430660DEST_PATH_IMAGE011
are designated as
Figure 40764DEST_PATH_IMAGE002
,
Figure 649600DEST_PATH_IMAGE003
...
Figure 2011102983263100002DEST_PATH_IMAGE020
, all high bass waves that decomposited
Figure 931415DEST_PATH_IMAGE002
,
Figure 942096DEST_PATH_IMAGE003
... The peak value of
Figure 972369DEST_PATH_IMAGE020
is respectively parameter
Figure 881550DEST_PATH_IMAGE021
,
Figure 2011102983263100002DEST_PATH_IMAGE022
... The initial value of
Figure 655471DEST_PATH_IMAGE023
.Because each target of decomposing all is the highest waveform of contained peak value in the current waveform, the initial value of parameter is at parameter
Figure 301926DEST_PATH_IMAGE021
,
Figure 619775DEST_PATH_IMAGE022
... Maximum in the initial value of
Figure 948120DEST_PATH_IMAGE023
.The initial value that utilizes parameter
Figure 995710DEST_PATH_IMAGE021
is as denominator, and the present invention has designed end condition (3) and can supply to select for use.The ratio of the peak value of the peak value of high bass wave
Figure 617053DEST_PATH_IMAGE011
and high bass wave
Figure 422198DEST_PATH_IMAGE002
stops during less than preset peakedness ratio threshold value, and the ratio of the peak value of the peak value of high bass wave and high bass wave
Figure 190751DEST_PATH_IMAGE002
is the continuation iteration during more than or equal to preset peakedness ratio threshold value.
During practical implementation, the peak value resolution of radar equipment is decided according to the radar equipment that produces original radar signal, and area can be set up on their own by those skilled in the art than threshold value, peakedness ratio threshold value and the decomposition waveform number upper limit as the case may be.For example area gets 0.02 than threshold value, and the peakedness ratio threshold value gets 0.1, decomposes the waveform number upper limit and gets 5.During practical implementation, one of above-mentioned four kinds of conditions be can select separately for use, two kinds or above combination also can be selected for use.Can be set in when satisfying any one condition and stop, also can be set in when satisfying two kinds or above condition simultaneously and stop.
The end condition of embodiment is made as the high bass wave number
Figure 733728DEST_PATH_IMAGE006
that has decomposited and reaches the preset decomposition waveform number upper limit, decomposes the waveform number upper limit and gets 5.For ease of implement with reference to for the purpose of, it is following to provide embodiment to estimate that Gauss decomposes the concrete iterative process of initial value:
Iteration for the first time: operating procedure 1 and step 2, step 3; Obtain the parameter initial value
Figure 271816DEST_PATH_IMAGE021
=52 of the maximum high bass wave
Figure 617981DEST_PATH_IMAGE002
of peak value;
Figure 2011102983263100002DEST_PATH_IMAGE024
=25;
Figure 44731DEST_PATH_IMAGE025
=7.207207; The result is as shown in Figure 2; Curve 1 is original radar signal waveform
Figure 758609DEST_PATH_IMAGE001
among the figure, curve 2 first high bass wave
Figure 849931DEST_PATH_IMAGE002
for decompositing.Curve 3 is the current waveform (this time for the raw radar signal waveform
Figure 245140DEST_PATH_IMAGE001
) minus
Figure 121829DEST_PATH_IMAGE002
results.
Iteration for the second time: operating procedure 4; Do not run into end condition; Current waveform is the result that original radar signal waveform
Figure 757341DEST_PATH_IMAGE001
deducts
Figure 24375DEST_PATH_IMAGE002
; Continue step 2,3; Obtain the parameter initial value
Figure 268984DEST_PATH_IMAGE022
=32.316383 of second high bass wave
Figure 285592DEST_PATH_IMAGE003
;
Figure 2011102983263100002DEST_PATH_IMAGE026
=40;
Figure 324665DEST_PATH_IMAGE027
=6.006006; Decomposition result is as shown in Figure 3; Curve 1 is original radar signal waveform among the figure, and curve 2 is first high bass wave
Figure 894635DEST_PATH_IMAGE002
that decomposites and the stack of second high bass wave
Figure 480337DEST_PATH_IMAGE003
.The result that curve 3 deducts
Figure 644602DEST_PATH_IMAGE002
for current waveform, the result of stack that promptly original radar signal waveform
Figure 197812DEST_PATH_IMAGE001
deducts
Figure 800832DEST_PATH_IMAGE002
with
Figure 178724DEST_PATH_IMAGE003
.
Iteration for the third time: operating procedure 4; Do not run into end condition, current waveform is that original radar signal waveform
Figure 592518DEST_PATH_IMAGE001
deducts the result that adds
Figure 462571DEST_PATH_IMAGE003
.Continue step 2,3; Obtain the parameter initial value
Figure 9484DEST_PATH_IMAGE029
=8.255985 of the 3rd Gaussian waveform
Figure 2011102983263100002DEST_PATH_IMAGE028
;
Figure 2011102983263100002DEST_PATH_IMAGE030
=35;
Figure 656497DEST_PATH_IMAGE031
=2.402402; Decomposition result is as shown in Figure 4; Curve 1 is original radar signal waveform
Figure 669452DEST_PATH_IMAGE001
among the figure, and curve 2 is the high bass wave
Figure 863542DEST_PATH_IMAGE002
that decomposites, the stack of
Figure 216026DEST_PATH_IMAGE003
and
Figure 220891DEST_PATH_IMAGE028
.The result that curve 3 deducts
Figure 471875DEST_PATH_IMAGE028
for current waveform, promptly original radar signal waveform deducts
Figure 427379DEST_PATH_IMAGE002
, the result of stack with
Figure 780573DEST_PATH_IMAGE028
.
The 4th iteration: operating procedure 4; Do not run into end condition, current waveform is that original radar signal waveform
Figure 67198DEST_PATH_IMAGE001
deducts , the result of
Figure 226095DEST_PATH_IMAGE003
stack with
Figure 700939DEST_PATH_IMAGE028
.Continue step 2 and step 3; Obtain the parameter initial value
Figure 40522DEST_PATH_IMAGE033
=6.577729 of the 4th Gaussian waveform
Figure 2011102983263100002DEST_PATH_IMAGE032
;
Figure 2011102983263100002DEST_PATH_IMAGE034
=49;
Figure 97471DEST_PATH_IMAGE035
=2.402402; Decomposition result is as shown in Figure 5; Curve 1 is original radar signal waveform among the figure, and curve 2 is the stack of the high bass wave that decomposites,
Figure 118189DEST_PATH_IMAGE003
, and
Figure 108328DEST_PATH_IMAGE032
.The result that curve 3 deducts
Figure 42917DEST_PATH_IMAGE032
for current waveform, promptly original radar signal waveform
Figure 475035DEST_PATH_IMAGE001
deducts the result of ,
Figure 739849DEST_PATH_IMAGE003
, stack with
Figure 318915DEST_PATH_IMAGE032
.
The 5th iteration: operating procedure 4; Do not run into end condition, current waveform is the result that original radar signal waveform
Figure 142646DEST_PATH_IMAGE001
deducts , , stack with .Continue step 2 and step 3; Obtain the parameter initial value
Figure 548482DEST_PATH_IMAGE037
=4.772379 of the 5th Gaussian waveform ;
Figure 2011102983263100002DEST_PATH_IMAGE038
=54; =3.603604; Decomposition result is as shown in Figure 6; Curve 1 is an original waveform among the figure, and curve 2 is for decompositing the stack of high bass wave
Figure 522571DEST_PATH_IMAGE002
,
Figure 304582DEST_PATH_IMAGE003
,
Figure 928855DEST_PATH_IMAGE028
,
Figure 999579DEST_PATH_IMAGE032
and
Figure 928352DEST_PATH_IMAGE036
.The result that curve 3 deducts
Figure 830449DEST_PATH_IMAGE036
for current waveform, promptly original radar signal waveform
Figure 622693DEST_PATH_IMAGE001
deducts the result of
Figure 243031DEST_PATH_IMAGE002
,
Figure 647598DEST_PATH_IMAGE003
,
Figure 341885DEST_PATH_IMAGE028
,
Figure 321342DEST_PATH_IMAGE032
and Die Jia.
Operating procedure 4 runs into end condition (decomposition come out waveform number 5 reach the predetermined upper limit), accomplishes parameter estimation.
Embodiment uses the result of L-BFGS algorithm completion parameter optimization as shown in Figure 7.L-BFGS belongs to the plan Newton's algorithm, and it finds the solution the inverse matrix of the matrix of second derivatives of objective function to be optimized through alternative manner, and then the direction of search of definite optimum solution, along direction of search search optimum solution.Parameter optimization result is following:
Figure 810803DEST_PATH_IMAGE021
=11.518133, =19.407736,
Figure 323004DEST_PATH_IMAGE025
=3.848433,
Figure 590037DEST_PATH_IMAGE022
=31.053619, =38.696826,
Figure 520133DEST_PATH_IMAGE027
=7.753454,
Figure 825081DEST_PATH_IMAGE029
=47.163967,
Figure 641727DEST_PATH_IMAGE030
=25.237025,
Figure 378739DEST_PATH_IMAGE031
=5.939188
Figure 715174DEST_PATH_IMAGE033
=1.584873,
Figure 207335DEST_PATH_IMAGE034
=51.384519,
Figure 448960DEST_PATH_IMAGE035
=3.262579
=2.075186,
Figure 756105DEST_PATH_IMAGE038
=47.053530,
Figure 904320DEST_PATH_IMAGE039
=33.409674
Curve 1 is original radar signal waveform
Figure 633242DEST_PATH_IMAGE001
among Fig. 7; Curve 2 is the curve of decomposition result
Figure 774373DEST_PATH_IMAGE014
; Curve 3 is a graph of errors, expresses the poor of original radar signal waveform and decomposition result.
Accompanying drawing 2～7th, the match oscillogram, the overlapping more expression precision of curve is high more.
Specific embodiment described herein only is that the present invention's spirit is illustrated.Person of ordinary skill in the field of the present invention can make various modifications or replenishes or adopt similar mode to substitute described specific embodiment, but can't depart from spirit of the present invention or surmount the defined scope of appended claims.

Claims (2)

1. small light spot radar signal decomposition method; It is characterized in that: original radar signal waveform
Figure 518810DEST_PATH_IMAGE001
is considered as a plurality of high bass waves
Figure 2011102983263100001DEST_PATH_IMAGE002
,
Figure 498268DEST_PATH_IMAGE003
... The stack of
Figure 2011102983263100001DEST_PATH_IMAGE004
Figure 91054DEST_PATH_IMAGE005
，
Figure 2011102983263100001DEST_PATH_IMAGE006
=1，2，…
Figure 548580DEST_PATH_IMAGE007
；
Wherein,
Figure 2011102983263100001DEST_PATH_IMAGE008
expression waveform point reflection interval position; The truth of a matter of expression natural logarithm; The amplitude of parameter expression high bass wave
Figure 387014DEST_PATH_IMAGE011
; The intermediate value of parameter
Figure 2011102983263100001DEST_PATH_IMAGE012
expression high bass wave , the variance of parameter
Figure 226849DEST_PATH_IMAGE013
expression high bass wave ;
By iterating on Gaussian
Figure 951409DEST_PATH_IMAGE011
number
Figure 705738DEST_PATH_IMAGE007
, and each Gaussian parameter
Figure 255799DEST_PATH_IMAGE010
, and for estimation, parameter ,
Figure 164926DEST_PATH_IMAGE012
and
Figure 621446DEST_PATH_IMAGE013
initial value, then the initial value for the optimized parameters
Figure 956612DEST_PATH_IMAGE010
, and
Figure 826665DEST_PATH_IMAGE013
The determined value, according to the Gaussian
Figure 698544DEST_PATH_IMAGE011
number
Figure 283240DEST_PATH_IMAGE007
and the parameter
Figure 233879DEST_PATH_IMAGE010
,
Figure 178701DEST_PATH_IMAGE012
and
Figure 845699DEST_PATH_IMAGE013
The determined value decomposition of raw radar signal waveform
Figure 850564DEST_PATH_IMAGE001
,
Figure 2011102983263100001DEST_PATH_IMAGE014
Estimate that Gauss decomposes initial value and may further comprise the steps,
Step 1; If
Figure 101548DEST_PATH_IMAGE006
=1, initial current waveform are original radar signal waveform
Figure 522165DEST_PATH_IMAGE001
;
Step 2; Extract the maximum point of peak value in the current waveform; The horizontal ordinate of this point
Figure 57052DEST_PATH_IMAGE015
is as the initial value of parameter
Figure 170501DEST_PATH_IMAGE012
; Peak value is as the initial value of parameter
Figure 407316DEST_PATH_IMAGE010
; Search peak is the horizontal ordinate of 0.5 *
Figure 631624DEST_PATH_IMAGE010
point respectively along this peak value left and right sides; Obtain
Figure 2011102983263100001DEST_PATH_IMAGE016
and ; Get the littler person of
Figure 915155DEST_PATH_IMAGE016
and
Figure 642196DEST_PATH_IMAGE017
middle distance and be designated as
Figure 2011102983263100001DEST_PATH_IMAGE018
, then the initial value of parameter
Figure 914094DEST_PATH_IMAGE013
is (
Figure 916817DEST_PATH_IMAGE018
-
Figure 878956DEST_PATH_IMAGE012
)/0.8326;
Step 3; Initial value according to step 2 gained parameter
Figure 445067DEST_PATH_IMAGE010
,
Figure 58320DEST_PATH_IMAGE012
and ; From current waveform, deduct the corresponding component of high bass wave
Figure 868330DEST_PATH_IMAGE005
, obtain a new Wave data;
Step 4; Judge whether to satisfy end condition; Satisfy end condition and then make n=
Figure 51181DEST_PATH_IMAGE006
; Withdraw from then; Otherwise make
Figure 941776DEST_PATH_IMAGE006
=
Figure 801148DEST_PATH_IMAGE006
+1; The Wave data that step 3 gained is new is as current waveform, repeating step 2 and step 3; The combination of one or more that said end condition is the following option,
(1) peak value of high bass wave
Figure 712779DEST_PATH_IMAGE011
is less than the peak value resolution of radar equipment;
(2) Gaussian to the area of the original radar signal waveform
Figure 959270DEST_PATH_IMAGE001
The area is less than a preset threshold area ratio;
(3) Gaussian
Figure 927226DEST_PATH_IMAGE011
The peak Gaussian
Figure 836407DEST_PATH_IMAGE002
The peak of the peak ratio is less than a preset threshold value;
(4) the high bass wave number that has decomposited
Figure 610328DEST_PATH_IMAGE006
reaches the preset decomposition waveform number upper limit.
2. small light spot radar signal decomposition method as claimed in claim 1 is characterized in that: the initial value to parameter
Figure 787100DEST_PATH_IMAGE010
,
Figure 925958DEST_PATH_IMAGE012
and
Figure 571703DEST_PATH_IMAGE013
is optimized employing memory-limited BFGS algorithm or expectation-maximization algorithm or non-linear least square algorithm.
CN 201110298326 2011-09-28 2011-09-28 Small spot radar signal decomposition method Expired - Fee Related CN102445685B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110298326 CN102445685B (en) 2011-09-28 2011-09-28 Small spot radar signal decomposition method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110298326 CN102445685B (en) 2011-09-28 2011-09-28 Small spot radar signal decomposition method

Publications (2)

Publication Number Publication Date
CN102445685A true CN102445685A (en) 2012-05-09
CN102445685B CN102445685B (en) 2013-05-08

Family

ID=46008368

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110298326 Expired - Fee Related CN102445685B (en) 2011-09-28 2011-09-28 Small spot radar signal decomposition method

Country Status (1)

Country Link
CN (1) CN102445685B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105044702A (en) * 2015-09-18 2015-11-11 宁波华仪宁创智能科技有限公司 Fitting method for pulse waveforms
CN104408018B (en) * 2014-11-19 2018-02-09 武汉大学 A kind of conformal LiDAR waveforms denoising method and system
CN111077532A (en) * 2019-11-22 2020-04-28 同济大学 Surface feature space information acquisition method based on deconvolution and Gaussian decomposition

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2914432A1 (en) * 2007-03-27 2008-10-03 Thales Sa Stationary signal i.e. synthetic aperture radar video signal, extrapolating method for remote or Doppler high resolution analysis of signal, involves performing spectral decomposition of sub-signals and concatenation of spectrums
CN101833087A (en) * 2010-05-13 2010-09-15 王成 Satellite-bone laser radar waveform data resolving method based on wavelet analysis
CN101923158A (en) * 2009-06-08 2010-12-22 霍尼韦尔国际公司 Be used for Gauss's decomposition of weather radar data for communication for the system and method that transmits
US8032319B1 (en) * 2007-06-29 2011-10-04 Hal Enterprises, LLC Methods for analyzing streaming composite waveforms

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2914432A1 (en) * 2007-03-27 2008-10-03 Thales Sa Stationary signal i.e. synthetic aperture radar video signal, extrapolating method for remote or Doppler high resolution analysis of signal, involves performing spectral decomposition of sub-signals and concatenation of spectrums
US8032319B1 (en) * 2007-06-29 2011-10-04 Hal Enterprises, LLC Methods for analyzing streaming composite waveforms
CN101923158A (en) * 2009-06-08 2010-12-22 霍尼韦尔国际公司 Be used for Gauss's decomposition of weather radar data for communication for the system and method that transmits
CN101833087A (en) * 2010-05-13 2010-09-15 王成 Satellite-bone laser radar waveform data resolving method based on wavelet analysis

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
黄涛 等: "激光雷达反遮蔽探测技术研究", 《光电工程》, vol. 38, no. 3, 31 March 2011 (2011-03-31) *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104408018B (en) * 2014-11-19 2018-02-09 武汉大学 A kind of conformal LiDAR waveforms denoising method and system
CN105044702A (en) * 2015-09-18 2015-11-11 宁波华仪宁创智能科技有限公司 Fitting method for pulse waveforms
CN111077532A (en) * 2019-11-22 2020-04-28 同济大学 Surface feature space information acquisition method based on deconvolution and Gaussian decomposition

Also Published As

Publication number Publication date
CN102445685B (en) 2013-05-08

Similar Documents

Publication Publication Date Title
Wang et al. Wind power forecasting based on singular spectrum analysis and a new hybrid Laguerre neural network
Hu et al. High resolution wind speed forecasting based on wavelet decomposed phase space reconstruction and self-organizing map
CN102967854B (en) Multi-fractal detection method of targets in FRFT (Fractional Fourier Transformation) domain sea clutter
CN102608620B (en) Laser scanning point cloud vegetation filtering method on basis of reflection strength and terrain
CN103729637A (en) Extended target probability hypothesis density filtering method based on cubature Kalman filtering
CN107942308A (en) Adaptive Rao detection methods based on gamma texture under complex Gaussian environment
CN104331583B (en) A kind of Multifractal Modeling method based on Observed sea clutter
CN103927556A (en) ECG signal classifying method based on wavelet packet and approximate entropy
CN102445685A (en) Small spot radar signal decomposition method
CN105527617A (en) Ground penetrating radar data background removing method based on robust principal component analysis
CN108804736A (en) A kind of method and apparatus multiple degrees of freedom FM signal design and optimized
CN105676205A (en) Airborne LiDAR waveform data Gaussian decomposition method
CN104268567A (en) Extended target tracking method using observation data clustering and dividing
Alexopoulos et al. Fractional order Pareto distributions with application to X‐band maritime radar clutter
CN110954860B (en) DOA and polarization parameter estimation method
CN104215939A (en) Knowledge assisted space-time adaptive processing method integrating generalized symmetrical structure information
CN106019257A (en) Interpolation method based on high-frequency ground-wave radar ocean current observation result space-time characteristics
CN101739718A (en) Parameter template-based corn leaf virtual simulation modeling method
CN103308829A (en) GIS (gas insulated switchgear) single-time partial discharge signal extraction and trigger moment adjusting method
CN106125039B (en) Improvement space-time adaptive Monopulse estimation method based on local Combined Treatment
CN103049653A (en) G0 distribution parameter maximum likelihood estimation method based on EM (Expectation-Maximization) algorithm
CN103218807A (en) Multispectral remote-sensing image change detection method
CN110658517B (en) Dimensionality reduction sparse STAP method and device based on uncertain priori knowledge
DE102020107429A1 (en) Determination of the temporal position of a signal
CN110516806A (en) The rarefaction method and apparatus of neural network parameter matrix

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130508

Termination date: 20190928

CF01 Termination of patent right due to non-payment of annual fee