JP2011203276A - Doppler measuring instrument and tidal current meter - Google Patents

Abstract

PROBLEM TO BE SOLVED: To provide a Doppler measuring instrument having superior noise-resistant performance, capable of measuring a Doppler shift amount with high accuracy.SOLUTION: Each gravity center calculation section Wt [1:n] is a section including a peak of a power spectrum Pt [fi], and a gravity center frequency fwt [1:n] is calculated in each section Wt [1:n]. A tentative Doppler shift amount is a frequency shift amount when an output in mutual correlation processing of the power spectra Pt [fi], Pr [fi] gets maximum. A gravity center calculation section Wr [1:n] correlated with the gravity center calculation section Wt [1:n] is determined based on the tentative Doppler shift amount, and a gravity center frequency fwr [1:n] of the power spectrum Pr [fi] is calculated in each section Wr [1:n]. Compensation of frequency dependency of a Doppler shift or the like is executed for a difference between the gravity center frequencies fwr [1:n], fwt [1:n] to determine the Doppler shift amount fd. A frequency range of the gravity center calculation section Wr [1:n] is corrected by the Doppler shift amount fd, and the Doppler shift amount fd is determined repeatedly until converged.

Description

  The present invention relates to a Doppler measuring instrument that measures the amount of Doppler shift of ultrasonic waves transmitted into water, and a tide meter that includes the Doppler measuring instrument.

Traditionally, tidal currents have been measured with a tide meter for fishing assistance and marine research. The tide meter is attached to the bottom of a ship, and includes a transducer that transmits and receives ultrasonic waves at a certain depression angle θ with respect to directions that are 120 degrees apart from each other in the horizontal direction, for example. Then, the ship speed (vs. water speed) with respect to the tidal current at that depth is obtained from the Doppler shift amount generated in the reflected wave (vs. water echo) coming from countless scatterers (plankton, etc.) in the sea located at the set depth. . Further, the speed of the ship (ground speed) with respect to the sea floor is obtained from the amount of Doppler shift generated in the reflected wave (ground echo) from the sea floor. Furthermore, the speed of the tidal current is obtained from the difference between the speed of the ground ship and the speed of the water ship. Assuming that the ship speed (or ship speed) is V, the underwater sound speed is c, the ultrasonic frequency is f0, and the Doppler shift amount is fd, c is sufficiently larger than V. Ship speed) V is approximated by the following equation (1).
V = fd · c / (2f0 · cosθ) (1)
Therefore, in order to accurately measure the speed V against the water (or the speed V against the ground), it is necessary to measure the Doppler shift amount fd with high accuracy.

  In Patent Document 1 below, a first pulse and a second pulse delayed by a predetermined time from the first pulse are transmitted as an ultrasonic signal in water, and a first quadrature phase sample and a second pulse of an echo signal by the first pulse are transmitted. It has been proposed to generate an autocorrelation of the echo signal with the second quadrature phase sample and to obtain a Doppler shift amount based on the phase shift obtained by the autocorrelation. It has also been proposed to use wideband coded pulses as the first and second pulses. In Patent Document 2 below, a Fourier transform is performed on a data string indicating a phase difference between a zero-cross signal (binary signal) generated from a received echo signal and a reference signal generated inside the apparatus, and obtained by Fourier transform. It has been proposed to obtain the Doppler shift amount based on the Fourier spectrum obtained.

