CN113640763A - Estimation method of lognormal texture sea clutter amplitude distribution shape parameters based on fractional order moment - Google Patents

Estimation method of lognormal texture sea clutter amplitude distribution shape parameters based on fractional order moment Download PDF

Info

Publication number
CN113640763A
CN113640763A CN202110865345.3A CN202110865345A CN113640763A CN 113640763 A CN113640763 A CN 113640763A CN 202110865345 A CN202110865345 A CN 202110865345A CN 113640763 A CN113640763 A CN 113640763A
Authority
CN
China
Prior art keywords
sea clutter
lognormal
texture
moment
amplitude
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
CN202110865345.3A
Other languages
Chinese (zh)
Other versions
CN113640763B (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.)
Xian University of Posts and Telecommunications
Original Assignee
Xian University of Posts and Telecommunications
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 Xian University of Posts and Telecommunications filed Critical Xian University of Posts and Telecommunications
Priority to CN202110865345.3A priority Critical patent/CN113640763B/en
Publication of CN113640763A publication Critical patent/CN113640763A/en
Application granted granted Critical
Publication of CN113640763B publication Critical patent/CN113640763B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details 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/418Theoretical aspects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details 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/414Discriminating targets with respect to background clutter
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

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 relates to the technical field of radar signal processing, in particular to a fractional order moment-based method for estimating the shape parameter of amplitude distribution of lognormal texture sea clutter, which adopts the negative fractional order moment and the positive fractional order moment of the amplitude of a sea clutter sample, reduces the calculation complexity and improves the estimation precision of the shape parameter compared with the existing estimation method based on zlogz and the estimation method based on the high order moment, and further obtains the logarithm normal texture sea clutter amplitude distribution shape parameter estimation value by the method
Figure DDA0003187347600000011
The method is more accurate in the value of the detection statistic in the radar target detection method obtained by the radar target detection method, so that the detection performance of the radar target is further improved.

Description

