CN107728121B - Local goodness-of-fit inspection method based on variable window - Google Patents

Local goodness-of-fit inspection method based on variable window Download PDF

Info

Publication number
CN107728121B
CN107728121B CN201710762845.8A CN201710762845A CN107728121B CN 107728121 B CN107728121 B CN 107728121B CN 201710762845 A CN201710762845 A CN 201710762845A CN 107728121 B CN107728121 B CN 107728121B
Authority
CN
China
Prior art keywords
sea clutter
local
goodness
amplitude
fitting
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.)
Active
Application number
CN201710762845.8A
Other languages
Chinese (zh)
Other versions
CN107728121A (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 Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
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 Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201710762845.8A priority Critical patent/CN107728121B/en
Publication of CN107728121A publication Critical patent/CN107728121A/en
Application granted granted Critical
Publication of CN107728121B publication Critical patent/CN107728121B/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/40Means for monitoring or calibrating
    • G01S7/4052Means for monitoring or calibrating by simulation of echoes
    • 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

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 a local goodness-of-fit inspection method based on a variable window. The goodness-of-fit test is a quantitative judgment on a sea clutter amplitude distribution model, the problem that amplitude observation intervals are single by mean square error test (MSD) and modified mean square test (MMSD) methods is solved, the weighted result of the goodness-of-fit test of the selected local intervals is calculated by selecting one or more radar-interested local amplitude observation intervals, the method enhances the comprehensiveness of the goodness-of-fit test result of the statistical model, and provides a more reliable basis for the judgment of a sea clutter amplitude sequence statistical model and the design of a sea surface target constant false alarm detection (CFAR) method.

Description

Local goodness-of-fit inspection method based on variable window
Technical Field
The invention relates to a local goodness-of-fit testing method for sea clutter amplitude statistical distribution based on a variable window, and belongs to the technical field of radars.
Background
The sea target detection is an important task of the radar, however, the sea clutter is an important factor influencing the radar target detection, and in order to effectively realize the target detection under the background of the sea clutter, a statistical distribution model of the sea clutter must be analyzed. For actually measured sea clutter data, it is inaccurate to select the best-fitting sea clutter amplitude model from various distribution models by visual observation only, so a statistic for quantitatively describing the fitting degree of the sea clutter amplitude model and the actually measured data, namely, a goodness-of-fit test, is needed. The traditional classical goodness-of-fit inspection method mainly comprises the following steps: MSD method and MMSD method. However, the MSD method only focuses on global statistical fitting of sea clutter amplitude statistics, and neglects the fitting effect of amplitude local intervals; the MMSD method only focuses on the fitting effect of the sea clutter large-amplitude sequence, namely the tailing fitting effect. The two typical goodness-of-fit inspection methods only consider the global or tail fitting result of the sea clutter amplitude singly, and neglect the fitting condition in the sea clutter local amplitude interval under different sea condition conditions. When the radar actually works, complex sea clutter echo signals under different sea conditions can be met, and the model fitting result of analyzing the sea clutter from the whole situation or the tail part only singly is incomplete and inaccurate.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method for testing the local goodness of fit based on the variable window overcomes the defects of the prior art, can test the fitting condition in the local amplitude interval of the sea clutter under different sea conditions, enhances the reliability and comprehensiveness of the goodness of fit test, and provides a more reliable basis for target constant false alarm detection.
The technical solution of the invention is as follows: a local goodness-of-fit inspection method based on a variable window comprises the following steps:
(1) extracting a sea clutter amplitude sequence X ═ X from sea clutter echo data measured by radarkAnd equally dividing the sea clutter amplitude sequence into L intervals to obtain the actual probability density distribution p of the sea clutter amplitude statisticse(xl) Wherein N is the length of the sea clutter amplitude sequence, L is 1,2,3, …, L;
(2) according to the sea clutter amplitude sequence, calculating the amplitude statistical distribution model parameter of the sequence by adopting a moment estimation method, and fitting the sea clutter theoretical probability distribution function p of the statistical model by utilizing a common fitting methodt(xl);
(3) Determining the number, position and length of local observation windows;
(4) sea clutter theoretical probability distribution function p fitted by aiming at each fitting methodt(xl) Calculating the local goodness-of-fit test result of each observation window, and obtaining the sea clutter theory fitted by each fitting method through weighting processingProbability distribution function pt(xl) And the actual probability density distribution pe(xl) The final fitting result of (1);
(5) the minimum value of the final fitting result corresponds to the theoretical probability distribution function p of the sea cluttert(xl) Namely the statistical distribution model of the amplitude sequence of the sea clutter with the best fitting.
In the step (1), the sea clutter amplitude sequence is equally divided into L intervals to obtain the actual probability density distribution p of sea clutter amplitude statisticse(xl) The method comprises the following steps:
equally dividing the sea clutter amplitude sequence into L intervals according to the maximum value and the minimum value in the sea clutter amplitude sequence, counting the number of the sea clutter amplitudes falling into each interval, and obtaining the actual probability density distribution function p of sea clutter amplitude statisticse(xl),pe(xl) The number of sea clutter amplitudes falling into the ith interval/N. Wherein, the interval L>20。
In the step (2), the number of the local observation windows is equal to the number of the radar-interested local amplitude sequence areas, and the position [ N ] of the ith local observation windowi1,Ni2]Coinciding with the i-th local amplitude sequence region of interest of the radar, Ni1,Ni2Respectively represents the starting position and the ending position of the ith local observation window, and the length len of the ith local observation windowi=Ni2-Ni1≥4,Ni2、Ni1Are all positive integers.
In the step (4), the goodness-of-fit test result of the ith local observation window is calculated by using the following formula:
Figure BDA0001393505210000021
in the step (4), the final fitting result D of the j fitting method is calculated by using the following formulalmsd
Dlmsd=α1Dlmsd,12Dlmsd,2+…+αiDlmsd,i+…+αMDlmsd,M
0<αi<1, and α12+…+αi+…+αM<And 1 and M are the total number of the local observation windows.
αiIs proportional to the sea state grade.
Compared with the prior art, the invention has the following beneficial effects:
the method solves the problem that the amplitude observation interval is single in the traditional method, provides the idea of windowing, can select the number, the position and the length of the local observation windows according to the actual observation requirement of the radar, matches the most suitable distribution model by calculating the inspection result of goodness of fit after windowing, and can inspect the fitting condition in the sea clutter local amplitude observation interval under different sea conditions. When the radar needs to take the fitting results of multiple groups of different intervals (multiple local intervals or multiple local intervals and the whole amplitude interval) into consideration, the combined fitting results of the multiple groups of intervals can be weighted and evaluated, the comprehensiveness of the test results of the fitting goodness of the statistical model is enhanced, and a powerful support is provided for the design of the sea surface target constant false alarm detection algorithm.
Drawings
FIG. 1 is a flow chart of a method for goodness of fit of variable window based local statistics;
FIG. 2 is a statistical fit result of an actually measured sea clutter amplitude sequence;
fig. 3 shows the goodness-of-fit test results of the actually measured sea clutter data, (a) shows MSD test results of four common distributions, (b) shows MMSD test results of four common distributions, and (c) shows local goodness-of-fit test results of four common distributions.
Detailed Description
In consideration of model fitting conditions under different sea conditions, the invention provides a variable window-based local goodness of fit (LMSD) method, which overcomes the problem that the existing method is single in amplitude observation area, and can calculate the goodness of fit inspection result of a local amplitude observation interval through windowing according to the observation requirement of a radar; when a plurality of observation intervals with different amplitudes are interested, the amplitude statistical goodness-of-fit test result under the combined action of the plurality of intervals can be calculated through weighting processing of the plurality of observation intervals, so that the reliability and the comprehensiveness of the goodness-of-fit test are enhanced, and a more reliable basis is provided for judgment of a sea clutter amplitude statistical model and design of a sea surface target constant false alarm detection method.
As shown in fig. 1, the main implementation steps of the present invention are as follows:
step 1, extracting a sea clutter amplitude sequence X ═ X from sea clutter echo data actually measured by a radarkAnd k is 1,2,3,.. N }, and N is the length of the sea clutter amplitude sequence. Equally dividing the sea clutter amplitude sequence into L intervals according to the maximum value and the minimum value in the sea clutter amplitude sequence, counting the occurrence frequency of the sea clutter amplitude in each interval, and obtaining the actual probability density distribution function p of sea clutter amplitude statisticse(xl). Wherein, L is not selected too small, and L is small, which may result in inaccurate fitting result of the amplitude distribution function, therefore, L is generally selected>20。l=1,2,3,…,L。
The interval between two adjacent intervals is (maximum value-minimum value in amplitude sequence)/L, pe(xl) Dividing the number of the sea clutter amplitudes falling into the ith interval by the total number N of the sea clutter amplitudes.
Step 2, sea clutter amplitude statistical fitting:
according to the sea clutter amplitude sequence, calculating the amplitude statistical distribution model parameter of the sequence by adopting a moment estimation method, and fitting the sea clutter theoretical probability distribution function p of the statistical model by utilizing a common fitting methodt(xl)。
Step 3, determining the number, position and length of the local observation windows:
3.1 select the number M of observation windows. If the radar is only interested in a single local region of the sea clutter amplitude statistics, then a single rectangular window is selected, i.e. M is 1. If the radar is interested in a plurality of local areas of sea clutter amplitude statistics, a rectangular window is selected for each area, namely M >1, and then weighting processing is carried out on a plurality of local area fitting results. The invention mainly takes the fitting result of one group or two groups of observation intervals as an example for analysis, if more observation intervals are interested, M can also be selected to be a larger value, and the analysis method is similar to the case that M is 2.
3.2 after the number of the observation windows is determined according to the number 3.1, the position of the observation windows can be selected according to the interested area of the radar. When the radar is interested in the statistical distribution of the local amplitude of a certain part of the sea clutter, a rectangular window can be added to any observation interval; a rectangular window may be added over the entire observation interval when the radar is interested in a global fitting statistical distribution of sea clutter.
3.3 selection of the local observation window length. Let the ith local observation window be at [ N ]i1,Ni2],Ni1,Ni2Respectively representing the start and end positions, N, of the ith local observation windowi2、Ni1Are all positive integers, and the length len of the ith local observation windowi=Ni2-Ni1And the amplitude range of the windowed sequence is more than 4N/L of the amplitude of the whole sequence, so that enough amplitude samples can participate in operation, and the random error influence caused by too few noise and samples is reduced.
Step 4, calculating the local goodness of fit (LMSD) result of each observation window
For the jth fitting method, the result D of goodness-of-fit test of the ith local observation windowlmsd,i
Figure BDA0001393505210000051
The final fitting result D of the jth fitting method is calculated by using the following formulalmsd
Dlmsd=α1Dlmsd,12Dlmsd,2+…+αiDlmsd,i+…+αMDlmsd,M
0<αi<1, and α12+…+αi+…+αM<And 1 and M are the total number of the local observation windows.
αiIs proportional to the level of the sea state, α is selected according to the actual observation requirementiThe higher the sea state grade is, the higher the weight of the large-amplitude sequence can be increased; sea stateThe lower the rank, the more equally the weight of the amplitude sequence can be reduced.
The choice of M weights is analyzed by taking M ═ 2 as an example, and so on for the case of M >2, that is:
Dlmsd=α1Dlmsd,12Dlmsd,2,0<α1<1,0<α2<1
α therein12Respectively represents the weight values of the two groups of goodness-of-fit results, and satisfies α12=1。
Pair weight α12Can be selected according to the actual radar observation needs. Under the condition of high sea, the weight of the large-amplitude observation interval can be increased; under the condition of low sea, the weight of the small-amplitude observation interval can be increased.
Suppose Dlmsd,1And Dlmsd,2Respectively representing the fitting test result of the global amplitude and the fitting test result of the local amplitude, and if the radar focuses on the result of the statistical combined action of the amplitude of the sea clutter and the fitting test result of the global amplitude and the fitting test result of the local amplitude, α can be selected1=α20.5, α may be selected if the radar is more focused on global fit test results1=0.7,α20.3, and α may be selected if the radar is more focused on local fit test results1=0.3,α20.7, wherein, α12The proportion of (A) can be properly adjusted according to actual needs.
Step 5, judging a statistical distribution model:
obtaining a local goodness-of-fit test result D of the sea clutter amplitude according to the step 4lmsdSimilar to conventional MSD and MMSD, DlmsdA smaller value of (a) indicates a better fit over this interval. By comparing different statistical distribution modelslmsdAs a result, the set of sea clutter data can be determined to fit the best statistical distribution model.
The effect of the invention is demonstrated by the test of measured data:
simulation conditions
Extracting sea clutter echo from sea clutter echo data actually measured by certain shore-based S-band radarWave amplitude sequence X ═ { X ═ Xk,k=1,2,3,...,N}。
(II) simulation content
FIG. 2 is a statistical fit result of the actually measured sea clutter amplitude sequence. The test divides the sea clutter amplitude sequence into L (L is 50) observation intervals, and the histogram is the actual probability density distribution p of the sea clutter amplitude statisticse(xl) The four curves respectively correspond to the sea clutter theoretical probability distribution function p obtained by fitting the four fitting methodst(xl)。
Aiming at the fitting result of the measured data in fig. 2, the MSD test, the MMSD test and the LMSD test proposed by the present invention are respectively adopted to evaluate the goodness and badness of the fitting result.
Fig. 3 (a) shows the results of four common MSD tests, and it can be seen that: for the global statistical fitting of the sea clutter amplitude, the fitting effect of the K distribution is the best, the Weibull distribution and the Rayleigh distribution are followed, and the fitting effect of the log-normal distribution is the worst.
FIG. 3 (b) shows the results of MMSD test for four common distributions, which can be seen: for the statistical fitting (trailing fitting result) of the sea clutter large amplitude sequence, the fitting effect of the Weibull distribution is the best, and the fitting effect of the K distribution, the Rayleigh distribution and the log-normal distribution is the worst.
FIG. 3(c) shows the local goodness-of-fit test results for four common distributions, where two observation windows are selected, the first window being [1, 50 ], taking the weighted fit result as an example]Calculating a global fitting test result of the amplitude; the second observation window is [40, 50 ]]In addition, the weight of the two groups of fitting test results is selected as α1=α20.5. As can be seen in fig. 3 (c): for the statistical fitting of the overall sea clutter amplitude and the combined action of the tail, the fitting effect of K distribution is the best, and then Weibull distribution and Rayleigh distribution are carried out, and the fitting effect of log-normal distribution is the worst.
MSD and MMSD only consider partial information of sea clutter amplitude, however, the local goodness-of-fit inspection method provided by the invention integrates global and local information, and is more comprehensive than the number of samples adopted from the whole and tail part intervals and other partial intervals, so that the conclusion obtained is more accurate and reliable.
Those skilled in the art will appreciate that the details of the invention not described in detail in the specification are within the skill of those skilled in the art.

