CN109506763B - Vector hydrophone sparse array arrangement method based on teaching and learning optimization - Google Patents

Vector hydrophone sparse array arrangement method based on teaching and learning optimization Download PDF

Info

Publication number
CN109506763B
CN109506763B CN201811451174.4A CN201811451174A CN109506763B CN 109506763 B CN109506763 B CN 109506763B CN 201811451174 A CN201811451174 A CN 201811451174A CN 109506763 B CN109506763 B CN 109506763B
Authority
CN
China
Prior art keywords
individual
teacher
class
teaching
vector hydrophone
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
CN201811451174.4A
Other languages
Chinese (zh)
Other versions
CN109506763A (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.)
National Defense Technology Innovation Institute PLA Academy of Military Science
Original Assignee
National Defense Technology Innovation Institute PLA Academy of Military Science
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 National Defense Technology Innovation Institute PLA Academy of Military Science filed Critical National Defense Technology Innovation Institute PLA Academy of Military Science
Priority to CN201811451174.4A priority Critical patent/CN109506763B/en
Publication of CN109506763A publication Critical patent/CN109506763A/en
Application granted granted Critical
Publication of CN109506763B publication Critical patent/CN109506763B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H3/00Measuring characteristics of vibrations by using a detector in a fluid

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

The invention provides a vector hydrophone sparse arraying method based on teaching and learning optimization, which comprises the following steps of: (1) the method comprises the steps of determining a vector hydrophone combined directivity and composite array directivity function, (2) initializing an optimization process, (3) selecting individual teachers, (4) independently learning teachers, (5) teaching teachers, (6) mutually learning students, and (7) terminating judgment. The invention combines the practical classroom teaching process, introduces the teacher self-learning stage, utilizes the working mechanism of the vector hydrophones in the sparse array to switch and optimize in the two states of activation and closing, enhances the local searching capability of teaching and learning optimization, can obtain better sparse array optimization results, does not need to adjust the specific parameter factors of the algorithm, saves a large amount of algorithm adjusting time, and better accords with the practical engineering application.

Description

