CN115774239B - Multi-target signal detection method based on non-uniform mutual mass array - Google Patents

Multi-target signal detection method based on non-uniform mutual mass array Download PDF

Info

Publication number
CN115774239B
CN115774239B CN202310103902.7A CN202310103902A CN115774239B CN 115774239 B CN115774239 B CN 115774239B CN 202310103902 A CN202310103902 A CN 202310103902A CN 115774239 B CN115774239 B CN 115774239B
Authority
CN
China
Prior art keywords
array
detection line
target signal
signal
target
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
CN202310103902.7A
Other languages
Chinese (zh)
Other versions
CN115774239A (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.)
Hangzhou Crysound Electronics Co Ltd
Original Assignee
Hangzhou Crysound Electronics Co Ltd
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 Hangzhou Crysound Electronics Co Ltd filed Critical Hangzhou Crysound Electronics Co Ltd
Priority to CN202310103902.7A priority Critical patent/CN115774239B/en
Publication of CN115774239A publication Critical patent/CN115774239A/en
Application granted granted Critical
Publication of CN115774239B publication Critical patent/CN115774239B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE 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/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Abstract

The invention relates to a multi-target signal detection method based on a non-uniform mutual mass array, which comprises the following steps: s1, detecting narrowband signals of different frequency bands in a sound source total signal; s2, if the number of the narrow-band signals is smaller than or equal to the number of the uniform linear arrays in the non-uniform mutual mass linear arrays, performing a step S3, otherwise, ending the detection; s3, adjusting the intervals among array elements in the plurality of detection line arrays based on the wavelength of the center frequency of each narrowband signal so that the main imaging frequency of the plurality of detection line arrays is consistent with the center frequency of the plurality of narrowband signals respectively; and S4, calculating one or more target signal angles detected by the detection line array at the main imaging frequency based on the sound source total signals received by the detection line array, the position coordinates of each array element in the detection line array, the main imaging frequency of the detection line array and the azimuth arrival angle estimation model. The invention introduces non-uniform mutual mass arrays capable of detecting a plurality of frequency points, and each uniform array can detect and obtain corresponding target signal angles at the main imaging frequency.

Description