Claims (5)

1. A local goodness-of-fit inspection method based on a variable window is characterized by comprising the following steps:
(1) extracting a sea clutter amplitude sequence X ═ X from sea clutter echo data measured by radarkAnd equally dividing the sea clutter amplitude sequence into L intervals to obtain the actual probability density distribution p of the sea clutter amplitude statisticse(xl) Wherein N is the length of the sea clutter amplitude sequence, L is 1,2,3, …, L;
equally dividing the sea clutter amplitude sequence into L intervals to obtain the actual probability density distribution p of sea clutter amplitude statisticse(xl) The method comprises the following steps:
equally dividing the sea clutter amplitude sequence into L intervals according to the maximum value and the minimum value in the sea clutter amplitude sequence, counting the number of the sea clutter amplitudes falling into each interval, and obtaining the actual probability density distribution function p of sea clutter amplitude statisticse(xl),pe(xl) Number of sea clutter amplitudes/N falling into the L-th interval, wherein the interval L>20;
(2) According to the sea clutter amplitude sequence, calculating the amplitude statistical distribution model parameter of the sequence by adopting a moment estimation method, and fitting the sea clutter theoretical probability distribution function p of the statistical model by utilizing a common fitting methodt(xl);
(3) Determining the number, position and length of local observation windows;
(4) sea clutter theoretical probability distribution function p fitted by aiming at each fitting methodt(xl) Calculating the local goodness-of-fit test result of each observation window, and obtaining the sea clutter theoretical probability distribution function p fitted by each fitting method through weighting processingt(xl) And the actual probability density distribution pe(xl) The final fitting result of (1);
(5) the minimum value of the final fitting result corresponds to the theoretical probability distribution function p of the sea cluttert(xl) Namely the statistical distribution model of the amplitude sequence of the sea clutter with the best fitting.
2. The method of claim 1, wherein the local goodness-of-fit test based on the variable window is as follows: in the step (2), the number of the local observation windows is equal to the number of the radar-interested local amplitude sequence areas, and the position [ N ] of the ith local observation windowi1,Ni2]Coinciding with the i-th local amplitude sequence region of interest of the radar, Ni1,Ni2Respectively represents the starting position and the ending position of the ith local observation window, and the length len of the ith local observation windowi=Ni2-Ni1≥4,Ni2、Ni1Are all positive integers.
3. The method of claim 2, wherein the local goodness-of-fit test based on the variable window is as follows: in the step (4), the goodness-of-fit test result of the ith local observation window is calculated by using the following formula:
Figure FDA0002166784820000021
4. the method of claim 3, wherein the local goodness-of-fit test based on the variable window is as follows: in the step (4), the final fitting result D of the j fitting method is calculated by using the following formulalmsd
Dlmsd=α1Dlmsd,12Dlmsd,2+…+αiDlmsd,i+…+αMDlmsd,M
0<αi<1, and α12+…+αi+…+αM<And 1 and M are the total number of the local observation windows.
5. According toThe method of claim 4, wherein α is a method for testing local goodness-of-fit based on variable windowiIs proportional to the sea state grade.
CN201710762845.8A 2017-08-30 2017-08-30 Local goodness-of-fit inspection method based on variable window Active CN107728121B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710762845.8A CN107728121B (en) 2017-08-30 2017-08-30 Local goodness-of-fit inspection method based on variable window

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710762845.8A CN107728121B (en) 2017-08-30 2017-08-30 Local goodness-of-fit inspection method based on variable window

Publications (2)

Publication Number Publication Date
CN107728121A CN107728121A (en) 2018-02-23
CN107728121B true CN107728121B (en) 2020-04-10

Family

ID=61205556

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710762845.8A Active CN107728121B (en) 2017-08-30 2017-08-30 Local goodness-of-fit inspection method based on variable window

Country Status (1)

Country Link
CN (1) CN107728121B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11347152B2 (en) 2018-01-17 2022-05-31 Asml Netherlands B.V. Scan signal characterization diagnostics
CN111580091B (en) * 2020-05-17 2022-08-05 西北工业大学 Sea surface weak target detection method based on AR spectrum singular intensity function

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104992159A (en) * 2015-07-13 2015-10-21 北京化工大学 SAR image ship detection method based on Gaussian function fitting of sea clutter
CN105046706A (en) * 2015-07-13 2015-11-11 北京化工大学 Rational polynomial function fitting sea clutter based SAR image ship detection method
CN105069775A (en) * 2015-07-13 2015-11-18 北京化工大学 SAR image ship detection method based on Fourier series

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104992159A (en) * 2015-07-13 2015-10-21 北京化工大学 SAR image ship detection method based on Gaussian function fitting of sea clutter
CN105046706A (en) * 2015-07-13 2015-11-11 北京化工大学 Rational polynomial function fitting sea clutter based SAR image ship detection method
CN105069775A (en) * 2015-07-13 2015-11-18 北京化工大学 SAR image ship detection method based on Fourier series

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"Weak target detection in sea clutter background using local-multifractal spectrum with adaptive window length";Yifei Fan et al.;《IET Radar, Sonar & Navigation》;20150930;第9卷(第7期);第835-842页 *