Japanese Patent No. 2695989 (Claim 7, page 6, left column, line 25 to page 11, right column, line 32) Japanese Patent No. 3028376 (Claim 1, page 5, right column, line 32 to page 6, right column, line 11)

  However, since the method as in Patent Document 1 uses all information of the received signal band, the Doppler shift amount can be measured with high accuracy when the noise included in the received signal is small. There is a problem that the measurement accuracy decreases. The noise includes noise caused by waves and screws, and noise caused by ultrasonic waves transmitted from other ships. Further, as a disadvantage of the method as in Patent Document 1, since the autocorrelation calculation cannot consider the frequency dependency of the Doppler shift amount, the frequency dependency of the Doppler shift amount is not particularly compensated for wideband signals. There is a problem. On the other hand, the method as disclosed in Patent Document 2 has a problem that the measurement accuracy cannot be increased because a received signal having one spectrum peak is an analysis target and the Doppler shift amount is measured from the single spectrum peak. In addition, in order to improve the measurement accuracy, ultrasonic waves were repeatedly transmitted, and the Doppler shift amount obtained from each echo signal was subjected to statistical processing such as averaging, and the final Doppler shift amount was output. Therefore, it takes time to measure the Doppler shift amount, and the throughput decreases.

  The present invention solves the above-described problems, and an object of the present invention is to provide a Doppler measuring instrument that is excellent in noise resistance and that can measure the Doppler shift amount with high accuracy.

  In the first invention, in the Doppler measuring device that measures the Doppler shift amount of the ultrasonic wave transmitted into the water, the cross-correlation process between the spectrum of the ultrasonic transmission waveform and the spectrum of the ultrasonic reception signal is performed, The Doppler shift amount is obtained based on the frequency at which the output of the correlation processing is maximized. Here, the frequency corresponds to the provisional Doppler shift amount shown in the embodiment.

  By doing so, the peak position information of the received signal (echo signal) (peak position information in the frequency domain) from the spectrum of the received signal including aperiodic noise using the spectrum of the transmission waveform and the spectrum of the received signal. Is extracted. That is, when there are a plurality of peaks, the Doppler shift amount of each peak is collected in one place and extracted. As a result, the above-mentioned frequency has high noise resistance performance, and therefore has little variation, and is also substantially equal to the true Doppler shift amount (or a value obtained by adding a constant to the true Doppler shift amount), so that high noise resistance performance is achieved. And a high measurement accuracy can be realized.

  In the second invention, the ground ship speed and the water ship speed are obtained by measuring the Doppler shift amount of the echo signal caused by the ultrasonic wave transmitted underwater, and the tidal current is determined from the ground ship speed and the water ship speed. A tidal meter for determining speed includes the Doppler shift device according to claim 1, which measures a Doppler shift amount of an echo signal based on an ultrasonic transmission waveform and an echo signal, and measured by the Doppler shift device. The ground ship speed and the water ship speed are obtained based on the quantity.

  By doing in this way, the tidal meter which can obtain | require the tidal velocity which is excellent in noise-proof performance and is highly accurate is obtained.

  ADVANTAGE OF THE INVENTION According to this invention, it is excellent in noise-proof performance, and also can implement | achieve the Doppler measuring device and tidal meter which can measure a Doppler shift amount with high precision.

It is a figure which shows the structure of the tide meter incorporating a Doppler measuring device. It is a figure which shows the structure of the Doppler measuring device which concerns on this invention. It is a figure which shows the waveform of a transmission signal. It is a figure which shows the power spectrum of a transmission waveform. It is a figure which shows the power spectrum of a received signal. It is a figure which shows the result of the cross correlation process of two power spectra. It is a figure which shows the aspect which narrows a gravity center calculation area. It is a figure which shows another wideband transmission signal.

  Hereinafter, embodiments of the present invention will be described with reference to the drawings. FIG. 1 shows a configuration of a tide meter incorporating a Doppler measuring instrument, FIG. 2 shows a configuration of a Doppler measuring instrument according to the present invention, and FIG. 3 shows a waveform of a transmission signal. In this tide meter 1, ultrasonic waves are transmitted from the three transducers 11 to the sea at a constant depression angle with respect to directions whose horizontal directions differ from each other by 120 °, and echo signals reflected by the seabed and countless scatterers in the sea. Is received by each transducer 11. The transducer 11 is driven by the transmission signal generated by the transmission waveform generator 14 and amplified by the transmission amplifier 13 via the transmission / reception switching circuit 12. The transmission waveform generation unit 14 also supplies the waveform (transmission waveform) of the transmission signal to the Doppler measuring device 3. This transmission waveform consists of a waveform data string equivalent to that obtained by A / D conversion of the transmission signal.

  The transmission signal (also referred to as a transmission pulse) shown in FIG. 3 is a transmission signal in which a plurality of wideband signals encoded by M-sequence BPSK (maximum sequence binary phase shift keying) are connected. It is transmitted from the transducer 11. The transmission signal is composed of four identical elements, and each element is composed of seven sub-pulses that are encoded as described above. In the figure, Ta is the time width of the transmission signal, Tb is the time width of the element, and Tc is the time width of the subpulse. The time width Ta is about 0.7 ms, for example, and the frequency of the subpulse (carrier frequency) is about 250 kHz, for example. If another code pattern (for example, + 1 + 1 + 1-1 + 1-1-1) in the M sequence is used, the envelope shape of the power spectrum of the transmission waveform shown in FIG. 4 changes.

  A reception signal including an echo signal received by each transducer 11 is amplified by a reception amplifier 15 and A / D converted by an A / D converter 16. The A / D converted sample data series is stored in the buffer memory 17. The Doppler measuring device 3 reads the sample data of the echo signal at one or a plurality of set depths in the received signal and the echo signal from the seabed from the buffer memory 17 and obtains the respective Doppler shift amounts fd based on the sample data. Output. The tidal current calculation unit 18 calculates the ground ship speed and the water ship speed by substituting the Doppler shift amounts obtained from the three received signals into the above equation (1), and further the ground ship speed and the water ship speed. Based on the above, the speed and direction of the tidal current are obtained. The display unit 19 displays, for example, arrows based on the navigation route at predetermined intervals, the direction of the tidal current is indicated by the direction of the arrow, and the speed of the tidal current is indicated by the length of the arrow.

  The control unit 20 includes a CPU 20a, a DSP (digital signal processor) 20b, a memory (program memory and data memory) 20c, and the like, and performs various calculations and control of each unit of the tide meter 1. Further, when a key of an operation unit (not shown) is pressed, the control unit 20 controls each unit according to the type of the pressed key. For example, when the set depth value is input from the key, the input value is stored in the memory 20c, and the tidal current at the set depth is measured.

  Next, the Doppler measuring device 3 for measuring the Doppler shift amount generated in the echo signal will be described with reference to FIG. A DFT (Discrete Fourier Transform) unit 31 performs discrete Fourier transform on the transmission waveform using a fast Fourier transform algorithm to convert the transmission waveform in the time domain into an amplitude spectrum in the frequency domain. Further, the power spectrum shown in FIG. 4 is generated by squaring the amplitude spectrum. The power spectrum Pt [fi] (fi is the frequency of the power spectrum Pt [i]) is generated with a resolution of 10 Hz, and is symmetric with respect to the peak 51, such as five peaks 52, 53, etc. There is a small peak 54. Since the spectrum is not an amplitude spectrum but a power spectrum, the peak shape is sharp. Further, due to the nature of the discrete Fourier transform and the transmission waveform, the power spectrum Pt [fi] is distributed in the range of 2 / Tc (FIG. 3), and the center frequency fc, that is, the frequency fc of the peak 51 is equal to the frequency of the subpulse. The frequency difference between adjacent peaks is 1 / Tb, and the width of each peak (zero cross width of each peak) is 2 / Ta.

