Sonar apparatus, target velocity display method, and non-transitory medium

Abstract

A processor of a sonar system extracts a signal in a predetermined time range from a signal received by an acoustic array, calculates a correlation value for a signal waveform extracted, detects a peak value of the correlation value, calculates an integrated value by performing time-integration of a square of the correlation value, obtains a target velocity based on the peak value and the integrated value, and displays the target velocity on a display apparatus.

Claims

1. A sonar apparatus comprising: a processor in circuit communication to an acoustic array and a display apparatus; and a memory storing program instructions executable by the processor, wherein the processor, when executing the program instructions stored in the memory: extracts one or more signal waveforms in a predetermined time range from a signal received by the acoustic array; calculates a correlation value of the one or more signal waveforms that has been extracted; detects a peak value of the correlation value; calculates an integrated value by performing time-integration of a square of the correlation value; derives a target velocity, based on the peak value and the integrated value; and displays the target velocity on a display apparatus.

2. The sonar apparatus according to claim 1, wherein the one or more signal waveforms comprise two signal waveforms, the signal received by the acoustic array is a pulse signal modulated by Linear Frequency Modulation (LFM), and wherein the processor, when executing the program instructions stored in the memory: calculates a cross-correlation value between the two signal waveforms extracted, as the correlation value; and calculates a Doppler coefficient, based on the integrated value and the peak value of the cross-correlation value, and an absolute value of a frequency change rate and a pulse length to derive the target velocity.

3. The sonar apparatus according to claim 1, wherein the one or more signal waveforms comprise two signal waveforms, the correlation value is a cross-correlation value between the two signal waveforms, the signal received by the acoustic array is a pulse signal modulated by Linear Frequency Modulation (LFM), and wherein the processor, when executing the program instructions stored in the memory: calculates a Doppler coefficient, based on a peak value R.sub.(0) of an autocorrelation R.sub.() of the cross-correlation value, a peak value R.sub.aa(0) of an autocorrelation R.sub.aa() of the autocorrelation R.sub.(), and an absolute value of a frequency change rate and a pulse length to derive the target velocity.

4. The sonar apparatus according to claim 1, wherein the processor, when executing the program instructions stored in the memory: calculates a cross-correlation value R.sub.SS(), an autocorrelation R.sub.() of the cross-correlation value, and an autocorrelation R.sub.aa() of the autocorrelation R.sub.() of the cross-correlation value, and performs frequency analysis of at least one of the cross-correlation value R.sub.SS(), the autocorrelation R.sub.() of the cross-correlation value, and the autocorrelation R.sub.aa() of the autocorrelation R.sub.() of the cross-correlation value to derive the target velocity.

5. The sonar apparatus according to claim 1, wherein the signal received by the acoustic array is a Pulse Continuous Wave (PCW) signal, and wherein the processor, when executing the program instructions stored in the memory: calculates a Doppler coefficient based on the peak value of the correlation value, the integrated value of an autocorrelation of the correlation value, and a pulse length, to derive the target velocity.

6. The sonar apparatus according to claim 1, wherein a relationship among at least the peak value, the integrated value, and a relative velocity between the sonar apparatus and a target having the target velocity are stored in advance in a storage apparatus, and wherein the processor, when executing the program instructions stored in the memory: derives the target velocity, based on at least the peak value and the integrated value obtained for the signal received by the acoustic array, with reference to the storage apparatus.

7. The sonar apparatus according to claim 1, wherein the signal received by the acoustic array is a single pulse, and wherein the processor, when executing the program instructions stored in the memory: calculates an autocorrelation value of the single pulse, as the correlation value; calculates the integrated value by performing time-integration of a square of the autocorrelation value; and derives the target velocity, based on the peak value of the autocorrelation value of the single pulse and an integrated value of the square of the autocorrelation value.

8. The sonar apparatus according to claim 1, wherein the one or more signal waveforms comprise three or more signal waveforms, and wherein the processor, when executing the program instructions stored in the memory: selects one or more pairs from among the three or more signal waveforms extracted; calculates a cross-correlation value between each of the one or more pairs; obtains one or more target velocities, for the one or more pairs, each of the one or more target velocities being derived based on the peak value and an integrated value of the square of the cross-correlation value between each of the one or more pairs to calculate, as the target velocity, an average value of the one or more target velocities; and displays the average value of the target velocity for each of the one or more pairs on the display.

9. The sonar apparatus according to claim 1, wherein the processor, when executing the program instructions stored in the memory: outputs a transmission signal at a predetermined transmission interval; and calculates the correlation value between two signal waveforms extracted, a time interval therebetween corresponding to the transmission interval for a plurality of the signals extracted.

10. The sonar apparatus according to claim 1, wherein the predetermined time range is longer than a pulse length of the received signal.

11. A computer-implemented target velocity display method, comprising: extracting, by a processor, one or more signal waveforms in a predetermined time range from a signal received by the acoustic array; calculating, by the processor, a correlation value of the one or more signal waveforms that has been extracted; detecting, by the processor, a peak value of the correlation value; calculating, by the processor, an integrated value by performing time-integration of a square of the correlation value; and deriving, by the processor, a target velocity based on the peak value and the integrated value, to display the target velocity on a display apparatus.

12. The computer-implemented target velocity display method according to claim 11, wherein the one or more signal waveforms comprise two signal waveforms, the signal received by the acoustic array is a pulse signal modulated by Linear Frequency Modulation (LFM), and wherein the method further comprises: calculating, by the processor, a cross-correlation value between the two signal waveforms extracted, as the correlation value; and calculating, by the processor, a Doppler coefficient, based on the integrated value and the peak value of the cross-correlation value, and an absolute value of a frequency change rate and a pulse length to derive the target velocity.

13. The computer-implemented target velocity display method according to claim 11, wherein the one or more signal waveforms comprise two signal waveforms, the correlation value is a cross-correlation value between the two signal waveforms, the signal received by the acoustic array is a pulse signal modulated by Linear Frequency Modulation (LFM), and wherein the method further comprises: calculating, by the processor, a Doppler coefficient, based on a peak value R.sub.(0) of an autocorrelation R.sub.() of the cross-correlation value between two signal waveforms extracted, a peak value R.sub.aa(0) of an autocorrelation R.sub.aa() of the autocorrelation R.sub.(), and an absolute value of a frequency change rate and a pulse length to derive the target velocity.

14. The computer-implemented target velocity display method according to claim 11, further comprising: calculating, by the processor, a cross-correlation value R.sub.SS(), an autocorrelation of the cross-correlation value, and an autocorrelation R.sub.aa() of the autocorrelation R.sub.() of the cross-correlation value; and performing, by the processor, frequency analysis of at least one of the cross-correlation value R.sub.SS(), the autocorrelation R.sub.() of the cross-correlation value, and the autocorrelation R.sub.aa() of the autocorrelation R.sub.() of the cross-correlation value to derive the target velocity.

15. The computer-implemented target velocity display method according to claim 11, wherein the signal received by the acoustic array is a Pulse Continuous Wave (PCW) signal, and wherein the method further comprises: calculating, by the processor, a Doppler coefficient based on the peak value of the correlation value, the integrated value of an autocorrelation of the correlation value, and a pulse length, to derive the target velocity.