Multi-target signal detection method based on non-uniform mutual mass array
Technical Field
The invention belongs to the technical field of multi-signal detection, and particularly relates to a multi-target signal detection method based on a non-uniform mutual mass array.
Background
In living and production environments, most of airtight pipelines or containers are leaked for various reasons, so that the leakage not only causes waste and economic loss of resources, but also causes serious safety problems due to leakage of harmful flammable and explosive gases or liquids, and therefore, many methods for effectively detecting the leakage have been applied for many years. When any pressure gas leaks or vacuum leaks, the gas flow forms turbulence, the turbulence emits sound wave energy, and sound waves cover the audible range and the ultrasonic frequency range. In a general environment, noisy ambient noise is very easy to mask the acoustic wave signal in the audible domain. The device will typically detect its ultrasonic frequency band.
Conventional detection devices are typically based on microphone arrays, detecting at ultrasonic frequency bands of no more than 50 kHz. Microphone arrays typically employ uniform line arrays and planar arrays, as well as non-uniform arrays, mostly sparse line arrays or planar arrays. The area array is generally a circular aperture array, a rectangular array of grids, a regular polygon array, or the like. Based on the microphone pickup array formed by the arrays, bandpass or high-pass filtering effects can be generated on specific frequency points or frequency bands, and the interference superposition is used for enhancing the directional and fixed-frequency pickup capability of the microphone pickup array, so that leakage abnormal sound is extracted. This has the advantage that a stronger and more accurate leak detection capability may be provided for a particular frequency bin of a single target. However, in the actual situation with multiple leakage targets and multiple leakage signal characteristic frequency points, the conventional detection means is used for detecting that the imaging of the system is interfered, and erroneous imaging is generated or misjudgment about whether leakage occurs is generated.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a multi-target signal detection method based on a non-uniform mutual mass array, which can realize multi-target leakage signal detection and imaging by utilizing incoherent multi-frequency points. The invention adopts the following technical scheme:
a multi-target signal detection method based on a non-uniform mutual mass array comprises the following steps:
s1, detecting a sound source total signal based on a non-uniform mutual linear array, and detecting narrowband signals of different frequency bands in the sound source total signal, wherein the non-uniform mutual linear array consists of a plurality of uniform linear arrays, and the number of array elements in each uniform linear array is recorded as
Figure SMS_1
,/>
Figure SMS_2
All of them are mutually equal in nature and are treated by>
Figure SMS_3
Representing the number of uniform linear arrays;
s2, if the number of the narrow-band signals is smaller than or equal to the number of the uniform linear arrays in the non-uniform mutual mass linear arrays, performing a step S3, otherwise ending detection;
s3, selecting a plurality of uniform linear arrays with the same quantity as the narrow-band signals from the non-uniform mutual linear arrays to serve as the linear arrays for detecting the multi-target signals, and adjusting each array element in the plurality of linear arrays based on the wavelength of the center frequency of each narrow-band signalThe distance between the two arrays is set so that the main imaging frequency of the arrays is consistent with the central frequency of the narrow-band signals, and the distance between array elements in each array is recorded as
Figure SMS_4
And->
Figure SMS_5
All of them are mutually equal in nature and are treated by>
Figure SMS_6
Representing the number of the detection line arrays; />
S4, establishing an azimuth arrival angle estimation model, and calculating one or more target signal angles detected by the detection line array at the main imaging frequency of the detection line array according to the total sound source signals received by the detection line array, the position coordinates of each array element in the detection line array, the main imaging frequency of the detection line array and the azimuth arrival angle estimation model;
s5, detecting the multi-target signal based on the target signal angles detected by the detection line arrays at the main imaging frequency of the detection line arrays.
In a preferred embodiment, in step S4, the azimuth angle of arrival estimation model is specifically:
suppose for the first
Figure SMS_7
The incidence angles of the target signals detected by the detection line arrays are respectively
Figure SMS_8
,/>
Figure SMS_9
Indicate->
Figure SMS_10
The detection of the individual detection lines +.>
Figure SMS_11
Incidence angle of the individual target signals, +.>
Figure SMS_12
The number of target signals is represented by:
Figure SMS_13
wherein ,
Figure SMS_14
indicate->
Figure SMS_15
The total signal of sound sources detected by the individual detection line arrays, < >>
Figure SMS_16
Vector matrix representing each target signal, +.>
Figure SMS_17
Indicate->
Figure SMS_18
The phase difference matrix of each target signal in each detection line array reaching each array element is +.>
Figure SMS_19
Indicate->
Figure SMS_20
Complex-vector gaussian white noise detected by the detection line arrays;
Figure SMS_21
wherein ,
Figure SMS_22
indicate->
Figure SMS_23
The first line of detection>
Figure SMS_24
A phase difference matrix of the target signals reaching each array element based on the +.>
Figure SMS_25
Position coordinates of each array element in each detection line array +.>
Figure SMS_26
Main imaging frequencies of the detection line arrays are obtained through calculation;
for the first
Figure SMS_27
Covariance matrix of cross matrix is calculated by sound source total signal detected by each detection line array>
Figure SMS_28
Figure SMS_29
wherein ,
Figure SMS_30
representing the desired value>
Figure SMS_31
Represents the conjugate transpose->
Figure SMS_32
,/>
Figure SMS_33
Representing noise signal power, < >>
Figure SMS_34
Representing the identity matrix;
for a pair of
Figure SMS_35
Performing eigenvalue decomposition to obtain:
Figure SMS_36
wherein ,
Figure SMS_37
representing feature vectors +_>
Figure SMS_38
Representing the feature vector;
will be
Figure SMS_39
The signal space and the noise space are decomposed to obtain:
Figure SMS_40
wherein ,
Figure SMS_41
representing feature vectors in signal space, +.>
Figure SMS_42
Representing feature vectors in signal space, +.>
Figure SMS_43
Representing feature vectors in noise space, +.>
Figure SMS_44
Representing feature vectors in noise space; />
Based on the formula
Figure SMS_45
Calculating to obtain the qualified +.>
Figure SMS_46
Eligible->
Figure SMS_47
Namely +.>
Figure SMS_48
The individual linear arrays detect the resulting target signal angle at their primary imaging frequency.
As a preferable scheme:
Figure SMS_49
wherein ,
Figure SMS_50
indicate->
Figure SMS_51
Vector of individual target signals,/>
Figure SMS_52
Representing the transpose.
As a preferable scheme:
Figure SMS_53
wherein ,
Figure SMS_56
indicate->
Figure SMS_60
The first line of detection>
Figure SMS_64
The target signal reaches the->
Figure SMS_55
Phase difference between each element and the first element is reached,/->
Figure SMS_59
Indicate->
Figure SMS_63
The>
Figure SMS_66
Position coordinates of individual array elements,/->
Figure SMS_54
Indicate->
Figure SMS_58
Personal examinationThe number of array elements in the linear array, +.>
Figure SMS_61
,/>
Figure SMS_65
Representing sound speed,/->
Figure SMS_57
Indicate->
Figure SMS_62
The main imaging frequency of each of the linear arrays.
In step S4, for any one of the detection line arrays, one or more target signal angles detected by the detection line array at a plurality of secondary primary imaging frequencies thereof are calculated based on the total sound source signal received by the detection line array, the position coordinates of each array element in the detection line array, a plurality of secondary primary imaging frequencies of the detection line array, and an azimuth arrival angle estimation model, wherein the plurality of secondary primary imaging frequencies are all within a preset range up and down with respect to the center frequency of the narrowband signal corresponding to the detection line array.
Preferably, between step S4 and step S5, step B is further included:
B. for the same detection line array, filtering target signals received at a plurality of imaging frequencies and at different target signal angles, wherein the filtering process comprises the following steps:
b1, calculating the intensity of target signals of the detection line array at one or more target signal angles received at a main imaging frequency, and the intensity of target signals of the detection line array at one or more target signal angles received at a secondary main imaging frequency;
and B2, if the difference value between the intensity of the target signal at the corresponding target signal angle received at the secondary main imaging frequency and the intensity of the strongest target signal received at the main imaging frequency is larger than a preset sound intensity difference threshold value, filtering the target signal at the corresponding target signal angle received at the secondary main imaging frequency.
Preferably, the intensity calculation formula of the target signal is:
Figure SMS_67
wherein ,
Figure SMS_68
is expressed in the imaging frequency +.>
Figure SMS_69
The angle of the target signal received at the receiver is +.>
Figure SMS_70
The strength of the target signal at that location.
Preferably, step a is further included between step S4 and step S5:
A. judging whether a difference value between two target signal angles is smaller than a preset threshold value degree or not among the plurality of target signal angles, if so, considering that the two target signal angles correspond to the same target signal, and calculating the target signal based on the average value of the two target signal angles to obtain the final corresponding target signal angle.
Preferably, the preset threshold degree is 5 °.
In step S3, the spacing between the array elements in the detection line array is half the wavelength of the center frequency of the narrowband signal corresponding to the spacing.
The beneficial effects of the invention are as follows:
and a non-uniform mutual mass array capable of detecting a plurality of frequency points is introduced, and compared with a traditional linear array and an area array of a microphone, each uniform array in the non-uniform mutual mass array can be detected at a main imaging frequency to obtain a corresponding target signal angle, so that estimation errors of arrival angles caused by a plurality of signal source gas leakage points can not occur.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a multi-target signal detection method based on a non-uniform mutual mass array according to the present invention;
FIG. 2 is a schematic diagram of a non-uniform mutual mass array according to the present invention.
Detailed Description
The following specific examples are presented to illustrate the present invention, and those skilled in the art will readily appreciate the additional advantages and capabilities of the present invention as disclosed herein. The invention may be practiced or carried out in other embodiments that depart from the specific details, and the details of the present description may be modified or varied from the spirit and scope of the present invention. It should be noted that the following embodiments and features in the embodiments may be combined with each other without conflict.
Referring to fig. 1, the embodiment provides a multi-target signal detection method based on a non-uniform mutual mass array, which includes the steps of:
s1, detecting a sound source total signal based on a non-uniform mutual linear array, and detecting narrowband signals of different frequency bands in the sound source total signal, wherein the non-uniform mutual linear array consists of a plurality of uniform linear arrays, and the number of array elements in each uniform linear array is recorded as
Figure SMS_71
,/>
Figure SMS_72
All of them are mutually equal in nature and are treated by>
Figure SMS_73
Representing the number of uniform linear arrays.
S2, if the number of the narrow-band signals is smaller than or equal to the number of the uniform linear arrays in the non-uniform mutual mass linear arrays, performing step S3, otherwise ending detection.
Because the number of the narrow-band signals is larger than the number of the uniform linear arrays in the non-uniform mutual linear array and exceeds the number of main imaging frequencies which can be provided by the non-uniform mutual linear array, the final multi-target signal detection precision can be reduced, and therefore the situation that the number of the narrow-band signals is larger than the number of the uniform linear arrays in the non-uniform mutual linear array is not considered in the embodiment.
S3, randomly selecting a plurality of uniform linear arrays with the same number as the narrow-band signals from the non-uniform mutual linear arrays to serve as the linear arrays for detecting the multi-target signals, and adjusting the intervals among the array elements in the plurality of the linear arrays based on the wavelength of the central frequency of each narrow-band signal (the adjustment of the positions of the array elements is needed to be explained here can be realized through a position adjusting device) so that the main imaging frequencies of the plurality of the linear arrays are respectively consistent with the central frequency of the plurality of narrow-band signals, wherein the intervals among the array elements in the plurality of the linear arrays are recorded as
Figure SMS_74
And->
Figure SMS_75
All of them are mutually equal in nature and are treated by>
Figure SMS_76
Representing the number of the detection lines.
The main imaging frequencies of the plurality of the detection line arrays are respectively consistent with the center frequencies of the plurality of the narrowband signals, and the intervals among the array elements in the detection line arrays are half of the wavelengths of the center frequencies of the narrowband signals corresponding to the intervals.
Referring to fig. 2, a non-uniform mutual linear array formed by two uniform linear arrays is shown, and because the number of array elements and the spacing between the array elements are all of the same mass, when the array elements are on a straight line, only the first array element of the two uniform linear arrays can be overlapped.
S4, an azimuth arrival angle estimation model is established, and for any detection line array, one or more target signal angles (to be noted, one or more target signal angles detected at the main imaging frequency may exist in the target signal angles detected at the main imaging frequency) of the detection line array are obtained through calculation based on the total sound source signals received by the detection line array, the position coordinates of each array element in the detection line array, the main imaging frequency of the detection line array and the azimuth arrival angle estimation model (to be explained specifically below).
S5, detecting the multi-target signal based on the target signal angles detected by the detection line arrays at the main imaging frequency, and performing final imaging.
It should be noted that, in this embodiment, the number of target signals is generally less than 4, otherwise, the scale of the heterogeneous linear array will be too large, and it is difficult to obtain a realistic array.
Therefore, the invention introduces the non-uniform mutual mass array capable of detecting a plurality of frequency points, and compared with the traditional linear array and the area array of the microphone, each uniform array in the non-uniform mutual mass array can detect and obtain a corresponding target signal angle at the main imaging frequency, and the estimation error of the arrival angle caused by a plurality of signal source gas leakage points can not occur.
Specifically:
when in far field
Figure SMS_77
The target signals have different frequency bands +.>
Figure SMS_78
The narrow-band signals are respectively->
Figure SMS_79
The receiving of the detection line array is realized by different line array gating when continuous signal sampling is carried out at equal time intervals, and the time intervals are +.>
Figure SMS_80
A certain group of uniform linear arrays is strobed in turn.
In step S4, the azimuth angle of arrival estimation model specifically includes:
suppose for the first
Figure SMS_81
The incidence angles of the target signals detected by the detection line arrays are respectively
Figure SMS_82
,/>
Figure SMS_83
Indicate->
Figure SMS_84
The detection of the individual detection lines +.>
Figure SMS_85
Incidence angle of the individual target signals, +.>
Figure SMS_86
The number of target signals is represented by:
Figure SMS_87
wherein ,
Figure SMS_88
indicate->
Figure SMS_89
The total signal of sound sources detected by the individual detection line arrays, < >>
Figure SMS_90
Vector matrix representing each target signal, +.>
Figure SMS_91
Indicate->
Figure SMS_92
The phase difference matrix of each target signal in each detection line array reaching each array element is +.>
Figure SMS_93
Indicate->
Figure SMS_94
Complex vector gaussian white noise detected by each detection line array.
Figure SMS_95
wherein ,
Figure SMS_96
indicate->
Figure SMS_97
Vector of individual target signals,/>
Figure SMS_98
Representing the transpose. />
Figure SMS_99
wherein ,
Figure SMS_100
indicate->
Figure SMS_101
The first line of detection>
Figure SMS_102
A phase difference matrix of the target signals reaching each array element based on the +.>
Figure SMS_103
Position coordinates of each array element in each detection line array +.>
Figure SMS_104
And calculating the main imaging frequency of each detection line array.
Figure SMS_105
wherein ,
Figure SMS_108
indicate->
Figure SMS_113
The first line of detection>
Figure SMS_116
The target signal reaches the->
Figure SMS_109
Phase difference between each element and the first element is reached,/->
Figure SMS_112
Indicate->
Figure SMS_115
The>
Figure SMS_118
Position coordinates of individual array elements,/->
Figure SMS_106
Indicate->
Figure SMS_110
The number of array elements in each detection line array, +.>
Figure SMS_114
,/>
Figure SMS_117
Representing the sound velocity, 340m/s in this embodiment, a more accurate value may be used in the actual value, and +.>
Figure SMS_107
Indicate->
Figure SMS_111
The main imaging frequency of each of the linear arrays.
For the first
Figure SMS_119
Covariance matrix of cross matrix is calculated by sound source total signal detected by each detection line array>
Figure SMS_120
Figure SMS_121
wherein ,
Figure SMS_122
representing the desired value>
Figure SMS_123
Represents the conjugate transpose->
Figure SMS_124
,/>
Figure SMS_125
Representing noise signal power, < >>
Figure SMS_126
Representing the identity matrix.
For a pair of
Figure SMS_127
Performing eigenvalue decomposition to obtain:
Figure SMS_128
wherein ,
Figure SMS_129
representing feature vectors +_>
Figure SMS_130
Representing the feature vector.
Will be
Figure SMS_131
The signal space and the noise space are decomposed to obtain:
Figure SMS_132
wherein ,
Figure SMS_133
representing feature vectors in signal space, +.>
Figure SMS_134
Representing feature vectors in signal space, +.>
Figure SMS_135
Representing feature vectors in noise space, +.>
Figure SMS_136
Representing the feature vector in noise space.
Based on the formula
Figure SMS_137
Calculating to obtain the qualified +.>
Figure SMS_138
Eligible->
Figure SMS_139
Namely +.>
Figure SMS_140
The individual linear arrays detect the resulting target signal angle at their primary imaging frequency.
Here, the formula is satisfied
Figure SMS_141
There may be a plurality of target signal angles and thus there may be a plurality of target signal angles detected at the primary imaging frequency.
In the actual operation process, the formula is as follows
Figure SMS_142
Can be adjusted to->
Figure SMS_143
Only the formula->
Figure SMS_144
As far as possible, 0, a range can be given in which the angles are all solved.
And calculating target signal angles obtained by detecting the other detection line arrays at the respective main imaging frequencies sequentially based on the azimuth arrival angle estimation model, thus finishing multi-target signal detection and carrying out final imaging.
In this embodiment, in step S4, for any one of the detection line arrays, one or more target signal angles (similarly, a plurality of target signal angles may be detected at a secondary primary imaging frequency) detected by the detection line array at a plurality of secondary primary imaging frequencies are also calculated based on the total sound source signal received by the detection line array, the position coordinates of each array element in the detection line array, the plurality of secondary primary imaging frequencies of the detection line array, and the azimuth arrival angle estimation model, where the plurality of secondary primary imaging frequencies are all within a preset range above and below the center frequency of the narrowband signal corresponding to the detection line array.
That is, in this embodiment, for any one of the detection line arrays, a plurality of secondary primary imaging frequencies are also taken within the upper and lower preset ranges of the primary imaging frequencies, and the calculation of the azimuth arrival angle estimation model is repeatedly performed to calculate the target signal angles detected by the detection line array at the plurality of secondary primary imaging frequencies (it should be noted that, when the secondary primary imaging frequencies are adopted for calculation, the above formula is given in the following manner
Figure SMS_145
And (3) replacing the value of the (c) and other calculation steps are consistent). The beneficial effects of this operation procedure are: image frequency points are added, and errors caused by physical positions of the array are eliminated.
After the step S4 is performed, a plurality of target signal angles may be obtained, and when the target signal angles corresponding to the plurality of target signals have a certain large angular resolution (generally greater than or equal to 5 °), multi-target multi-frequency point imaging may be achieved. If the signal is smaller than 5 °, the two target signals are considered to be the same target signal, so step a is further included between step S4 and step S5:
A. judging whether a difference value between two target signal angles is smaller than a preset threshold value degree or not among the plurality of target signal angles, if so, considering that the two target signal angles correspond to the same target signal, and calculating the target signal based on the average value of the two target signal angles to obtain the final corresponding target signal angle.
Finally, between step S4 and step S5, step B is further included:
B. for the same detection line array, filtering target signals received at a plurality of imaging frequencies and at different target signal angles, wherein the filtering process comprises the following steps:
b1, calculating the intensity of target signals of the detection line array at one or more target signal angles received at a main imaging frequency, and the intensity of target signals of the detection line array at one or more target signal angles received at a secondary main imaging frequency;
and B2, if the difference value between the intensity of the target signal at the corresponding target signal angle received at the secondary main imaging frequency and the intensity of the strongest target signal received at the main imaging frequency is larger than a preset sound intensity difference threshold (the sound intensity difference threshold is set to be 12dB in the embodiment), filtering the target signal at the corresponding target signal angle received at the secondary main imaging frequency.
The beneficial effects of this operation procedure are: interference noise is eliminated, so that the actual sound source localization is more obvious.
The intensity calculation formula of the target signal is:
Figure SMS_146
,/>
wherein ,
Figure SMS_147
is expressed in the imaging frequency +.>
Figure SMS_148
The angle of the target signal received at the receiver is +.>
Figure SMS_149
The strength of the target signal at that location.
The above examples are merely illustrative of the preferred embodiments of the present invention and are not intended to limit the scope of the present invention, and various modifications and improvements made by those skilled in the art to the technical solution of the present invention should fall within the protection scope of the present invention without departing from the design spirit of the present invention.