Vector hydrophone sparse array arrangement method based on teaching and learning optimization
Technical Field
The invention belongs to the technical field of sonar, and particularly relates to a sparse array arrangement method of vector hydrophones based on teaching and learning optimization.
Background
The vector hydrophone is used as a novel underwater sound measuring sensor, is structurally formed by compounding a traditional nondirectional sound pressure sensor and a dipole directional vibration velocity sensor, can synchronously and jointly measure sound pressure and vibration velocity information, fundamentally solves the problem of port and starboard fuzzy, and is widely applied to the fields of underwater sound warning sonar, towed line array sonar, shipboard array conformal sonar, multi-base sonar and the like.
The sparse array is an array with array elements distributed in a sparse form obtained by removing part of array elements from a traditional densely-arranged full array on the basis of meeting array performance constraint, and can reduce the number of sensors and reduce the system hardware cost on the premise of obtaining higher resolution. The vector hydrophone sparse array distribution technology is beneficial to solving the problems of invalid array element repair and distributed multi-base sonar array distribution in the actual sonar array, and has important engineering value.
R.V Rao et al propose a teaching and learning optimization algorithm based on the process of sharing and acquiring knowledge of teachers and students in the course of classroom teaching. In the teaching stage, class members learn according to the difference between the teacher matrix and the average level of students; in the learning stage, the members in the class are compared with each other, and the member with low fitness learns from the member with high fitness. The knowledge sharing mode omits the autonomous learning stage of the teacher, and generally speaking, the teacher is the member with the highest fitness in the whole population, and needs to spend a large amount of time for autonomous learning to prepare teaching contents before giving lessons in a classroom so as to improve the classroom teaching quality.
The existing sparse array arrangement technology generally adopts a random optimization algorithm to search for a global optimal solution of the problems, such as a simulated annealing algorithm, a genetic algorithm, a particle swarm algorithm and the like, and the random optimization algorithm has a common problem that the adjustment of specific parameter factors of the algorithm is very complicated, such as a cooling coefficient used for controlling the temperature drop speed in the simulated annealing algorithm, a variation factor and a cross factor in the genetic algorithm, and a learning factor and an inertia weight in the particle swarm algorithm. Such parameter factors require a lot of time to adjust in the algorithm implementation process, and once the problem model or application occasion changes, the algorithm parameters need to be adjusted again to adapt to a new problem model or application occasion, which is very unfavorable for practical engineering application.
Disclosure of Invention
Aiming at the defects or improvement requirements in the prior art, the invention provides a vector hydrophone sparse array arrangement method based on teaching and learning optimization, which mainly comprises the following steps:
(1) determining a vector hydrophone combination directivity and composite array directivity function:
the sound pressure and vibration speed channel output signals of a single vector hydrophone are utilized to obtain and combine to form a directivity function R (f, theta), wherein f is the center frequency of the hydrophone, and theta is the horizontal azimuth angle. The sound pressure signal is recorded as p, the vibration velocity signal is recorded as v, and the vibration velocity sensors can form dipole directivity in a three-dimensional space. In the ocean waveguide, since a vertical direction is a standing wave, two-dimensional directivity in a horizontal direction is generally considered.
Determining a composite array directivity function F consisting of a plurality of vector hydrophones according to the combined directivity function R (F, theta) of the single vector hydrophonem(f, θ) expressed by
Fm(f,θ)=R(f,θ)·F(f,θ)
Where F (F, θ) is the directivity function of the scalar array of sound pressures.
(2) Initializing an optimization process:
control parameters in the teaching and learning optimization algorithm are initialized, and the control parameters comprise the number Pz of the class individuals and the maximum iteration number Gen. And initializing and randomly assigning excitation weight values w of the vector hydrophone array, wherein the random excitation weight values are random real numbers between 0 and 1, the matrix scale is PzxM, and M represents the number of vector hydrophones in the array. The excitation weight value is indiscriminately acted on all sound pressure and vibration velocity channels of the vector hydrophone, namely the sound pressure and the excitation weight value on the vibration velocity channel in a single vector hydrophone are the same.
Defining an ith class individual matrix wiCorresponding objective function f (w)i) Is expressed as
Figure GSB0000188253080000031
||·||0Is represented by 10Norm ofOperation for measuring the sparsity, k, of the matrix1The penalty coefficient is used for controlling the peak level of the side lobe, and the value size of the penalty coefficient is related to the number of the beams in the whole visible area and the resolution ratio of the beams; s denotes the region where the beam energy value is higher than the desired sidelobe peak.
(3) Selecting individual teacher step
Evaluating the individual fitness of the class, selecting the individual with the best fitness as the individual of the teacher, wherein the fitness function of the class can be expressed as
Figure GSB0000188253080000032
Where g represents the current corresponding number of iterations. Selecting the member with the highest fitness in the class individuals as the teacher
Figure GSB0000188253080000033
The remaining class members are students.
(4) Teacher self-learning step
The teacher prepares classroom teaching contents through autonomous learning, the autonomous learning is developed by adopting a local search method, and the autonomous learning stage of the teacher comprises the following steps:
(4.1) selecting the mth vector hydrophone according to the random traversal sequence, and judging the excitation weight w of the sound pressure and vibration velocity channels of the vector hydrophonest,mIf the value is 0, starting an 'activation' process if the value is 0, and otherwise, starting a 'shutdown' process.
The "activation" process includes the steps of:
(4.1.1) assigning a value to the vector hydrophone by using a random excitation weight, wherein the random excitation weight is a real number between 0 and 1, and the excitation weight acts on all sound pressure and vibration velocity channels of the vector hydrophone.
(4.1.2) judging whether the new teacher matrix can obtain better fitness after being assigned, if so, receiving a new excitation weight, and otherwise, resetting all the excitation weights of the sound pressure and vibration velocity channels of the vector hydrophone to be 0.
The "shutdown" procedure includes the steps of:
and (4.2.1) setting the excitation weight of all sound pressure and vibration velocity channels of the vector hydrophone to be 0.
(4.2.2) judging whether the new teacher matrix can obtain better fitness after being assigned, if so, receiving a new excitation weight, and otherwise, returning all the excitation weights of the sound pressure and vibration velocity channels of the vector hydrophone to the original state.
(4.2) judging whether all vector hydrophones in the teacher individual traverse, if so, ending the self-learning step of the teacher, and otherwise, turning to the step (4.1).
(5) Teaching stage steps of teachers
In the stage, the teacher gives lessons according to the overall knowledge level of the class members, and the aim is to enable the class average level to be closer to the knowledge level of the teacher. The difference between the teacher and the average knowledge level of the class is defined as follows:
Figure GSB0000188253080000041
wherein rand represents [0, 1 ]]Random number of (a), representing the average knowledge level of the class, TfAnd expressing a teaching factor, wherein the factor is used for expressing a teaching effect and takes the value of 1 or 2 with equal probability. The class member knowledge updating mode at the stage is
Figure GSB0000188253080000042
If the updated value of the member individual
Figure GSB0000188253080000043
Better results than originally achieved are retained, otherwise the original stimulus state is returned.
(6) The step of mutual learning stage for students
(6.1) randomly selecting two student individual weight matrixes
Figure GSB0000188253080000044
And
Figure GSB0000188253080000045
the corresponding objective function values are respectively
Figure GSB0000188253080000046
And
Figure GSB0000188253080000047
if it is
Figure GSB0000188253080000048
I.e. the fitness of the individual i is higher than that of the individual j, at the moment, the individual j learns from the individual i,
Figure GSB0000188253080000049
otherwise the individual i learns from the individual j,
Figure GSB00001882530800000410
(6.2) if the updated weight matrix
Figure GSB0000188253080000051
If the corresponding objective function value is lower than the original value, then order
Figure GSB0000188253080000052
Otherwise, not processing; if the updated weight matrix
Figure GSB0000188253080000053
If the corresponding objective function value is lower than the original value, then order
Figure GSB0000188253080000054
Otherwise, not processing;
(6.3) repeating the steps (6.1) - (6.2) to make the multiple class individuals have an opportunity to learn each other, wherein the number of times of repeated execution is determined by the number of the class individuals.
(7) Termination judgment step
If the current iteration times g reach the maximum iteration times Gen, terminating the optimization process, outputting a final class excitation weight matrix, and selecting the optimal individual as a final vector hydrophone sparse array scheme; otherwise, go to step (3).
Compared with the prior art, the invention has the following advantages:
(1) the vector sparse array arrangement method provided by the invention does not need to adjust specific parameter factors of the algorithm, and can execute the sparse array arrangement optimization process by only setting the population scale and the maximum iteration number of the algorithm, thereby saving a large amount of algorithm adjustment time and being more in line with the practical application of engineering.
(2) The invention combines the practical classroom teaching process, introduces the teacher self-learning stage, utilizes the working mechanism of the vector hydrophones in the sparse array to switch and optimize in the states of activation and closing, enhances the local searching capability of teaching and learning optimization, and can obtain better sparse array optimization results.
Drawings
FIG. 1 is a flow chart of the sparse arraying method of vector hydrophones according to the present invention;
FIG. 2 is a diagram showing a vector hydrophone array excitation weight distribution after sparse arrangement;
fig. 3 is a diagram of a sparse vector hydrophone array versus the original array beam.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Fig. 1 is a flow chart of a sparse arraying method for vector hydrophones based on teaching and learning optimization, based on which, the sparse arraying is performed on the most common linear arrays in the sonar field, and the main steps are as follows:
(1) determining a vector hydrophone combination directivity and composite array directivity function:
obtaining a directional function R (f, theta) formed by combining sound pressure and vibration speed channel output signals of a single vector hydrophone, wherein f is the center frequency of the vector hydrophone, theta is the horizontal azimuth angle, the sound pressure signal is recorded as p, the vibration speed signal is recorded as v, and the combined directivity is
Figure GSB0000188253080000061
In the formula
Figure GSB0000188253080000062
To guide orientation, vcFor combining the vibration velocities, from two orthogonal vibration velocity components vxAnd vyIs formed by combining the following formula
Figure GSB0000188253080000063
In the example, the most common half-wavelength uniform linear array in the sonar field is considered, the number of array elements is 32, the central frequency of the vector hydrophone array is 8kHz, and the focusing direction of the main beam is 90 degrees.
Determining a composite array directivity function F consisting of a plurality of vector hydrophones according to the combined directivity function R (F, theta) of the single vector hydrophonem(f, θ) expressed by
Fm(f,θ)=R(f,θ)·F(f,θ)
Wherein F (F, theta) is a directivity function of a sound pressure scalar array and is expressed as
Figure GSB0000188253080000064
(2) Initializing an optimization process:
control parameters in the teaching and learning optimization algorithm are initialized, and the control parameters comprise the number Pz of the class individuals and the maximum iteration number Gen. The number of individual classes in this example is set to 50 and the maximum number of iterations is 100.
Initializing and randomly assigning an excitation weight w of the vector hydrophone array, wherein the random excitation weight is a random real number between 0 and 1, the matrix scale is 50 multiplied by 32, and the excitation weight is indiscriminately acted on all sound pressure and vibration speed channels of the vector hydrophone, namely the sound pressure in a single vector hydrophone and the excitation weight on the vibration speed channel are the same.
Defining an ith class individual matrix wiCorresponding objective function f (w)i) Is expressed as
Figure GSB0000188253080000071
k1The penalty coefficient is used for controlling the peak level of the side lobe, and the value size of the penalty coefficient is related to the number of the beams in the whole visible area and the resolution ratio of the beams; s denotes the region where the beam energy value is higher than the desired sidelobe peak. The penalty factor is set to k in this example10.5, the side lobe region is set to [0, 86 °]∪[94°,360°]。
(3) Selecting individual teacher step
Evaluating the individual fitness of the class, selecting the individual with the best fitness as the individual of the teacher, wherein the fitness function of the class can be expressed as
Figure GSB0000188253080000072
Where g represents the current corresponding number of iterations. Selecting the member with the highest fitness in the class individuals as the teacher
Figure GSB0000188253080000073
The remaining class members are students.
(4) Teacher self-learning step
The teacher prepares classroom teaching contents through autonomous learning, the autonomous learning is developed by adopting a local search method, and the autonomous learning stage of the teacher comprises the following steps:
(4.1) selecting the mth vector hydrophone according to the random traversal sequence, and judging the excitation weight w of the sound pressure and vibration velocity channels of the vector hydrophonest,mIf the value is 0, starting an 'activation' process if the value is 0, and otherwise, starting a 'shutdown' process.
The "activation" process includes the steps of:
(4.1.1) assigning a value to the vector hydrophone by using a random excitation weight, wherein the random excitation weight is a real number between 0 and 1, and the excitation weight acts on all sound pressure and vibration velocity channels of the vector hydrophone.
(4.1.2) judging whether the new teacher matrix can obtain better fitness after being assigned, if so, receiving a new excitation weight, and otherwise, resetting the excitation weights of all the sound pressure and vibration velocity channels of the vector hydrophone to be 0.
The "shutdown" procedure includes the steps of:
and (4.2.1) setting the excitation weight of all sound pressure and vibration velocity channels of the vector hydrophone to be 0.
(4.2.2) judging whether the new teacher matrix can obtain better fitness after being assigned, if so, receiving a new excitation weight, and otherwise, returning all the excitation weights of the sound pressure and vibration velocity channels of the vector hydrophone to the original state.
(4.2) judging whether all vector hydrophones in the teacher individual traverse, if so, ending the self-learning step of the teacher, and otherwise, turning to the step (4.1).
(5) Teaching stage steps of teachers
In the stage, the teacher gives lessons according to the overall knowledge level of the class members, and the aim is to enable the class average level to be closer to the knowledge level of the teacher. The difference between the teacher and the average knowledge level of the class is defined as follows:
Figure GSB0000188253080000081
wherein rand represents [0, 1 ]]Random number of (d), wgRepresenting the average knowledge level, T, of a classfAnd expressing a teaching factor, wherein the factor is used for expressing a teaching effect and takes the value of 1 or 2 with equal probability. The class member knowledge updating mode at the stage is
Figure GSB0000188253080000082
If the updated value of the member individual
Figure GSB0000188253080000083
Compared with the original method, the method obtains better resultsWill be retained otherwise the original stimulus state is returned.
(6) The step of mutual learning stage for students
(6.1) randomly selecting two student individual weight matrixes
Figure GSB0000188253080000091
And
Figure GSB0000188253080000092
the corresponding objective function values are respectively
Figure GSB0000188253080000093
And
Figure GSB0000188253080000094
if it is
Figure GSB0000188253080000095
I.e. the fitness of the individual i is higher than that of the individual j, at the moment, the individual j learns from the individual i,
Figure GSB0000188253080000096
otherwise the individual i learns from the individual j,
Figure GSB0000188253080000097
(6.2) if the updated weight matrix
Figure GSB0000188253080000098
If the corresponding objective function value is lower than the original value, then order
Figure GSB0000188253080000099
Otherwise, not processing; if the updated weight matrix
Figure GSB00001882530800000910
If the corresponding objective function value is lower than the original value, then order
Figure GSB00001882530800000911
Otherwise, not processing;
(6.3) repeating the steps (6.1) - (6.2) to make the multiple class individuals have an opportunity to learn each other, wherein the number of times of repeated execution is determined by the number of the class individuals.
(7) Termination judgment step
If the current iteration times g reach the maximum iteration times 100, terminating the optimization process, outputting a final class excitation weight matrix, and selecting the optimal individual as a final vector hydrophone sparse array scheme; otherwise, go to step (3).
The number of hydrophones in the sparse vector array finally obtained in the present example is 24, the excitation weight distribution is shown in fig. 2, and fig. 3 is a spatial beam contrast diagram of the sparse vector array and the original array. According to the final sparse array arrangement result, the sparse vector array has the advantages that the number of hydrophones is reduced by 25% compared with the original array, the maximum side lobe peak value is reduced to-18 dB from the original-13.3 dB, and the better side lobe suppression effect is achieved.
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and that any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (2)

