CN1110701C - Method for processing sonar wave beams and its data processing system - Google Patents

Method for processing sonar wave beams and its data processing system Download PDF

Info

Publication number
CN1110701C
CN1110701C CN 98112959 CN98112959A CN1110701C CN 1110701 C CN1110701 C CN 1110701C CN 98112959 CN98112959 CN 98112959 CN 98112959 A CN98112959 A CN 98112959A CN 1110701 C CN1110701 C CN 1110701C
Authority
CN
China
Prior art keywords
data
array
centerdot
sigma
opt
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
Application number
CN 98112959
Other languages
Chinese (zh)
Other versions
CN1246615A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN 98112959 priority Critical patent/CN1110701C/en
Publication of CN1246615A publication Critical patent/CN1246615A/en
Application granted granted Critical
Publication of CN1110701C publication Critical patent/CN1110701C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

The present invention relates to a method for processing sonar wave beams and a data processing system thereof. A groove noise field method is adopted, and an array manifold alpha(that<i>) which is estimated from actually measured array data is introduced to determine an optimal weight coefficient vector W<OPT>. The array data processing system which realizes the wave beam optimization method is composed of a computer and a data acquisition card, wherein the data acquisition card is a high-speed array data acquisition card which is composed of a data buffer isolator, a handshake circuit, a CPU, a program memory, a data memory and a decoding circuit. The present invention solves the problem of optimization design of determining sensor array directional patterns in optional structural shapes by academic modeling, and has the advantages of reasonable system design and convenient use.

Description