Estimation method of lognormal texture sea clutter amplitude distribution shape parameters based on fractional order moment
Technical Field
The invention relates to the technical field of radar signal processing, in particular to a method for estimating distribution shape parameters of lognormal texture sea clutter amplitude based on fractional order moment.
Background
The sea surveillance radar inevitably receives, in operation, backscattered signals of electromagnetic waves transmitted by the sea surface to the radar, which are commonly referred to as sea clutter. The radar target detection method in the sea clutter background is usually related to the shape parameter of the sea clutter amplitude distribution, so in order to better detect the radar target, the shape parameter of the sea clutter amplitude distribution must be accurately estimated. When the radar resolution unit length is much larger than the wavelength of the swell or the ground-scraping angle observed by the radar is larger than 10 degrees, a Gaussian model can be used for modeling the sea clutter. However, as the radar resolution increases or the observed ground angle decreases, the sea clutter will appear significantly non-gaussian, and the gaussian model will be severely mismatched. At present, a model capable of effectively modeling non-Gaussian sea clutter is a composite Gaussian model, the model is based on a sea clutter physical formation mechanism, and theoretical verification and experimental inspection are carried out. The complex gaussian model describes the sea clutter by using the product of the slowly varying texture component and the rapidly varying speckle component. The speckle component of sea clutter is an electromagnetic wave speckle component caused by small-scale capillary waves on the sea surface, which is a complex gaussian random variable obeying zero-mean unit power. The texture component of the sea clutter is generated by large-scale gravity waves of the sea clutter, and is a random variable. The amplitude characteristic of the composite gaussian sea clutter is determined by the probability density function of the texture component. Under the composite Gaussian model, according to the probability density function of the texture component and the speckle component, the amplitude distribution of the sea clutter can be obtained based on a total probability formula. The texture component of the sea clutter is modeled as a random variable obeying gamma distribution, and K distribution which is widely used in the field of radar target detection and used for describing the amplitude of the sea clutter can be obtained. However, as radar resolution increases, the K-distribution has not been able to accurately describe non-gaussian sea clutter. To more accurately model the amplitude of non-Gaussian sea clutter, a lognormal texture complex Gaussian distribution is proposed. The lognormal texture composite gaussian distribution models the texture component of the sea clutter using lognormal distribution, which is more suitable for describing the amplitude characteristics of non-gaussian sea clutter than K distribution.
The shape parameters contained in the lognormal texture complex gaussian distribution control the non-gaussian nature of sea clutter: the larger the shape parameter is, the more serious the non-gaussian property of the sea clutter is; the smaller the shape parameter, the weaker the non-gaussian nature of the sea clutter. The shape parameter occurs in a radar target detection method under a lognormal texture sea clutter background, so the shape parameter must be estimated by using sea clutter data received by a radar before detecting a target. At present, the estimation methods of the lognormal texture sea clutter amplitude distribution shape parameters include a high-order moment estimation method and a zlogz method, and the estimation methods are shown in documents of I.Chalabi and A.Mezache, Estimators of compound Gaussian distributor with log-normal texture, Remote Sensing Letters, vol.10, No.7, pp.709-716, Jul.2019, doi:10.1080/2150704X.2019.1601275. However, when the sea clutter is not gaussian (i.e. when the shape parameter is large), the estimation error of the amplitude distribution shape parameter estimation method of the lognormal texture sea clutter based on zlogz is too large, which affects the detection performance of the radar target under the background of the lognormal texture sea clutter.
Disclosure of Invention
The invention aims to solve the technical problem that the existing lognormal texture sea clutter amplitude distribution shape parameter estimation method is insufficient in precision, and provides a shape parameter estimation method of lognormal texture sea clutter amplitude distribution based on fractional order moment aiming at the defects in the prior art so as to improve the estimation precision of shape parameters in the lognormal texture sea clutter amplitude distribution.
The technical problem solved by the invention can be realized by adopting the following technical scheme:
the method for estimating the distribution shape parameters of the amplitude of the sea clutter based on the lognormal texture of the fractional order moment comprises the following steps
S1: calculating n-order theoretical moment m of lognormal texture sea clutter amplitude distributionnWherein m represents the theoretical moment of the sea clutter amplitude distribution, and n represents the order of the sea clutter amplitude distribution moment;
s2: according to the n-order theoretical moment m obtained in S1nDeducing a mathematical equation of a lognormal texture sea clutter amplitude distribution shape parameter gamma based on the fractional order moment;
s3: selecting N sea clutter amplitude samples from sea clutter data received by a radar, and respectively calculating statistical moments of the N sea clutter amplitude samples based on fractional order;
s4: substituting the fractional order statistical moment obtained in the step S3 into the mathematical equation in the step S2 to obtain the lognormal texture sea clutter based on the fractional order momentEstimation of amplitude profile shape parameters
Figure BDA0003187347580000031
S5: according to the obtained lognormal texture sea clutter amplitude distribution shape parameter estimation value based on fractional order moment
Figure BDA0003187347580000032
The parameter estimation value is applied to a radar target detection method under a lognormal texture sea clutter background, and test statistics of the radar target detection method is calculated and used for improving detection performance of the radar target.
Further, the mathematical equation of the shape parameter γ of the distribution of the amplitude of the sea clutter based on the lognormal texture of the fractional order moment in S2 is obtained by derivation
Figure BDA0003187347580000033
And
Figure BDA0003187347580000034
and (3) a mathematical equation of a shape parameter gamma of the amplitude distribution of the time lognormal texture sea clutter.
Further, the step of calculating the fractional order-based statistical moments of the N sea clutter amplitude samples in S3 is to calculate the N sea clutter amplitude samples based on
Figure BDA0003187347580000035
Order sum
Figure BDA0003187347580000036
Statistical moments of order
Figure BDA0003187347580000037
And
Figure BDA0003187347580000038
further, the nth order theoretical moment m of the log-normal texture sea clutter amplitude distribution is calculated in S1nThe method also comprises the following steps:
s01: calculating a lognormal texture sea clutter amplitude distribution probability density function f (r) according to the following formula:
Figure BDA0003187347580000041
wherein r represents the amplitude of the sea clutter, τ represents the texture component of the sea clutter, γ represents the shape parameter of the amplitude distribution of the lognormal texture sea clutter, δ represents the scale parameter of the amplitude distribution of the lognormal texture sea clutter, and ln (·) represents a natural logarithm function; e.g. of the type(·)Representing a natural exponential function;
s02: calculating the n-order theoretical moment m of the lognormal texture sea clutter amplitude distribution according to the lognormal texture sea clutter amplitude distribution probability density function f (r) obtained in the S01nThe calculation formula is as follows:
Figure BDA0003187347580000042
where E (-) represents the statistical averaging and Γ (-) represents the gamma function.
Further, the derivation in S2 is based on fractional order moment
Figure BDA0003187347580000043
And
Figure BDA0003187347580000044
the mathematical equation of the shape parameter gamma of the distribution of the amplitude of the time lognormal texture sea clutter comprises the following methods:
s001: n-order theoretical moment m of sea clutter amplitude distribution using lognormal texturenTo obtain
Figure BDA0003187347580000045
The equation for the order moment is as follows:
Figure BDA0003187347580000046
s002: n-order theoretical moment m of sea clutter amplitude distribution using lognormal texturenTo obtain
Figure BDA0003187347580000047
The equation for the order moment is as follows:
Figure BDA0003187347580000048
s003: and solving a mathematical equation of a lognormal texture sea clutter amplitude distribution shape parameter gamma based on the fractional order moment by using the equations in S001 and S002, wherein the mathematical equation comprises the following steps:
Figure BDA0003187347580000051
further, the N sea clutter amplitude samples in S3 are r1,r2,...,rNRespectively calculating N sea clutter amplitude samples r1,r2,...,rNIs/are as follows
Figure BDA0003187347580000052
Order sum
Figure BDA0003187347580000053
Statistical moments of order
Figure BDA0003187347580000054
And
Figure BDA0003187347580000055
wherein N represents the number of sea clutter samples, and r represents the amplitude of the sea clutter; the calculation formula is as follows:
Figure BDA0003187347580000056
further, in the step S4, the order is
Figure BDA0003187347580000057
Is taken as
Figure BDA0003187347580000058
Is taken as
Figure BDA0003187347580000059
Will be provided with
Figure BDA00031873475800000510
And
Figure BDA00031873475800000511
substituting into the mathematical equation in S2 to obtain the estimation value of the distribution shape parameter of the amplitude of the lognormal texture sea clutter based on the fractional order moment
Figure BDA00031873475800000512
The formula is as follows:
Figure BDA00031873475800000513
the invention has the beneficial effects that:
compared with the prior art, the method for estimating the shape parameter of the lognormal texture sea clutter amplitude distribution based on the fractional order moment adopts the negative fractional order moment and the positive fractional order moment of the sea clutter sample amplitude, reduces the calculation complexity and improves the estimation precision of the shape parameter compared with the existing estimation method based on zlogz and the estimation method based on the high order moment, and further obtains the shape parameter estimation value of the lognormal texture sea clutter amplitude distribution by the method
Figure BDA00031873475800000514
The method is more accurate in the value of the detection statistic in the radar target detection method obtained by the radar target detection method, so that the detection performance of the radar target is further improved.
Drawings
The invention is further illustrated with reference to the following figures and examples.
Fig. 1 is a schematic flow chart of the implementation of the invention.
FIG. 2 is a diagram of a relative root mean square error curve of the lognormal texture sea clutter amplitude distribution shape parameter estimation obtained by the present invention and the existing method.
Detailed Description
In the following, a scheme of the method for estimating the shape parameter of the sea clutter amplitude distribution based on the lognormal texture with fractional order moments according to the present invention will be described in detail through several specific embodiments.
It should be noted that the same symbols are used to represent the same meanings in all the formulas of the present invention.
Wherein m represents the theoretical moment of the sea clutter amplitude distribution;
n represents the order of the sea clutter amplitude distribution moment;
r represents the amplitude of the sea clutter;
τ represents the texture component of the sea clutter;
gamma represents the shape parameter of the amplitude distribution of the lognormal texture sea clutter;
δ represents a scale parameter of the lognormal texture sea clutter amplitude distribution.
Referring to fig. 1, the invention relates to a method for estimating the distribution shape parameters of the amplitude of the lognormal texture sea clutter based on fractional order moment, which comprises the following steps
S1: calculating n-order theoretical moment m of lognormal texture sea clutter amplitude distributionn
Calculating n-order theoretical moment m of lognormal texture sea clutter amplitude distributionnThe method specifically comprises the following steps:
s01: calculating a lognormal texture sea clutter amplitude distribution probability density function f (r) according to a lognormal distribution probability density function obeyed by a texture component tau of the sea clutter and a conditional Rayleigh distribution probability density function obeyed by a sea clutter amplitude r, wherein the calculation formula is as follows:
Figure BDA0003187347580000071
wherein ln (·) represents a natural logarithmic function; e.g. of the type(·)Representing a natural exponential function;
s02: calculating the n-order theoretical moment m of the lognormal texture sea clutter amplitude distribution according to the lognormal texture sea clutter amplitude distribution probability density function f (r) obtained in the S01nThe calculation formula is as follows:
Figure BDA0003187347580000072
wherein, E (-) represents the statistical average, and Γ (-) represents the gamma function;
s2: according to the n-order theoretical moment m obtained in S1nDeducing a mathematical equation of a lognormal texture sea clutter amplitude distribution shape parameter gamma based on the fractional order moment;
further, in the step S2, the mathematical equation of the lognormal texture sea clutter amplitude distribution shape parameter γ based on the fractional order moment is obtained by derivation
Figure BDA0003187347580000073
And
Figure BDA0003187347580000074
a mathematical equation of a time lognormal texture sea clutter amplitude distribution shape parameter gamma;
further, derived to obtain
Figure BDA0003187347580000075
And
Figure BDA0003187347580000076
the mathematical equation of the shape parameter gamma of the distribution of the amplitude of the time lognormal texture sea clutter comprises the following methods:
s001: n-order theoretical moment m of sea clutter amplitude distribution using lognormal texturenTo obtain
Figure BDA0003187347580000077
The equation for the order moment is as follows:
Figure BDA0003187347580000078
s002: n-order theoretical moment m of sea clutter amplitude distribution using lognormal texturenTo obtain
Figure BDA0003187347580000079
The equation for the order moment is as follows:
Figure BDA0003187347580000081
s003: and solving a mathematical equation of a lognormal texture sea clutter amplitude distribution shape parameter gamma based on the fractional order moment by using the equations in S001 and S002, wherein the mathematical equation comprises the following steps:
Figure BDA0003187347580000082
s3: selecting N sea clutter amplitude samples from sea clutter data received by a radar, and respectively calculating statistical moments of the N sea clutter amplitude samples based on fractional order;
further, the N sea clutter amplitude samples in S3 are r1,r2,...,rN
Furthermore, N sea clutter amplitude samples r are respectively selected from the sea clutter data received by the radar1,r2,...,rNThe acquisition method comprises the following steps:
a1, representing echo data received after a radar transmits pulse signals to the sea surface as a three-dimensional matrix Z, wherein Z is a P multiplied by L multiplied by Q three-dimensional matrix, P represents the azimuth number of the echo data matrix, L represents the distance unit number of the echo data matrix, and Q represents the pulse number of the echo data matrix;
a2 at echo numberSelecting a sea clutter data area G from a matrix Z, wherein G is a three-dimensional matrix P multiplied by L multiplied by Q, P is more than or equal to 1 and less than or equal to P, and L is more than or equal to 1 and less than or equal to L, and the amplitude of the sea clutter data contained in G is represented as r1,r2,...,rN,N=p×l×Q;
A3 calculating N sea clutter amplitude samples r respectively1,r2,...,rNIs/are as follows
Figure BDA0003187347580000083
Order sum
Figure BDA0003187347580000084
Statistical moments of order
Figure BDA0003187347580000085
And
Figure BDA0003187347580000086
wherein N represents the number of sea clutter samples, and r represents the amplitude of the sea clutter;
the formula is as follows:
Figure BDA0003187347580000091
s4: substituting the fractional order statistical moment obtained in the step S3 into the mathematical equation in the step S2 to obtain an estimated value of the log-normal texture sea clutter amplitude distribution shape parameter based on the fractional order moment
Figure BDA0003187347580000092
Further, in S4 is order
Figure BDA0003187347580000093
Is taken as
Figure BDA0003187347580000094
Is taken as
Figure BDA0003187347580000095
Will be provided with
Figure BDA0003187347580000096
And
Figure BDA0003187347580000097
substituting into the mathematical equation of the shape parameter gamma in S2, namely formula 5, to obtain the estimation value of the shape parameter of the lognormal texture sea clutter amplitude distribution based on the fractional order moment
Figure BDA0003187347580000098
The formula is as follows:
Figure BDA0003187347580000099
s5: according to the obtained lognormal texture sea clutter amplitude distribution shape parameter estimation value based on fractional order moment
Figure BDA00031873475800000910
The parameter estimation value is applied to a radar target detection method under a lognormal texture sea clutter background, and test statistics of the radar target detection method is calculated and used for improving detection performance of the radar target.
The radar target detection method is as follows:
Figure BDA00031873475800000911
where alpha represents the test statistic for the radar target detection method, M represents the number of accumulated coherent pulses,
Figure BDA00031873475800000912
representing an estimate of the texture of sea clutter in the absence of targets,
Figure BDA00031873475800000913
an estimate representing the texture of the sea clutter in the presence of a target, q0Indicating the power after whitening of the radar echo in the absence of a target,
Figure BDA00031873475800000914
representing the power of the whitened sea clutter under the condition of a target, wherein gamma represents the shape parameter of the amplitude distribution of the lognormal texture sea clutter, and delta represents the scale parameter of the amplitude distribution of the lognormal texture sea clutter;
the obtained lognormal texture sea clutter amplitude distribution shape parameter estimation value based on the fractional order moment
Figure BDA00031873475800000915
Introducing the equation to obtain test statistic of radar target detection method, and obtaining lognormal texture sea clutter amplitude distribution shape parameter estimation value obtained by the method
Figure BDA0003187347580000101
The method is more accurate in test statistic obtained by the method applied to radar target detection, and therefore the detection performance of the radar target is further improved.
The effect of the present invention will be further explained with the simulation experiment.
1. Simulation parameters
And simulating lognormal texture sea clutter amplitude samples by utilizing matlab software, namely selecting N sea clutter amplitude samples, setting experimental simulation parameters as sample number N being 10000, setting scale parameter delta being 1, increasing shape parameter gamma from 0.1 to 10, and increasing the step length to 0.1. The number of independent simulation experiments under each shape parameter was set to 10000. The Error of the shape parameter estimation is estimated by using a Relative Root Mean Square Error (RRMSE) with a calculation formula of RRMSE
Figure BDA0003187347580000102
2. Content of simulation experiment
In a simulation experiment, a high-order moment estimation method, a zlogz estimation method and the estimation method are used for estimating the shape parameters of the sea clutter data respectively, and an RRMSE curve graph of an estimation result is drawn.
The experimental result is shown in fig. 2, wherein the horizontal axis represents the true shape parameter γ of the sea clutter data, and the vertical axis represents the RRMSE corresponding to the estimation result. In fig. 2, the "- -" marked curve represents the RRMSE curve corresponding to the high-order moment estimation method, the "- - -" marked curve represents the RRMSE curve corresponding to the zlogz estimation method, and the "-" marked curve represents the RRMSE curve corresponding to the invention.
It can be seen from fig. 2 that the RRMSE of the present invention is smaller than the high-order moment estimation method and zlogz estimation method when the shape parameter is greater than 2.5. The larger the shape parameter is, the stronger the non-gaussian property of the sea clutter is, and the result of fig. 2 shows that the estimation error of the shape parameter of the invention is smaller than that of the existing method under the non-gaussian sea clutter environment.
In summary, the invention provides a lognormal texture sea clutter amplitude distribution shape parameter estimation method based on fractional order moment, and the method has high estimation precision on shape parameters under the non-Gaussian sea clutter background.
Therefore, the lognormal texture sea clutter amplitude distribution shape parameter estimation value with high estimation precision obtained by the invention
Figure BDA0003187347580000111
The method is applied to radar target detection under the background of the lognormal texture sea clutter, and the detection performance of the radar target is improved.
While the embodiments of the present invention have been described in detail with reference to the drawings, the present invention is not limited to the above embodiments, and various changes can be made without departing from the spirit of the present invention within the knowledge of those skilled in the art, and the scope of the present invention is within the scope of the claims.
Technical solutions between various embodiments may be combined with each other, but must be realized by those skilled in the art, and when the technical solutions are contradictory or cannot be realized, such a combination should not be considered to exist, and is not within the protection scope of the present invention.