1. A vector hydrophone sparse arraying method based on teaching and learning optimization is characterized by comprising the following steps:
(1) determining a vector hydrophone combination directivity and composite array directivity function:
utilizing the sound pressure and vibration speed channel output signals of a single vector hydrophone to obtain and combine to form a directivity function R (f, theta), wherein f is the hydrophone center frequency, and theta is a horizontal azimuth angle; the sound pressure signal is recorded as p, the vibration velocity signal is recorded as v, and the vibration velocity sensor can form dipole directivity in a three-dimensional space; in the ocean waveguide, since the vertical direction is a standing wave, the two-dimensional directivity in the horizontal direction is generally considered;
determining a composite array directivity function composed of a plurality of vector hydrophones according to the combined directivity function R (f, theta) of the single vector hydrophoneNumber Fm(f, θ) expressed by
Fm(f,θ)=R(f,θ)·F(f,θ)
Wherein F (F, θ) is a directivity function of the sound pressure scalar array;
(2) initializing an optimization process:
initializing control parameters in a teaching and learning optimization algorithm, wherein the control parameters comprise the number Pz of class individuals and the maximum iteration number Gen; initializing random assignment to an excitation weight w of the vector hydrophone array, wherein the random excitation weight is a random real number between 0 and 1, and the matrix scale is PzxM, wherein M represents the number of vector hydrophones in the array; the excitation weight value is indiscriminately acted on all sound pressure and vibration velocity channels of the vector hydrophone, namely the sound pressure and the excitation weight value on the vibration velocity channel in a single vector hydrophone are the same;
defining an ith class individual matrix wiCorresponding objective function f (w)i) Is expressed as
Figure FSB0000188253070000011
||·||0Is represented by0Norm operation for measuring the sparsity, k, of matrix1The penalty coefficient is used for controlling the peak level of the side lobe, and the value size of the penalty coefficient is related to the number of the beams in the whole visible area and the resolution ratio of the beams; s represents the region where the beam energy value is higher than the desired sidelobe peak;
(3) selecting individual teacher step
Evaluating the individual fitness of the class, selecting the individual with the best fitness as the individual of the teacher, wherein the fitness function of the class can be expressed as
Figure FSB0000188253070000021
Wherein g represents the current corresponding iteration number; selecting the member with the highest fitness in the class individuals as the teacher
Figure FSB0000188253070000022
The rest members of the class are students;
(4) teacher self-learning step
The teacher prepares classroom teaching contents through autonomous learning, the autonomous learning is developed by adopting a local search method, and the autonomous learning stage of the teacher comprises the following steps:
(4.1) selecting the mth vector hydrophone according to the random traversal sequence, and judging the excitation weight w of the sound pressure and vibration velocity channels of the vector hydrophonest,mIf the value is 0, starting an 'activation' process if the value is 0, otherwise, starting a 'closing' process;
(4.2) judging whether all vector hydrophones in the teacher individual traverse, if so, ending the self-learning step of the teacher, otherwise, turning to the step (4.1);
(5) teaching stage steps of teachers
In the stage, a teacher gives lessons according to the overall knowledge level of class members, and the aim is to enable the class average level to be closer to the knowledge level of the teacher; the difference between the teacher and the average knowledge level of the class is defined as follows:
Figure FSB0000188253070000023
wherein rand represents [0, 1 ]]Random number of (d), wgRepresenting the average knowledge level, T, of a classfExpressing a teaching factor, wherein the factor is used for expressing a teaching effect and takes the value of 1 or 2 with equal probability; the class member knowledge updating mode at the stage is
Figure FSB0000188253070000031
If the updated value of the member individual
Figure FSB0000188253070000032
Compared with the original method, the method obtains better results, and is retained, otherwise, the method returns to the original excitation state;
(6) the step of mutual learning stage for students
(6.1) random selection of two schoolsWeight matrix of the living entity
Figure FSB0000188253070000033
And
Figure FSB0000188253070000034
the corresponding objective function values are respectively
Figure FSB0000188253070000035
And
Figure FSB0000188253070000036
if it is
Figure FSB0000188253070000037
I.e. the fitness of the individual i is higher than that of the individual j, at the moment, the individual j learns from the individual i,
Figure FSB0000188253070000038
otherwise the individual i learns from the individual j,
Figure FSB0000188253070000039
(6.2) if the updated weight matrix
Figure FSB00001882530700000310
If the corresponding objective function value is lower than the original value, then order
Figure FSB00001882530700000311
Otherwise, not processing; if the updated weight matrix
Figure FSB00001882530700000312
If the corresponding objective function value is lower than the original value, then order
Figure FSB00001882530700000313
Otherwise, not processing;
(6.3) repeatedly executing the steps (6.1) - (6.2) to enable a plurality of class individuals to have an opportunity to learn each other, wherein the repeated execution times are determined by the number of the class individuals;
(7) termination judgment step
If the current iteration times g reach the maximum iteration times Gen, terminating the optimization process, outputting a final class excitation weight matrix, and selecting the optimal individual as a final vector hydrophone sparse array scheme; otherwise, go to step (3).
2. The sparse vector hydrophone arraying method based on teaching and learning optimization as claimed in claim 1, wherein the activating process comprises the following steps:
(4.1.1) assigning a value to the vector hydrophone by using a random excitation weight, wherein the random excitation weight is a real number between 0 and 1, and the excitation weight acts on all sound pressure and vibration velocity channels of the vector hydrophone;
(4.1.2) judging whether the new teacher matrix can obtain better fitness after assignment, if so, receiving a new excitation weight, and otherwise, resetting the excitation weights of all the sound pressure and vibration velocity channels of the vector hydrophone to be 0;
the "shutdown" procedure includes the steps of:
(4.2.1) setting excitation weights of all sound pressure and vibration velocity channels of the vector hydrophone to be 0;
(4.2.2) judging whether the new teacher matrix can obtain better fitness after being assigned, if so, receiving a new excitation weight, and otherwise, returning all the excitation weights of the sound pressure and vibration velocity channels of the vector hydrophone to the original state.
CN201811451174.4A 2018-11-30 2018-11-30 Vector hydrophone sparse array arrangement method based on teaching and learning optimization Active CN109506763B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811451174.4A CN109506763B (en) 2018-11-30 2018-11-30 Vector hydrophone sparse array arrangement method based on teaching and learning optimization

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811451174.4A CN109506763B (en) 2018-11-30 2018-11-30 Vector hydrophone sparse array arrangement method based on teaching and learning optimization