16. The computer-implemented target velocity display method according to claim 11, further comprising: storing, by the processor, a relationship among at least the peak value, the integrated value, and a relative velocity between the sonar apparatus and a target having the target velocity in advance in a storage apparatus; and deriving, by the processor, the target velocity, based on at least the peak value and the integrated value obtained for the signal received by the acoustic array, with reference to the storage apparatus.

17. The computer-implemented target velocity display method according to claim 11, wherein the signal received by the acoustic array is a single pulse, and wherein the method further comprises: calculating, by the processor, an autocorrelation value of the single pulse, as the correlation value; calculating, by the processor, the integrated value by performing time-integration of a square of the autocorrelation value; and deriving, by the processor, the target velocity, based on the peak value of the autocorrelation value of the single pulse and an integrated value of the square of the autocorrelation value.

18. A non-transitory computer-readable medium storing thereon a program that when executed by a processor causes the processor to execute processing comprising: extracting one or more signal waveforms in a predetermined time range from a signal received by the acoustic array; calculating a correlation value of the one or more signal waveforms that has been extracted; detecting a peak value of the correlation value; calculating an integrated value by performing time-integration of a square of the correlation value; and deriving a target velocity based on the peak value and the integrated value, to display the target velocity on a display apparatus.

19. The non-transitory computer-readable medium according to claim 18, wherein the one or more signal waveforms comprise two signal waveforms, the signal received by the acoustic array is a pulse signal modulated by Linear Frequency Modulation (LFM), and wherein the program when executed by the processor causes the processor to execute the processing further comprising: calculating a cross-correlation value between the two signal waveforms extracted; and calculating a Doppler coefficient, based on the integrated value and the peak value of the cross-correlation value, and an absolute value of a frequency change rate and a pulse length to derive the target velocity.

20. The non-transitory computer-readable medium according to claim 18, wherein the one or more signal waveforms comprise two signal waveforms, the correlation value is a cross-correlation value between the two signal waveforms, the signal received by the acoustic array is a pulse signal modulated by Linear Frequency Modulation (LFM), and wherein the program when executed by the processor causes the processor to execute the processing further comprising: calculating a Doppler coefficient, based on a peak value R.sub.(0) of an autocorrelation R.sub.() of the cross-correlation value, a peak value R.sub.aa(0) of an autocorrelation R.sub.aa() of the autocorrelation R.sub.(), and an absolute value of a frequency change rate and a pulse length to derive the target velocity.

Description

BRIEF DESCRIPTION OF DRAWINGS

(1) FIG. 1 is a diagram illustrating an example embodiment according to the present invention.

(2) FIG. 2 is a diagram illustrating an example of a configuration of an apparatus of an example embodiment according to the present invention.

(3) FIG. 3 is a diagram illustrating a processing procedure of an apparatus of an example embodiment according to the present invention.

(4) FIG. 4 is a diagram illustrating an example embodiment according to the present invention.

(5) FIG. 5 is a diagram illustrating an example embodiment according to the present invention.

(6) FIG. 6 is a diagram illustrating TDI.

(7) FIG. 7 is a diagram illustrating a tandem transmission waveform I.

(8) FIG. 8 is a diagram illustrating correlation by a plurality of continuous HFMs.

DETAILED DESCRIPTION

(9) An example embodiment of the present invention will be described. FIG. 1 illustrates an example embodiment according to the present invention, and schematically illustrates target echoes when an LFM signal is transmitted twice. Referring to FIG. 1, a sonar system (sonar apparatus) is configured to transmit, multiple times at intervals, a signal waveform, such as an LFM signal or the like, which can provide high distance resolution by taking correlation thereof and performs correlation processing among a plurality of target echoes received. When transmitting an LFM signal twice, as illustrated in FIG. 1, the active sonar apparatus extracts two signal waveforms, each in a predetermined time range (time window) from the received signal, and calculates a cross-correlation between the extracted signal waveforms. Furthermore, the sonar apparatus performs time-integration of the correlation processing result (a square of an absolute value of the cross-correlation value) and estimates a target velocity, based on a relationship between the integration result and a peak of the cross-correlation value. As illustrated in FIG. 1, the time window (predetermined time range) used to sequentially extract a received signal is preferably longer than an LFM waveform (a pulse length). Furthermore, the time interval between a plurality of signals to be extracted corresponds to a time interval of the transmission waveform (LFM).

(10) According to the example embodiment of the present invention, a computational load is not increased even in a low SNR environment and is not affected by change of a transmission path in time, and hence it is possible to provide a display which increases target distance accuracy and enhances target visibility.

(11) FIG. 2 illustrates an example of a configuration of a sonar system (sonar apparatus) 100 (target velocity display system) of an example embodiment according to the present invention. In the sonar system 100, an acoustic array 101 includes a plurality of acoustic elements which transmit an acoustic wave into a water by converting an electrical signal outputted from a transmission controller (transmission control part) 102 to an acoustic signal and converts a received acoustic signal to an electrical signal. The transmission controller 102 generates an electrical signal with a waveform specified. An AD converter (analog-to-digital converter) 103 converts an analog electrical signal supplied from each acoustic element to a digital signal. A phase adjuster (phaser) (phase adjusting part) 104 performs phase adjustment (phasing) for each azimuth direction on digital signals, each received from each acoustic element to perform beamforming and outputs a beam for each azimuth direction. A correlation processor (correlation processing part) 105 performs a correlation processing for each beam and calculates a target velocity.

(12) A display apparatus 106 displays a target velocity on a screen (not illustrated). The acoustic array 101, the AD converter 103, the phase adjuster 104, and the display apparatus 106 themselves are apparatuses generally available, and those employed in various existing sonar systems may be used.

(13) The following describes the transmission controller 102 and the correlation processor 105, focusing operation thereof. The transmission controller 102 generates an electrical signal that enables a pulse having a predetermined waveform to be transmitted multiple times. On this occasion, the transmission controller 102 shares information on an interval between pulses in addition to each pulse waveform with the correlation processor 105. Transmission times of transmission waveforms, for example, twice or more may be sufficient, and the number of times is not particularly limited.

(14) The correlation processor 105 operates according to a flow as shown in FIG. 3.

(15) The correlation processor 105 extracts a signal in a predetermined time range from a received signal 1051 for each beam outputted from the phase adjuster 104 and stores the signal in a storage part (not shown) (S101). This time range preferably exceeds one pulse length.

(16) Next, the correlation processor 105 selects, from a plurality of extracted signal waveforms, a combination of two extracted signal waveforms, with time interval therebetween corresponding to a time interval of the transmitted pulse and calculates a cross-correlation value between the two extracted signal waveforms (S102).

(17) Next, the correlation processor 105 detects a peak of the cross-correlation value (S103) and stores the peak value in the storage part 1052.

(18) The correlation processor 105 integrates a cross-correlation result (time integration result of a square of the cross-correlation) and stores an integrated value in the storage part 1053.