A kind of sonar wave beams disposal route and data handling system thereof
Technical field: the present invention relates to a kind of sonar wave beams disposal route and data handling system thereof, be mainly used in the data acquisition of arbitrary shape sensor array under water and the optimal design of beam pattern.Sonar wave beams figure optimization method adopts " groove noise field method " and adaptive array treatment technology to finish, for the data handling system of supporting sonar wave beams disposal route proposed by the invention to design is made up of data acquisition and data processor.
Background technology: sonar wave beams has clear and definite requirement to main lobe shape, beam position error, beamwidth, secondary lobe level etc.For solving the optimal design of sonar wave beams figure, desirable modeling optimization method was proposed once in theory.As " groove noise field method ", be a kind of pure theory modeling beam optimization method.This method has proposed the notion of " groove noise field method " and " design of arbitrary shape sensor array directional diagram optimum optimization ", and its main points have three: all primitives of (1), composition basic matrix, its roomage response is a uniformity.This only is theoretic desirable imagination, and may there be the roomage response of uniformity in the actual conditions on the engineering; (2), the main lobe direction of actual beam, and make that its noise intensity is zero, all the other directions are even noise intensity, this is a groove noise field method.This also is desirable modeling; (3) signal response that will retrain the wave beam main lobe direction is made as 1, also is with the Utopian processing of engineering problem.Then, the adaptive array treatment theory is applied to the sensor basic matrix of arbitrary shape.But, the result who relies on idealized modeling to obtain, when it was applied to actual sonar system, measured beam pattern quality was not good.This is because theoretical and in esse distance.In practical engineering application, because complicated formation structure can cause that sound covers, scattering and transmission effect, and the reason of manufacturing process can cause the shape and size error of basic matrix, makes that array element response (array manifold) is inconsistent; Everything makes the beam pattern can't be met actual needs, the phenomenon of measured beam pattern performance substantial deviation expected results occurred.So in engineering is used, wish to find a kind of new weight coefficient optimization method, always so that all can obtain practical beam pattern to the sensor array of the arbitrary shape of any reality.
For fear of the deficiencies in the prior art part, the present invention proposes a kind of sonar wave beams disposal route and data handling system thereof.This method is based upon on the actual measurement array data basis fully based on the adaptive beam-forming algorithm of array data restructuring array stream shape.The array data disposal system is a kind of integrated digital display circuit of function that is independent of signal processor, not only can realize array data collection and beam pattern optimal design, and can do beam pattern and measure.
Summary of the invention: thought of the present invention is to solve optimal design and the measurement of the random geometry sensor array beam figure in the actual engineering.Consideration is for the sensor array of a given arbitrary structures shape, adopts " groove noise field method ", keeping under main lobe direction and beamwidth is certain and output noise power the is minimum condition, and the application self-adapting wave beam forms theory and finds the solution-organize the best weights coefficient.The setting array signal phasor is C, and noise vector is X, and W is the best weights coefficient, tries to achieve beamformer output signal Y=W T(X+C); It is characterized in that: the best weights coefficient amount of owing W OPTAsk for not
Figure C9811295900031
It is resulting to be that theory is found the solution again W opt = R XX - 1 C * C T R XX - 1 C * , But introduce the array manifold α (θ that estimates from the actual measurement array data i) back and the optimum weight vector W that tries to achieve OPTThe optimum weight vector W that the present invention proposes OPTTry to achieve by following method: at first groove noise field vector is introduced actual measurement array manifold α (θ i) revised, and establish by J respectively from θ 1, θ 2, θ 3... θ JThe zero-mean narrow band noise σ of direction iN i(t), then X ( t ) = &Sigma; i = 1 J &sigma; i &CenterDot; N i ( t ) &CenterDot; &alpha; ( &theta; i ) , Then according to adaptive modeling beam optimization method desired output signal d ( t ) = &Sigma; i = 1 J &sigma; i &CenterDot; N i ( t ) &CenterDot; D ( &theta; i ) ; Promptly get the covariance matrix of array data R = E [ X ( t ) &CenterDot; X H ( t ) ] = &Sigma; i = 1 J &sigma; i 2 &CenterDot; &alpha; ( &theta; i ) &CenterDot; &alpha; H ( &theta; i ) , The covariance matrix of data and wanted signal P = E [ d * ( t ) &CenterDot; X ( t ) ] = &Sigma; i = 1 J &sigma; i 2 &CenterDot; D * ( &theta; i ) &CenterDot; &alpha; ( &theta; i ) , Try to achieve optimum weight vector W according to the covariance matrix of above-mentioned array data and the covariance matrix of data and wanted signal Opt=R -1P.Because the effect of having implied objective factor in the array manifold, so optimum weight vector W OPTPlayed the caused deviation of the revised theory value of finding the solution.This improved beam optimization method can obtain beam pattern desirable, that satisfy actual needs.
The present invention proposes a kind of array data disposal system that realizes above-mentioned beam optimization method, is made up of computing machine 2 and data collecting card 1, can finish array data collection, array beams figure optimal design and array beams figure and measure.It is characterized in that: data collecting card 1 is a kind of high-speed array data collecting card, is made up of data buffering isolator 3, handshake circuit 4, CPU5, program storage 6 and data-carrier store 7 and decoding scheme 8.Handshake circuit 4 receives the handshake of sonars, and CPU5 control data buffer isolator 3 reads in data from sonar, after CPU5 searches synchronous head, just by the sensor sequence number array data is deposited in the data-carrier store 7.When a collection of data acquisition finishes, again the content in the data-carrier store 7 is deposited in the hard disc of computing machine 2.Program storage 6 is placed the executive routine of CPU5, the chip selection signal of decoding scheme 8 generating routine storeies 6 and data-carrier store 7.
Description of drawings:
Fig. 1: array data disposal system theory diagram
Fig. 2: data collecting card theory diagram
Fig. 3: embodiment block diagram 1-data collecting card 2-computing machine 3-data buffering isolator 4-handshake circuit 5-CPU 6-data-carrier store 7-program storage 8-decoding scheme 9-volume array 10-Signal Pretreatment 11-concentric cable 12-demodulating equipment
Embodiment: now in conjunction with the embodiments the present invention is further described:
Present embodiment is to utilize a kind of array data disposal system of the present invention, measures with beam pattern at the beam optimization design that a volume array of being made up of 64 array elements is finished.The computing machine 2 of array data disposal system is selected commercial PC for use, the data buffering isolator 3 of data collecting card 1 adopts ternary output bus to transmit receiver, handshake circuit 4 adopts monostalbe trigger, CPU5 adopts the high-speed dsp chip, program storage 6 adopts high-speed RAM, data-carrier store 7 adopts multiport memory, and decoding scheme 8 adopts the GAL circuit.
Front end in system, signal source are the volume arrays 9 of one 64 array element composition, and array signal is through sending into demodulating equipment 12 with the Signal Pretreatment 10 and the concentric cable 11 of volume array 9 connections.The signal of demodulation is sent into the array data disposal system.The array data disposal system is estimated array manifold α (θ according to the array data of the reality of gathering i), calculate R and P, try to achieve optimum weight vector W OPT, basis again Y = W opt T &CenterDot; X + C Calculate beamformer output signal, thereby survey and draw out the panorama multi-beam figure of volume array.
Beneficial effect: the present invention's advantage compared to existing technology is:
Estimate the beam optimization method of array manifold based on array data, beam pattern this difficult problem that seriously degenerates when having solved optimization power that the pure theory modeling finds the solution and being applied to actual sonar system.The inventive method can be applied to the optimal design of arbitrary structures shape sensor array pattern, and this is that additive method can't be finished.Be the array data disposal system that realizes that this method designs, can finish the optimal design and the engineering beam pattern of complex array beam pattern and measure.System design is reasonable, and is easy to use.