Publications (2)

Publication Number Publication Date
CN109506763A CN109506763A (en) 2019-03-22
CN109506763B true CN109506763B (en) 2020-09-25

Family

ID=65751437

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811451174.4A Active CN109506763B (en) 2018-11-30 2018-11-30 Vector hydrophone sparse array arrangement method based on teaching and learning optimization

Country Status (1)

Country Link
CN (1) CN109506763B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110298109B (en) * 2019-06-28 2023-07-14 中国人民解放军军事科学院国防科技创新研究院 Hydrophone array beam restoration method under failure array element condition
CN113671507B (en) * 2021-07-14 2024-01-05 中国人民解放军军事科学院国防科技创新研究院 Waveguide invariant estimation method based on deep sea vertical array

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101420068A (en) * 2008-11-25 2009-04-29 电子科技大学 Distribution method for sensor antenna array
CN101625408A (en) * 2009-08-03 2010-01-13 浙江大学 Method for rarefying and optimizing energy transducer array of three-dimensional imaging sonar system
CN103675799A (en) * 2013-10-24 2014-03-26 华中科技大学 Sparse planar array optimizing method for energy transducers of phased array sonar system
CN103744072A (en) * 2013-12-24 2014-04-23 浙江大学 Sparse array optimization method based on simulated annealing algorithm and distributed parallel subarray beam forming algorithm
CN105354351A (en) * 2015-09-25 2016-02-24 浙江大学 Improved simulated annealing algorithm based limit sparse array optimization method
CN105842702A (en) * 2016-05-13 2016-08-10 南京航空航天大学 Multi-beam imaging sonar sidelobe suppression method and array sparsity method
CN106099393A (en) * 2016-08-25 2016-11-09 电子科技大学 Submatrix bare cloth based on genetic algorithm and array element sparse antenna array are comprehensively structured the formation method
CN106596724A (en) * 2016-12-08 2017-04-26 清华大学 Method for sparse compression optimization reconstruction of narrowband Lamb wave detection data

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040158154A1 (en) * 2003-02-06 2004-08-12 Siemens Medical Solutions Usa, Inc. Portable three dimensional diagnostic ultrasound imaging methods and systems
CN107145943A (en) * 2017-04-25 2017-09-08 南京信息工程大学 Method based on the detection echo state network small-signal for improving teaching optimized algorithm in Chaotic Background

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101420068A (en) * 2008-11-25 2009-04-29 电子科技大学 Distribution method for sensor antenna array
CN101625408A (en) * 2009-08-03 2010-01-13 浙江大学 Method for rarefying and optimizing energy transducer array of three-dimensional imaging sonar system
CN103675799A (en) * 2013-10-24 2014-03-26 华中科技大学 Sparse planar array optimizing method for energy transducers of phased array sonar system
CN103744072A (en) * 2013-12-24 2014-04-23 浙江大学 Sparse array optimization method based on simulated annealing algorithm and distributed parallel subarray beam forming algorithm
CN105354351A (en) * 2015-09-25 2016-02-24 浙江大学 Improved simulated annealing algorithm based limit sparse array optimization method
CN105842702A (en) * 2016-05-13 2016-08-10 南京航空航天大学 Multi-beam imaging sonar sidelobe suppression method and array sparsity method
CN106099393A (en) * 2016-08-25 2016-11-09 电子科技大学 Submatrix bare cloth based on genetic algorithm and array element sparse antenna array are comprehensively structured the formation method
CN106596724A (en) * 2016-12-08 2017-04-26 清华大学 Method for sparse compression optimization reconstruction of narrowband Lamb wave detection data

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"教与学"优化算法研究综述;拓守恒等;《计算机应用研究》;20130730;第30卷(第7期);第1933-1938页 *
Synthesis of Thinned Square Array Using Modified Teaching-Learning-Based Optimization Algorithm;Luo Zailei et al.;《Proceedings of ISAP 2014》;20141205;第145-146页 *
基于教与学优化算法的PID控制器参数寻优;何学明等;《计算机工程》;20150830;第41卷(第8期);第313-316页 *
相控阵三维摄像声纳系统的稀疏阵列优化设计;袁龙涛等;《华南理工大学学报(自然科学版)》;20130130;第41卷(第l期);第29-37页 *