(19) Then, the correlation processor 105 calculates a target velocity by supplying, as inputs, the peak value and the integrated value respectively obtained in steps S103 and S104 to a function which represents relationship between a peak value and an integrated value stored in the storage part 1054 in advance (S105). Since the relationship between a peak value and an integrated value depends on Doppler effect, it is possible to derive a target velocity based on the relationship between a peak value and an integrated value.

(20) The correlation processor 105 supplies the target velocity to the display apparatus 106 (S106).

(21) An example where a transmission waveform is LFM will be described below.

(22) Let an acoustic velocity be c, let a velocity component of the sonar in a direction from the sonar to the target be v.sub.s under a condition where a direction from the sonar to a target is positive and let a velocity component of the target in a direction from the sonar to the target be v.sub.o (that is, v.sub.s is positive when the sonar has a velocity directed to the target and v.sub.o is positive when the target has a velocity component moving away from the sonar). Let a relative velocity between the sonar system 100 and the target (not shown) be v=v.sub.sv.sub.o, a Doppler coefficient (and its square) is given (approximated) by an expression (1a) (an expression (1b)).

(23) $\begin{array}{cc}1+\frac{2\left(v_s-v_o\right)}{c}=1+\frac{2v}{c}& \left(1a\right)\end{array}$$\begin{array}{cc}{}^{2}1+\frac{4v}{c}& \left(1b\right)\end{array}$

(24) That is, let a frequency of a sound source be f.sub.0, a frequency f.sub.1 of an acoustic wave from the sound source to a target is given as f.sub.1=f.sub.0c/(cv) and a frequency f.sub.2 of a reflected wave from the target is given as f.sub.2=f.sub.1(c+v)/c. Hence, we have the following expression.

(25) $\begin{array}{cc}=\frac{f_2}{f_0}=\frac{\left(c+v\right)}{\left(c-v\right)}=\frac{\left(1+\frac{v}{c}\right)}{\left(1-\frac{v}{c}\right)}{\left(1+\frac{v}{c}\right)}^{2}1+2\frac{v}{c}\left(vc\right)& \left(2\right)\end{array}$

(26) In a case where the pulse has a single frequency, the frequency f.sub.2 from the target becomes f.sub.0 by Doppler effect.

(27) Here, a transmission waveform (complex representation) (replica) is assumed to be LFM as represented in the following expression (3).