First, the center-of-gravity frequency calculation unit 32 is set so that each peak 51 and the like are positioned at the center of each center-of-gravity calculation section Wt [k] (k = 1 to n) (refer to FIG. 4, also expressed as Wt [1: n]). The center of gravity calculation section Wt [k] is determined. This center-of-gravity calculation section Wt [k] is defined by a lower limit value and an upper limit value (or by a lower limit value and a frequency width) of the frequency, and the frequency width is equal to the frequency difference (1 / Tb) between adjacent peaks. Next, using the power spectrum Pt [fi] as a weight, the centroid frequency fwt [k] (k = 1 to n) is calculated for each centroid calculation section Wt [k] by the following equation.
fwt [k] = Σ (Pt [fj] · fj) / ΣPt [fj]
Here, Pt [fj] is a power spectrum belonging to the centroid calculation section Wt [k].

  Note that the time width Tb of the element of the transmission signal is determined so that the width of the center-of-gravity calculation section Wt [k] is at least twice the assumed maximum value of the Doppler shift amount. This is because the center of gravity calculation section Wt [k] of the transmission waveform can be associated with the center of gravity calculation section Wr [k] of the received signal described later. Further, every time an ultrasonic wave is transmitted, the power spectrum Pt [fi], the centroid calculation section Wt [k] and the centroid frequency fwt [k] do not change, so the power spectrum Pt [fi] and the like are obtained in advance. The value stored in the memory 20c is used when the Doppler shift amount is obtained.

  Next, the receiving system will be described. The DFT unit 33 performs a discrete Fourier transform using a fast Fourier transform algorithm on a sample data sequence corresponding to a time width Ta corresponding to an echo signal from a set depth, which is stored in the buffer memory 17, and The received signal (echo signal) is converted into an amplitude spectrum in the frequency domain. Further, the power spectrum shown in FIG. 5 is generated by squaring the amplitude spectrum. This power spectrum Pr [fi] (fi is the frequency of the power spectrum Pr [i]) is also generated with a resolution of 10 Hz, and the peak of the power spectrum Pt [fi] of the transmission waveform is respectively to the left and right of the central peak 61. There are five peaks 62, 63, etc. corresponding to 52, 53, etc. However, the thing corresponding to the peak 54 is buried in noise and cannot be discriminated.

  Further, it can be seen from FIG. 5 that the center frequency of the power spectrum Pr [fi], that is, the frequency of the peak 61 is shifted by Δfa from the center frequency fc of the power spectrum Pt [fi] of the transmission waveform. This Δfa is an approximate Doppler shift amount. The frequency difference between adjacent peaks is substantially equal to the frequency difference (1 / Tb) between the peaks in FIG. 4, but is slightly larger or smaller than 1 / Tb due to the frequency dependence of the Doppler shift. Further, the power spectrum Pr [fi] obtained here is stored in the memory 20c so that it can be used by the center-of-gravity frequency calculation unit 42 and the weighting factor calculation unit 46 described later.

  The provisional Doppler shift amount detection unit 34 includes a cross correlation processing unit 35 and a peak detection unit 36, and calculates a provisional Doppler shift amount from the power spectrum Pt [fi] of the transmission waveform and the power spectrum Pr [fi] of the received signal. Find and output. The cross-correlation processing unit 35 performs a cross-correlation process between the power spectra Pt [fi] and Pr [fi]. That is, while shifting the power spectrum Pt [fi] by the resolution (10 Hz) with respect to the power spectrum Pr [fi], the product-sum operation of both power spectra is performed in each shift state, and the calculation result is output. FIG. 6 shows the cross-correlation output as the calculation result. The maximum peak 71 located at the center is when the peak 51 of the power spectrum Pt [fi] coincides with the peak 61 of the power spectrum Pr [fi], and is shifted from the center frequency fc of the transmission waveform by Δfb. . The peak 72 is obtained when the peak 51 coincides with the peak 62, and the peak 73 is obtained when the peak 51 coincides with the peak 63. As a cross-correlation process, the power spectra to be correlated are set to the same number of discrete points by padding with zeros, etc., and then the discrete Fourier transform is performed. One of them takes a complex conjugate to perform product operation, and the operation result is subjected to inverse discrete Fourier transform. Also good.

