CN103713289B - Based on the object detection method of distributed Phased-MIMO Combined Treatment - Google Patents
Based on the object detection method of distributed Phased-MIMO Combined Treatment Download PDFInfo
- Publication number
- CN103713289B CN103713289B CN201310695778.4A CN201310695778A CN103713289B CN 103713289 B CN103713289 B CN 103713289B CN 201310695778 A CN201310695778 A CN 201310695778A CN 103713289 B CN103713289 B CN 103713289B
- Authority
- CN
- China
- Prior art keywords
- transmitting
- subarray
- array
- signal
- receiving
- 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.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 10
- 238000003491 array Methods 0.000 claims description 19
- 238000012545 processing Methods 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 5
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 3
- 230000001427 coherent effect Effects 0.000 abstract description 3
- 230000015572 biosynthetic process Effects 0.000 abstract 1
- 238000001914 filtration Methods 0.000 abstract 1
- 230000001629 suppression Effects 0.000 abstract 1
- 230000000875 corresponding effect Effects 0.000 description 17
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/04—Systems determining presence of a target
-
- 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/52003—Techniques for enhancing spatial resolution of targets
-
- 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/52015—Diversity systems
-
- 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/521—Constructional features
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a kind of object detection method based on distributed Phased-MIMO Combined Treatment.The method is by being arranged in different positions in a distributed manner by multiple submatrix, the inner formation by launching beam of submatrix launches coherent signal acquisition coherent gain, between submatrix, the mutually orthogonal target diverse location that is irradiated to of signal obtains target diversity, and receiving array is formed by received beam and detects realization of goal with matched filtering.The inventive method can obtain relevant and incoherent gain, and eliminates interference, suppression reverberation, anti-target glint characteristic.In distributed object, when the beam angle of submatrix is greater than target size, detection perform is better than traditional phased array.
Description
Technical Field
The invention relates to an underwater target detection technology based on distributed Phased-MIMO combined processing, and belongs to the technical field of underwater sound target detection.
Background
When shallow sea target detection is performed, the dominant disturbance affecting active detection is reverberation. Both reverberation and target echo are caused by the transmitted signal, whose spectral characteristics are strongly correlated with the small doppler shift. The characteristics of reverberation are also related to propagation channels, and the reverberation has space-time variation characteristics and presents strong non-stationarity due to the influence of shallow sea underwater acoustic channels. Therefore, anti-reverberation is a difficult research point for detecting underwater targets. On the other hand, the spatial nonuniformity of the target also causes the flicker instability of the active sonar detection performance, and the interference of the emission array of the single base on the distributed target and the signal hardly ensures the stability of the detected target.
In the last decade, the research and development of Multiple Input Multiple Output (MIMO) technology in the radar field has been fast. The traditional phased array can acquire the diversity gain of the transmitting array by driving and transmitting signals, and effective focusing irradiation on a target is realized, so that the purposes of reverberation resistance and interference influence reduction are achieved. Meanwhile, the distributed MIMO system can obtain the transmitting waveform diversity and the target diversity, and is widely applied to target detection and positioning. However, due to the delay and doppler double spreading of the underwater acoustic channel, the MIMO technology develops slowly in the sonar field. In 2006, the i.bekkerman and j.tabrikaan first proposed a unified MIMO processing framework for radar and sonar. Considering the double extension of the underwater acoustic channel, and combining Phased array and MIMO technology, a distributed Phased-MIMO combined processing technology based on a broadband signal model is proposed.
Disclosure of Invention
The invention aims to provide a target detection method based on distributed Phased-MIMO combined processing, aiming at the problems of the existing single-base active sonar detection technology and considering reverberation interference and the spatial characteristics of distributed targets.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows: the invention relates to a target detection method based on distributed phase-MIMO combined processing, which comprises the following steps:
1) arranging M transmitting sub-arrays at a transmitting end, wherein each transmitting sub-array comprises MkA plurality of transmit transducer elements; each transmitting transducer array element in each transmitting subarray is connected with the output end of a power amplifier, and the input end of the power amplifier is connected with the output end of corresponding signal generating equipment; the transmitting transducer array elements in each transmitting sub-array are arranged into a uniform linear array, and the distance between the adjacent transmitting transducer array elements is smaller than or equal to the minimum value of half-wavelength of the transmitting signal of the transmitting sub-array; the spacing between each transmitting subarray meets the requirement shown in formula (1):
arranging N receiving sub-arrays at a receiving end, wherein each receiving sub-array comprises NlEach hydrophone array element; the hydrophone array elements in each receiving subarray are connected with corresponding signal acquisition equipment, the hydrophone array elements in each receiving subarray are arranged into uniform linear arrays, and the distance between every two adjacent hydrophone array elements is smaller than or equal to the minimum value of half wavelength of corresponding transmitting signals; the spacing between each receiving subarray satisfies the requirement shown in formula (2):
wherein M is a positive integer greater than or equal to 1, MkIs a positive integer of 1 or more andRtrepresenting the distance of the corresponding transmitting subarray of the transmitting end to the center of the target, D representing the size of the target, DtRepresenting the spacing between adjacent transmit sub-arrays, drRepresenting the spacing, λ, between adjacent receiving sub-arraysmaxRepresents the maximum of the wavelengths of all the transmitted signals; n represents a positive integer of 1 or more, NlRepresents a positive integer of 1 or more andRrrepresenting the distance from the corresponding receiving subarray of the receiving end to the target center;
2) each of the signal generating devices generates a transmission signalWhere E represents the total energy of the transmitted signal,
m denotes the number of transmit sub-arrays,
{Sk(fm) Denotes the transmit signal of the transmit sub-array k, | Sk(fm)||2=1/L,j ∈ {1, 2.,. M } and j ≠ k.
Emission signal of each of the signal generating devicesThe signals are amplified by a power amplifier and then transmitted to each transmitting transducer array element in a corresponding transmitting subarray to be converted into sound wave signals to be transmitted to a water area to be detected;
3) each receiving sub-array receives the transmitting signal S of the corresponding transmitting sub-arrayk(fm) Echo R oflk(fm) To a corresponding signal acquisition device which processes the received echo Rlk(fm) Performing receive beamforming as shown in equation (3):
in formula (3), (…)HIs a conjugate of the transposed symbol and,
will Ylk(fm) Matching is carried out according to the method shown in the formula (4):
in the formula (4), ∑ represents a summation symbol, | … | | non-volatile phosphor2Indicating the modulo of a vector, Sk(fm) Representing the transmitted signal of the transmitting subarray k, N being the number of receiving subarrays, M being the number of transmitting subarrays, and L being the number of frequencies of the signal transmitted by each transmitting subarray;
distance R from target center to all receiving sub-arrayslThe corresponding T value is plotted as Rl-T two-dimensional images; then observe Rl-T two-dimensional image if at the center of the objectA peak exists between the receiving subarrays indicating the presence of a target.
Compared with the prior art, the invention has the beneficial effects that: the traditional phased array is combined with a distributed MIMO system, so that the target diversity information is obtained while certain resolution is obtained to inhibit reverberation interference, the target flicker resistance is achieved, and the active sonar detection performance is stabilized. Numerical simulation and experiment verification show that the phase-MIMO active target detection technology is the traditional Phased array detection technology under the conditions that the target is a distributed target and the signal-to-noise ratio is high.
Drawings
FIG. 1 is a schematic layout of a transmit array and a receive array relative to a target;
FIG. 2(a) is a graph of sound speed as a function of depth; FIG. 2(b) is an environmental schematic of one embodiment of the present invention;
FIG. 3(a) is the detection output result of the Phased-MIMO active sonar system corresponding to the emission of the PCW signal; FIG. 3(b) is the detection output result of the Phased-Array active sonar system corresponding to the emission of the PCW signal;
FIG. 4(a) is the detection output result of the Phased-MIMO active sonar system corresponding to the emission of LFM signals; fig. 4(b) is the detection output result of the Phased-Array active sonar system corresponding to the emission of the LFM signal.
Detailed Description
The method comprises the following steps:
1) m transmitting sub-arrays are arranged at a transmitting end, wherein the serial number is 1. The k (k) ═ 1, 2.., M) th transmit subarray contains MkA transmitting transducer array, numbered 1,2k,MkIs a positive integer of 1 or more andRtdenote the distance from the emitting end to the center of the target, D denotes the size of the target, one would have M × MkThe output of the signal generating device of the circuit is connected with a circuit with M × MkThe input end of the power amplifier of the circuit and the output end of the power amplifier are connected to M × MkAnd a plurality of transmitting transducer elements. Spacing d between adjacent transmit transducer elements within the kth transmit sub-arraykThe minimum value of half wavelength of the transmitting signal of the kth transmitting sub-array is less than or equal to, and the distance between the adjacent transmitting sub-arrays is dtSatisfies the requirement shown in formula (1):
n receiving sub-arrays are arranged at a receiving end, wherein the number is 1, and N represents a positive integer which is larger than or equal to 1. The l (l ═ 1, 2.., N) th receiving subarray includes NlArray element of individual hydrophone, NlRepresents a positive integer of 1 or more andRrindicating the distance of the receiving end to the center of the target. The distance d between adjacent hydrophone array elements in the first receiving subarraylLess than or equal to the minimum of half wavelength of the signal received by the first receiving subarray, and the spacing between adjacent receiving subarrays is drSatisfies the requirement expressed by the formula (2):
the qth (q 1, 2.., M) of the k (k 1, 2., M) transmit subarraysk) The mth frequency point signal emitted by each emitting transducer array element isWherein, [ … ]]TIs a transposed symbol, θkIs the angle of arrival from the target center to the transmit subarray k, k being 1,2mIs the mth frequency of the signal transmitted by the transmitting subarray k, m is 1, 2.. L, L is the number of frequencies of the signal transmitted by each transmitting subarray, c is the propagation speed of the transmitted signal, d is the propagation speed of the transmitted signalkIs the spacing between adjacent transmit transducer elements within the transmit subarray k, | Sk(fm)||2=1/L,m=1,2,...L,j ∈ {1, 2.,. M } and j ≠ kWhereinFor the k-th transmitting sub-array at frequency fmA driving vector pointing to the center of the target.
The frequency of the signal transmitted by the kth transmitting subarray and received by the pth hydrophone array element in the ith receiving subarray is fmThe echo of (c) can be expressed as:
wherein,the equivalent scattering coefficient is the m frequency point emitted by the kth emission subarray after being irradiated on a scatterer;representing the propagation loss of the transmit sub-array to the target,respectively representing the transmission loss from the transmitting subarray to the target and the propagation loss from the target to the receiving subarray; written in vector form:
at frequency point f for target center to the l receiving subarraymThe response vector of (2);zero mean gaussian complex noise.
3) For each receiving subarray, the received transmission signal S of the corresponding transmitting subarrayk(fm) Echo R oflk(fm) To a corresponding signal acquisition device which processes the received echo Rlk(fm) Performing receive beamforming as shown in equation (3):
in formula (3), (…)HIs a conjugate of the transposed symbol and,
4) the processing data obtained in step 3) can be expressed as:
definition ofRepresenting the equivalent attenuation coefficient from the mth transmit sub-array to the jth receive sub-array, equation (7) can be rewritten as:
that is, the matrix form of equation (8) is as shown in equation (9):
wherein, MNL is multiplied by MNL dimension diagonal matrix Q is shown as formula (10):
Q=diag[M1N1s1(f1),...,M1N1s1(fL),
M2N1s2(fL+1),...,M2N1s2(fL+L),…,
MMN1sM(f(M-1)L+1),...,MMN1s2(fML),(10)
M1N2s1(f1),...,M1N2s1(fL),…,
MMN2sM(f(M-1)L+1),...,MMN2sM(fML),…,
MMNNsM(f(M-1)L+1),...,MMNNsM(f2L)]
the MNLx 1-dimensional column vector λ is as shown in equation (11):
MNL × 1-dimensional column vectorTo obey CN (0, sigma)2IMNL) Distributed complex gaussian noise.
4) Will Ylk(fm) Matching is carried out according to the method shown in the formula (4):
in the formula (4), ∑ represents a summation symbol, | … | | non-volatile phosphor2Indicating the modulo of a vector, Sk(fm) Representing the transmitted signal of transmit subarray k, N being the number of receive subarrays, M being the number of transmit subarrays, and L being the number of frequencies of the signal transmitted by each transmit subarray. When the received signal is judged to be present or absent, firstly, a threshold is determined, and if the received signal is matched with output energy Y2(namely the value T in the formula (4)) is greater than the threshold, judging that the target exists in the detection area; otherwise, judging that no target exists.
The MIMO detector based on NP criterion is given by the following equation (12):
according to the complex Gaussian noise and the statistical characteristics of the equivalent scattering coefficient, the statistic | Y | non-calculation amount is tested2Compliance
Wherein M is0=Mk,k=1,2,...,M,N0=Nl1, 2. Threshold according to false alarm probability Pfa(i.e., the probability of a false positive being present in the absence of a target) is selected, i.e., according to equation (14),
to give formula (15):
the detection probability P that can be obtained at this timed(i.e., the probability of correctly deciding that a target exists in the presence of a target) is as shown in equation (16):
the method of the present invention is specifically described below by taking the environment shown in fig. 2 as an example:
as shown in fig. 1,2 transmitting sub-arrays are arranged at the transmitting end, and are numbered as 1 and 2; each transmitting subarray comprises 3 transmitting transducer array elements which are numbered as 1,2 and 3; the distance between adjacent transmitting transducer array elements in each transmitting subarray is dk0.075m, the spacing between the transmitting subarrays is dt4 m. The output of an NI device with 6 channels is connected to a 6-channel power amplifier, the 6 output channels of which are connected to the inputs of 6 transmit transducer elements.
A linear target 1.5m long is located at 22m, satisfyingThe distributed condition of (2).
1 receiving subarray is arranged at a receiving end, the receiving subarray comprises 14 hydrophone array elements which are numbered as 1,2l=0.075m。
Signal sampling frequency fs50kHz, the signal duration T is 10ms, i.e. the signal time T is 0, 1/fs.
Scenario 1: signal frequency of f18kHz and f2In the single frequency Pulse (PCW) signal of 6kHz, 3 array elements in the transmitting sub-array 1 respectively transmit signals
3 array elements in the transmitting subarray 2 respectively transmit signals
Scenario 2: linear Frequency Modulation (LFM) signals with signal frequencies of 6-8 kHz and 8-10 kHz are transmitted, and 3 array elements in the sub-array 1 are transmitted at time toTransmitting signals separately
3 array elements in the transmitting subarray 2 respectively transmit signals
Wherein,representing taking the real part of the complex number, cos () being a cosine function, d1=0.075m,c≈1450m/s,θ1≈150I represents an imaginary number; d2=0.075m,θ2≈-150
The environment as shown in fig. 2 consists of three layers: a water medium layer, a sediment layer and a water bottom layer. Wherein the water depth is 21.53m, and the density is rho 1.0g/cm3The sound velocity profile is shown in fig. 2(a), and it can be seen that the sound velocity changes rapidly on the water surface, and has a relatively large negative gradient, and can be basically regarded as an equal sound velocity environment below 10m in water depth, and the energy of the transmitted signal irradiated on the target is relatively strong in the equal sound velocity environment, which is beneficial to improving the detection performance of the target. The deposited layer had a thickness of 7.343m and a density of 2.2g/cm3The speed of sound is cs1624.3 m/s. The water bottom is a uniform half space, and the density is rho ═ 3.0g/cm3The speed of sound is cs=1792.4m/s。
And the 14 channels of hydrophones convert the received echoes into electric signals and transmit the electric signals to signal acquisition equipment for acquisition and storage, and the electric signals are recorded as r (t). The received signal firstly obtains the receiving array gain through the beam former, so as to obtain x (t) ═ alpha (theta) r (t), and for the stationary target (namely neglecting the Doppler effect) due to the existence of the time delay Doppler double extension in the underwater acoustic channel, a copy coherent integration detector (RCI) is adopted to comprehensively consider the target echo returned due to multipath. Namely, it is
The presence or absence of the target is judged by comparing the output energy value of the corresponding position with the output energy value of no target around. Under the condition that the target position is unknown and the interference energy is small, the position corresponding to the maximum energy output value can be found through the technology to realize the positioning function of the target.
Conventional phased arrays, i.e., directional transmission to a target at one location and directional reception to a target at one location. Transmitting with one of 3 transmitting transducer elementsFor example, the subarray transmits 10ms of 8kHz single-frequency Pulse (PCW) signals or 10ms of 6-8 kHz Linear Frequency Modulation (LFM) signals, and the amplitude of the signals is the original amplitudeAnd (4) doubling. The receiving array is a single horizontal array with 14 paths, and the processing method is consistent with the Phased-MIMO array processing method.
As can be seen from fig. 3 and 4, two sonar systems can estimate the position of the target: the target real position is 22m, and the phase-MIMO system estimated value using the PCW signal isThe phase-Array system estimates asThe phase-MIMO system estimation value using LFM signal isThe phase-Array system estimates asUnder the condition of limited maximum instantaneous power, the pulse width of the transmitting signal is increased to improve the energy of the transmitting signal and the input signal-to-noise ratio of the receiver, thereby improving the detection performance of the receiver, and the increase of the pulse width simultaneously increases the ambiguity of distance estimation.LFM signal can better solve the contradiction, so that comparing fig. 3 and fig. 4, the matching result of LFM signal to the target position is shown to be larger than the resolution of the result of PCW signal, and simultaneously comparing the correlation matching output peak values of the two systems, the output peak value of Phased-MIMO system adopting PCW signal at the target position is shown to be 2.56 × 106Greater than the output peak of the Phased-Array system at the target position (1.034 × 10)6) The difference is obvious because the Phased-MIMO system acquires the diversity gain of the target.
Claims (1)
1. A target detection method based on distributed phase-MIMO combined processing is characterized by comprising the following steps:
1) arranging M transmitting sub-arrays at a transmitting end, wherein each transmitting sub-array comprises MkA plurality of transmit transducer elements; each transmitting transducer array element in each transmitting subarray is connected with the output end of a power amplifier, and the input end of the power amplifier is connected with the output end of corresponding signal generating equipment; the transmitting transducer array elements in each transmitting sub-array are arranged into a uniform linear array, and the distance between the adjacent transmitting transducer array elements is smaller than or equal to the transmitting sub-array elementsThe minimum value of half wavelength of the transmitting signal of the transmitting subarray; the spacing between each transmitting subarray meets the requirement shown in formula (1):
arranging N receiving sub-arrays at a receiving end, wherein each receiving sub-array comprises NlEach hydrophone array element; the hydrophone array elements in each receiving subarray are connected with corresponding signal acquisition equipment, the hydrophone array elements in each receiving subarray are arranged into uniform linear arrays, and the distance between every two adjacent hydrophone array elements is smaller than or equal to the minimum value of half wavelength of corresponding transmitting signals; the spacing between each receiving subarray satisfies the requirement shown in formula (2):
wherein M is a positive integer greater than or equal to 1, MkIs a positive integer of 1 or more andRtrepresenting the distance of the corresponding transmitting subarray of the transmitting end to the center of the target, D representing the size of the target, DtRepresenting the spacing between adjacent transmit sub-arrays, drRepresenting the spacing, λ, between adjacent receiving sub-arraysmaxRepresents the maximum of the wavelengths of all the transmitted signals; n represents a positive integer of 1 or more, NlRepresents a positive integer of 1 or more andRrrepresenting the distance from the corresponding receiving subarray of the receiving end to the target center;
2) each of the signal generating devices generates a transmission signalWhere E represents the total energy of the transmitted signal,
m denotes the number of transmit sub-arrays,
wherein, [ … ]]TIs a transposed symbol, θkIs the angle of arrival from the target center to the transmit subarray k, k being 1,2mIs the mth frequency of the signal transmitted by the transmitting subarray k, m is 1, 2.. L, L is the number of frequencies of the signal transmitted by each transmitting subarray, c is the propagation speed of the transmitted signal, d is the propagation speed of the transmitted signalkIs the spacing, M, between adjacent transmit transducer elements within a transmit sub-array kkThe number of transmitting transducer array elements contained in a transmitting subarray k, and i represents an imaginary number;
{Sk(fm) Denotes the transmit signal of the transmit sub-array k, | Sk(fm)||2=1/L,j ∈ {1, 2., M } with j ≠ k;
emission signal of each of the signal generating devicesThe signals are amplified by a power amplifier and then transmitted to each transmitting transducer array element in a corresponding transmitting subarray to be converted into sound wave signals to be transmitted to a water area to be detected;
3) each receiving sub-array receives the transmitting signal S of the corresponding transmitting sub-arrayk(fm) Echo R oflk(fm) To a corresponding signal acquisition device which processes the received echo Rlk(fm) Performing receive beamforming as shown in equation (3):
in formula (3), (…)HIs a conjugate of the transposed symbol and,
wherein, thetalIs the angle of arrival from the target center to the receiving subarray l, l is 1,2, …, N, dlIs the spacing, N, between adjacent hydrophone elements within the receiving sub-array llThe number of hydrophone array elements contained in a receiving subarray l, and i represents an imaginary number;
will Ylk(fm) Matching is carried out according to the method shown in the formula (4):
in the formula (4), ∑ represents a summation symbol, | … | | non-volatile phosphor2Indicating the modulo of a vector, Sk(fm) Representing the transmitted signal of the transmitting subarray k, N being the number of receiving subarrays, M being the number of transmitting subarrays, and L being the number of frequencies of the signal transmitted by each transmitting subarray;
distance R from target center to all receiving sub-arrayslThe corresponding T value is plotted as Rl-T two-dimensional images; then observe Rl-T two-dimensional image, indicating the presence of the target if there is a peak between the center of the target and the receiving subarray.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310695778.4A CN103713289B (en) | 2013-12-15 | 2013-12-15 | Based on the object detection method of distributed Phased-MIMO Combined Treatment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310695778.4A CN103713289B (en) | 2013-12-15 | 2013-12-15 | Based on the object detection method of distributed Phased-MIMO Combined Treatment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103713289A CN103713289A (en) | 2014-04-09 |
CN103713289B true CN103713289B (en) | 2016-03-23 |
Family
ID=50406414
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310695778.4A Expired - Fee Related CN103713289B (en) | 2013-12-15 | 2013-12-15 | Based on the object detection method of distributed Phased-MIMO Combined Treatment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103713289B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104535169B (en) * | 2014-12-03 | 2017-05-24 | 北京神州普惠科技股份有限公司 | Noise measurement apparatus based on fiber optic hydrophone array and measurement method thereof |
CN106324602B (en) * | 2016-07-29 | 2019-05-10 | 中国科学院声学研究所 | A kind of MIMO Sonar system |
CN106772326B (en) * | 2016-12-05 | 2019-04-30 | 中国人民解放军91388部队 | A kind of multiple submatrixes synthetic aperture sonar phase error analysis method |
CN110082741B (en) * | 2019-03-14 | 2022-07-15 | 哈尔滨工程大学 | Super-resolution angle of arrival estimation calculation method based on pseudo data reconstruction |
CN111505581B (en) * | 2020-04-16 | 2022-04-29 | 浙江大学 | Passive target detection method based on distributed sensor nodes |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288966A (en) * | 2011-05-12 | 2011-12-21 | 浙江大学 | MIMO (Multiple Input Multiple Output) detection method for forward small target |
CN102608609A (en) * | 2011-12-19 | 2012-07-25 | 浙江大学 | MIMO (multiple input multiple output) based high-resolution underwater target detection device and method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0719995D0 (en) * | 2007-10-12 | 2007-11-21 | Qinetiq Ltd | Radar method and apparatus suitable for use in multipath environments |
-
2013
- 2013-12-15 CN CN201310695778.4A patent/CN103713289B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288966A (en) * | 2011-05-12 | 2011-12-21 | 浙江大学 | MIMO (Multiple Input Multiple Output) detection method for forward small target |
CN102608609A (en) * | 2011-12-19 | 2012-07-25 | 浙江大学 | MIMO (multiple input multiple output) based high-resolution underwater target detection device and method |
Non-Patent Citations (1)
Title |
---|
多输入多输出阵列目标方位估计方法;蒋敏等;《***工程与电子技术》;20090731;第31卷(第07期);1602-1605 * |
Also Published As
Publication number | Publication date |
---|---|
CN103713289A (en) | 2014-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang | Phased-MIMO radar with frequency diversity for range-dependent beamforming | |
CN111198374B (en) | Doppler sensitive signal moving target underwater sound detection method based on space-time-frequency joint interference suppression | |
CN103713289B (en) | Based on the object detection method of distributed Phased-MIMO Combined Treatment | |
Frazer et al. | Multiple-input multiple-output over-the-horizon radar: experimental results | |
CN107505596B (en) | MIMO active detection signal design and detection system and method based on dual-extension underwater acoustic channel environment | |
CN103913725B (en) | Airborne radar Ground moving targets detection method under intensive repeating jamming environment | |
CN109597041B (en) | Segmented linear frequency modulation waveform design method based on coherent FDA | |
CN105301580A (en) | Passive detection method based on split array cross-spectrum phase difference variance weighing | |
Bekkerman et al. | Spatially coded signal model for active arrays | |
CN109061654B (en) | Single-circular-ring-array active three-dimensional positioning method in deep sea environment | |
Tabrikian et al. | Transmission diversity smoothing for multi-target localization [radar/sonar systems] | |
Ma et al. | Underwater broadband source detection using an acoustic vector sensor with an adaptive passive matched filter | |
Tabrikian | Barankin bounds for target localization by MIMO radars | |
Frazer et al. | Recent results in MIMO over-the-horizon radar | |
Rubinstein et al. | Frequency diverse array signal optimization: From non-cognitive to cognitive radar | |
CN108896974B (en) | Improved MIMO array high-resolution spatial spectrum estimation method | |
US8116169B2 (en) | Active sonar system and active sonar method using noise reduction techniques and advanced signal processing techniques | |
Wang et al. | Snapshot-deficient active target localization in beam-time domain using multi-frequency expectation-maximization algorithm | |
Pan et al. | TR-MIMO detection of a small target in a shallow water waveguide environment | |
Pan et al. | Distributed MIMO sonar for detection of moving targets in shallow sea environments | |
Pan et al. | Robust time-reversal is combined with distributed multiple-input multiple-output sonar for detection of small targets in shallow water environments | |
Jiang et al. | Research on target DOA estimation method using MIMO sonar | |
Das et al. | Analysis of effective signal design for active sensing of undersea objects/bottoms in tropical shallow waters | |
Xenaki et al. | Sparse MIMO synthetic aperture sonar processing with distributed optimization | |
Yan et al. | Bistatic MIMO sonar space-time adaptive processing based on knowledge-aided transform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160323 |