Claims (3)

1, a kind of sonar wave beams disposal route, optimal design and measurement for random geometry sensor array beam figure, consideration is for the sensor array of a given arbitrary structures shape, adopt " groove noise field method ", keeping under main lobe direction and beamwidth is certain and output noise power the is minimum condition, the application self-adapting wave beam forms theory and finds the solution one group of best weights coefficient, the setting array signal phasor is C, noise vector is X, and W is the best weights coefficient, tries to achieve beamformer output signal Y=W T(X+C); It is characterized in that: optimum weight vector W OPTAsk for that no longer to be that theory is found the solution resulting W opt = R XX - 1 C * C T R XX - 1 C * , But introduce the array manifold α (θ that estimates from the actual measurement array data i) back and the optimum weight vector W that tries to achieve OPTOptimum weight vector W OPTTry to achieve by following method: at first groove noise field vector is introduced actual measurement array manifold α (θ i) revised, and establish by J respectively from θ 1, θ 2, θ 3... θ JThe zero-mean narrow band noise σ of direction iN i(t), then X ( t ) = &Sigma; i = 1 J &sigma; i &CenterDot; N i ( t ) &CenterDot; &alpha; ( &theta; i ) , Then according to adaptive modeling beam optimization method desired output signal d ( t ) = &Sigma; i = 1 J &sigma; i &CenterDot; N i ( t ) &CenterDot; D ( &theta; i ) ; Promptly get the covariance matrix of array data R = E [ X ( t ) &CenterDot; X H ( t ) ] = &Sigma; i = 1 J &sigma; i 2 &CenterDot; &alpha; ( &theta; i ) &CenterDot; &alpha; H ( &theta; i ) ; The covariance matrix of data and wanted signal P = E [ d * ( t ) &CenterDot; X ( t ) ] = &Sigma; i = 1 J &sigma; i 2 &CenterDot; D * ( &theta; i ) &CenterDot; &alpha; ( &theta; i ) ; Try to achieve optimum weight vector according to the covariance matrix of above-mentioned array data and the covariance matrix of data and wanted signal
W OPT=R -1·P。
2, a kind of data handling system that realizes the sonar wave beams disposal route of claim 1, form by computing machine (2) and data collecting card (1), it is characterized in that: data collecting card (1) is a kind of high-speed array data collecting card, it consists of the handshake circuit (4) of accepting the sonar synchronizing signal and connects CPU (5), and the data buffering isolator (3) of preposition receiving array data, program storage (6), data-carrier store (7) and decoding scheme (8) connect on the bus of CPU (5).
3, a kind of data handling system as claimed in claim 2, it is characterized in that: the data buffering isolator (3) of data collecting card (1) adopts ternary output bus to transmit receiver, handshake circuit (4) adopts monostalbe trigger, CPU (5) adopts the high-speed dsp chip, program storage (6) adopts high-speed RAM, data-carrier store (7) adopts multiport memory, and decoding scheme (8) adopts the GAL circuit.
CN 98112959 1998-09-02 1998-09-02 Method for processing sonar wave beams and its data processing system Expired - Fee Related CN1110701C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 98112959 CN1110701C (en) 1998-09-02 1998-09-02 Method for processing sonar wave beams and its data processing system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 98112959 CN1110701C (en) 1998-09-02 1998-09-02 Method for processing sonar wave beams and its data processing system