Also Published As

Publication number Publication date
CN107728121A (en) 2018-02-23

Similar Documents

Publication Publication Date Title
CN111722199B (en) Radar signal detection method based on convolutional neural network
CN103995259B (en) Radar target adaptive-filtering fusion detection method under intensive interference environment
CN110208806B (en) Marine radar image rainfall identification method
CN109143195B (en) Radar target detection method based on full KL divergence
CN110658508B (en) K distribution sea clutter parameter estimation method based on characteristic quantity
CN104898103B (en) Low velocity target detection method based on multichannel clutter map
CN103760542A (en) MMVI-CFAR target detection method
CN112859012B (en) Radar spoofing interference identification method based on cascade convolution neural network
CN109212503B (en) Offshore target detection method and device based on AR spectrum fractal correlation characteristics
Newman et al. Range-correcting azimuthal shear in Doppler radar data
CN109444869A (en) A kind of radar extension target component tunable detector for SLM Signal Label Mismatch
CN111381216B (en) Hybrid distribution radar sea clutter analysis method and device
CN107728121B (en) Local goodness-of-fit inspection method based on variable window
CN104198998B (en) Clustering treatment based CFAR (Constant False Alarm Rate) detection method under non-uniform background
Zou et al. A method of radar echo extrapolation based on TREC and Barnes filter
CN105116412A (en) Broadband radar ISAR image performance quantitative discrimination method
CN106526545A (en) Method for detecting robust CFAR detector
CN114563771A (en) Double-threshold laser radar cloud layer detection algorithm based on cluster analysis
CN105548986B (en) Object detection method based on prewhitening median ratio detector under sea clutter background
CN105528657A (en) Building earthquake damage prediction method based on Beidou and vector machine
Palencia et al. Dent overlap in hailpads: Error estimation and measurement correction
CN112731388A (en) Target detection method based on effective scattering point energy accumulation
CN106353752A (en) Method for measuring soil water content based on ultra wide band and non-singleton fuzzy logic
CN114578333B (en) Active sonar target dynamic and static identification method
CN110780275B (en) Method for removing batch increase by signal sorting

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