Claims (7)

1. The method for estimating the distribution shape parameters of the lognormal texture sea clutter amplitude based on the fractional order moment is characterized by comprising the following steps of: comprises the following steps
S1: calculating n-order theoretical moment m of lognormal texture sea clutter amplitude distributionnWherein m represents the theoretical moment of the sea clutter amplitude distribution, and n represents the order of the sea clutter amplitude distribution moment;
s2: according to the n-order theoretical moment m obtained in S1nDeducing a mathematical equation of a lognormal texture sea clutter amplitude distribution shape parameter gamma based on the fractional order moment;
s3: selecting N sea clutter amplitude samples from sea clutter data received by a radar, and respectively calculating statistical moments of the N sea clutter amplitude samples based on fractional order;
s4: substituting the fractional order statistical moment obtained in the step S3 into the mathematical equation in the step S2 to obtain an estimated value of the log-normal texture sea clutter amplitude distribution shape parameter based on the fractional order moment
Figure FDA0003187347570000011
S5: according to the obtained lognormal texture sea clutter amplitude distribution shape parameter estimation value based on fractional order moment
Figure FDA0003187347570000012
The parameter estimation value is applied to a radar target detection method under a lognormal texture sea clutter background, and test statistics of the radar target detection method is calculated and used for improving detection performance of the radar target.
2. The method for estimating the shape parameter of the distribution of the lognormal texture sea clutter amplitude based on the fractional order moment as claimed in claim 1, wherein: the mathematical equation of the lognormal texture sea clutter amplitude distribution shape parameter gamma in the S2 based on the fractional order moment is obtained by derivation
Figure FDA0003187347570000013
And
Figure FDA0003187347570000014
and (3) a mathematical equation of a shape parameter gamma of the amplitude distribution of the time lognormal texture sea clutter.
3. The method for estimating the shape parameter of the distribution of the lognormal texture sea clutter amplitude based on the fractional order moment as claimed in claim 2, wherein: the step of S3 is that the fractional order-based statistical moments of the N sea clutter amplitude samples are respectively calculated
Figure FDA0003187347570000015
Order sum
Figure FDA0003187347570000016
Statistical moments of order
Figure FDA0003187347570000017
And
Figure FDA0003187347570000018
4. the method for estimating the shape parameter of the distribution of the lognormal texture sea clutter amplitude based on the fractional order moment as claimed in claim 1, wherein: in S1, an nth order theoretical moment m of the logarithmic normal texture sea clutter amplitude distribution is calculatednThe method also comprises the following steps:
s01: calculating a lognormal texture sea clutter amplitude distribution probability density function f (r) according to the following formula:
Figure FDA0003187347570000021
wherein r represents the amplitude of the sea clutter, τ represents the texture component of the sea clutter, γ represents the shape parameter of the amplitude distribution of the lognormal texture sea clutter, δ represents the scale parameter of the amplitude distribution of the lognormal texture sea clutter, and ln (·) represents a natural logarithm function; e.g. of the type(·)Representing natural exponential functions;
S02: calculating the n-order theoretical moment m of the lognormal texture sea clutter amplitude distribution according to the lognormal texture sea clutter amplitude distribution probability density function f (r) obtained in the S01nThe calculation formula is as follows:
Figure FDA0003187347570000022
where E (-) represents the statistical averaging and Γ (-) represents the gamma function.
5. The method for estimating the shape parameter of the distribution of the lognormal texture sea clutter amplitude based on the fractional order moment as claimed in claim 2, wherein: the derivation of S2 is based on fractional order moment
Figure FDA0003187347570000023
And
Figure FDA0003187347570000024
the mathematical equation of the shape parameter gamma of the distribution of the amplitude of the time lognormal texture sea clutter comprises the following methods:
s001: n-order theoretical moment m of sea clutter amplitude distribution using lognormal texturenTo obtain
Figure FDA0003187347570000025
The equation for the order moment is as follows:
Figure FDA0003187347570000026
s002: n-order theoretical moment m of sea clutter amplitude distribution using lognormal texturenTo obtain
Figure FDA0003187347570000027
The equation for the order moment is as follows:
Figure FDA0003187347570000031
s003: and solving a mathematical equation of a lognormal texture sea clutter amplitude distribution shape parameter gamma based on the fractional order moment by using the equations in S001 and S002, wherein the mathematical equation comprises the following steps:
Figure FDA0003187347570000032
6. the method according to claim 3, wherein the method comprises: the N sea clutter amplitude samples in the S3 are r1,r2,...,rNRespectively calculating N sea clutter amplitude samples r1,r2,...,rNIs/are as follows
Figure FDA0003187347570000033
Order sum
Figure FDA0003187347570000034
Statistical moments of order
Figure FDA0003187347570000035
And
Figure FDA0003187347570000036
wherein N represents the number of sea clutter samples, and r represents the amplitude of the sea clutter; the calculation formula is as follows:
Figure FDA0003187347570000037
7. the method according to claim 3, wherein the fractional order moment-based method for estimating the shape parameter of the amplitude distribution of the lognormal texture sea clutter is characterized in thatCharacterized in that: in the step S4, the order is
Figure FDA0003187347570000038
Is taken as
Figure FDA0003187347570000039
Is taken as
Figure FDA00031873475700000310
Will be provided with
Figure FDA00031873475700000311
And
Figure FDA00031873475700000312
substituting into the mathematical equation in S2 to obtain the estimation value of the distribution shape parameter of the amplitude of the lognormal texture sea clutter based on the fractional order moment
Figure FDA00031873475700000313
The formula is as follows:
Figure FDA00031873475700000314
CN202110865345.3A 2021-07-29 2021-07-29 Method for estimating distribution shape parameters of sea clutter amplitude of lognormal texture based on fractional order moment Active CN113640763B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110865345.3A CN113640763B (en) 2021-07-29 2021-07-29 Method for estimating distribution shape parameters of sea clutter amplitude of lognormal texture based on fractional order moment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110865345.3A CN113640763B (en) 2021-07-29 2021-07-29 Method for estimating distribution shape parameters of sea clutter amplitude of lognormal texture based on fractional order moment