Publications (2)

Publication Number Publication Date
CN1246615A CN1246615A (en) 2000-03-08
CN1110701C true CN1110701C (en) 2003-06-04

Family

ID=5222746

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 98112959 Expired - Fee Related CN1110701C (en) 1998-09-02 1998-09-02 Method for processing sonar wave beams and its data processing system

Country Status (1)

Country Link
CN (1) CN1110701C (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100334464C (en) * 2003-06-19 2007-08-29 中国科学院声学研究所 Method for sonar array signal simulation by utilizing interpolation filter
CN101256234B (en) * 2008-04-15 2010-12-01 哈尔滨工程大学 Method for modeling lens sound field
CN104765037A (en) * 2015-04-22 2015-07-08 国家深海基地管理中心 Underwater target positioning robust method based on short perpendicular array
CN104765039B (en) * 2015-04-22 2017-06-13 国家深海基地管理中心 The method that passive sonar operating distance is improved using shallow sound field Spatial characteristic
RU2723145C1 (en) * 2019-11-18 2020-06-09 Российская Федерация, от имени которой выступает ФОНД ПЕРСПЕКТИВНЫХ ИССЛЕДОВАНИЙ Method and device for detecting noisy objects in the sea with onboard antenna

Also Published As

Publication number Publication date
CN1246615A (en) 2000-03-08

Similar Documents

Publication Publication Date Title
US11733381B2 (en) Sound velocity profile inversion method based on inverted multi-beam echo sounder
CN103292792B (en) Actual measurement SVP reconstruction method suitable for submarine detection and pseudo-landform processing
CN108540248B (en) Dynamic multi-path channel model for offshore wireless communication and construction method
CN101573624A (en) Received field intensity estimating device, and received field intensity estimating program
CN103308904B (en) A kind of sparse optimization method taking into account the three-dimensional camera shooting sonar system transducer array of near field and far-field performance
CN101504458B (en) Phase filtering based beam forming method
CN111679245B (en) Sub-array coincidence split beam orientation method based on uniform circular array
CN1110701C (en) Method for processing sonar wave beams and its data processing system
CN108519602A (en) Cold seepage fast imaging system and its array signal processing method
CN112130112B (en) Information source number estimation method based on acoustic vector array joint information processing
CN105445718A (en) DOA estimation method for distributed multi-carrier-warship BVR (beyond visual range) radar based on array reconfiguration
CN107729627B (en) Quick selection method for unit of opportunistic array antenna
CN112098938B (en) Six-element cone vector array-based underwater acoustic target dimension reduction matching sound field positioning method
CN113534161B (en) Beam mirror image focusing method for remotely positioning underwater sound source
CN110390073B (en) Multi-channel space synthesis azimuth filtering method for vector sensing
CN111381212A (en) Virtual ultrashort baseline positioning method based on subarray division
CN101398482B (en) Noise field numerical computation method in passiveness wideband detection of sound reception array
CN107402371A (en) A kind of smart antenna array
CN113126029B (en) Multi-sensor pulse sound source positioning method suitable for deep sea reliable acoustic path environment
CN109799477A (en) A kind of sequential vehicle fingerprint localization method and device towards millimeter wave car networking
CN106821418B (en) A kind of method and device reducing beam synthesizer power consumption
CN110244288A (en) A kind of sonar array signal processing method based on focusing playback principle
CN103792528B (en) A kind of method that underwater sound array Bartlett wave beam based on diagonal angle off-load forms
CN116299435B (en) Method and system for checking and converting echo data of coal-rock interface recognition radar
CN114791992B (en) Deep sea target forward scattering sound field calculation method based on ray theory

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C06 Publication
PB01 Publication
C14 Grant of patent or utility model
GR01 Patent grant
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee