METHOD AND APPARATUS FOR EXTRACTING PHYSICAL PARAMETERS FROM AN ACOUSTIC SIGNAL

Abstract: Method and apparatus for extracting at least one physical parameter from a signal. The signal, preferably an acoustic signal, is received by at least two transducers, which are located above a water surface or above land, is converted into an electric signal. The electric signal is then digitised and transformed to a complex digital signal. From the phase history and the amplitude history of the complex digital signal the physical parameter or parameters are determined.

Description: METHOD AND APPARATUS FOR EXTRACTING PHYSICAL PARAMETERS FROM AN ACOUSTIC SIGNAL Method and apparatus for measuring sound The invention relates to a method for extracting at least one physical parameter from at least one first signal, where the signal is received by at least two first transducers which are located above a water surface or above land, the signal is converted into an electric signal, and the electric signal is applied to a signal processor, in which signal processor the electric signal is converted by an analog-to-digital converter into at least one first digital signal. The invention further relates to an apparatus for extracting at least one physical parameter from an signal, comprising at least two receiver means which are located above a water surface or above land for converting the signal into an electric signal, a signal processor and a digital signal processor, connected to the receiver means and provided with an analog-to-digital converter. In general, such a method and apparatus are known. These methods and apparatus have the disadvantage, however, that for a low signal-to-noise ratio it is difficult, if not impossible, to derive information from a digital signal. It has been tried to overcome these problems by using an analog phase network. This method entails analog and digital circuitry enlarging the relative phase errors. Especially the absolute phase error will be sensitive for operational conditions. Thus, the use of an analog phase network has the disadvantage that it can only be used on electric signals with a small frequency bandwidth. The method and the apparatus according to the invention obviate these drawbacks. The method of the known type mentioned in the opening paragraph is characterised, according to an aspect of the invention, in that the digital signal is transformed into a complex digital signal, and physical parameters are determined on the basis of a phase history and an amplitude history of the complex digital signal. The apparatus of the known type mentioned in the second paragraph according to the invention is characterised in that the digital signal processor is provided with a means for deriving a first complex digital signal, a converter means connected to the means for deriving a first complex digital signal for converting the complex digital signal into an amplitude signal and a phase signal, an amplitude memory connected to the converter means for storing the amplitude signal and a phase memory means for storing the phase signal. The invention provides a method and an apparatus which is capable of extracting physical parameters from a signal having a low signal-to-noise ratio. In fact, the method and apparatus are capable of extracting physical parameters from signals having a signal-to-noise ratio of 20 dB, and even lower than OdB. Furthermore, the method and apparatus are capable of doing so with signals having a broad frequency bandwidth. In general this means that the bandwidth is at least 25% of the mean frequency, but depends also on the required accuracy in the phase calculation. Furthermore, the method can be used and the apparatus can be build using very simple components and a standard computer means. The transducers are located above a water surface or above land, meaning that the transducers are not located below a water surface, i.e. in the water, like for example in sonar applications. The transducers can for example be land-based or air-based, or be mounted on a vehicle, in motion or stationary. The signal can be an acoustic signal or an electromagnetic signal. Preferably, the signal is an acoustic signal, more preferably a sound signal. The method and apparatus according to the current invention preferably handles signals having a rather low signal-to-noise ratio. Noise from aeroplanes is an example of a signal which is processed using the method and apparatus according to the current invention. The invention, however, is not limited to such signals. The signal can be converted into a complex digital signal using mathematical methods which are, as such, well known in the art. A preferred method is a Hubert filter, although other procedures can be used as well. The advantage of using a Hubert filter is that it is very fast when implemented on a computer and needs a relative small amount of memory. An other mathematical method which is, as such, well known in the art and which can be used for the same purpose is a Fast Fourier Transform (FFT) . The signals are transformed using a real FFT and transformed back using a complex FFT. It should be noted that the Hubert filter is known per se from, for example, "Theory and application of digital signal processing, L.R. Rabiner and B. Gold, Prentice-Hall, Inc, 1 975, page. 70-73. A further favourable embodiment of the method according to the invention, which can be applied in situations where the signal-to-noise ratio of a received signal is relatively small, is characterised in that the digital signal processor transforms the complex digital signal into a digital amplitude signal and a corresponding digital phase signal. According to the underlying inventive principle, it is essential that as much information available in the signal as possible be used while trying to determine a physical parameter from it. According to a further aspect of the invention the method may be used for an even smaller signal-to-noise ratio if the history of the amplitude of the received signal is taken into account. This embodiment of the method is accordingly characterised in that the digital signal processor comprises an amplitude memory, in which a filtered digital amplitude signal, as obtained from previous measurements, is stored and that every time a new filtered digital amplitude signal is determined from the digital amplitude signal and from the stored filtered digital amplitude signal. According to a further aspect of the invention the method may be used with an even smaller signal-to-noise ratio if the history of the phase of the received signal is taken into account. This embodiment is characterised in that the digital signal processor comprises a frequency memory, in which an estimated average frequency of the signal, obtained from previously obtained digital phase signals, is stored. Subsequently, the physical parameters are estimated using an estimation technique like least squares, maximum likelihood or momentum. Preferably, the weighted least square method is used to estimate a phase history from the estimated average frequency and from the digital phase signal. In fact, a straight line is fitted to the series of phases of which the 2ττ ambiguity is removed by using the previous estimated frequency, with weighting factors as determined by the filtered digital amplitude signals. The history is also used in the fitting procedure depending on the time variance of the mean frequency. By using a least square method a number, the so-called x 2 (chi square) is obtained, which is the residue of the least square procedure and which, apart from the system error, is a measure for the statistical error or noise. The physical parameter may subsequently be determined on the basis of a zero crossing in the estimated phase history that is closest to a maximum in the filtered digital amplitude signal and a new estimated frequency. Subsequently the error in the estimated mean frequency and the zero crossing error are derived quantities from the χ 2 error. A favourable embodiment of the inventive method, is characterised in that the digital signal processor comprises a crosstalk memory, in which a first complex digital signal is stored, representing crosstalk signals received by other receiving means and calculated using the characteristics of the specific crosstalk signals, and that the first complex digital signal is subtracted from a complex digital signal before the complex digital signal is further processed. It is also possible to subtract the crosstalk signal in the time-domain. This can be done using analog circuitry, i.e. before the analog-to-digital converter, or using the digital signals. The transducer according to the invention is capable of converting signals into electric signals. The signals preferably are acoustic signals, which are to be converted into electric signals. Acoustic signals suitable to be used in the invention can vary in a very broad range. It is possible to use 20 Hz signals, but also 1 2 kHz, up to 1 00 kHz and even MHz range signals can be applied. Preferably, the method of the current invention is used for acoustic signals relating to acoustic signals audible to the human ear, i.e. 20 Hz up to 20 KHz. Any signal may be used for the technique. In case of correlation of the receiving signals from two transducers located at dif- ferent positions the signal due to correlation is symmetric and can be used by this technique. The length of procedure is reciprocal with the bandwidth. The method and apparatus according to the present invention are specifically suited for application in the field of measurement of sound and specifically determination of sound levels, tracking of moving objects producing sound, like aeroplanes, cars, trains etcetera. The signal-to-noise ratio of these signals often can be very low, resulting in large errors in determined physical parameters, if these parameters can be determined altogether. In this field, the method and apparatus according to the current invention can, with considerable success, be used to detect moving objects, like aeroplanes. Even better results in the field of measurement of sound and tracking of moving sound sources in the air, like for example vehicles and aeroplanes, are obtained in the following way. In this case the transducer comprises at least two microphones ('phones') . If more than two phones are used, it is advantageous if the phones are not placed in one line if using three phones, and not placed in one plane when using four phones. It is thus advantageous if the phone are placed in such a way that three phones define a plane or four or more phones define a space. The distances between the phones are such that the angle sensitivity suffices, taking the turbulence of the medium into account. In case of two phones one directional angle is obtained, in case of three phones two angles are obtained with an ambiguity in the sign of the angle perpendicular to, or outside the plane defined by the three phones. In case of four phones, preferably the configuration is defined so that no phone is in the plane defined by the other phones. Preferably, the phones are placed on the vertices of a tetrahedron. Thus, the distances between the phones are the same, thus simplifying the mathematical procedures. By combining the signals from the different phones, the horizontal angle and the vertical angle are obtained. The phones convert the acoustic signals into electric signals, which, using cross correlation, are combined such, that complex correlation time signals are generated from the real signals. These correlated signals are averaged so that less frequency components are missing and thus the averaged correlation signals are more wide frequency banded, lowering the ambiguity in the correlation signals. In case of three phones, the number of correlation signals and averaged signals are three, and in case of four phones the number of correlation signals are six. Detection and tracking motion analysis is applied on these sources via the procedures described above in order to obtain accurate time differences relating to the directional information of these sources. The parameters of the motion analysis are used for the enhancement of the correlation signals. In case of four phones the ambiguity is fully resolved by combining the time differences in the averaged correlation signals. The acoustic signal from the moving or stationary sources is converted by the phones into electrical signals and further filtered/amplified by the preamplifier and converted by the analog to digital converter into series of numbers, called digital signals. By correlating these signals via a real Fast Fourier Transform (FFT) and transforming back with a complex FFT, the complex correlation signal is obtained, without using a Hubert filter, with that reducing the sample frequency by a factor 2. However, a Hubert filter can also be used. More generally, a digital filter is used to convolute the two signals and calculate the complex signal. The correlation signal is averaged, reducing the intrinsic ambiguity for noise signals. The procedure is performed for all two input signal combinations. By analysing the averaged correlation signal(s), using the noisiness or variation in the signals as a measure, one or more peaks are obtained which are related to the sound sources. By using the method described above, very accurate time differences are obtained for these peaks. Motion analysis predicts the places of the peaks in the correlation signals. By taking the motion into account, the source signal information is enhanced with regard to noise of other sources. The source information is analyzed by the algorithm described above and thus very accurate time differences are obtained for these source signals. By combining the time differences of the peaks in the correlation source signals, the three phones configuration gives both angles with a sign ambiguity for the vertical angle. For the four phones configuration the sign ambiguity is removed and the timing and angle errors are reduced. An even better embodiment of the procedure for the spatial configuration is to introduce the vertical elevation and the horizontal azimuth angles as variables in the least square method using all the correlated signals. By using the rate of change of the angles the relative motion or velocity vector can be calculated, which is a uniform motion for passing aeroplanes, for example. The relative motion is used to update the relative position and to calculate, using this relative position, both angles for the next measurement cycle. The variance of the measurements determines the gain or weighting of the history. Actually, this sequence of steps is well described by the term Kaiman Filtering. In this application of the Kalman Filtering, the phases of the correlated acoustic signals are used as input data, and the equation of uniform motion is used as the model. This type of filtering is optimal towards the accuracy of the relative position of the source and thus towards the sound measurement. Relative means in this context that two of the three parameters describing the motion relate to the third parameter. For example, the flight hight of an aeroplane is most likely to be taken as the third parameter. An even better embodiment is to omit the angles and calculate directly from the input data the relative position and motion of the source. To start tracking a source, peaks in the correlation signals have to be related to each other by a coincidence filter so that it corresponds with both angles, within an error depending on the variability of the environment. This relation is a function of time delay between the different signals and the strength of the signals. By combining two configurations on a large base, the range is obtained from these measurements with an error, depending on the atmospheric conditions and even more so on the distance between the source and the phone configurations and its phone base. As mentioned above, the method is specifically suited for signals having a large bandwidth. The noise of aeroplanes, for example, is a source of sound which is very well suited for analysis using the method and apparatus of the present invention. Using more apparatus as described, the distance of the sources can be determined. It is also possible to determine the position of a certain object in rest producing sound using the apparatus or method according to the present invention, mounted on a moving vehicle. If the position and speed of the vehicle are known, for example using GPS, Global positioning system, the exact location of a certain sound producing object can be determined. It is also possible to determine the exact position of the vehicle using a source having a known position. In case all data are related to position of the vehicle as parameter a technique can be applied more or less equivalent with the technique well known under the name synthetic aperture. Disturbing factors like the sound velocity and wind velocity can easily be estimated by modifying the time differences between the signals received by the different microphones due to the windvector, soundvelocity and position relative to the microphone configuration. The method and apparatus described above relating to acoustic signals can also be applied to eg. electromagnetic signals, preferably signals travelling through the air. The various aspects of the invention will now be further explained with reference to the following figures showing specific embodiments of the invention, in which: Fig. 1 shows a flow chart of the method and the apparatus according to the invention; Fig. 2A represents a possible wave shape of a received signal; Fig. 2B represents a frequency spectrum of this wave shape; Fig. 3 represents a block diagram of a possible embodiment of an apparatus for measuring the sound and location of objects; Fig. 4A presents a schematic diagram of the set-up for measuring sound properties using two phones; Fig. 4B presents a schematic diagram of the set-up for measuring sound properties using three phones; Fig. 4C presents a schematic diagram of the set-up for measuring sound properties using four phones; Fig. 5 shows a more detailed diagram according to elucidating component 60 of figures 4A, B and C; Fig. 6 presents a more detailed scheme of part 1 00 of figure 5; Fig. 7 presents a more detailed flow chart of part 1 20 and 140 of figure 5; Fig. 8 Presents a flow chart of an alternative embodiment of the data processor; Fig. 8A Presents a detailed flow chart of part 1 70 of figure 8; Fig. 8B Presents a detailed flow chart of part 1 71 of figure 8B; Fig. 9 presents an example of measurements of moving sound producing object; Fig. 10 presents an example of measurements of moving sound producing object. Fig. 1 shows a generalised, possible flow chart of the method and apparatus. A signal is detected by a first transducer 000 and converted into an electric signal. The electric signal is further processed by circuitry 1000. The circuitry 1 000 is connected to an analog-to-digital converter (ADC) 2000. The ADC 2000 converts the electric signal into a digital signal. If crosstalk or background noise removal circuitry is added, the signal is fed to a subtracter 3000 to subtract this crosstalk or background noise 360 in crosstalk memory 350 from the digital signal. The digital signal without crosstalk is subsequently fed to a digital filter 400, which converts the digital signal into a complex digital signal. This filter can be a combination of Fourier transforms or Fast Fourier Transforms (FFT) . Preferably, a Hubert filter is used. The Hubert filter converts a flow of single digital measurements into a flow of complex digital measurements with half the sample rate, in a manner known in the art. The complex digital signal is subsequently supplied to converter 500. The converter converts the complex digital signal to a digital amplitude signal and a digital phase signal. This conversion as such is well known in the art. If the complex digital signal is for example represented by a series: s (t, t n } then converter 500 produces from this the digital amplitude signal and the digital phase signal: A(t„ ... , t n ) and ø(t 1 f ... , t n ) The digital amplitude is stored in amplitude memory 560; the digital phase is stored in phase memory 570. The digital amplitude signal is supplied to a procedure for averaging the amplitude 600 in such a way that the history of the amplitude is taken into account, using for example a so-called Infinite Impulse Response (MR) filter like the analog equivalent, an RC filter, or other means known in the art. The amplitude signal and the history averaged amplitude signal and peak position of a previous peak are sup- plied to peak labelling procedure 700. In this peak labelling procedure a rough estimate of the peak position is estimated. The linearised phase data is unwrapped in order to remove the 2rt ambiguity resulting in datapoints on a continuous curve, and using a linear least squares method 800 a straight line is fitted through these phase data. In this least squares fit, the arrival time of the pulse is determined from the abscissa, the average frequency of the signal is calculated from the tangent, and the noise is determined from / 2 (chi-square) . This average frequency and arrival time of the detected pulse are calculated in 900. It is clear to a man skilled in the art that it might also be possible to perform a least squares fit, linear or non-linear, to an curve or function other than a straight line. It is preferred, however, to first convert the data to a smooth function, preferably a straight line, and then perform a least squares fit to fit the data points to a smooth function. Fig. 2A shows a possible wave shape of a received signal in the time domain. Usually a very short, unipolar electric signal is applied to the transmitting second transducer, comparable to a Dirac-pulse in response to which the transducer outputs a short, damped vibration. The frequency of the vibration is thereby determined by the resonance frequency of the transducer while the length of the signal is substantially determined by the internal damping of the transducer, coupling of the medium, and the frequency dependant damping caused by the medium in which the signal is transmitted, for example air, the object, in case the method of the invention is used for measurements above in the air. The receiving transducer is also set responding; for example vibrating upon the reception of one or more reflected signals, and then outputs an electric signal as shown. In Fig. 2A the signal-to-noise ratio is about 20 dB; consequently it is a strong signal. In an apparatus according to the invention, signals with a far smaller signal-to-noise ratio may be processed. In Fig. 2A the horizontal scale is in microseconds, the measured values have been sampled with a frequency of 400 Kc/s. The vertical scale only shows a ratio with respect to the maximum signal. Fig. 2B shows the spectrum of the received signal as shown in Fig. 2A in the frequency domain. The horizontal scale is in Kc/s, the vertical scale in dB. It clearly shows that the signal is a wide band signal. Fig. 3 represents a block diagram of a possible embodiment of apparatus according to the present invention. As shown, the transducer 31 0 is connected to receiver means 320, comprising a preamplifier and a bandpass filter. The receiver means 320 is connected to a computer 330 comprising an analog-to-digital converter 340, memory means 350 and software 360. The computer 330 is further connected to well known input, output and storage devices like for example a display 365. The digital signal is stored in the memory means 350 and led to the software 360 where it is processed using a Hubert filter subroutine. The Hubert filter converts a flow of single digital measurements into a flow of complex digital measurements with half the sample rate, in a manner known in the art. This data is also stored in the memory means 350. The software will be further elucidated. The digital signals, representing the received electric signals, are directed to a Hubert filter routine 370, supplied to a subtracter 380 in which, if desired, a crosstalk or background signal, is removed from the complex digital signal. The crosstalk can also be subtracted in the time domain, or before the Hilbert-filter step. Before the actual measurement of the physical parameter is started, crosstalk signals are determined and stored in a crosstalk memory 350. This is for example easily effected by measuring in a direction wherein no variation of the signal, due to the physical parameter to be measured, is present, it is also possible to determine the crosstalk signal by measuring a number of signals under different circumstances and by averaging the complex digital signals, in which averaging process a changing component in these complex digital signals disappears and the crosstalk signal is retained. Once the crosstalk signal is known, it may simply be subtracted from the complex digital signal in subtracter 380 for obtaining a substantially crosstalk free complex digital signal. The crosstalk-free complex digital signal thus obtained is subsequently supplied to a converter 390, arranged for converting this signal to a digital amplitude signal and a digital phase signal, which conversion as such is well known in the art and described above. The digital amplitude signal is supplied for further processing to an amplitude memory 405 in which a filtered digital amplitude signal, obtained from preceding measurements, is stored. The digital amplitude signal is supplied to an amplitude filter 41 0 together with a digital amplitude signal. Amplitude filter 410 determines a new filtered digital amplitude signal from a times the filtered digital amplitude signal and ( 1 -σ) times the digital amplitude signal, with a = 0.6-0.95. The new filtered amplitude signal and the maximum in the new filtered amplitude signal are then supplied to weighting means 420, which determines a weighting factor. The digital phase signal is processed by applying a weighted least square method 430. According to this method, a line is fitted through the successive unwrapped phases of the digital phase signal, in which process the new estimated average frequency determines the direction of the line and the new filtered amplitude signal supplies the weighting factors, because the phase information emanating from a small measured value is relatively unreliable. Finally a physical parameter within the signal is determined from the zero crossing of the estimated phase history that is nearest to the maximum of the new filtered amplitude signal. If the signal-to-noise ratio becomes very small, the weighting means 420 could incidentally point to a faulty zero crossing in the estimated phase history, causing a step-like change in the physical parameter as determined. In order to prevent this, a filter is added. The value of the parameter as measured is supplied to this filter. This value is then compared with values from other transducers. If there is a step-difference between these transducers, history of the parameters is taken into account to determine which values are most likely, and if this happens more then a specific number of times, then this type of filtering is not applied and a possible faulty condition is prevented. After an other specific number of times continued correspondence between parameter values, obtained from the data from all transducers, the filter is again used. If there is just one transducer present, stepwise changes are validated only by the amplitude history. This method and apparatus is preferably used with at least two transducers for receiving signals in order to track moving objects. It is possible to improve the signal-to-noise ratio by arranging transmitter means 31 0 such that a relatively long, coded pulse is transmitted by the transmitting transducer, and by providing receiver means 320 with a correlator for the coded pulse. The correlation peak, as produced by the correlator, has substantially the same shape as the pulse shown in Fig. 2A, which implicates that the processing in fact can be used unchanged for the main lob. Between adjacent lobs 1 80 degree phase-jumps occur. Since the information content in the side lobs is much less, the information may be discarded otherwise, the 1 80 degree phase-jumps have to be removed in the unwrapping process of the phase data. In figure 4A a set-up using two phones 50 is shown. The signals from the phones are directed to a data processor 60. Furthermore, parameters about the environment such as temperature of the air and humidity are fed from 70 into the data processor 60. Using two phones 50, information about the direction of the sources can be calculated. From these parameters other parameters such as the total amount of sound energy the source emits, etc. can be calculated. The distance 80 between the two phones 50 preferably is about 20 cm to 2 meter. The actual distance has to be accurately known. It is for example possible to determine the distance using a calibrated source. In figure 4B the same set-up as in figure 4A is presented when using three phones 50. Using three phones 50 in, for example, a triangular setting with more or less equal distances between the phones 50, information can be obtained about the direction, but no information disclosing at which side of the plane defined by the three phones 50 a source is located. In figure 4C a set-up similar to figures 4A and 4B is presented, using four phones 50. Using four phones 50 in, for example, a tetrahedron configuration with about the same distance 80 between each microphone, the direction of a source can be calculated. From this information, other parameters can be calculated. In figure 5 the data processor 60 of figures 4A-4C is elucidated, for four phones. The electric signals S enter the data processor 60 at circuitry 90 comprising a pre-amplifier and ADC. The four digital signals are input into six correlators 1 00 where each signal is correlated with each one of the other signals. Thus, four input signals result in 6 correlated signals. These correlated signals are each input a separate source detector/tracker circuitry 1 1 0. There, the arrival time and the time difference between each phone are calculated using source detector/tracker 1 1 0. Each of the six correlated signals enters averager 1 20, which is connected to the input of detector/tracker 1 30. The output of detector/tracker 1 30 is connected to tracking filter 140. The input of the source detector/tracker is also connected to the input of tracking filter 140. The output from the tracking filter 140 is looped back to the input of detector/tracker 1 30. The six signals output from the six source detector/tracker circuits are input into one coincidence filter 1 50. From the coincidence filter 1 50 several parameters 1 60 are obtained such as the sound level (dB) of the source, the sound level (dB) of the environment, the spectrum of the sound of the source, the spectrum of the sound of the environment and the direction of the source with respect to the detectors. In figure 6 the circuitry of the correlator 1 00 is elucidated. Two signals from different phones enter real Fast Fourier Transforming (FFT) means 102. The two Fourier transformed signals are subsequently input into a correlator 1 04 where the two signals are combined. The output side of the correlator 104 is connected to the input side of a FFT means where the complex inverse FFT of the input signal is calculated. Thus, the signal is converted into a complex signal. In figure 7 the averager 1 20 and the tracking filter 140 of figure 6 are further elucidated. The correlated signal from correlator 100 is input in averager 1 20. The input is connected to averager 1 21 and envelope averager 1 25. The averager 1 21 is connected to a converter, which converts the complex signal into an amplitude signal and a phase signal. The phase data and amplitude information is input in a phase unwrapper 1 23. The envelope averager 1 25 is further connected to peak detector 1 26 which is connected to phase unwrapper 123. The phase unwrapper 1 23 is connected to least squares estimator 1 24. The correlated signal from correlator 1 00 and the output from detector/tracker 1 30 are input to tracking filter 140. The input of tracking filter 140 is input in tracking filter 141 , the output of tracking filter 141 is input in converter 142 which converts the complex signal into an amplitude signal and a phase signal. The phase and amplitude data from converter 142 are input to phase unwrapper 143. The unwrapped phase data and amplitude history is input to least squares estimator 144 where time delay, frequency and strength are calculated. The result from least squares estimator 144 is input to source motion analyzer 145. The output from the source motion analyzer is also coupled back to tracking filter 141 and unwrapper 1 43. In figure 8, an alternative embodiment of the invention is shown, resulting in even better tracking of sound sources. The signals are conditioned and correlated as in figure 7. The source detector locates peaks in all correlation signals taking into account correlation peaks of already detected sources. Peaks are found by averaging the envelopes obtained from the complex correlation signal. From the averaged envelope signal the mean and variance is obtained. The peaks are found as those above a certain value defined by the mean and the variance of the envelope and the bandwidth of the correlation signal. The coincidence filter calculates the square error of possible azimuth and elevation angles (σ,e) j of source j. There are a limited possible directions for time delays r 3 ^. The correct direction and combination has the lowest square error. These angles are used for source tracking. The source tracker uses a steady state model for the motion of the plane. It predicts the angles of source j from the previous time and the history of the relative coordinates. From the complex correlation the envelope A a b k and the unwrapped phases ø a b k are calculated taking into account the time delays. These phases are used to calculate by a least square method the updates in angles and relative coordinates, motion of the source and the auxiliary phase rate parameter omega, which is equivalent with the mean frequency of the correlated radial frequency. Since the strength can vary 20 dB due to ground reaction, different weighting schemes can be applied to reduce these effects. By optimizing the filtering the Kalman filter on angles/relative positions is obtained using the phase signals. It should be noted that the Kalman filter is known per se from, for example, "Modern Control Theory", William T. Brogran, Prentice-Hall, 1 974. The sound velocity in air, a function of temperature and humidity, is supplied by a subsystem. The windvector, which manifests itself in the propagation of sound by changing the arrival times of the signals at the different transducers, can be determined by measuring directly or as auxiliary parameters of the least squares method. Since the sound velocity changes very slowly the sound velocity (scalar) may also be obtained via the least square method. A constraint is that the number of directions should be present in the data of a substantial time interval otherwise a singularity may occur due to the dependency of sound propagation vectors and wind velocity vectors. In case the errors are too large, the tracking is lost. In the source data extractor the source is parametrized according to the user specifications. For example, the angles are used to obtain the correct sound intensity of the source expressed in units like dB(A) or user defined units. Beamforming can be applied to obtain, from all phones used in the source detector and sound tracker, highly directional sound data. Due to the known coordinates of the source a limited number of calculations has to be made to obtain an acoustical time record of the source via the beamformer. The following components are depicted in figure 8. The electric signals S(t) enter the data processor 60 of figure 4A-4C at circuitry 90 comprising a pre-amplifier and ADC. The four digital signals are input into six correlators 1 00 where each signal is correlated with one of the other signals. The signals are convoluted through a real Fourier transform and transformed back using a complex inverse Fourier transform, resulting in a complex signal S x ' v (x = I, II, III, y = II, III, IV). Thus, four input signals result in 6 correlated complex signals. The correlator is already explained in figure 5. The correlated signals and the original signals are input in a data bus means or memory means 1 60. From the memory or data bus means 1 60 the signals are input in a source detector 1 70, a source tracker 1 90 and a source extractor 200. The source detector 1 70 is further elucidated in figure 8A and 8B. In the source detector 1 70 and source tracker 1 90 the arrival time and the time difference between each phone are calculated. In the source data extractor 200, finally, the physical parameters, in this case the relative position x and y relative to the height in case of aeroplanes, strength of the source in dB(A) and the wind vector in the horizontal plane in m/s are determined in the least square method in 1 90. In figure 8A, the source detector is further elucidated. The correlated complex signals S are entered into peak detector 1 71 . This peak detector is equivalent to the peak detector described in figure 7. From the peak detector 1 71 the time differences r and the mean amplitudes of the correlation signal at time differences τ are entered into the coincidence filter 1 80. The coincidence filter equals the coincidence filter described in figure 5. In figure 8B, the peak detector 1 71 is further explained. The signal S is directed to transformation means 1 72 for transforming the complex correlated signals into amplitude and phase signals. The signal S is also input into history means 1 75, where a predetermined amount of signal history is added to the current signal. From the transformation means 1 72 the amplitude and phase signals are input into a history filter 1 73, where a predetermined amount of history is added to the current signal. The signal is from history means 173 input into an average variance calculator 1 74. From average variance calculator 1 74 the signals are input into peak allocator 1 76. Into the peak allocator, also the signals with history from history means 175 is input. From the peak allocator, the signals with added history and the time differences between two microphones, r is output. In figure 9, a three dimensional plot is presented of measurements of two moving objects, i.e. two twin-engined jets. The y-axis represents the angle in the horizontal plane with regard to the measuring system. From the plot it is clear that the position of the jets can be followed very accurately. In figure 1 0, a plot is presented showing the intensity of the sound of the moving objects of figure 8 as function of time. The second source in the time interval 50-70 s is a Learjet. The plot shows that using the measuring system according to the invention for measuring sound properties, the sound intensity of moving objects as well as the direction of the objects with regard to the measuring system can be accurately determined. PvE