Claims (10)

1. A multi-target signal detection method based on a non-uniform mutual mass array is characterized by comprising the following steps:
s1, detecting a sound source total signal based on a non-uniform mutual linear array, and detecting narrowband signals of different frequency bands in the sound source total signal, wherein the non-uniform mutual linear array consists of a plurality of uniform linear arrays, and the number of array elements in each uniform linear array is recorded as
Figure QLYQS_1
,/>
Figure QLYQS_2
All of them are mutually equal in nature and are treated by>
Figure QLYQS_3
Representing the number of uniform linear arrays;
s2, if the number of the narrow-band signals is smaller than or equal to the number of the uniform linear arrays in the non-uniform mutual mass linear arrays, performing a step S3, otherwise ending detection;
s3, selecting a plurality of uniform linear arrays with the same number as the narrow-band signals from the non-uniform mutual linear arrays to serve as the linear arrays for detecting the multi-target signals, and adjusting the intervals among the array elements in the plurality of the linear arrays based on the wavelength of the central frequency of each narrow-band signal, so that the main imaging frequency of the plurality of the linear arrays is respectively consistent with the central frequency of the plurality of the narrow-band signals, and the array element intervals in the plurality of the linear arrays are recorded as
Figure QLYQS_4
And->
Figure QLYQS_5
All betweenMutual mass (I/II)>
Figure QLYQS_6
Representing the number of the detection line arrays;
s4, establishing an azimuth arrival angle estimation model, and calculating one or more target signal angles detected by the detection line array at the main imaging frequency of the detection line array according to the total sound source signals received by the detection line array, the position coordinates of each array element in the detection line array, the main imaging frequency of the detection line array and the azimuth arrival angle estimation model;
s5, detecting the multi-target signal based on the target signal angles detected by the detection line arrays at the main imaging frequency of the detection line arrays.
2. The method for detecting multi-target signals based on the non-uniform mutual mass array according to claim 1, wherein in step S4, the azimuth angle of arrival estimation model is specifically:
suppose for the first
Figure QLYQS_7
The incidence angles of the target signals detected by the detection line arrays are respectively
Figure QLYQS_8
,/>
Figure QLYQS_9
Indicate->
Figure QLYQS_10
The detection of the individual detection lines +.>
Figure QLYQS_11
Incidence angle of the individual target signals, +.>
Figure QLYQS_12
The number of target signals is represented by:
Figure QLYQS_13
wherein ,
Figure QLYQS_14
indicate->
Figure QLYQS_15
The total signal of sound sources detected by the individual detection line arrays, < >>
Figure QLYQS_16
Vector matrix representing each target signal, +.>
Figure QLYQS_17
Indicate->
Figure QLYQS_18
The phase difference matrix of each target signal in each detection line array reaching each array element is +.>
Figure QLYQS_19
Represent the first
Figure QLYQS_20
Complex-vector gaussian white noise detected by the detection line arrays;
Figure QLYQS_21
wherein ,
Figure QLYQS_22
indicate->
Figure QLYQS_23
The first line of detection>
Figure QLYQS_24
A phase difference matrix of the target signals reaching each array element based on the +.>
Figure QLYQS_25
Position coordinates of each array element in each detection line array +.>
Figure QLYQS_26
Main imaging frequencies of the detection line arrays are obtained through calculation;
for the first
Figure QLYQS_27
Covariance matrix of cross matrix is calculated by sound source total signal detected by each detection line array>
Figure QLYQS_28
Figure QLYQS_29
,/>
wherein ,
Figure QLYQS_30
representing the desired value>
Figure QLYQS_31
Represents the conjugate transpose->
Figure QLYQS_32
,/>
Figure QLYQS_33
Representing noise signal power, < >>
Figure QLYQS_34
Representing the identity matrix;
for a pair of
Figure QLYQS_35
Performing eigenvalue decomposition to obtain:
Figure QLYQS_36
wherein ,
Figure QLYQS_37
feature vector representing total signal of sound source, +.>
Figure QLYQS_38
A diagonal matrix representing the composition of the eigenvalues;
will be
Figure QLYQS_39
The signal space and the noise space are decomposed to obtain:
Figure QLYQS_40
wherein ,
Figure QLYQS_41
representing feature vectors in signal space, +.>
Figure QLYQS_42
Representing a diagonal matrix under signal space, +.>
Figure QLYQS_43
Representing feature vectors in noise space, +.>
Figure QLYQS_44
Representing a diagonal matrix under noise space;
based on the formula
Figure QLYQS_45
Calculating to obtain the qualified +.>
Figure QLYQS_46
Eligible->
Figure QLYQS_47
Namely +.>
Figure QLYQS_48
The individual linear arrays detect the resulting target signal angle at their primary imaging frequency.
3. The multi-target signal detection method based on the non-uniform mutual mass array according to claim 2, wherein the method is characterized in that:
Figure QLYQS_49
wherein ,
Figure QLYQS_50
indicate->
Figure QLYQS_51
Vector of individual target signals,/>
Figure QLYQS_52
Representing the transpose.
4. The multi-target signal detection method based on the non-uniform mutual mass array according to claim 2, wherein the method is characterized in that:
Figure QLYQS_53
wherein ,
Figure QLYQS_57
indicate->
Figure QLYQS_59
The first line of detection>
Figure QLYQS_63
The target signal reaches the->
Figure QLYQS_56
Phase difference between each element and the first element is reached,/->
Figure QLYQS_60
Indicate->
Figure QLYQS_64
The>
Figure QLYQS_66
Position coordinates of individual array elements,/->
Figure QLYQS_54
Indicate->
Figure QLYQS_58
The number of array elements in each detection line array, +.>
Figure QLYQS_62
,/>
Figure QLYQS_65
Representing sound speed,/->
Figure QLYQS_55
Indicate->
Figure QLYQS_61
The principal imaging frequency of each of the linear arrays, j, represents an imaginary unit.
5. The method for detecting multiple target signals based on non-uniform mutual mass arrays according to claim 4, wherein in step S4, for any one of the detection line arrays, one or more target signal angles detected by the detection line array at a plurality of secondary primary imaging frequencies thereof are calculated based on a sound source total signal received by the detection line array, position coordinates of each array element in the detection line array, a plurality of secondary primary imaging frequencies of the detection line array and an azimuth arrival angle estimation model, and the plurality of secondary primary imaging frequencies are all within a preset range up and down with respect to a center frequency of a narrowband signal corresponding to the detection line array.
6. The method for detecting multiple target signals based on non-uniform mutual mass array according to claim 5, wherein between step S4 and step S5, further comprising step B:
B. for the same detection line array, filtering target signals received at a plurality of imaging frequencies and at different target signal angles, wherein the filtering process comprises the following steps:
b1, calculating the intensity of target signals of the detection line array at one or more target signal angles received at a main imaging frequency, and the intensity of target signals of the detection line array at one or more target signal angles received at a secondary main imaging frequency;
and B2, if the difference value between the intensity of the target signal at the corresponding target signal angle received at the secondary main imaging frequency and the intensity of the strongest target signal received at the main imaging frequency is larger than a preset sound intensity difference threshold value, filtering the target signal at the corresponding target signal angle received at the secondary main imaging frequency.
7. The method for detecting multiple target signals based on the non-uniform mutual mass array according to claim 6, wherein the intensity calculation formula of the target signals is:
Figure QLYQS_67
wherein ,
Figure QLYQS_68
is expressed in the imaging frequency +.>
Figure QLYQS_69
The angle of the target signal received at the receiver is +.>
Figure QLYQS_70
The strength of the target signal at that location.
8. The method for detecting multiple target signals based on non-uniform mutual mass array according to claim 1, wherein step a is further included between step S4 and step S5:
A. judging whether a difference value between two target signal angles is smaller than a preset threshold value degree or not among the plurality of target signal angles, if so, considering that the two target signal angles correspond to the same target signal, and calculating the target signal based on the average value of the two target signal angles to obtain the final corresponding target signal angle.
9. The method for detecting multiple target signals based on a non-uniform mutual mass array according to claim 8, wherein the preset threshold degree is 5 °.
10. The method for detecting multiple target signals based on non-uniform mutual mass arrays according to claim 1, wherein in step S3, the spacing between each array element in the detection line array is half the wavelength of the center frequency of the narrowband signal corresponding to the spacing.
CN202310103902.7A 2023-02-13 2023-02-13 Multi-target signal detection method based on non-uniform mutual mass array Active CN115774239B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310103902.7A CN115774239B (en) 2023-02-13 2023-02-13 Multi-target signal detection method based on non-uniform mutual mass array

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310103902.7A CN115774239B (en) 2023-02-13 2023-02-13 Multi-target signal detection method based on non-uniform mutual mass array