Also Published As

Publication number Publication date
CN109506763A (en) 2019-03-22

Similar Documents

Publication Publication Date Title
CN103675799B (en) A kind of phased array sonar system transducer sparse battle array optimization method
CN105572658B (en) The a burst of first sparse optimization method of three-dimensional imaging sonar receiving plane based on improved adaptive GA-IAGA
CN108828603B (en) Cross-based sparse optimization method for three-dimensional imaging sonar array
CN111983619B (en) Underwater acoustic target forward scattering acoustic disturbance positioning method based on transfer learning
CN110068795A (en) A kind of indoor microphone array sound localization method based on convolutional neural networks
CN108318862A (en) A kind of sound localization method based on neural network
CN109557526B (en) Vector hydrophone sparse array arrangement method based on compressed sensing theory
CN109506763B (en) Vector hydrophone sparse array arrangement method based on teaching and learning optimization
US5537511A (en) Neural network based data fusion system for source localization
CN109489796A (en) A kind of underwater complex structural radiation noise source fixation and recognition based on unit radiation method and acoustic radiation forecasting procedure
CN109635486A (en) A kind of high resolution three-dimensional imaging sonar transducer array sparse optimization method
CN113607447A (en) Acoustic-optical combined fan fault positioning device and method
CN109709510A (en) A kind of estimation method and system of coherent 2-d direction finding
CN108614235B (en) Single-snapshot direction finding method for information interaction of multiple pigeon groups
US5488589A (en) Neural network based three dimensional ocean modeler
CN110736976A (en) sonar beam former performance estimation method of arbitrary array
CN111352075B (en) Underwater multi-sound-source positioning method and system based on deep learning
CN110864802A (en) Self-noise forecasting method for ship-shell sonar platform area based on virtual sound source wave superposition
CN112511199B (en) Low-complexity 3D-MIMO related channel fast generation method
Frank et al. Comparing vision-based to sonar-based 3D reconstruction
CN111458698B (en) Passive sonar sparse bit optimization method
Wang et al. Maximum likelihood DOA estimation based on improved invasive weed optimization algorithm and application of MEMS vector hydrophone array
CN111965601A (en) Underwater sound source passive positioning method based on nuclear extreme learning machine
CN116466334B (en) Two-dimensional sonar array sparsification method based on longhorn beetle whisker search algorithm
CN114429154B (en) Underwater sound array orientation method for multipath assisted deep learning

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