CN113267767A - Multi-base sonar buoy detection efficiency analysis method and storage medium - Google Patents
Multi-base sonar buoy detection efficiency analysis method and storage medium Download PDFInfo
- Publication number
- CN113267767A CN113267767A CN202110430578.0A CN202110430578A CN113267767A CN 113267767 A CN113267767 A CN 113267767A CN 202110430578 A CN202110430578 A CN 202110430578A CN 113267767 A CN113267767 A CN 113267767A
- Authority
- CN
- China
- Prior art keywords
- sonobuoy
- field
- analyzed
- probability
- density
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52004—Means for monitoring or calibrating
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a multi-base sonobuoy detection efficiency analysis method and a storage medium, the invention identifies the type of a sonobuoy field to be analyzed, then calculates the density of the sonobuoy field based on the type of the sonobuoy field, further determines the detection probability of the sonobuoy field to be analyzed according to the density of the sonobuoy field, finally adjusts the type and the layout of the sonobuoys in the sonobuoy field to be analyzed based on the detection probability and the comprehensive performance of the sonobuoys, so as to reduce the layout cost of the sonobuoys in the sonobuoy field under the condition of reaching the preset detection probability, thereby realizing the analysis and the adjustment of the detection efficiency of the multi-base sonobuoys with any flexible quantity and combination, providing an analysis means for the multi-base sonobuoy task sea area, resource configuration planning and the like, and further effectively solving the problems of the prior art that a transmitter is used for the multi-, The number and combination of the receivers are arbitrary and flexible, and the detection efficiency of the multi-base sonar buoy cannot be effectively analyzed.
Description
Technical Field
The invention relates to the technical field of sonobuoy detection, in particular to a multi-base sonobuoy detection efficiency analysis method and a computer readable storage medium.
Background
With the development of submarine damping and noise reduction technology and the use of a large number of silencing tiles, the level of submarine radiation noise is gradually reduced. The traditional single-base sonar is more and more difficult to meet the underwater target detection requirement, and the multi-base sonar system is valued by various countries due to the advantages of strong anti-interference capability, better detection performance, tactical use concealment and the like. However, in the current multi-base sonar system, the multi-base sonar buoy array detection efficiency can be improved by optimizing the layout of the sonar buoys only under the condition that the types of the transmitters and the receivers are uniform, the array type and the number of the transmitters and the receivers are known, but the multi-base sonar buoy detection efficiency which is arbitrary and flexible in number and combination of the transmitters and the receivers cannot be solved, so that the task sea area and resource allocation planning are influenced.
Disclosure of Invention
The invention provides a multi-base sonobuoy detection efficiency analysis method and a computer readable storage medium, which aim to solve the problem that the multi-base sonobuoy detection efficiency with any flexible combination of the number of transmitters and receivers cannot be solved in the prior art.
In a first aspect, the present invention provides a method for analyzing detection efficiency of a multi-base sonobuoy, the method comprising: identifying the type of sonobuoy field to be analyzed; calculating the density of a sonobuoy field of the sonobuoy field to be analyzed based on the type of the sonobuoy field, determining the detection probability of the sonobuoy field to be analyzed according to the density of the sonobuoy field, and adjusting the type and layout of the sonobuoys in the sonobuoy field to be analyzed based on the detection probability and the comprehensive performance of the sonobuoys so as to reduce the layout cost of the sonobuoys in the sonobuoy field under the condition of reaching the preset detection probability; the sonobuoy field to be analyzed is a multi-base sonobuoy, and the sonobuoy comprises a transmitter SBuoy, a receiver RBuoy and a transmitting and receiving common-position machine C _ SRBuoy which are independently and randomly distributed.
Optionally, before the identifying the sonobuoy field type to be analyzed, the method further comprises:
presetting an I-type sonobuoy field and a II-type sonobuoy field;
the I-type sonobuoy field is provided with an independent randomly distributed transmitter SBuoy and an independent randomly distributed receiver RBuoy;
the II-type sonobuoy field is an independent and randomly distributed transmitting and receiving common-position machine C _ SRBuoy and an independent and randomly distributed receiver RBuoy.
Optionally, the identifying a sonobuoy field type to be analyzed includes: and identifying the type of the sonobuoy field to which the sonobuoy field to be analyzed belongs based on the composition among the transmitter, the receiver and the transmitting and receiving co-location machine distributed in the sonobuoy field to be analyzed.
Optionally, when the sonobuoy field to be analyzed is a type I sonobuoy field, the calculating the sonobuoy field density of the sonobuoy field to be analyzed based on the type of the sonobuoy field includes: calculating the area A formed by all the sonobuoys in the sonobuoy field to be analyzed1Random distribution Density of SBuoy
Random placement density of RBuoy
The density of the sonobuoy field to be analyzed is: λ ═ λs,λr](ii) a Wherein Ns is the number of SBuoy, and Nr is the number of RBuoy.
Optionally, the area a formed by all the sonobuoys in the sonobuoy field to be analyzed is obtained1The method comprises the following steps: solving the area of the convex polygon sonar buoy field according to a convex polygon area formula:
Optionally, determining a detection probability of the sonobuoy field according to the density of the sonobuoy field to be analyzed includes: to sonar buoy yardThe probability that there is no sonobuoy within an arbitrary radius R, i.e. the probability that the distance R of a sonobuoy point from the random target T is greater than R, is the probability that a poisson random variable τ is 0:
where λ is the random deployment density of the sonobuoy, then the probability of no SBuoy sonobuoy within radius r of a random target T
Probability of no RBuoy sonobuoy within radius r of random target T
According to probability density function
To obtain
According to
Calculating the detection failure probability of the sonobuoy field:
further according to P (T) ═ 1-QICalculating the detection probability of the sonar buoy field;
wherein R issDistance, λ, of SBuoy from random target TsRandom distribution density for SBuoy, RrDistance, λ, of RBuoy from random target TrAnd randomly distributing density for RBuoy.
Optionally, when the sonobuoy field to be analyzed is a type II sonobuoy field, the calculating the sonobuoy field density of the sonobuoy field to be analyzed based on the type of the sonobuoy field includes:
calculating the area A formed by all the sonobuoys in the sonobuoy field to be analyzed2Random distribution density of C _ SRBuoy
Random placement density of RBuoy
The density of the sonobuoy field to be analyzed is: λ ═ λsr,λr](ii) a Wherein Nsr is the number of C _ SRBuoy and Nr is the number of RBuoy.
Optionally, the area a formed by all the sonobuoys in the sonobuoy field to be analyzed is obtained2The method comprises the following steps: solving the area of the convex polygon sonar buoy field according to a convex polygon area formula:
wherein, the first sonobuoy (x21, y21), the second sonobuoy (x22, y22), the first sonobuoy (x23, y23) … … nth sonobuoy (x2n, y2n) are the coordinates of each sonobuoy in the type II sonobuoy field.
Optionally, determining a detection probability of the sonobuoy field according to the density of the sonobuoy field to be analyzed includes: for a random target T in the sonobuoy field, the probability that there is no sonobuoy within any radius R, i.e. the probability that the distance R of a sonobuoy point from the random target T is greater than R, is the probability that a poisson random variable τ is 0:
wherein lambda is the random deployment density of the sonobuoy, the probability that no RBuoy sonobuoy exists within the radius r of the random target T
Probability of no C _ SRBuoy sonobuoy within radius r of random target T
Because of the probability density function
Then
Will be provided with
Substitution into QIIQ can be calculatedII
Further according to P (T) ═ 1-QIICalculating the detection probability of the sonar buoy field;
wherein R issrIs the distance, λ, of C _ SRBuoy from the random target TsrIs the random distribution density of C _ SRBuoy, RrDistance, λ, of RBuoy from random target TrThe density is randomly laid for the RBuoy,
is the bistatic buoy system equivalent tactical action radius to the random target T.
In a second aspect, the present invention provides a computer-readable storage medium, which stores a signal-mapped computer program, and when the computer program is executed by at least one processor, the computer program implements any one of the multi-base sonar buoy detection performance analysis methods described above.
The invention has the following beneficial effects:
the invention identifies the type of the sonobuoy field to be analyzed, then calculates the density of the sonobuoy field based on the type of the sonobuoy field, further determines the detection probability of the sonobuoy field to be analyzed according to the density of the sonobuoy field, finally adjusts the type and the layout of the sonobuoy in the sonobuoy field to be analyzed based on the detection probability and the comprehensive performance of the sonobuoy, so as to reduce the layout cost of the sonobuoy in the sonobuoy field under the condition of reaching the preset detection probability, thereby realizing the analysis and adjustment of the detection efficiency of the multi-base sonar buoy with any flexible quantity and combination, providing an analysis means for the multi-base sonar buoy task sea area, resource allocation planning and the like, and further effectively solves the problem that the detection efficiency of the existing multi-base sonar buoy with any flexible combination of the number of transmitters and receivers cannot be effectively analyzed.
The foregoing description is only an overview of the technical solutions of the present invention, and the embodiments of the present invention are described below in order to make the technical means of the present invention more clearly understood and to make the above and other objects, features, and advantages of the present invention more clearly understandable.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
FIG. 1 is a schematic view of a multi-base sonobuoy detection performance analysis process according to a first embodiment of the present invention;
FIG. 2 is a schematic diagram of the detection principle of the present double-base sonar buoy;
FIG. 3 is a schematic diagram of the field type of two types of multi-base sonar buoys provided by the first embodiment of the present invention;
FIG. 4 is a schematic view of an alternative multi-base sonar buoy detection efficiency analysis process according to the first embodiment of the present invention;
fig. 5 is a schematic layout view of a type I sonobuoy provided in a first embodiment of the present invention;
fig. 6 is a schematic layout view of a type II sonobuoy provided in the first embodiment of the present invention.
Detailed Description
Aiming at the problem that the prior multi-base sonar system can only improve the detection efficiency of a multi-base sonar buoy array by optimizing the layout of a sonar buoy under the condition that the types of transmitters and receivers are uniform, the array type and the quantity are known, but cannot solve the detection efficiency of any flexible multi-base sonar buoy formed by combining the numbers of the transmitters and the receivers, the invention identifies the type of a sonar buoy field to be analyzed, calculates the density of the sonar buoy field based on the type of the sonar buoy field, further determines the detection probability of the sonar buoy field to be analyzed according to the density of the sonar buoy field, and finally adjusts the type and the layout of the sonar buoy in the sonar buoy field to be analyzed based on the detection probability and the comprehensive performance of the sonar buoy so as to reduce the layout cost of the sonar buoy in the sonar buoy field under the condition of reaching the preset detection probability, therefore, the multi-base sonar buoy detection efficiency with any flexible combination and quantity can be analyzed and adjusted, an analysis means is provided for multi-base sonar buoy task sea areas, resource allocation planning and the like, and the problem that the existing multi-base sonar buoy detection efficiency with any flexible combination and quantity of transmitters and receivers cannot be effectively analyzed is effectively solved. The present invention will be described in further detail below with reference to the drawings and examples. It should be understood that the specific embodiments described herein are merely illustrative of the invention and do not limit the invention.
The embodiment of the invention provides a multi-base sonar buoy detection efficiency analysis method, and referring to fig. 1, the method comprises the following steps:
s101, identifying the type of a sonobuoy field to be analyzed;
it should be noted that the embodiment of the present invention is directed to a multi-base sonobuoy, and the sonobuoy includes a transmitter SBuoy, a receiver RBuoy, and a transmitting-receiving co-located computer C _ srbuy, which are independently and randomly deployed.
In specific implementation, before step S101, the method according to the embodiment of the present invention further includes: presetting an I-type sonobuoy field and a II-type sonobuoy field; the type I sonobuoy field is a transmitter SBuoy and a receiver RBuoy which are independently and randomly distributed; the II-type sonobuoy field is an independent and randomly distributed transmitting and receiving common-position machine C _ SRBuoy and a receiver RBuoy.
Specifically, the identifying a type of a sonobuoy field to be analyzed according to the embodiment of the present invention includes: and identifying the type of the sonobuoy field to which the sonobuoy field to be analyzed belongs based on the composition among the transmitter, the receiver and the transmitting and receiving co-location machine distributed in the sonobuoy field to be analyzed.
That is, in the embodiment of the present invention, different sonobuoy field types are set based on a combined layout among a transmitter, a receiver, and a transmitting-receiving co-located computer, and in order to facilitate implementation, two types of sonobuoy fields are set in the embodiment of the present invention, which are an I-type sonobuoy field and an II-type sonobuoy field, respectively, where the I-type sonobuoy field is a transmitter SBuoy and a receiver RBuoy that are independently and randomly distributed; the type II sonobuoy field is an independent and randomly distributed transmitting and receiving common-position machine C _ SRBuoy and a receiver RBuoy.
S102, calculating the density of the sonobuoy field to be analyzed based on the type of the sonobuoy field;
when the sonobuoy field to be analyzed is a type I sonobuoy field, the calculating the sonobuoy field density of the sonobuoy field to be analyzed based on the type of the sonobuoy field includes: calculating the area A formed by all the sonobuoys in the sonobuoy field to be analyzed1Random distribution Density of SBuoy
Random placement density of RBuoy
The density of the sonobuoy field to be analyzed is: λ ═ λs,λr](ii) a Wherein Ns is the number of SBuoy, and Nr is the number of RBuoy.
When the sonobuoy field to be analyzed is of type IIA sonobuoy field, said calculating a sonobuoy field density of said sonobuoy field to be analyzed based on said sonobuoy field type, comprising: calculating the area A formed by all the sonobuoys in the sonobuoy field to be analyzed2Random distribution density of C _ SRBuoy
Random placement density of RBuoy
The density of the sonobuoy field to be analyzed is: λ ═ λsr,λr](ii) a Wherein Nsr is the number of C _ SRBuoy and Nr is the number of RBuoy.
S103, determining the detection probability of the sonobuoy field to be analyzed according to the density of the sonobuoy field;
when the sonobuoy field to be analyzed is a type I sonobuoy field: for a random target T in the sonobuoy field, the probability that there is no sonobuoy within any radius R, i.e. the probability that the distance R of a sonobuoy point from the random target T is greater than R, is the probability that a poisson random variable τ is 0:
where lambda is the random deployment density of the sonobuoy, then,
probability of no SBuoy sonobuoy within radius r of random target T
Probability of no RBuoy sonobuoy within radius r of random target T
According to probability density function
To obtain
According to
Calculating the detection failure probability of the sonobuoy field:
further according to P (T) ═ 1-QICalculating the detection probability of the sonar buoy field;
wherein R issDistance, λ, of SBuoy from random target TsRandomly laying density, R, for RBuoyrDistance, λ, of SBuoy from random target TrAnd randomly distributing density for RBuoy.
When the sonobuoy field to be analyzed is a type II sonobuoy field, for a random target T in the sonobuoy field, the probability that there is no sonobuoy within any radius R, that is, the probability that the distance R of a sonobuoy point from the random target T is greater than R is the probability that a poisson random variable τ is 0:
where lambda is the random deployment density of the sonobuoy, then,
probability of no RBuoy sonobuoy within radius r of random target T
Probability of no C _ SRBuoy sonobuoy within radius r of random target T
Because of the probability density function
Then
Will be provided with
Substitute QIIQ can be calculatedII
Further according to P (T) ═ 1-QIICalculating the detection probability of the sonar buoy field;
wherein R issrIs the distance, λ, of C _ SRBuoy from the random target TsrIs the random distribution density of C _ SRBuoy, RrDistance, λ, of RBuoy from random target TrAnd randomly distributing density for RBuoy.
S104, adjusting the type and layout of the sonobuoy in the sonobuoy field to be analyzed based on the calculated detection probability, the preset detection probability and the comprehensive performance of the sonobuoy so as to reduce the layout cost of the sonobuoy in the sonobuoy field under the condition of reaching the preset detection probability.
It should be noted that the preset detection probability in the embodiment of the present invention is a detection probability that needs to be achieved by a preset operation of an operator, and a person skilled in the art may set the preset detection probability according to actual needs in specific implementation, which is not specifically limited in the present invention.
Specifically, the embodiment of the invention performs comprehensive analysis and adjustment based on the detection probability and the preset detection probability obtained by calculation, the performance of the transmitting and receiving signals of each sonobuoy in the to-be-analyzed sonobuoy field, the cost of each sonobuoy, the planned payment cost, the area of the sonobuoy field and the like, so that the layout cost of the sonobuoys in the sonobuoy field is reduced as much as possible under the condition that the type and the layout of the sonobuoys in the to-be-analyzed sonobuoy field reach the preset detection probability, and finally, the sonobuoy combination with the minimum resource quantity, the minimum expenditure and the optimal detection probability is given.
That is, in the embodiment of the present invention, by adjusting the area a and the resource number N (the resource data is the number of the above-mentioned transmitter SBuoy, receiver RBuoy, and transmitting/receiving co-location machine C _ SRBuoy), the sonar buoy field density λ under different resource numbers and different combinations is calculated, and the multi-base sonar buoy detection efficiency p (t) with any flexible resource number and combination is given.
That is to say, in the embodiment of the present invention, the density of the sonobuoy field in the sonobuoy field is calculated based on the type of the sonobuoy field to be analyzed, the detection probability of the sonobuoy field to be analyzed is further determined according to the density of the sonobuoy field, and finally the type and layout of the sonobuoys in the sonobuoy field to be analyzed are adjusted based on the detection probability and the comprehensive performance of the sonobuoys, so that the layout cost of the sonobuoys in the sonobuoy field is reduced when the preset detection probability is reached, thereby realizing the analysis and adjustment of the detection efficiency of the multi-base sonobuoys in any flexible quantity and combination, and providing an analysis means for the multi-base sonobuoy mission sea area, resource allocation planning, and the like.
In specific implementation, 3 types of sonobuoy application types and 2 types of sonobuoy fields are divided from an application angle so as to be convenient for selection in planning application of a mission sea area, resource allocation and the like, wherein the 3 types of sonobuoy application types are a transmitter SBuoy, a receiver RBuoy and a transmitting-receiving common computer C _ srbuy which are independently and randomly distributed, and the 2 types of sonobuoy fields are an I type sonobuoy field and a II type sonobuoy field. The I-type sonobuoy field is a transmitter SBuoy and a receiver RBuoy which are independently and randomly distributed, and the II-type sonobuoy field is a transmitting and receiving common-position machine C _ SRBuoy and a receiver RBuoy which are independently and randomly distributed.
In general, the core idea of the method according to the embodiment of the present invention is: the detection probability of the known buoy field can be determined, namely the detection effect of the known sonar buoy field can be determined by the method, and the best effect can be achieved by the method under the buoy field with the precondition of detection probability requirement, field area requirement or layout requirement and the like.
In brief, the embodiment of the invention is based on the basic principle of sonar equation, and the equivalent tactical action radius gamma of the bistatic sonar buoy system in transmitting and receiving is derived according to the spherical wave propagation loss formula TL. And representing the probability of the target T not being detected by using the probability Q that the Poisson random number increment is 0, and giving a detection probability calculation formula P (T) 1-Q of the multi-base sonobuoy to the specific target T.
For a type II sonobuoy field, when the target T distance C _ SRBuoy is greater than gamma, the target T cannot be detected by the C _ SRBuoy. When the distance product of the target T and the distances C _ SRBuoy and RBuoy is more than gamma2In time, the sonobuoy field composed of C _ SRBuoy and RBuoy cannot detect the target. The type II sonobuoy field cannot detect the target only when the target T cannot be detected by the C _ SRBuoy and the sonobuoy field consisting of the C _ SRBuoy and the RBuoy, so that the target T cannot be detected by the type II sonobuoy field
Because C _ SRBuoy and RBuoy are independent of each other, the embodiment of the present invention obtains:
By adjusting the proportion of the area A to the resource quantity N, calculating the density lambda of the sonar buoy field formed by different quantities and combined transmitters and receivers, giving out the multi-base sonar buoy detection efficiency P (T) with uniform types, quantity and combination of the transmitters and the receivers, and providing an analysis means for the multi-base sonar buoy task sea area and resource allocation planning. And (4) selecting the optimal sonar buoy field density lambda by traversing the detection efficiency curve, thereby providing reference for the optimization of the multi-base sonar buoy task sea area and resource allocation.
As shown in fig. 4, the method for detecting the efficiency of the multi-base sonobuoy provided by the present invention comprises the following steps:
the method comprises the following steps: the sonobuoy type is initialized. In the present invention, 3 unified sonobuoy types are defined according to the deployment method of the transmitter and the receiver, as shown in fig. 2:
independent transmitter SBuoy: the transmitters are independently and randomly distributed;
independent receiver RBuoy: the receivers are independently and randomly distributed;
co-located transceiver C _ srbus: co-location and random arrangement of a transmitter and a receiver;
step two: a sonobuoy field is defined as in fig. 3.
(1) Type I sonobuoy field: SBuoy and RBuoy covered task area A1；
As shown in fig. 5, in general, the sonobuoy is a convex polygon, and the area of the convex polygon sonobuoy field is solved according to the convex polygon area formula:
the first sonobuoy (x11, y11), the second sonobuoy (x12, y12), the first sonobuoy (x13, y13) … … nth sonobuoy (x1n, y1n) are the coordinates of each sonobuoy in the I-type sonobuoy field.
(2) Type II sonobuoy field: task area A covered by C _ SRBuoy and RBuoy2。
As shown in fig. 6, the convex polygon sonobuoy field area is solved according to the convex polygon area formula:
wherein, the first sonobuoy (x21, y21), the second sonobuoy (x22, y22), the first sonobuoy (x23, y23) … … nth sonobuoy (x2n, y2n) are the coordinates of each sonobuoy in the type II sonobuoy field.
Step three: definition of detection probability formula of sonobuoy field to target
According to the sonar equation basic principle of 'signal level-background interference ═ detection threshold', the signal level of the bistatic sonar buoy is SL-TLT+TSWhere SL is the sound source level of the transmitter, TLTFor propagation loss of acoustic waves from transmitter to target, TSIs the target intensity. Background interference is NL + TLRDI, where NL is the ocean ambient noise level, TLRIs the propagation loss of the acoustic wave from the target to the receiver, DI directionalAnd (4) sexual energy.
The bistatic sonar equation is (SL + T)S-TLT)–(NL+TLRDI) ═ DT, where DT is the detection threshold.
Since the multi-base transmitter is non-directional, no consideration is given to the absorption attenuation of seawater medium, according to the formula TL (propagation loss) n × 10lg (r), n is equal to 2 when the acoustic energy is calculated as the spherical wave attenuation, and r is the diffusion distance. Meanwhile, the assumption that NL, DI and DT are constant in a certain range basically holds, the transmitter and the target are also determined, so SL and TSIs determined. Obtaining the following data according to a bistatic sonar equation: SL + TS-(NL-DI+DT)＝TLT+TLR＝20lg(Rs)+20lg(Rr)＝20lg(Rs×Rr) Is a determined value.
Definition of the invention
Representing the radius of equivalent tactical action of the bistatic buoy system on a specific target T, it is clear that the probability of detection for a specific target T
As the sonobuoys are randomly distributed, the detection frequency of the target T by the sonobuoy field is a Poisson random variable obeying the distribution density lambda of the sonobuoys, the probability Q that the random number increment in the Poisson random process is 0 (the probability that the target T is not detected by the sonobuoy field) is firstly calculated, and then the detection probability P (T) of the target T is calculated to be 1-Q.
I-type sonar buoy field
For a type II sonobuoy field, when the target T distance C _ SRBuoy is greater than gamma, the target T cannot be detected by the C _ SRBuoy. When the distance product of the target T and the distances C _ SRBuoy and RBuoy is more than gamma2In time, the sonobuoy field composed of C _ SRBuoy and RBuoy cannot detect the target. Type II only when the target T can not be detected by the C _ SRBuoy nor the sonobuoy field consisting of the C _ SRBuoy and the RBuoyThe sonobuoy field cannot detect the target, so
Since C _ SRBuoy and RBuoy are independent from each other,
step four: calculating the detection probability of the sonobuoy field to the target
Calculating the density of SBuoy, RBuoy and C _ SRBuoy according to the number of SBuoy, RBuoy and C _ SRBuoy
Wherein N iss、Nr、NsrThe distribution quantities of SBuoy, RBuoy and C _ SRBuoy are respectively. By using the mathematical property of independent stationary increment of the Poisson process, the lambda is converteds、λr、λsrSubstituting the parameters into (I) and (II) as poisson strength parameters, solving the detection failure probability of the I-type sonobuoy field and the II-type sonobuoy field to the target T:
type I sonobuoy field detection failure probability
Type II sonobuoy field failure probability
I sonar buoy field detection probability P (T) ═ 1-QI。
II sonobuoy field detection probability P (T) ═ 1-QII。
Compared with a Monte Carlo simulation method, the multi-base sonar buoy detection efficiency analysis method provided by the invention has the advantages of easiness in computer programming realization, quick calculation time, clear input and output relation, support for detection efficiency optimization, comparison efficiency and cost ratio and the like.
The invention calculates the detection probability of the target T by assuming that gamma is 5km, the area A of the target T is possible to be 500km2, 10 transmitters and 20 receivers, and the result shows that the proportion (namely the detection probability) of the optimal layout area A' of the sonobuoy to the target area A presents a convex function relationship, the detection probability is increased along with the increase of the average layout distance of the transmitters and the receivers, and when the average layout distance of the transmitters and the receivers is continuously increased, the bistatic sonobuoy gradually fails, and the detection probability begins to decrease. If the cost difference between the transmitter and the receiver is large, the number of the transmitters is reduced and the number of the receivers is increased in the sonobuoy field, the average distribution density of the sonobuoys can be increased under the condition that the cost is not increased or even reduced, the detection probability is not reduced or even increased, and the cost-effectiveness ratio is good.
A second embodiment of the present invention provides a computer-readable storage medium, which stores a signal-mapped computer program, and when the computer program is executed by at least one processor, the computer program implements the method for analyzing detection efficiency of a multi-base sonobuoy according to any one of the first embodiments of the present invention. The relevant content of the embodiments of the present invention can be understood by referring to the first embodiment of the present invention, and will not be discussed in detail herein.
Although the preferred embodiments of the present invention have been disclosed for illustrative purposes, those skilled in the art will appreciate that various modifications, additions and substitutions are possible, and the scope of the invention should not be limited to the embodiments described above.
Claims (10)
1. A multi-base sonar buoy detection efficiency analysis method is characterized by comprising the following steps:
identifying the type of sonobuoy field to be analyzed;
calculating the density of the sonobuoy field to be analyzed based on the type of the sonobuoy field;
determining the detection probability of the sonobuoy field to be analyzed according to the density of the sonobuoy field;
adjusting the type and layout of the sonobuoy in the sonobuoy field to be analyzed based on the calculated detection probability, the preset detection probability and the comprehensive performance of the sonobuoy so as to reduce the layout cost of the sonobuoy in the sonobuoy field under the condition of reaching the preset detection probability;
the sonobuoy field to be analyzed is a multi-base sonobuoy, and the sonobuoy comprises a transmitter SBuoy, a receiver RBuoy and a transmitting and receiving common-position machine C _ SRBuoy which are independently and randomly distributed.
2. The method of claim 1, wherein prior to identifying the sonobuoy field type to be analyzed, the method further comprises:
presetting an I-type sonobuoy field and a II-type sonobuoy field;
the I-type sonobuoy field is provided with an independent randomly distributed transmitter SBuoy and an independent randomly distributed receiver RBuoy;
the II-type sonobuoy field is an independent and randomly distributed transmitting and receiving common-position machine C _ SRBuoy and an independent and randomly distributed receiver RBuoy.
3. The method of claim 2, wherein the identifying a sonobuoy field type to be analyzed comprises:
and identifying the type of the sonobuoy field to which the sonobuoy field to be analyzed belongs based on the composition among the transmitter, the receiver and the transmitting and receiving co-location machine distributed in the sonobuoy field to be analyzed.
4. The method of claim 2, wherein when the sonobuoy field to be analyzed is a type I sonobuoy field, said calculating a sonobuoy field density of the sonobuoy field to be analyzed based on the sonobuoy field type comprises:
calculating the area A formed by all the sonobuoys in the sonobuoy field to be analyzed1Random distribution Density of SBuoy
Random placement density of RBuoy
The density of the sonobuoy field to be analyzed is: λ ═ λs,λr]；
Wherein Ns is the number of SBuoy, and Nr is the number of RBuoy.
5. The method of claim 4, wherein the area A of all sonobuoys in the sonobuoy field to be analyzed is determined1The method comprises the following steps:
solving the area of the convex polygon sonar buoy field according to a convex polygon area formula:
the first sonobuoy (x11, y11), the second sonobuoy (x12, y12), the first sonobuoy (x13, y13) … … nth sonobuoy (x1n, y1n) are the coordinates of each sonobuoy in the I-type sonobuoy field.
6. The method of claim 4, wherein determining a probability of detection of a sonobuoy field from a sonobuoy field density of the sonobuoy field to be analyzed comprises:
for a random target T in the sonobuoy field, the probability that there is no sonobuoy within any radius R, i.e. the probability that the distance R of a sonobuoy point from the random target T is greater than R, is the probability that a poisson random variable τ is 0:
where lambda is the random deployment density of the sonobuoy, then,
probability of no SBuoy sonobuoy within radius r of the random target T
The random meshProbability of no RBuoy sonobuoy within a standard T radius r
According to probability density function
To obtain
According to
Calculating the detection failure probability of the sonobuoy field:
further according to P (T) ═ 1-QICalculating the detection probability of the sonar buoy field;
wherein R issIs the distance, λ, of SBuoy from said random target TsRandomly laying density, R, for RBuoyrDistance, λ, of SBuoy from random target TrThe density is randomly laid for the RBuoy,
is the bistatic buoy system equivalent tactical action radius to the random target T.
7. The method of claim 2, wherein when the sonobuoy field to be analyzed is a type II sonobuoy field, said calculating a sonobuoy field density of the sonobuoy field to be analyzed based on the sonobuoy field type comprises:
calculating the area A formed by all the sonobuoys in the sonobuoy field to be analyzed2Random distribution density of C _ SRBuoy
Random placement density of RBuoy
The density of the sonobuoy field to be analyzed is: λ ═ λsr,λr]；
Wherein Nsr is the number of C _ SRBuoy and Nr is the number of RBuoy.
8. The method of claim 7, wherein the area A of all sonobuoys in the sonobuoy field to be analyzed is determined2The method comprises the following steps:
solving the area of the convex polygon sonar buoy field according to a convex polygon area formula:
wherein, the first sonobuoy (x21, y21), the second sonobuoy (x22, y22), the first sonobuoy (x23, y23) … … nth sonobuoy (x2n, y2n) are the coordinates of each sonobuoy in the type II sonobuoy field.
9. The method of claim 7, wherein determining a probability of detection of a sonobuoy field from a sonobuoy field density of the sonobuoy field to be analyzed comprises:
for a random target T in the sonobuoy field, the probability that there is no sonobuoy within any radius R, i.e. the probability that the distance R of a sonobuoy point from the random target T is greater than R, is the probability that a poisson random variable τ is 0:
where lambda is the random deployment density of the sonobuoy, then,
no RBuoy sonar within radius r of random target TProbability of buoy
Probability of no C _ SRBuoy sonobuoy within radius r of random target T
Because of the probability density function
Then
Will be provided with
Substitute QIIQ can be calculatedII
Further according to P (T) ═ 1-Q2Calculating the detection probability of the sonar buoy field;
wherein R issrIs the distance, λ, of C _ SRBuoy from the random target TsrIs the random distribution density of C _ SRBuoy, RrDistance, λ, of RBuoy from random target TrThe density is randomly laid for the RBuoy,
is the bistatic buoy system equivalent tactical action radius to the random target T.
10. A computer-readable storage medium storing a signal-mapped computer program which, when executed by at least one processor, implements the multi-base sonobuoy detection performance analysis method of any one of claims 1 to 9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110430578.0A CN113267767B (en) | 2021-04-21 | 2021-04-21 | Multi-base sonar buoy detection efficiency analysis method and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110430578.0A CN113267767B (en) | 2021-04-21 | 2021-04-21 | Multi-base sonar buoy detection efficiency analysis method and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113267767A true CN113267767A (en) | 2021-08-17 |
| CN113267767B CN113267767B (en) | 2023-04-07 |
Family
ID=77229116
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110430578.0A Active CN113267767B (en) | 2021-04-21 | 2021-04-21 | Multi-base sonar buoy detection efficiency analysis method and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113267767B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117849807A (en) * | 2024-03-06 | 2024-04-09 | 西北工业大学青岛研究院 | Method for optimizing tripwire sonar node layout of forward scattering detection |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN205608179U (en) * | 2016-05-13 | 2016-09-28 | 武汉大学 | A two standing posture UHF doppler radar systems for hydrologic monitoring |
| RU2713005C1 (en) * | 2019-04-24 | 2020-02-03 | Акционерное общество "Концерн "Центральный научно-исследовательский институт "Электроприбор" | Multi-static underwater surveillance system |
| CN112287547A (en) * | 2020-10-29 | 2021-01-29 | 西北工业大学 | Passive buoy array optimization method based on NSGA-II |
-
2021
- 2021-04-21 CN CN202110430578.0A patent/CN113267767B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN205608179U (en) * | 2016-05-13 | 2016-09-28 | 武汉大学 | A two standing posture UHF doppler radar systems for hydrologic monitoring |
| RU2713005C1 (en) * | 2019-04-24 | 2020-02-03 | Акционерное общество "Концерн "Центральный научно-исследовательский институт "Электроприбор" | Multi-static underwater surveillance system |
| CN112287547A (en) * | 2020-10-29 | 2021-01-29 | 西北工业大学 | Passive buoy array optimization method based on NSGA-II |
Non-Patent Citations (3)
| Title |
|---|
| 凌青等: "双(多)基地声纳浮标系统在反潜中的应用研究", 《海军工程大学学报》 * |
| 周旭等: "基于遗传算法的被动浮标阵优化布放技术研究", 《电子与信息学报》 * |
| 鞠建波等: "舰机协同下多基地声纳阵应召搜潜效能研究", 《计算机仿真》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117849807A (en) * | 2024-03-06 | 2024-04-09 | 西北工业大学青岛研究院 | Method for optimizing tripwire sonar node layout of forward scattering detection |
| CN117849807B (en) * | 2024-03-06 | 2024-05-10 | 西北工业大学青岛研究院 | Method for optimizing tripwire sonar node layout of forward scattering detection |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113267767B (en) | 2023-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH10503915A (en) | Method for wireless communication system planning | |
| RU2528556C1 (en) | Method of processing sonar echo signal | |
| KR102161246B1 (en) | Bistatic sonar system having bottom mounted hydrophones and method for source deployment based on ocean environment using the same | |
| CN105158756A (en) | Centralized MIMO radar radio frequency stealth multi-target tracking wave beam pointing method | |
| CN112162266B (en) | Conformal array two-dimensional beam optimization method based on convex optimization theory | |
| CN114781190B (en) | Radar detection capability simulation method and device | |
| CN117849807B (en) | Method for optimizing tripwire sonar node layout of forward scattering detection | |
| CN109655834A (en) | Multibeam sonar sounding method and system based on CFAR detection | |
| GB2496030A (en) | Monopulse radar | |
| CN110118956B (en) | Bistatic SAR anti-interference performance evaluation method for immune ground jammer area | |
| CN113267767B (en) | Multi-base sonar buoy detection efficiency analysis method and storage medium | |
| Chi | Underwater Real-Time 3D Acoustical Imaging: Theory, Algorithm and System Design | |
| RU2156984C1 (en) | Process of generation of information on noisy object at sea and process of obtainment of color scales for it | |
| CN113376612B (en) | Radar clutter generation method based on terrain matrixing and detection | |
| JP2009002921A (en) | Radar device | |
| US6861974B1 (en) | Clutter normalization by adaptation of transmit waveform | |
| CN109901131B (en) | Multipath utilization coherent beam forming method based on oblique projection | |
| Kerstens et al. | An optimized planar MIMO array approach to in-air synthetic aperture sonar | |
| Malyshkin | Experimental testing of the efficiency of fast projective adaptive algorithms | |
| Wang et al. | Passive distributed arrays cross-location algorithm with windowed deconvolution | |
| CN110244288A (en) | A kind of sonar array signal processing method based on focusing playback principle | |
| JP2008180541A (en) | Radar system | |
| CN118050690B (en) | Robust wide main lobe transmitting beam forming method | |
| Langli et al. | The first results with a new multibeam subbottom profiler | |
| CN109031316B (en) | Optimal detection frequency graphical estimation method under array size constraint |
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 |