Publications (2)

Publication Number Publication Date
CN115774239A CN115774239A (en) 2023-03-10
CN115774239B true CN115774239B (en) 2023-05-12

Family

ID=85393596

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310103902.7A Active CN115774239B (en) 2023-02-13 2023-02-13 Multi-target signal detection method based on non-uniform mutual mass array

Country Status (1)

Country Link
CN (1) CN115774239B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106324558A (en) * 2016-08-30 2017-01-11 东北大学秦皇岛分校 Broadband signal DOA estimation method based on co-prime array

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106324558A (en) * 2016-08-30 2017-01-11 东北大学秦皇岛分校 Broadband signal DOA estimation method based on co-prime array

Also Published As

Publication number Publication date
CN115774239A (en) 2023-03-10

Similar Documents

Publication Publication Date Title
AU2022259835B2 (en) Direction of arrival estimation
CN110764053B (en) Multi-target passive positioning method based on underwater sensor network
CN107272005B (en) Active positioning method based on target echo arrival time delay and arrival angle under reliable acoustic path
CN108872970B (en) Grating lobe discrimination method suitable for general equidistant sparse array single-frequency signal beam forming
CN113884986B (en) Beam focusing enhanced strong impact signal space-time domain joint detection method and system
CN109581317B (en) Corner target positioning method based on echo peak matching
CN113868583B (en) Method and system for calculating sound source distance focused by subarray wave beams
CN104502904A (en) Torpedo homing beam sharpening method
US7952513B2 (en) Counter target acquisition radar and acoustic adjunct for classification
CN109001297A (en) Large sample underwater acoustic materials acoustical reflection factor measurement method based on single vector hydrophone
CN115774239B (en) Multi-target signal detection method based on non-uniform mutual mass array
Ma et al. An underwater bistatic positioning system based on an acoustic vector sensor and experimental investigation
CN109100711A (en) Active sonar low operand 3-D positioning method in single base under a kind of deep-marine-environment
CN109061654A (en) Single ring array active 3-D positioning method under a kind of deep-marine-environment
RU2723145C1 (en) Method and device for detecting noisy objects in the sea with onboard antenna
EP2317335A1 (en) Improved beamforming method for analysing signals received by a transducer arrray, and relative detection system
CN108088547A (en) A kind of weak target passive detection method based on small-bore two-dimensional vector hydrophone battle array
RU2305297C2 (en) Device with correlation shaper of directivity characteristics for detection of signals and determination of direction to their source
CN102121982A (en) Orientation method based on sound level difference
US20220412835A1 (en) Ultrasonic detector
CN115656994A (en) Real-time calibration method for double-base active detection towed array formation
JP2010133907A (en) Angle measurement processing apparatus
CN111077515B (en) Target detection method based on analog television external radiation source radar
Kaouri Left-right ambiguity resolution of a towed array sonar
US11863280B2 (en) Method and apparatus for determining the directional frequency response of an arrangement of transducer elements

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