The peak detection unit 36 detects a peak 71 having the maximum value among the cross-correlation outputs output from the cross-correlation processing unit 35, and temporarily calculates a frequency difference Δfb obtained by subtracting the center frequency fc of the transmission waveform from the frequency of the peak 71. Output as Doppler shift amount. Although the center frequency of the cross-correlation output is fc here, if the cross-correlation processing is performed so that the center frequency of the cross-correlation output becomes 0 Hz, the frequency of the peak 71 itself becomes the provisional Doppler shift amount. Further, the temporary Doppler shift amount may be obtained by another method instead of the above method. For example, the centroid frequency of the cross-correlation output is calculated, and the provisional Doppler shift amount is obtained from the centroid frequency. In this case, if the cross-correlation output is L [fi], the provisional Doppler shift amount is given by the following equation.
Provisional Doppler shift amount = (Σ (L [fi] · fi) / ΣL [fi]) − fc
In the present invention, when calculating the provisional Doppler shift amount, the frequency of an arbitrary peak at which the output of the cross-correlation processing is maximized can be used, and the frequency of the peak 71 located in the center of FIG. There is no. Therefore, the temporary Doppler shift amount may be calculated based on the frequencies of the peak 72 and the peak 73.

  In the cross-correlation processing described above, the peak position of the periodic echo signal is derived from the power spectrum Pr [fi] including aperiodic noise using the periodicity of the peaks in the power spectra Pt [fi] and Pr [fi]. Information (peak position information in the frequency domain) is extracted. That is, when there are a plurality of peaks, the Doppler shift amount of each peak is collected in one place and extracted. As a result, the noise resistance performance is improved by the cross-correlation processing, and the provisional Doppler shift amount obtained as described above has a small variation and is substantially equal to the true Doppler shift amount. Therefore, the temporary Doppler shift amount can also be used as the final Doppler shift amount. Although the frequency dependence of the Doppler shift is not considered here, the provisional Doppler shift amount described above is a Doppler shift amount that can withstand practical use.

  As described above, the provisional Doppler shift amount can be set as the final Doppler shift amount. However, in the present embodiment, the provisional Doppler shift amount is given as an initial value to the center-of-gravity frequency processing unit 40, which will be described later. A highly accurate Doppler shift amount fd is obtained. The initial value is preferably a value close to the true Doppler shift amount (for example, the provisional Doppler shift amount). However, since the Doppler shift amount fd is increased in accuracy by the center-of-gravity frequency processing, a rough value such as a conventional value is used. It is also possible to set the Doppler shift amount or 0 Hz obtained by the autocorrelation processing of the echo signal as an initial value. Further, the shift amount at the frequency position where the spectrum becomes the maximum value may be used instead of the above-described maximum value of the cross-correlation processing output. However, when the rough value is set as the initial value, the amount of deviation between the initial value and the true Doppler shift amount is large, so that the above Doppler shift amount fd converges compared to the case where the temporary Doppler shift amount is set as the initial value. The number of repetitions of the center-of-gravity frequency processing until this is increased. Further, when the S / N of the received signal is poor and the Doppler shift amount obtained by autocorrelation processing or the like is greatly deviated from the true Doppler shift amount due to noise, the Doppler shift amount fd may not converge.

  Next, the gravity center frequency processing unit 40 will be described. The center-of-gravity frequency processing unit 40 includes a center-of-gravity calculation section determination unit 41 and the like, and is output from the weighted average processing unit 44 using the provisional Doppler shift amount as an initial value after the cross-correlation processing described above is performed once. The center-of-gravity frequency processing is repeated until the Doppler shift amount fd is converged. Convergence refers to a state in which the difference between the previous value and the current value of the Doppler shift amount fd is equal to or less than a predetermined value (including zero). The center-of-gravity frequency process is a process performed by the center-of-gravity calculation section determination unit 41 or the like. Hereinafter, the center-of-gravity frequency processing performed repeatedly is referred to as “repetitive processing”. The advantage of cross-correlation processing is that it is resistant to noise, and the advantage of center-of-gravity frequency processing is that measurement accuracy is high. Therefore, by first performing cross-correlation processing and eliminating the influence of noise, calculating the Doppler shift amount by performing centroid frequency processing on the range near the peak of the spectrum with good S / N in the frequency domain, In addition to improving noise resistance, measurement accuracy can be improved.

  The center-of-gravity calculation section determination unit 41 uses the provisional Doppler shift amount output from the peak detection unit 36 to calculate the center-of-gravity calculation section Wr [k] (k = 1 to n) (FIG. 5) is determined such that the center of the center of gravity calculation section substantially coincides with the peak 61. If expressed accurately, the center-of-gravity calculation section Wr [k] is on the right side (high frequency side) of the above-described center-of-gravity calculation section Wt [k] (temporary Doppler shift amount (Δfb (FIG. 6) instead of Δfa (FIG. 5)). It is only shifted. Here, the width of each centroid calculation section Wr [k] is made equal to the width of the centroid calculation section Wt [k], but may be somewhat narrower or wider than the width of the centroid calculation section Wt [k]. When the Doppler shift amount fd is output from the convergence determination unit 45 described later, the center of gravity calculation section Wr [k] is similarly determined using the Doppler shift amount fd. That is, the center-of-gravity calculation section Wr [k] is corrected by the Doppler shift amount fd.

The center-of-gravity frequency calculation unit 42 calculates the center-of-gravity frequency fwr [k] (k = 1 to n) for each center-of-gravity calculation section Wr [k] using the power spectrum Pr [fi] as a weight.
fwr [k] = Σ (Pr [fj] · fj) / ΣPr [fj]
Here, Pr [fj] is a power spectrum belonging to the centroid calculation section Wr [k]. At this time, if the peak 62 or the like of the power spectrum Pr [fi] is positioned at the approximate center of the centroid calculation section Wr [k], noise and side lobes that appear before and after the peak 62 and the like are canceled out, and the calculated centroid frequency fwr [k] is approximately equal to the true frequency after receiving the Doppler effect. On the other hand, when the peak 62 and the like are shifted to the left (right) from the center of the centroid calculation section Wr [k], the calculated centroid frequency fwr [k] is higher than the true frequency after receiving the Doppler effect ( Lower). The center of gravity calculation section Wr [k] is corrected by the Doppler shift amount fd output from the convergence determination unit 45, whereby the center of the center of gravity calculation section Wr [k] approaches the peak 62 of the power spectrum Pr [fi]. The center-of-gravity frequency fwr [k] approaches the true frequency after receiving the Doppler effect to be obtained.

The Doppler shift correction unit 43 determines the difference (Doppler shift amount) between the centroid frequency fwr [k] calculated by the centroid frequency calculation unit 42 and the centroid frequency fwt [k] calculated by the centroid frequency calculation unit 32 as a transmission waveform. Is converted into a Doppler shift amount fd [k] (k = 1 to n) corresponding to a Doppler shift at the center frequency (carrier frequency) fc. This conversion is expressed by the following formula.
fd [k] = (fwr [k] −fwt [k]) · fc / fwt [k]
The Doppler shift amount fd [k] is obtained by correcting the Doppler shift amount (fwr [k] −fwt [k]) in which the frequency dependency of the Doppler shift is not considered, that is, the frequency dependency is compensated. is there. Thus, since the frequency dependence of the Doppler shift is compensated, the measurement accuracy of the Doppler shift amount fd described later can be improved.

  The weighting factor calculation unit 46 calculates the total value s [k] (k = 1 to n) of the power spectrum Pr [fi] for each centroid calculation section Wr [k], that is, substantially the integral value, and calculates the total. A weight coefficient w [k] (k = 1 to n) of each centroid calculation section Wr [k] is calculated from the value s [k]. When the power spectrum in each centroid calculation section Wr [k] is Pr [fj], the total value s [k] is represented by ΣPr [fj], and the weight coefficient w [k] is s [k] / Σs [k It is represented by Here, the weighting coefficient w [k] is calculated from the power spectrum Pr [fi], but it is also possible to calculate the weighting coefficient w [k] from the power spectrum Pt [fi] of the transmission waveform not including noise. is there. However, it has been experimentally confirmed that the measurement accuracy of the Doppler shift amount fd is somewhat lowered in this way.

The weighted average processing unit 44 performs weighted average processing of the Doppler shift amount fd [k] using the weighting factor w [k] calculated by the weighting factor calculating unit 46 as a weight to obtain the Doppler shift amount fd. The weighted average process is expressed by the following equation.
fd = Σ (fd [k] · w [k]) / Σw [k]
By this weighted average processing, the highly accurate Doppler shift amount fd with reduced influence of noise is obtained according to the magnitude and distribution of the power spectrum of the received signal, that is, according to the accuracy of information characterizing the received signal (echo signal). Desired. Instead of the above-described weighting factor w [k], an integral value for each centroid calculation section Wt [k] of the power spectrum Pt [fi] of the transmission waveform can be used as the weighting factor. In the above example, the weighted average process is performed using the integral value of the spectrum for each centroid calculation section as a weight, but the weighted average process may be performed using the maximum value of the spectrum for each centroid calculation section as a weight. . Further, the spectrum for each centroid calculation section may be simply averaged without weighting.

  The convergence determination unit 45 determines whether or not the Doppler shift amount fd output from the weighted average processing unit 44 has converged. When it determines with having converged, the Doppler shift amount fd is passed to the tidal current calculation part 18 (FIG. 1). On the other hand, when it is determined that it has not converged, the Doppler shift amount fd is passed to the center-of-gravity calculation section determination unit 41. Correction of the center-of-gravity calculation section Wr [k] by the Doppler shift amount fd, calculation of the center-of-gravity frequency fwr [k] of the corrected center-of-gravity calculation section Wr [k] is performed, and 2 until the Doppler shift amount fd converges. The Doppler shift amount fd for the third time, the third time,... Is calculated. That is, the iterative process is performed until the Doppler shift amount fd converges.

  Specifically, the convergence determining unit 45 performs the iterative process for the convergence of the Doppler shift amount fd a predetermined number of times (for example, twice), and then the difference between the previous and current Doppler shift amount fd is a predetermined value (for example, 1 Hz) or less, it is determined that the calculated Doppler shift amount fd has converged, and the current Doppler shift amount fd is output as the final value. On the other hand, if the Doppler shift amount fd does not converge even after a predetermined number of repetitions (for example, 30), error processing is performed. Further, as described above, since the temporary Doppler shift amount is a value close to the true Doppler shift amount, the temporary Doppler shift amount is output as the final Doppler shift amount fd without passing through the centroid frequency processing unit 40. May be.

  Next, the reason for performing the repetition process will be described. The first reason is that the frequency dependence of the Doppler shift is not considered in the provisional Doppler shift amount, and the center frequency of the centroid calculation section Wr [k] including the peak 62 other than the peak 61 is the power spectrum Pr [fi. This is because the center-of-gravity frequency fwr [k] is affected by the noise component and a slight error occurs in the calculated Doppler shift amount fd. As described above, this error is obtained by correcting the centroid calculation section Wr [k] by the calculated Doppler shift amount fd, and performing the Doppler shift based on the centroid frequency fwr [k] of the corrected centroid calculation section Wr [k]. The amount fd is gradually reduced by calculating the amount fd. That is, the measurement accuracy of the Doppler shift amount fd is improved by iterative processing.

  The second reason is that an error occurs in the calculated Doppler shift amount due to a frequency error caused by the number of discrete data when the received signal is sampled in the provisional Doppler shift amount. This error is also sufficiently reduced by performing the above-described repetitive processing. In addition, if the centroid frequency of the cross-correlation output is set as the provisional Doppler shift amount as described above instead of using peak detection, the computation processing time increases to obtain the provisional Doppler shift amount, but the number of discrete data at the time of sampling The error due to is not generated.

  As described above, the peak 62 of the power spectrum Pr [fi] of the received signal (echo signal) corresponding to the peak 52 of the power spectrum Pt [fi] included in the centroid calculation section Wt [k] is the centroid calculation section. It is included in the center-of-gravity calculation section Wr [k] associated with Wt [k]. Moreover, both peaks 52, 62, etc. are information indicating the characteristics of the transmission signal and the reception signal having substantially the same waveform. Therefore, a practically accurate Doppler shift amount fd is obtained based on the difference between the centroid frequency fwt [k] for each centroid calculation section Wt [k] and the centroid frequency fwr [k] for each centroid calculation section Wr [k]. be able to.

  In addition, since the power spectra Pr [fi] and Pt [fi] have a plurality of peaks, the amount of information for obtaining the Doppler shift amount fd increases, and the measurement accuracy of the Doppler shift amount fd increases. Furthermore, the measurement accuracy of the Doppler shift amount fd can be further improved by the above-described repetitive processing, the compensation of the frequency dependence of the Doppler shift performed by the Doppler shift correction unit 43, and the weighted average processing performed by the weighted average processing unit 44. it can. In the above embodiment, the configuration of the Doppler measuring instrument 3 has been described using functional blocks. However, the Doppler measuring instrument 3 is actually realized by software by the CPU 20a and the DSP 20b.

  In the embodiment described above, the width (frequency width) of the center-of-gravity calculation section Wr [k] is constant even when the number of iterations increases. However, the width is gradually narrowed as the number increases. (For example, the current width may be 90% of the previous width, and the width may be narrowed within a range that does not fall below the peak width (2 / Ta)). FIG. 7 is a diagram for explaining the above operation when there are three centroid calculation sections. The initial centroid calculation section of the power spectrum of the received signal shown in (b) has the same width as the centroid calculation section of the power spectrum of the transmission waveform shown in (a), and is further to the right by the provisional Doppler shift amount DA. Has moved. The upward arrow indicates the initial center-of-gravity frequency. Then, the width of the second center of gravity calculation section shown in (c) is made narrower than the first time, and is obtained based on the first center of gravity calculation section (b) instead of the center of gravity calculation section (a) of the power spectrum of the transmission waveform. The center of the center-of-gravity calculation section is moved rightward by the amount of Doppler shift DB.

  In this way, the target range when calculating the center-of-gravity frequency is narrowed to the vicinity of the peak of the spectrum, and the influence of noise of frequency components away from the peak of the power spectrum Pr [fi] can be further reduced. The measurement accuracy of the Doppler shift amount can be further increased while maintaining the noise resistance performance. Here, since the Doppler shift is originally expansion and contraction in the frequency axis direction of the spectrum waveform with 0 Hz as a fixed point, the width of the center of gravity calculation section is changed according to the center frequency of each section, and the center of gravity calculation section is changed on the high frequency side. It is preferable that the center-of-gravity calculation section be narrow on the low frequency side. Instead of moving the second center of gravity calculation section by the Doppler shift amount DB in a lump as described above, the second time so that the first calculated center of gravity frequency is located at the center of the second center of gravity calculation section. Although it is possible to individually determine the center value of the centroid calculation section for each centroid calculation section, in the case where noise appears only in a specific centroid calculation section, there is a risk of being susceptible to noise.

  In the above-described embodiment, a transmission signal in which a plurality of M-sequence BPSK signals that are one of wideband signals that cause a plurality of peaks in the power spectrum Pt [fi] is used as a transmission signal. A signal in which a plurality of BPSK signals of the above sequence are connected, or another wideband signal, for example, a linear FM signal whose frequency continuously changes from fmax to fmin (or from fmin to fmax) over time within a certain time T 8 (FIG. 8 (a)), or signals (FIG. 8 (b)) in which sine wave signals having different frequencies (f1-f3) are connected for a certain period of time (T1-T3) (2) having different frequencies. The present invention can also be implemented by using a signal in which two or more sine wave signals are superimposed as a transmission signal. That is, a wideband signal that generates a plurality of discrete peaks in the power spectrum of a transmission waveform obtained by discrete Fourier transform can be used as a transmission signal.

  Furthermore, in the above-described embodiment, a transmission signal in which a plurality of M-sequence BPSK signals that generate a plurality of peaks in the power spectrum Pt [fi] is used as a transmission signal. However, a narrow signal that generates only one peak in the power spectrum is used. The present invention can also be implemented using a band signal, for example, a sine wave signal as a transmission signal. In this case, the center of gravity calculation interval Wt [k], Wr [k] of the transmission waveform and the reception signal is one, but the Doppler shift amount fd can be obtained by performing the same processing as in the case of the wideband signal. However, since the power spectra Pt [fi] and Pr [fi] have one peak and the amount of information related to the frequency of the transmission waveform and the received signal (echo signal) that can be used for frequency detection is reduced, the Doppler shift amount fd The accuracy of the measured value (calculated value) is somewhat inferior to the case where a broadband signal is used.

  Further, in the above-described embodiment, the Doppler shift amount fd is obtained using the power spectra Pt [fi] and Pr [fi], but a spectrum such as an amplitude spectrum or a power spectrum of 3 or more can also be used. If the amplitude spectrum is used, the peak width (width at a level lower by 3 dB from the maximum value of the main lobe, for example) becomes wider, so the accuracy of the calculated Doppler shift amount may be somewhat lowered. On the other hand, in the power spectrum of 3 or more, since the peak width due to the echo signal becomes narrower and the peak level due to the noise component becomes higher, the calculated center-of-gravity frequency fwr [k] is easily affected by noise and Noise performance may be somewhat degraded.

  Moreover, although the said embodiment demonstrated about the example which calculates a spectrum using a discrete Fourier transform, in this invention, various systems other than a discrete Fourier transform can be used as a spectrum calculation means. For example, there are a Welch method as a nonparametric method, a Yule-Waker AR method as a parametric method, and a MUSIC method as a subspace method. Even if these methods are adopted, it is possible to obtain a measurement accuracy substantially equivalent to that in the case of discrete Fourier transform.

  Furthermore, although the case where the Doppler measuring instrument 3 was incorporated in the tidal current meter 1 was demonstrated in the said embodiment, the Doppler measuring instrument 3 of this invention can be used also in the underwater communication system etc. which consist of a transmitter and a receiver. . In this underwater communication system, a transmitter transmits, for example, a reference signal whose transmission waveform is known and an information signal whose frequency changes depending on transmission information. The Doppler measuring instrument 3 is incorporated in the receiver, and the Doppler measuring instrument 3 receives the power spectrum, the centroid calculation section and the centroid frequency of the transmission waveform of the reference signal stored in the storage unit, and the received reference signal. And the Doppler shift amount is obtained as described above. Then, the frequency of the received information signal is corrected by the amount of Doppler shift, and transmission information is obtained from the corrected frequency. As an example of the above-described underwater communication system, there is a fishnet depth meter including a transmitter attached to a fishnet and a receiver attached to the bottom of a ship. For example, a sine wave signal whose frequency changes according to the measured water pressure Is transmitted as an information signal.

1 Tidal current meter 3 Doppler measuring instrument 31 DFT section (for transmission waveform)
32 Center-of-gravity frequency calculator (for transmission waveform)
33 DFT section (for received signal)
34 Temporary Doppler shift amount detection unit 35 Cross correlation processing unit 36 Peak detection unit 40 Center of gravity frequency processing unit 41 Center of gravity calculation section determination unit 42 Center of gravity frequency calculation unit (for reception signal)
43 Doppler shift correction unit 44 Weighted average processing unit 45 Convergence determination unit 46 Weight coefficient calculation unit Pt, Pr Power spectrum Wt [1: n] Center of gravity calculation section (first center of gravity calculation section)
Wr [1: n] centroid calculation section (second centroid calculation section)
fwt [1: n], fwr [1: n] center of gravity frequency w [1: n] weighting factor fd Doppler shift amount

Claims (2)

In a Doppler measuring instrument that measures the amount of Doppler shift of ultrasonic waves transmitted into the water,
Perform cross-correlation processing between the spectrum of the transmission waveform of the ultrasonic wave and the spectrum of the reception signal of the ultrasonic wave,
A Doppler measuring instrument characterized in that a Doppler shift amount is obtained based on a frequency at which the output of cross-correlation processing becomes maximum. In a tide meter that determines the speed of the ground ship and the speed of the water ship by measuring the amount of Doppler shift of the echo signal caused by the ultrasonic wave transmitted into the water, and determines the speed of the tidal current from the speed of the ground ship and the speed of the water ship ,
The Doppler measuring instrument according to claim 1, wherein the Doppler shift amount of the echo signal is measured based on the transmission waveform of the ultrasonic wave and the echo signal,
A tidal current meter characterized in that the ground ship speed and the water ship speed are obtained based on a Doppler shift amount measured by the Doppler measuring instrument.

Images