(28) $\begin{array}{cc}s\left(t\right)=\{\begin{array}{cc}B\exp \left(jt+\frac{j{t}^2}{2}\right)& .Math.t.Math.\frac{T_0}{2}\\ 0& \mathrm{other}\mathrm{than}\mathrm{avobe}\end{array}& \left(3\right)\end{array}$

(29) In the expression (3), B is an amplitude, (=2f.sub.0: f.sub.0 is a start frequency) is an angular frequency, t is time, and is a frequency change rate (chirp rate), where j.sup.2=1, and T.sub.0 is a pulse length.

(30) An instantaneous frequency f(t) of the transmission waveform is given by time-differentiating a phase component of the transmission waveform of the expression (3) =t+t.sup.2/2, as below.

(31) $\begin{array}{cc}f\left(t\right)=\left(\frac{1}{2}\right)\frac{d}{\mathrm{dt}}\left\{t+\left(\frac{1}{2}\right)t^2\right\}=f_0+\left(\frac{1}{2}\right)t& \left(4\right)\end{array}$

(32) The received waveform subjected to Doppler effect is given by the following expression (5).

(33) $\begin{array}{cc}S\left(t\right)=\{\begin{array}{cc}A\sqrt{}\exp \left\{j.Math.\left(t-t_0\right)+\frac{j.Math.{{}^2\left(t-t_0\right)}^2}{2}\right\}& .Math.\left(t-t_0\right).Math.\frac{T_0}{2}\\ 0& \mathrm{other}\mathrm{than}\mathrm{above}\end{array}& \left(5\right)\end{array}$

(34) It is assumed that, by Doppler effect, the instantaneous frequency of the received waveform (LFM) is given by f(t) for an instantaneous frequency f(t) (the expression (4)) of a transmission waveform as with a single pulse. It can also be regarded that a time is set from t to t. That is, a phase component of the received waveform may be obtained by time-integrating the instantaneous frequency f(t) of the above expression (4) for a range from 0 to t.

(35) $\begin{array}{cc}=2{}_o^t\left\{f_0+\left(\frac{1}{2}\right)t^{}\right\}{\mathrm{dt}}^{}=2\left\{f_{0}t+\left(\frac{1}{2}\right)\frac{{\left(t\right)}^2}{2}\right\}=t+\frac{{\left(t\right)}^2}{2}& \left(6\right)\end{array}$

(36) When t.sub.0=0 in the expression (5) only for the sake of calculation simplification of, the received waveform S(t) is represented as below.

(37) $\begin{array}{cc}S\left(t\right)=\{\begin{array}{cc}A\sqrt{}\exp \left\{j.Math.t+\frac{j.Math.{}^2{t}^2}{2}\right\}& .Math.t.Math.\frac{T_0}{2}\\ 0& \mathrm{other}\mathrm{than}\mathrm{above}\end{array}& \left(7\right)\end{array}$

(38) For example, a correlation (cross-correlation) between two received waveforms (target echoes) which are waveforms returned by reflection from a target when the same waveform is transmitted twice can be regarded as an autocorrelation of the received waveform.

(39) $$\begin{gathered} \begin{array}{cc}\mathrm{For}=0, \\ {R}_{\mathrm{SS}}\left(0\right)={}_{-}^{}{S}^{*}\left(t\right)S\left(t\right)\mathrm{dt}={.Math.A.Math.}^2{}_{-T_0/2}^{T_0/2}\exp \left(-jt-j{}^2{t}^2/2\right)\exp \left(jt+j{}^2{t}^2/2\right)\mathrm{dt}={.Math.A.Math.}^2{}_{-T_0/2}^{T_0/2}\mathrm{dt}={.Math.A.Math.}^2{T}_0& \left(8\right)\end{array} \end{gathered}$$
(where * denotes complex conjugate operator)

(40) Next, for 0<<T.sub.0,

(41) $\begin{array}{cc}{R}_{\mathrm{SS}}\left(\right)={}_{-}^{}{S}^{*}\left(t\right)S\left(t+\right)\mathrm{dt}={.Math.A.Math.}^2{}_{\frac{-T_0}{2}}^{\frac{T_0}{2}-}\exp \left\{-j.Math.t-\frac{j.Math.{}^2{t}^2}{2}\right\}\exp \left(jt+\frac{j{}^2{t}^2}{2}+j+\frac{j{}^2{}^2}{2}+j{}^{2}t\right)\mathrm{dt}={.Math.A.Math.}^2{}_{-T_0/2}^{T_0/2-}\exp \left(j+\frac{j{}^2{}^2}{2}+j{}^{2}t\right)\mathrm{dt}={.Math.A.Math.}^2\exp \left(j+\frac{j{}^2{}^2}{2}\right){}_{-T_0/2}^{T_0/2-}\exp \left(j{}^{2}t\right)\mathrm{dt}& \left(9\right)\end{array}$

(42) A definite integral of a function exp(j.sup.2t) of the above expression (9) with an integration section of [T.sub.0/(2),T.sub.0/(2)] is given by the following expression (10).

(43) 0$\begin{array}{cc}{}_{-T_0/2}^{T_0/2-}\exp \left(j{}^{2}t\right)\mathrm{dt}={\left[\frac{1}{j{}^2}\exp \left(j{}^{2}t\right)\right]}_{-T_0/2}^{T_0/2-}=\frac{1}{j{}^2}\left\{\exp \left(\frac{j{T}_0}{2}-j{}^2{}^2\right)-\exp \left(-\frac{j{T}_0}{2}\right)\right\}=\frac{\exp \left(-\frac{j{}^2{}^2}{2}\right)}{j{}^2}\left\{\exp \left(\frac{j{T}_0}{2}-\frac{j{}^2{}^2}{2}\right)-\exp \left(-\frac{j{T}_0}{2}+\frac{j{}^2{}^2}{2}\right)\right\}=\frac{2\exp \left(-\frac{j{}^2{}^2}{2}\right)}{{}^2}\sin \left(\frac{T_0}{2}-\frac{{}^2{}^2}{2}\right)=\frac{2\exp \left(-\frac{j{}^2{}^2}{2}\right)}{{}^2}\sin \left\{\frac{}{2}\left(T_0-\right)\right\}& \left(10\right)\end{array}$

(44) Accordingly, the expression (9) is given as follows.

(45) $\begin{array}{cc}{R}_{SS}\left(\right)={.Math.A.Math.}^2\exp \left(j+\frac{j{}^2{}^2}{2}\right)\frac{2\exp \left(-\frac{j{}^2{}^2}{2}\right)}{{}^2}\sin \left\{\frac{}{2}\left(T_0-\right)\right\}=\frac{2{.Math.A.Math.}^2}{}\exp \left(j\right)\sin \left\{\frac{}{2}\left(T_0-\right)\right\}=\frac{{.Math.A.Math.}^2\left(T_0-\right)\exp \left(j\right)\sin \left\{\frac{}{2}\left(T_0-\right)\right\}}{\frac{}{2}\left(T_0-\right)}={.Math.A.Math.}^2\left(T_0-\right)\exp \left(j\right)\mathrm{sinc}\left\{\frac{}{2}\left(T_0-\right)\right\}& \left(11\right)\end{array}$
where sin c(x) is a denormalized sin c function (sin x)/x.

(46) With such approximation as ||<<T.sub.0, the expression (11) is approximated as below.

(47) $\begin{array}{cc}{R}_{\mathrm{SS}}\left(\right){.Math.A.Math.}^2{T}_0\exp \left(j\right)\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(12\right)\end{array}$

(48) Form the above, the following is derived with respect to R.sub.SS().

(49) $\begin{array}{cc}{R}_{\mathrm{SS}}\left(\right)=\{\begin{array}{cc}\frac{2{.Math.A.Math.}^2}{}\exp \left(j\right)\mathrm{sinc}\left\{\frac{}{2}\left(T_0-\right)\right\}& 0<T_0\\ {.Math.A.Math.}^2{T}_0& =0\\ \frac{2{.Math.A.Math.}^2}{}\exp \left(j\right)\mathrm{sinc}\left\{\frac{}{2}\left(T_0+\right)\right\}& -T_0<0\\ 0& <-T_0\mathrm{or}>T_0\end{array}& \left(13\right)\end{array}$

(50) The above expression (13), when approximated by ||<<T.sub.0, is given as below.

(51) $\begin{array}{cc}{R}_{\mathrm{SS}}\left(\right)\{\begin{array}{cc}\frac{2{.Math.A.Math.}^2}{}\exp \left(j\right)\sin \left(\frac{T_0}{2}\right)=& 0<T_0\\ {.Math.A.Math.}^2{T}_0\exp \left(j\right)\mathrm{sinc}\left(\frac{T_0}{2}\right)& \\ {.Math.A.Math.}^2{T}_0& =0\\ \frac{2{.Math.A.Math.}^2}{}\exp \left(j\right)\sin \left(\frac{T_0}{2}\right)=& -T_0<0\\ {.Math.A.Math.}^2{T}_0\exp \left(j\right)\mathrm{sinc}\left(\frac{T_0}{2}\right)& \\ 0& <-T_0\mathrm{or}>T_0\end{array}& \left(14\right)\end{array}$

(52) Since R.sub.SS() is intensively concentrated at =0 and sharply drops to zero at .fwdarw., R.sub.SS () can be regarded to continue infinitely. That is, R.sub.SS() can be approximated by the following expression (15) for all .

(53) $\begin{array}{cc}{R}_{\mathrm{SS}}\left(\right){.Math.A.Math.}^2{T}_0\exp \left(j\right)\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(15\right)\end{array}$

(54) The correlation value R.sub.SS() in the above expression (15) takes its peak value at =0. Since exp(0)=1 and sinc(0)=1, the peak value P is given by the following equation (16).
P=R.sub.SS(0)=|A|.sup.2T.sub.0 (16)

(55) The following describes an autocorrelation of R.sub.SS(), wherein R.sub.SS(t) in the expression (15) is denoted by R.sub.() as in the following equation (17) and let its autocorrelation be denoted by R.sub.() (autocorrelation of R.sub.SS()).

(56) $\begin{array}{cc}{R}_{}\left(t\right)={.Math.A.Math.}^2{T}_0\exp \left(jt\right)\mathrm{sinc}\left(\frac{T_{0}t}{2}\right)& \left(17\right)\end{array}$$\begin{array}{cc}{R}_{}\left(\right)={}_{-}^{}{R}_{}^{*}\left(t\right)R_{}\left(t+\right)\mathrm{dt}={.Math.A.Math.}^4{T}_0^2{}_{-}^{}\exp \left(-jt\right)\mathrm{sinc}\left(\frac{T_{0}t}{2}\right)\exp \left(jt+j\right)\mathrm{sinc}\left(\frac{T_{0}t}{2}+\frac{T_0}{2}\right)\mathrm{dt}={.Math.A.Math.}^4{T}_0^2\exp \left(j\right){}_{-}^{}\mathrm{sinc}\left(\frac{T_{0}t}{2}\right)\mathrm{sinc}\left(\frac{T_{0}t}{2}+\frac{T_0}{2}\right)\mathrm{dt}={.Math.A.Math.}^4{T}_0^2\exp \left(j\right){}_{-}^{}\mathrm{sinc}\left(\frac{T_{0}t}{2}\right)\mathrm{sinc}(\frac{T_0}{2}\left(t+\right)\}\mathrm{dt}& \left(18\right)\end{array}$

(57) An outline of derivation of an integral term in the above expression (18) will be described below. A well-known mathematical formula

(58) $\begin{array}{cc}{}_{-}^{}\frac{\sin \left\{\left(x-a\right)\right\}}{x-a}\frac{\sin \left\{\left(x-b\right)\right\}}{x-b}dx=\frac{.Math.\sin \left\{\left(a-b\right)\right\}}{a-b}& \left(19\right)\end{array}$
is transformed as below.

(59) $\begin{array}{cc}{}_{-}^{}\frac{.Math.\sin \left\{\left(x-a\right)\right\}}{\left(x-a\right)}\frac{.Math.\sin \left\{\left(x-b\right)\right\}}{\left(x-b\right)}dx=\frac{.Math.\sin \left\{\left(a-b\right)\right\}}{\left(a-b\right)}& \left(20\right)\end{array}$

(60) Accordingly, we have the following equation.

(61) $\begin{array}{cc}{}_{-}^{}.Math.\mathrm{sinc}\left\{\left(x-a\right)\right\}.Math.\sin c\left\{\left(x-b\right)\right\}\mathrm{dx}=.Math.\mathrm{sinc}\left\{\left(a-b\right)\right\}& \left(21\right)\end{array}$$\begin{array}{cc}{}_{-}^{}\mathrm{sinc}\left\{\left(x-a\right)\right\}.Math.\sin c\left\{\left(x-b\right)\right\}\mathrm{dx}=\frac{}{}\mathrm{sinc}\left\{\left(a-b\right)\right\}& \left(22\right)\end{array}$

(62) In the above equation (22), when letting

(63) 0$==\frac{T_0}{2},x=t,a=0,b=-$

(64) we have the following equation.

(65) $\begin{array}{cc}{}_{-}^{}\sin c\left(\frac{T_0}{2}r\right).Math.\mathrm{sinc}\left\{\frac{T_0}{2}\left(t+\right)\right\}\mathrm{dt}=\frac{2}{T_0}\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(23\right)\end{array}$

(66) If >0,

(67) $\begin{array}{cc}{}_{-}^{}\mathrm{sinc}\left(\mathrm{ax}\right).Math.\sin c\left\{a\left(x+b\right)\right\}\mathrm{dx}=\frac{}{a}\mathrm{sinc}\left(\mathrm{ab}\right)& \left(24\right)\end{array}$

(68) If <0,

(69) $\begin{array}{cc}{}_{-}^{}\mathrm{sinc}\left(\mathrm{ax}\right).Math.\sin c\left\{a\left(x+b\right)\right\}\mathrm{dx}=-\frac{}{a}\mathrm{sinc}\left(\mathrm{ab}\right)& \left(25\right)\end{array}$

(70) If =0,
.sub..sup. sin c(ax).Math.sin c{(x+b)}dx=0 (26)

(71) Hence, if 0, the following holds.

(72) $\begin{array}{cc}{}_{-}^{}\mathrm{sinc}\left(\mathrm{ax}\right).Math.\mathrm{sinc}\left\{a\left(x+b\right)\right\}\mathrm{dx}=\frac{}{.Math.a.Math.}\mathrm{sinc}\left(\mathrm{ab}\right)& \left(27\right)\end{array}$

(73) Since , T.sub.0>0 always holds, the above equation (23) can be represented as below, if 0.

(74) $\begin{array}{cc}{}_{-}^{}\mathrm{sinc}\left(\frac{T_0}{2}t\right).Math.\mathrm{sinc}\left\{\frac{T_0}{2}\left(t+\right)\right\}dt=\frac{2}{.Math..Math.T_0}\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(28\right)\end{array}$

(75) Hence, if 0, the expression (18) can be given as the following expression (29).

(76) $\begin{array}{cc}{R}_{}\left(\right)={}_{-}^{}{R}_{}^{*}\left(t\right)R_{}\left(t+\right)\mathrm{dt}={.Math.A.Math.}^4{T}_0^2\exp \left(j\right){}_{-}^{}\mathrm{sinc}\left(\frac{T_{0}t}{2}\right)\mathrm{sinc}\left\{\frac{T_0}{2}\left(t+\right)\right\}\mathrm{dt}={.Math.A.Math.}^4{T}_0^2\exp \left(j\right)\frac{2}{.Math..Math.T_0}\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(29\right)\end{array}$

(77) R.sub.SS () is given by following expression (where above expression (13) is cited again).

(78) $\begin{array}{cc}{R}_{\mathrm{SS}}\left(\right)=\{\begin{array}{cc}\frac{2{.Math.A.Math.}^2}{}\exp \left(j\right)\sin \left\{\frac{}{2}\left(T_0-\right)\right\}& 0<T_0\\ {.Math.A.Math.}^2{T}_0& =0\\ \frac{2{.Math.A.Math.}^2}{}\exp \left(j\right)\sin \left\{\frac{}{2}\left(T_0+\right)\right\}& -T_0<0\\ 0& <-T_0\mathrm{or}>T_0\end{array}& \left(13\right)\end{array}$

(79) Hence, from

(80) $\begin{array}{cc}{R}_{\mathrm{SS}}\left(\right){.Math.A.Math.}^2{T}_0\exp \left(j\right)\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(30\right)\end{array}$

(81) we have the following equation.
R.sub.SS(0)=|A|.sup.2T.sub.0 (31)

(82) From the equation (29), we have

(83) $\begin{array}{cc}{R}_{}\left(0\right)={}_{-}^{}{R}_{}^{*}\left(t\right)R_{}\left(t\right)\mathrm{dt}={.Math.A.Math.}^4{T}_0^2\frac{2}{.Math..Math.T_0}& \left(32\right)\end{array}$

(84) From the equations (31) and (32), we have

(85) 0$\begin{array}{cc}{R}_{}\left(0\right)={\left\{R_{SS}\left(0\right)\right\}}^2\frac{2}{.Math..Math.T_0}& \left(33\right)\end{array}$

(86) Accordingly, from the above equation (33), the Doppler coefficient can be derived by following equation (34).

(87) $\begin{array}{cc}=\frac{2{\left\{R_{ss}\left(0\right)\right\}}^2}{.Math..Math.T_0{R}_{}\left(0\right)}& \left(34\right)\end{array}$

(88) In the equation (34), R.sub.SS(0) is a peak P of the cross-correlation value R.sub.SS(), and thus the equation (34) can be represented as follows.

(89) $\begin{array}{cc}=\frac{2P^2}{.Math..Math.T_0{R}_{}\left(0\right)}& \left(35\right)\end{array}$

(90) where R.sub.(0) is a time-integrated value of a square of an absolute value of a cross-correlation result R.sub.SS(). If the integrated value R.sub.(0) is obtained, it is not necessary to calculate an autocorrelation function form R.sub.(). The autocorrelation function form R.sub.() has been derived in order to obtain an autocorrelation of an autocorrelation described later.

(91) In the right side expression of the equation (35), since a frequency change rate is a value set in advance, it is possible to obtain a Doppler coefficient . Hence, using the above expression (1a), a relative velocity v can be obtained by the following expression (36).

(92) $\begin{array}{cc}v=\frac{c\left(-1\right)}{2}=\frac{c}{2}\left(\frac{2P^2}{.Math..Math.T_0{R}_{}\left(0\right)}-1\right)& \left(36\right)\end{array}$

(93) The display apparatus 106 receives the relative velocity v calculated using the above expression (36), from the correlation processor 105 and displays the relative velocity v, for example, as illustrated in FIG. 4. The display apparatus 106 may receive the peak value P of the correlation value, the integrated value S of the square of the correlation value, the frequency change rate , and the pulse length T.sub.0 from the correlation processor 105, perform calculation of the expression (36) to obtain the relative velocity v, and displays the relative velocity v on the screen in association with a target distance.

(94) By setting a range of integration of the square of a correlation value to, for example, a pulse interval (signal interval) or more of an LFM signal, a process gain which is a ratio of SNRs before and after the correlation can be the same or more than that obtainable when a replica correlation is performed. It is possible to raise visibility compared with a display based on PCW (Pulse Continuous Wave) with no correlation processing.

(95) For example, when substituting a pulse length T=1 s and an analysis width f=1 s in the following expression (37), a process gain for PCW is PG=0 dB

(96) $\begin{array}{cc}PG=5\log \left(\frac{T}{f}\right)\left[\mathrm{dB}\right]& \left(37\right)\end{array}$

(97) In contrast, a process gain for the replica correlation of LFM becomes significantly large to be PG=23 dB, where a band range width B=100 Hz and a pulse length T=1 s are substituted in the following expression.
PG=10 log(2BT) [dB](38)

(98) The expressions (37) and (38) are cited from expressions (11.84) and (11.65) in NPL 1.

(99) According to the present example embodiment, although a plurality of pulses are used, it is made possible to estimate a target velocity without being affected by variation of pulse intervals of a propagation path of an acoustic wave.

(100) As illustrated in FIG. 4, it is possible to grasp a position of a target with a high accuracy because a target position is indicated by a peak of correlation. Furthermore, because a transmission waveform of FH (Frequency Hopping) or the like can be used, a distance resolution can be raised higher than that in case of using a transmission waveform of HFM, visibility of a target can be increased.

(101) The following considers calculating an autocorrelation of an autocorrelation result R.sub.() of the expression (29).

(102) The autocorrelation result is given below.

(103) $\begin{array}{cc}{R}_{}\left(t\right)={.Math.A.Math.}^4{T}_0^2\exp \left(jt\right)\frac{2}{.Math..Math.T_0}\mathrm{sinc}\left(\frac{T_0}{2}t\right)={.Math.A.Math.}^4{T}_0\frac{2}{.Math..Math.}\exp \left(jt\right)\mathrm{sinc}\left(\frac{T_0}{2}t\right)& \left(39\right)\end{array}$

(104) The autocorrelation of the autocorrelation result is given below.

(105) $\begin{array}{cc}{R}_{\mathrm{aa}}\left(\right)={}_{-}^{}{R}_{}^{*}\left(t\right)R_{}\left(t+\right)\mathrm{dt}={.Math.A.Math.}^8{T}_0^2\frac{4{}^2}{{.Math..Math.}^2{}^2}{}_{-}^{}\exp \left(-jt\right)\mathrm{sinc}\left(\frac{T_{0}t}{2}\right)\exp \left(jt+jt\right)\mathrm{sinc}\left(\frac{T_{0}t}{2}+\frac{T_{0}t}{2}\right)\mathrm{dt}={.Math.A.Math.}^8{T}_0^2\frac{4{}^2}{{.Math..Math.}^2{}^2}\exp \left(jt\right){}_{-}^{}\mathrm{sinc}\left(\frac{T_{0}t}{2}\right)\mathrm{sinc}\left(\frac{T_{0}t}{2}+\frac{T_{0}t}{2}\right)\mathrm{dt}& \left(40\right)\end{array}$

(106) If 0, the following equation holds.

(107) $\begin{array}{cc}{}_{-}^{}\mathrm{sinc}\left(\frac{T_0}{2}t\right).Math.\mathrm{sinc}\left\{\frac{T_0}{2}\left(t+\right)\right\}\mathrm{dt}=\frac{2}{.Math..Math.T_0}\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(41\right)\end{array}$

(108) Therefore, the expression (40) is given as below.

(109) $\begin{array}{cc}{R}_{aa}\left(\right)={.Math.A.Math.}^8{T}_0^2\frac{4{}^2}{{.Math..Math.}^2{}^2}\exp \left(j\right)\frac{2}{.Math..Math.T_0}\mathrm{sinc}\left(\frac{T_0}{2}\right)={.Math.A.Math.}^8{T}_0\frac{8{}^3}{{.Math..Math.}^3{}^3}\exp \left(j\right)\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(42\right)\end{array}$

(110) From the above expression (39), we have

(111) $\begin{array}{cc}{R}_{}\left(0\right)={.Math.A.Math.}^4{T}_0\frac{2}{.Math..Math.}& \left(43\right)\end{array}$

(112) Furthermore, from the above expression (42), we have the following equation.

(113) 0$\begin{array}{cc}{R}_{aa}\left(0\right)={}_{-}^{}{R}_{}^{*}\left(t\right)R_{}\left(t\right)\mathrm{dt}={.Math.A.Math.}^8{T}_0\frac{8{}^3}{{.Math..Math.}^3{}^3}& \left(44\right)\end{array}$

(114) Thus, the following equation holds.

(115) $\begin{array}{cc}\frac{R_{aa}\left(0\right)}{{\left\{R_{}\left(0\right)\right\}}^2}=\frac{{.Math.A.Math.}^8{T}_0\frac{8{}^3}{{.Math..Math.}^3{}^3}}{{.Math.A.Math.}^8{T}_0^2\frac{4{}^2}{{.Math..Math.}^2{}^2}}=\frac{2}{T_0.Math..Math.}& \left(45\right)\end{array}$

(116) From the equation (45), the Doppler coefficient is derived by the following equation (46).

(117) $\begin{array}{cc}=\frac{2{\left\{R_{}\left(0\right)\right\}}^2}{T_0.Math..Math.R_{aa}\left(0\right)}& \left(46\right)\end{array}$

(118) In this manner, it is also possible to obtain the Doppler coefficient from the peak value R.sub.(0) of the autocorrelation R.sub.() of the cross-correlation value R.sub.SS() and the peak value R.sub.aa(0) of the autocorrelation of the autocorrelation R.sub.(0). The relative velocity v can be derived based on the above equation (36) from the obtained Doppler coefficient .

(119) Here, among cross-correlation in the expression (15) which reads:

(120) $\begin{array}{cc}{R}_{SS}\left(\right){.Math.A.Math.}^2{T}_0\exp \left(j\right)\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(47\right)\end{array}$

(121) In the autocorrelation of the cross-correlation in the above expression (29) which reads:

(122) $\begin{array}{cc}{R}_{}\left(\right)={.Math.A.Math.}^4{}_0^2\exp \left(j\right)\frac{2}{.Math..Math.T_0}\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(48\right)\end{array}$

(123) and the autocorrelation of the autocorrelation in the above expression (42) which reads:

(124) $\begin{array}{cc}{R}_{aa}\left(\right)={.Math.A.Math.}^8{T}_0\frac{8{}^3}{{.Math..Math.}^3{}^3}\exp \left(j\right)\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(49\right)\end{array}$$\begin{array}{cc}\exp \left(j\right)\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(50\right)\end{array}$
is a common term, wherein it can be regarded that a signal (exp(j)) whose angular frequency is modulated by

(125) $\begin{array}{cc}\mathrm{sinc}\left(\frac{T_0}{2}\right)& \left(51\right)\end{array}$

(126) That is, a frequency analysis is performed for correlation values of the expressions (47) to (49) (at least one thereof). For example, regarding a waveform which is modulated by,

(127) $\mathrm{sinc}\left(\frac{T_0}{2}\right)$
it is possible to obtain the Doppler coefficient , by obtaining an average frequency of the modulation.

(128) A time domain waveform sin c(Wx) (where W=T.sub.0/2), when Fourier transformed, becomes a rectangular wave in a frequency domain.

(129) $\begin{array}{cc}{}_{-}^{}\frac{\sin \left(\mathrm{Wt}\right)}{t}{e}^{-it}\mathrm{dt}=\{\begin{array}{cc}& \mathrm{if}.Math..Math.W\\ 0& \mathrm{if}.Math..Math.>W\end{array}& \left(52\right)\end{array}$

(130) Fourier transform of
exp(j)
becomes function in a frequency domain, and a rectangular wave is obtained as a result of a convolution operation of the function and the rectangular wave of the expression (52) in the frequency domain. From information of a width W of the rectangular wave, the Doppler coefficient may be obtained.

(131) An example of PCW will be described below, where it is assumed that a transmission waveform (replica) is given by following expression.

(132) $\begin{array}{cc}s\left(t\right)=\{\begin{array}{cc}B\exp \left(jt\right)& .Math.t.Math.\frac{T_0}{2}\\ 0& \mathrm{Other}\mathrm{than}\mathrm{above}\end{array}& \left(53\right)\end{array}$

(133) A received waveform is given below, with t.sub.0=0 for the sake of simplicity,

(134) 0$\begin{array}{cc}S\left(t\right)=\{\begin{array}{cc}A\sqrt{}\exp \left\{j.Math.t\right\}& .Math.t.Math.\frac{T_0}{2}\\ 0& \mathrm{other}\mathrm{than}\mathrm{above}\end{array}& \left(54\right)\end{array}$

(135) A transmission waveform is transmitted twice and autocorrelation of two received signals reflected by a target is calculated.

(136) For =0, the autocorrelation is given as below.

(137) $\begin{array}{cc}{R}_{\mathrm{SS}}\left(0\right)={}_{-}^{}{S}^{*}\left(t\right)S\left(t\right)\mathrm{dt}={.Math.A.Math.}^2{}_{-T_0/2}^{T_0/2}\exp \left\{-jt\right\}\exp \left(jt\right)\mathrm{dt}={.Math.A.Math.}^2{}_{-T_0/2}^{T_0/2}={.Math.A.Math.}^2{T}_0& \left(55\right)\end{array}$

(138) For 0<T.sub.0, the autocorrelation is given as below.

(139) $\begin{array}{cc}{R}_{\mathrm{SS}}\left(\right)={}_{-}^{}{S}^{*}\left(t\right)S\left(t+\right)\mathrm{dt}={.Math.A.Math.}^2{}_{-T_0/2}^{T_0/2-}\exp \left\{-j.Math.t\right\}\exp \left(jt+j\right)\mathrm{dt}={.Math.A.Math.}^2{}_{-T_0/2}^{T_0/2-}\exp \left(j\right)\mathrm{dt}={.Math.A.Math.}^2\exp \left(j\right){}_{-T_0/2}^{T_0/2-}\mathrm{dt}={.Math.A.Math.}^2\exp \left(j\right)\left(\frac{T_0}{}-\right)={.Math.A.Math.}^2\exp \left(j\right)\left(T_0-\right)& \left(56\right)\end{array}$

(140) On the other hand, for T.sub.0<0, the autocorrelation is given as below.

(141) $\begin{array}{cc}{R}_{\mathrm{SS}}\left(\right)={}_{-}^{}{S}^{*}\left(t\right)S\left(t+\right)\mathrm{dt}={.Math.A.Math.}^2{}_{-T_0/2-}^{T_0/2}\exp \left\{-j.Math.t\right\}\exp \left(jt+j\right)\mathrm{dt}={.Math.A.Math.}^2{}_{-T_0/2-}^{T_0/2}\exp \left(j\right)\mathrm{dt}={.Math.A.Math.}^2\exp \left(j\right){}_{-T_0/2-}^{T_0/2}\mathrm{dt}={.Math.A.Math.}^2\exp \left(j\right)\left(\frac{T_0}{}-\right)={.Math.A.Math.}^2\exp \left(j\right)\left(T_0-\right)& \left(57\right)\end{array}$

(142) The above is summarized below.

(143) For T.sub.0tT.sub.0, the autocorrelation is given as below.
R.sub.(t)R.sub.SS(t)=|A|.sup.2 exp(jt)(T.sub.0|t|) (58)

(144) An autocorrelation of the autocorrelation given in the equation (58) is derived for =0.

(145) $\begin{array}{cc}{R}_{}\left(0\right)={}_{-}^{}{R}_{}^{*}\left(t\right)R_{}\left(t\right)\mathrm{dt}={.Math.A.Math.}^4{}_{-}^{}{\left(T_0-.Math.t.Math.\right)}^2\mathrm{dt}={.Math.A.Math.}^4{}_{-T_0/}^{T_0/}{\left(T_0-.Math.t.Math.\right)}^2\mathrm{dt}=2{.Math.A.Math.}^4{}_0^{T_0/}{\left(T_0-t\right)}^2\mathrm{dt}=2{.Math.A.Math.}^4{}_0^{T_0/}\left(T_0^2-2T_{0}t+{}^2{t}^2\right)\mathrm{dt}=2{{.Math.A.Math.}^4\left[T_0^{2}t-T_0{t}^2+\frac{{}^2{t}^3}{3}\right]}_0^{T_0/}=2{.Math.A.Math.}^4\left(\frac{T_0^3}{}-\frac{T_0^3}{}+\frac{T_0^3}{3}\right)=\frac{2{.Math.A.Math.}^4{T}_0^3}{3}& \left(59\right)\end{array}$

(146) Therefore, we have the following equation.

(147) $\begin{array}{cc}{R}_{}\left(0\right)=\frac{2{\left\{R_{SS}\left(0\right)\right\}}^2{T}_0}{3}& \left(60\right)\end{array}$

(148) The Doppler coefficient is given as below from the equation (60).

(149) $\begin{array}{cc}=\frac{2{\left\{R_{\mathrm{SS}}\left(0\right)\right\}}^2{T}_0}{3R_{}\left(0\right)}& \left(61\right)\end{array}$

(150) If a pulse which has a low distance resolution and may cause a mistake between a noise peak and a correlation peak under a low SNR environment be assumed to be a PCW, the Doppler coefficient is given as following Equation.

(151) $\begin{array}{cc}=\frac{2P^2{T}_0}{3S}& \left(62\right)\end{array}$

(152) In the present example embodiment, instead of LFM, it is possible to use other waveforms such as FH (Frequency Hopping) and so on, with distance resolution higher than that of LFM. However, in this case, it is not possible to represent a target relative velocity using simple expressions as in the case of LFM or PCW as described above. Therefore, a relationship(s) among a peak value of correlation, an integrated value and a relative velocity may be obtained in advance, and then, a relative velocity may be estimated from the integrated value and the peak value of correlation using the relationship(s).

(153) In the present example embodiment, it is possible to increase estimation accuracy of a velocity in such manner that, when three or more transmission waveforms are used, for example, a plurality of pairs each having two waveforms are extracted, a correlation value between each pair is calculated, a target velocity is estimated based on a peak value of the correlation value and an integrated value of the correlation value for each pair, and the estimated velocities are averaged among the pairs.

(154) As for a received signal with a high SNR, the target velocity may be obtained from a relationship between a peak of the autocorrelation and an integrated value of the square of an absolute value of an autocorrelation obtained by autocorrelating a single pulse, without transmitting a plurality of pulses.

(155) In the above example embodiments, although described in view of a monostatic sonar that a transmitter and a receiver are arranged in the same place (transmission and reception are performed by a common acoustic array), it is possible to apply to a bistatic sonar (assembly) in which a transmitter and a receiver are arranged at a distance between them and a multi-static sonar (assembly) in which a plurality of receivers receive a signal(s) from one transmitter.

(156) FIG. 5 illustrates a configuration of an example embodiment according to the present invention, in which a sonar apparatus is implemented on a computer apparatus 200. Referring to FIG. 5, the computer apparatus 200 includes a processor 201, a memory 202 which may be a semiconductor memory, such as RAM (Random Access Memory), ROM (Read Only Memory), EEPROM (Electrically Erasable Programmable Read-Only Memory) and so on (otherwise, it may be HDD (Hard Disk Drive) or SSD (Solid State Drive) and so on), a display apparatus 203, and an interface 204 (bus interface) which connects to an acoustic array 101 and AD converter 103 as shown in FIG. 2 and so on. The processor 201 may be a DSP (Digital Signal Processor). The processor 201 performs processing of the transmission controller 102, the phase adjuster 104 and the correlation processor 105 by executing a program 205 stored in the memory 202.

(157) According to the present example embodiment, even in a low SNR situation, if there are two echoes from a target, it is possible to estimate a relative velocity in a line-of-sight direction of the target. Furthermore, in a high SNR situation, even if there is one echo from a target, it is possible to estimate a relative velocity in a line-of-sight direction of the target and no transmission waveform (replica signal) is necessary.

(158) Although the above example embodiments are described by taking a sonar system for example, application of the present invention is not limited to a sonar system and application to a radar system, a laser system or the like can be conceivable and possible.

(159) Each disclosure of the above PTL 1 and NPL 1 is incorporated herein by reference thereto. Variations and adjustments of the example embodiments and examples are possible within the scope of the overall disclosure (including the claims) of the present invention and based on the basic technical concept of the present invention. Furthermore, various combinations and selections of various disclosed elements (including the elements in each of the claims, example embodiments, examples, drawings, etc.) are possible within the scope of the claims. Namely, the present invention of course includes various variations and modifications that could be made by those skilled in the art according to the overall disclosure including the claims and the technical concept.

(160) The example embodiments described above can be described as the following supplementary notes (abbreviated as note).

(161) (Note 1)

(162) An apparatus (reception apparatus), comprising: a processor connected to a reception array and a display apparatus, by performing time-integration of a square of the correlation value wherein the processor extracts a signal in a predetermined time range from a signal received by the reception array, calculates a correlation value of a signal waveform extracted, detects a peak value of the correlation value, calculates an integrated value by performing time-integration of the square of the correlation value, derives a target velocity, based on the peak value and the integrated value, and displays the target velocity on a display apparatus.
(Note 2)

(163) The sonar apparatus (reception apparatus) according to note 1, wherein the signal received by the acoustic array is a pulse signal modulated by Linear Frequency Modulation (LFM), and wherein the processor calculates a cross-correlation value between two signal waveforms extracted, and calculates a Doppler coefficient, based on the integrated value and the peak value of the cross-correlation value, and an absolute value of a frequency change rate and a pulse length to derive the target velocity.
(Note 3)

(164) The apparatus (reception apparatus) according to note 1, wherein the signal received by the acoustic array is a pulse signal modulated by Linear Frequency Modulation (LFM), and wherein the processor calculates a Doppler coefficient based on a peak value R.sub.(0) of an autocorrelation R.sub.() of the correlation value, a peak value R.sub.aa(0) of an autocorrelation R.sub.aa() of the autocorrelation R.sub.(), and an absolute value of a frequency change rate and a pulse length to derive the target velocity.
(Note 4)

(165) The apparatus (reception apparatus) according to note 1, wherein the processor calculates a Doppler coefficient, based on a peak value R.sub.(0) of an autocorrelation R.sub.() of a cross-correlation value between two signal waveforms extracted, a peak value R.sub.aa(0) of an autocorrelation R.sub.aa() of the autocorrelation R.sub.(), and an absolute value of a frequency change rate and a pulse length to derive the target velocity.

(166) (Note 5)

(167) The apparatus (reception apparatus) according to note 1, wherein the signal received by the acoustic array is a Pulse Continuous Wave (PCW) signal, and wherein the processor, calculates a Doppler coefficient based on a peak value of the correlation value, the integrated value of an autocorrelation of the correlation value, and a pulse length, and obtains the target velocity.
(Note 6)

(168) The apparatus (reception apparatus) according to note 1, wherein relationship among at least, the peak value, the integrated value, and a relative velocity between the sonar apparatus and the target are stored in advance in a storage part, and wherein the processor, obtains the target velocity based on at least the peak value and the integrated value obtained for the signal received by the acoustic array, with reference to the storage apparatus.
(Note 7)

(169) The apparatus (reception apparatus) according to note 1, wherein the signal received by the reception array is a single pulse, and wherein the processor, calculates an autocorrelation value of the single pulse, calculates the integrated value by performing time-integration of a square of the autocorrelation value, and derives the target velocity, based on a peak value of the autocorrelation value of the single pulse and an integrated value of the square of the autocorrelation value.
(Note 8)

(170) The apparatus (reception apparatus) according to note 1, wherein the processor extracts a set of two pairs with respect to three or more extracted signal waveforms; calculates a correlation value between each of the pairs and obtains the target velocity based on the peak value and an integrated value of the square of the correlation value for each of the pairs; and performs averaging of the target velocities among the pairs
(Note 9)

(171) The apparatus (reception apparatus) according to any one of notes 1 to 3, wherein the processor controls to output a transmission signal at a predetermined transmission interval, and calculates the correlation value between two signal waveforms extracted, a time interval therebetween corresponding to the transmission interval for a plurality of extracted signals.

(172) (Note 10)

(173) The apparatus (reception apparatus) according to note 2 or 3, wherein the predetermined time range is longer than a pulse length of the received signal.