Publications (2)

Publication Number Publication Date
CN113640763A true CN113640763A (en) 2021-11-12
CN113640763B CN113640763B (en) 2023-08-04

Family

ID=78418900

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110865345.3A Active CN113640763B (en) 2021-07-29 2021-07-29 Method for estimating distribution shape parameters of sea clutter amplitude of lognormal texture based on fractional order moment

Country Status (1)

Country Link
CN (1) CN113640763B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6040797A (en) * 1997-09-02 2000-03-21 Daimlerchrysler Ag Method for detecting a radar target
US20060072816A1 (en) * 2004-09-23 2006-04-06 Mitsubishi Denki Kabushiki Kaisha Determining statistical descriptors of a signal from a set of its samples
US20100073218A1 (en) * 2008-09-24 2010-03-25 Lockheed Martin Corpration Method and apparatus for radar surveillance and detection of sea targets
US20120154208A1 (en) * 2010-12-17 2012-06-21 Raytheon Company Methods and apparatus for automatic stc from sea state measurement via radar sea clutter eccentricity
US20170139038A1 (en) * 2015-11-13 2017-05-18 Rohde & Schwarz Gmbh & Co. Kg Method and signal generator for simulation of sea clutter
CN107167781A (en) * 2017-03-31 2017-09-15 西安电子科技大学 The quantile method of estimation of sea clutter amplitude lognormal distribution parameter
US20180052224A1 (en) * 2014-12-15 2018-02-22 Airbus Singapore Private Limited Automated method for selecting training areas of sea clutter and detecting ship targets in polarimetric synthetic aperture radar imagery
CN110658508A (en) * 2019-10-17 2020-01-07 中国人民解放军火箭军工程大学 K distribution sea clutter parameter estimation method based on characteristic quantity
CN112764000A (en) * 2020-12-24 2021-05-07 西安电子科技大学 Log-normal texture sea clutter amplitude model parameter estimation method based on quantile

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6040797A (en) * 1997-09-02 2000-03-21 Daimlerchrysler Ag Method for detecting a radar target
US20060072816A1 (en) * 2004-09-23 2006-04-06 Mitsubishi Denki Kabushiki Kaisha Determining statistical descriptors of a signal from a set of its samples
US20100073218A1 (en) * 2008-09-24 2010-03-25 Lockheed Martin Corpration Method and apparatus for radar surveillance and detection of sea targets
US20120154208A1 (en) * 2010-12-17 2012-06-21 Raytheon Company Methods and apparatus for automatic stc from sea state measurement via radar sea clutter eccentricity
US20180052224A1 (en) * 2014-12-15 2018-02-22 Airbus Singapore Private Limited Automated method for selecting training areas of sea clutter and detecting ship targets in polarimetric synthetic aperture radar imagery
US20170139038A1 (en) * 2015-11-13 2017-05-18 Rohde & Schwarz Gmbh & Co. Kg Method and signal generator for simulation of sea clutter
CN107167781A (en) * 2017-03-31 2017-09-15 西安电子科技大学 The quantile method of estimation of sea clutter amplitude lognormal distribution parameter
CN110658508A (en) * 2019-10-17 2020-01-07 中国人民解放军火箭军工程大学 K distribution sea clutter parameter estimation method based on characteristic quantity
CN112764000A (en) * 2020-12-24 2021-05-07 西安电子科技大学 Log-normal texture sea clutter amplitude model parameter estimation method based on quantile

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
旷仁杰: "基于海杂波的雷达信号传输特性研究", 《CNKI优秀硕士学位论文》 *
李发宗: ""基于极化信息的高频地波雷达目标检测技术研究"", 《CNKI优秀硕士学位论文》 *
杨超 等: "雷达海杂波反演大气波导的改进回溯搜索算法", 《系统工程与电子技术》 *
石志广: "K分布海杂波参数估计方法研究", 《信号处理》, vol. 23, no. 3 *

Also Published As

Publication number Publication date
CN113640763B (en) 2023-08-04

Similar Documents

Publication Publication Date Title
CN107167781B (en) Quantile estimation method for sea clutter amplitude log-normal distribution parameter
CN110275148B (en) Sea clutter amplitude distribution parameter estimation method and system
CN108508439B (en) Method for three-dimensional positioning of target collaborative imaging by double airborne SAR
CN105699952B (en) Double quantile methods of estimation of sea clutter K profile shape parameters
CN111381216B (en) Hybrid distribution radar sea clutter analysis method and device
CN106872958B (en) Radar target self-adapting detecting method based on linear fusion
CN110503071A (en) Multi-object tracking method based on the more Bernoulli Jacob's Additive Models of variation Bayes's label
CN110058222B (en) Double-layer particle filter tracking-before-detection method based on sensor selection
CN110376582B (en) Maneuvering target tracking method of self-adaptive GM-PHD
CN110865343A (en) LMB-based particle filter tracking-before-detection method and system
CN109991580A (en) A kind of satellite radar heighting wave reset method based on waveform derivative
CN107255805A (en) The Forecasting Methodology of radar target RCS based on weighted least-squares
CN110889862B (en) Combined measurement method for multi-target tracking in network transmission attack environment
CN116500575A (en) Extended target tracking method and device based on variable decibel leaf theory
CN107064893B (en) Pareto distribution with wide scope method for parameter estimation based on logarithmic moment
CN109709526B (en) Knowledge-assisted grouping generalized likelihood ratio detection method
CN114646935A (en) Sea clutter distribution parameter estimation method
CN111830481B (en) Radar echo single-component amplitude distribution model parameter estimation method and device
CN113640763B (en) Method for estimating distribution shape parameters of sea clutter amplitude of lognormal texture based on fractional order moment
CN109190647B (en) Active and passive data fusion method
CN115469314A (en) Uniform circular array steady underwater target azimuth tracking method and system
CN113466812B (en) Three-point estimation method for complex Gaussian sea clutter model parameters of inverse Gaussian texture
CN111025280B (en) Moving target speed measurement method based on distributed minimum total error entropy
CN107315169B (en) Clutter covariance matrix estimation method based on second-order statistic similarity
CN107179531B (en) Modified sample covariance matrix estimation algorithm based on maximum posterior

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