Claims: CLAIMS 1 . Method for extracting at least one physical parameter from at least one first signal, where the signal is received by at least two transducers which are located above a water surface or above land, the signal is converted into an electric signal, and the electric signal is applied to a signal processor, in which signal processor the electric signal is converted by an analog- to-digital converter into at least one first digital signal, characterised in that the digital signal is transformed into a complex digital signal, and physical parameters are determined on the basis of a phase history and an amplitude history of the complex digital signal. 2. Method according to claim 1 , characterised in that the digital signal is transformed into a complex digital signal by a Hubert filter or by a real Fourier transformation followed by a complex inverse Fourier transformation. 3. Method of claim 1 or 2, characterised in that the signal processor transforms the complex digital signal into a digital amplitude signal and a corresponding digital phase signal. 4. Method according to any of the preceding claims, characterised in that the digital signal processor comprises an amplitude memory, in which a filtered digital amplitude signal as obtained from previous measurements is stored, and that every time a new filtered digital amplitude signal is determined from the digital amplitude signal and from the stored filtered digital amplitude signal. 5. Method according to claim 4, characterised in that the digital signal processor comprises a frequency memory, in which an estimated average frequency of the signal, obtained from previously obtained digital phase signals, is stored. 6. Method according to claim 5, characterised in that using the weighted least squares method a phase history is estimated from the estimated frequency and from the digital phase signal. 7. Method according to claim 6, characterised in that the filtered digital amplitude signal supplies the weighting factors. 8. Method according to claim 7, characterised in that the physical parameter is determined based on a zero crossing in the estimated phase history that is closest to a maximum in the filtered digital amplitude signal. 9. Method according to any of the preceding claims, characterised in that a crosstalk signal not comprising a physical parameter to be measured is subtracted from the signal, digital signal or complex digital signal before the complex digital signal is further processed. 10. Method according to any of the preceding claims, characterised in that signals from the at least two different transducers are correlated in such a way that the signal of each transducer is correlated with each one of the other signals, resulting in correlation of all different combinations of two signals. 1 1 . Method according to any of the preceding claims, characterised in that the history is taken into account through the physical model used. 1 2. Method according to any of the preceding claims, characterised in that the proportion of the variance of a current measurement and the variance of previous measurements determine the history to be taken into account. 1 3. Method according to any of the preceding claims, characterised in that the measured data is unwrapped in order to remove discontinuities. 1 4. Method according to any of the preceding claims, characterised in that all correlated signals from all the transducers are input in a single optimisation procedure, preferably a least squares procedure. 1 5. Method according to any of the preceding claims, characterised in that at least the relative position of a source and the wind velocity vector are determined. 1 6. Method according to any of the preceding claims, characterised in that physical parameters of moving objects located above a water surface or above land are measured. 1 7. Use of the method according to any of the preceding claims 1 -1 6 for determining the relative speed, relative trajectory and strength of sound of moving objects above a water surface or above land. 1 8. Apparatus for extracting at least one physical parameter from an signal, comprising at least two receiver means located above a water surface or above land for converting the signal into an electric signal, a signal processor and a digital signal processor, connected to the receiver means and provided with an analog-to-digital converter, characterised in that the digital signal processor is provided with a means for deriving a first complex digital signal, a converter means connected to the means for deriving a first complex digital signal for converting the complex digital signal into an amplitude signal and a phase signal, an amplitude memory connected to the converter means for storing the amplitude signal and a phase memory means for storing the phase signal. 1 9. Apparatus according to claim 1 8, characterised in that means for deriving a first complex digital signal comprises a Hubert filter or a filter comprising a real Fourier transform and a complex inverse Fourier transform. 20. Apparatus according to claim 1 8 or 1 9, characterised in that the digital signal processor is arranged for converting a first complex digital signal into a digital amplitude signal and a corresponding digital phase signal. 21 . Apparatus according to claim 20, characterised in that the digital signal processor is provided with an amplitude memory, for storing a filtered digital amplitude signal, derived from previous measurements, and a frequency memory, for storing an estimated average frequency, derived from digital phase signals from previous measurements. 22. Apparatus according to claim 21 , characterised in that the digital signal processor is provided with filter means for deriving a new filtered digital amplitude signal from a digital amplitude signal and a filtered digital amplitude signal stored in the amplitude memory and for deriving a new average frequency of the received signal from a digital phase signal and from an average frequency stored in the frequency memory. 23. Apparatus according to claim 22, characterised in that the digital signal processor is provided with weighting means, for estimating a phase history using a weighted least square method, from the estimated frequency of the signal and the series of phase estimations, where the filtered digital amplitude signal provides the weighting factors. 24. Apparatus according to any of the preceding claims 1 8-23, characterised in that the digital signal processor is arranged for deriving physical parameters from the signal from the phase history. 25. Apparatus according to any of the preceding claims 1 8-24, charac- terised in that the digital signal processor is moreover provided with a crosstalk memory for storing a crosstalk representing second complex digital signal, derived from received complex digital signals, and with a subtracter, for subtracting the second complex digital signal from a first complex digital signal before this first complex digital signal is further processed. 26. Vehicle with means for propagation or movement equipped for measuring physical parameters of objects moving above a water surface or above land comprising an apparatus according to any one of the preceding claims 1 8-25. 27. Data carrier loaded with software for performing the method or used in the apparatus according to any one the preceding claims, wherein the software comprises the steps of - converting digital signals into complex digital signals, - converting the complex digital signals into digital phase signals and digital amplitude signals, - storing the digital phase signals into phase memory means, - storing the digital amplitude signals into amplitude memory means, - finding a first estimate of the maximum amplitude of the current amplitude signal using previous amplitude data, - linearising the digital phase data as function of the associated time base using the first estimate of the maximum amplitude, - calculating the tangent and abscissa of the phase data as a function of the associated time basis. - calculating a physical parameter using the tangent, the abscissa and a previously calculated physical parameter. PvE

Inventors: ROOSNEK NICO, ROOSNEK, NICO
Applicants: ROOSNEK NICO, ROOSNEK, NICO
Application year: 2000
Main classification: G01S11
IPC codes: G01V1