A METHOD AND A SYSTEM FOR DECOMPOSITION OF ACOUSTIC SIGNAL INTO SOUND OBJECTS, A SOUND OBJECT AND ITS USE
A method and a system for decomposition of acoustic signal into sound objects, a sound object and its use FIELD OF INVENTION [0001] The object of the invention is a method and a system for decomposition of acoustic signal into sound objects having the form of signals with slowly varying amplitude and frequency, and sound objects and their use. The invention is applicable in the field of analysis and synthesis of acoustic signals, e.g. in particular to speech signal synthesis. STATE OF THE ART [0002] For a dozen of years the progress in analysis of sound signals has been inconsiderable. Still well-known methods are used such as neural networks, wavelet analysis or fuzzy logic. Beside these methods rather widespread is the use of the classic Fast Fourier Transform (FFT) algorithm for signal filtering, which allows to analyse the frequency of components with the use of relatively low computational power . [0003] One of the most difficult areas, but also the one being of the greatest interest within the analysis of sound signals is the analysis and synthesis of speech. [0004] Despite huge progress being observed in the development of digital technology, the progress in sound signal processing systems in this field is not significant. During the last few years, multiple applications have appeared, which attempt to fill the niche related to recognition of speech, but their common origin (mainly the analysis in the frequency domain with the use of Fourier transform) and the limitations related to it cause that they do not respond to the demand of the market. [0005] The main drawbacks of these systems are: 1) Vulnerability to external interference The existing sound analysis systems operate satisfactorily in conditions ensuring one source of the signal. If additional sources of sound appear, such as interference, ambient sounds or consonant sounds of multiple instruments, their spectrum overlap, causing the mathematical models being applied to fail. 2) Relative variation of spectral parameters Methods for calculating a sound signal ' s parameters which are currently used originate in the Fourier transformation. It
assumes a linear variation of analysed frequencies, meaning that a relative variation of two adjacent frequencies is not constant. For example, if a window of 1024 (210) data of a signal sampled with the rate of 44100 samples per second (SPS) is analysed with the use of the FFT algorithm, then the subsequent frequencies of the spectrum differ by 43.07 Hz. The first non-zero frequency is Fl = 43.07 Hz, the next one F2 = 86.13 Hz. The last frequencies are F510 = 21963.9 Hz, F511 = 22 006,9 Hz. At the beginning of the range a relative variation of the spectral frequency is 100% and leaves no opportunity to identify the sounds being closer. At the end of the range a relative variation of the spectral parameter is 0.0019% and is undetectable for human ear. 3) Limitation of parameters to spectral amplitude characteristics The algorithms based on the Fourier transformation use the amplitude characteristic for the analysis, and in particular the maximum of the amplitude of the spectrum. In the case of sounds with different frequencies close to each other this parameter will be strongly distorted. In this case, additional information could be obtained from the phase characteristic, analysing the signal's phase. However, since the spectrum is analysed in windows shifted e.g. by 256 samples, there is nothing to relate the calculated phase to. [0006] This problem has been partly solved by the speech information extraction system, described in patent US5214708. Disclosed therein is a bank of filters having central frequencies logarithmically spaced in relation to each other, according to the model of human ear perception. Due to an assumption that within a band of any of these filter banks there is only one tone, the problem of the uncertainty principle in the field of signal processing has been partially evaded. According to the solution disclosed in US5214708, information about modulation on each of harmonics, including frequency and time-domain waveform information, can be extracted based on the measurement of the logarithm of each harmonics' power. Logarithms of signal's amplitude in adjacent filters are obtained with the use of Gaussian filters and logarithmic amplifiers. However, the drawback of this solution is that the function FM(t) used for speech analysis does not effectively extract essential characteristic parameters of a single speech signal. The next much more significant drawback of this solution is an assumption that the audio signal comprises a signal from only one source, such simplification reducing significantly practical possibilities of using such system for decomposition . [0007] Therefore, an object of the invention is to provide a method and a system for decomposition of acoustic signal,
which would make possible an effective analysis of acoustic signal perceived as a signal incoming simultaneously from a number of sources, while maintaining a very good resolution in time and frequency. More generally, an object of the invention is to improve the reliability and to enhance the possibilities of sound signals' processing systems, including those for analysis and synthesis of speech. SUMMARY OF THE INVENTION [0008] The essence of the invention is that a method for decomposition of acoustic signal into sound objects having the form of sinusoidal wave with slowly-varying amplitude and frequency, comprising conversion of the signal from analogue to digital form, then splitting the signal into adjacent sub- bands with central frequencies distributed according to logarithmic scale by feeding samples of the acoustic signal to the digital filter bank's input, is characterised in that - at each filter's output for each sample of the input signal PIN the real value FC (n) and the imaginary value FS (n) of the filtered signal is determined, and then based on this - for each sample of the input signal PIN the frequency, the amplitude and the phase of all detected constituent elements of said acoustic signal are determined, for each sample of the input signal PIN an operation improving the frequency-domain resolution of said filtered signal is executed - for each detected element of said acoustic signal an active object in an active objects database is created for its tracking - for each sample of the input signal PIN subsequent detected elements of said acoustic signal are associated with at least selected active objects in said active objects database to create a new active object or to append said detected element to an active object, or to close an active object for each active object in the database values of the envelope of amplitude and values of frequency and their corresponding time instants are determined not less frequently than once per period of duration of a given filter's window W(n) so as to create characteristic points describing slowly-varying sinusoidal waveform of said sound obj ect for each sample of the input signal PIN at least one selected closed active object is transferred to a database of sound objects to obtain at least one decomposed sound object, defined by a set of characteristic points with coordinates in time-frequency-amplitude space. [0009] The essence of the invention is also that a system for decomposition of acoustic signal into sound objects in the form of sinusoidal waveforms with slowly-varying amplitude and frequency, comprising a converter system for digital-to- analog conversion from which said signal being input to a
filter bank having central frequencies of filters distributed according to geometrical distribution is characterised in that each filter is adapted to determine a real value FC (n) and an imaginary FS (n) of said filtered signal, said filter bank being connected to a system for tracking objects, wherein the system for tracking objects comprises a spectral analysis system, a voting system adapted to calculate the value of angular frequency weighted by amplitude, a resolution improving system adapted to adjust the filter's window and/or to determine a difference between spectrum between the input signal and at the filter bank's output, and/or to determine a difference between the input signal and the bank filter's output signal, a system for associating objects, a shape forming system adapted to determine characteristic points describing slowly-varying sinusoidal waveforms, a database of active objects (34) and a database of sound objects (35) . [0010] The essence of the invention is also that a sound object being a signal having slowly-varying amplitude and frequency is characterised in that it is obtained by the method according to any of claims 1 to 7. [0011] Furthermore, the essence of the invention is also that a sound object being a signal having slowly-varying amplitude and frequency is characterised in that it is defined by characteristic points having three coordinates in the time- amplitude-frequency space, wherein each characteristic point is distant from the next one by a value not greater than 20 periods of the object's frequency. [0012] The main advantage of the method and the system for decomposition of signal according to the invention is that it is suitable for effective analysis of a real acoustic signal, which usually is composed of signals incoming from a few different sources, e.g. a number of various instruments or a number of talking or singing persons. [0013] The method and the system according to the invention allow to decompose a sound signal into sinusoidal components having slow variation of amplitude and frequency of the components. Such process can be referred to as a vectorization of a sound signal, wherein vectors calculated as a result of the vectorization process can be referred to as sound objects. In the method and the system according to the invention a primary objective of decomposition is to extract at first all the signal's components (sound objects), next to group them according to a determined criterion, and afterwards to determine the information contained therein. [0014] In the method and the system according to the invention a signal is analysed both in the time domain and in the frequency domain sample by sample. Of course this
increases the demand for computational power. As already mentioned, the technologies applied so far, including the Fourier transformation with its implementation as the fast transform FFT and SFT, have played a very important role in the past when the computational power of computers was not high. However, during the last 20 years the computational power of computers has increased 100000 times. Therefore, the invention reaches for tools which are more laborious, but which offer improved accuracy and are better suited to the human hearing model. [0015] Due to the use of a filter bank having a very large number of filters (over 300 for the audible band) with logarithmically spaced central frequencies, and due to applied operations increasing the frequency-domain resolution, one obtains a system capable of extracting two simultaneous sources of sound separated from each other even by half a tone. [0016] A spectrum of the audio signal obtained at said filter bank's output comprises information about the current location and variations in the sound objects' signal. The task of the system and the method according to the invention is to precisely associate a variation of these parameters with existing objects, to create a new object, if the parameters do not fit to any of the existing objects, or to terminate an object if there are no further parameters for it. [0017] In order to precisely determine the parameters of an audio signal, which are intended to be associated with existing sound objects, the number of considered filters is increased and a voting system is used, allowing to more precisely localize frequencies of the present sounds. If close frequencies appear, the length of said filters is increased for example to improve the frequency-domain resolution or techniques for suppressing the already recognized sounds are applied so as to better extract newly appearing sound objects. [0018] The key point is that the method and the system according to the invention track objects having a frequency variable in time. This means that the system will analyse real phenomena, correctly identifying an object with a new frequency as an already existing object or an object belonging to the same group associated with the same source of signal. Precise localization of the objects' parameters in amplitude and frequency domain allows to group objects in order to identify their source. Assignment to a given group of objects is possible due to the use of specific relations between the fundamental frequency and its harmonics, determining the timbre of the sound.
A precise separation of objects makes a chance of further analysis for each group of objects, without interference, by means of already existing systems, which obtain good results for a clean signal (without interference) . Possessing precise information about sound objects which are present in a signal makes it possible to use them in completely new applications such as, for example, automatic generation of musical notation of individual instruments from an audio signal or voice control of devices even with high ambient interference. BRIEF DESCRIPTION OF DRAWINGS [0020] The invention has been depicted in an embodiment with reference to the drawings, wherein: [0021] FIG.l is a block diagram of a system for decomposition of audio signal into sound objects, [0022] FIG.2a is a parallel structure of a filter bank according to the first embodiment of the invention, [0023] FIG.2b is a tree structure of the filter bank according to the second embodiment of the invention, [0024] FIG.3 shows a general principle of operation of a passive filter bank system, [0025] FIG.4 shows exemplary parameters of filters, [0026] FIG.5 is the impulse response of a filter F(n) having the Blackman window, [0027] FIG.6 is a flowchart of a single filter, [0028] FIGS .7a and 7c show a part of a spectrum of the filter bank output signal, comprising the real component FC(n), the imaginary component FS (n) and the resulting amplitude of the spectrum FA(n) and the phase FF(n) [0029] FIGS .7b and 7d show the nominal angular frequency F# (n) of a corresponding filter group and angular frequency of the spectrum FQ(n) . [0030] FIG .8 is a block diagram of a system for tracking sound objects, [0031] FIGS .9a and 9b show exemplary results of operation of a voting system, [0032] FIG.10 is a flowchart of a sound system for associating objects,
FIG.11 shows the operation of a frequency resolution improvement system according to an embodiment, [0034] FIG.12 shows the operation of a frequency resolution improvement system according to another embodiment, [0035] FIG.13 shows the operation of a frequency resolution improvement system according to yet another embodiment, [0036] FIGS.14a, 14b, 14c, 14d show examples of representation of sound objects, [0037] FIG.15 shows an exemplary format of notation of information about sound objects, [0038] FIG.16 shows a first example of a sound object requiring correction, [0039] FIG.17 shows a second example of a sound object requiring correction, [0040] FIG.18 shows further examples of sound objects requiring correction; [0041] FIGS.19a, 19b, 19c, 19d, 19e, 19f, 19g, 19h show the process of extracting sound objects from an audio signal and synthesis of an audio signal from sound objects. DETAILED DESCRIPTION OF EMBODIMENTS [0042] In the present patent application the term "connected", in the context of a connection between any two systems, should be understood in the broadest possible sense as any possible single or multipath, as well as direct or indirect physical or operational connection. [0043] A system 1 for decomposition of acoustic signal into sound objects according to the invention is shown schematically in FIG.l. An audio signal in digital form is fed to its input. A digital form of said audio signal is obtained as a result of the application of typical and known A/D conversion techniques. The elements used to convert the acoustic signal from analogue to digital form have not been shown herein. The system 1 comprises a filter bank 2 with an output connected to a system for tracking objects 3, which is further connected with a correcting system 4. Between the system for tracking objects 3 and the filter bank there exists a feedback connection, used to control the parameters of the filter bank 2. Furthermore, the system for tracking objects 3 is connected to the input of the filter bank 2 via a differential system 5, which is an integral component of a frequency resolution improvement system 36 in FIG.8.
In order to extract sound objects from an acoustic signal, a time-domain and frequency-domain signal analysis has been used. Said digital input signal is input to the filter bank 2 sample by sample. Preferably, said filters are SOI filters. It is shown in FIG.2a a typical structure of the filter bank 2, in which individual filters 20 process in parallel the same signal with a given sampling rate. Typically, the sampling rate is at least two times higher than the highest expected audio signal's component, preferably 44.1 kHz. Since such a number of samples to be processed per 1 second requires large computational expense, preferably a filter bank tree structure of FIG.2b can be used. In the filter bank tree structure 2 the filters 20 are grouped according to the input signal sampling rate. For example, the splitting in the tree structure can be done at first for the whole octaves, For individual sub-bands with lower frequencies it is possible to cut off high frequency components using a low-pass filter and to sample them with a smaller rate. As a consequence, due to reduction of the number of samples a significant increase in processing speed is achieved. Preferably, for the interval up to 300 Hz the signal is sampled with fp=600 Hz, up to 2.5 kHz with fp=5 kHz. [0045] Since the main task of the method and the system according to the invention is to localize all sound objects in the spectrum, an important issue is possible accuracy of determination of signal's parameters and a resolution of simultaneously appearing sounds. The filter bank should provide a high frequency-domain resolution, i.e. greater than 2 filters per semitone, making it possible to separate two adjacent semitones. In the presented examples 4 filters per semitone are used. [0046] Preferably, in the method and the system according to the invention a scale corresponding to human ear's parameters has been adopted, with logarithmic distribution, however a person skilled in the art will know that other distributions of filters' central frequencies are allowed within the scope of the invention. Preferably, a pattern for the distribution of filters' central frequencies is the musical scale, wherein the subsequent octaves begin with a tone 2 times higher than the previous octave. Each octave is divided into 12 semitones, i.e. the frequency of two adjacent semitones differs by 5.94% (e.g. el=329.62 Hz, fl=349.20 Hz). To increase accuracy, there are four filters for each semitone in the method and the system according to the invention, wherein each filter listens to its own frequency, differing from an adjacent frequency by 1,45%. It has been assumed that the lowest audible frequency is C2=16,35 Hz. Preferably, the number of filters is greater than 300. A particular number of filters for a given embodiment depends on the sampling rate. With sampling at 22050 samples per second the highest
frequency is e6 = 10548 Hz, 450 filters being in this range. With sampling at 44100 samples per second the highest frequency is e7 = 21096 Hz, 498 filters being in this range. [0047] A general principle of operation of a passive filter bank is shown in FIG.3. The input signal which is fed to each filter 20 of the filter bank 2 is transformed as a result of relevant mathematical operations from the time domain into the frequency domain. In practice, a response to an excitation signal appears at the output of each filter 20, and the signal's spectrum jointly appears at the filter bank's output. [0048] FIG.4 shows exemplary parameters of selected filters 20 in the filter bank 2. As can be seen in the table, central frequencies correspond to tones to which a particular music note symbol can be attributed. The window width of each filter 20 is given by the relation: W(n) = K * fp / FN(n) (1) where: W(n) - window width of a filter n fp - sampling rate ( e.g. 44100 Hz ) FN(n)- nominal (central) frequency of a filter n K - window width coefficient ( e.g. 16 ) Since a higher frequency-domain resolution is necessary in the lower range of the musical scale, therefore for this range of frequencies the filter windows will be the widest. Thanks to an introduction of coefficient K and a normalization to the filter nominal frequency FN there is provided an identical amplitude and phase characteristic for all the filters. [0049] With regard to the implementation of said filter bank - a skilled person will know that one of possible ways of obtaining the coefficients of a SOI type band-pass filter is to determine the impulse response of the filter. An exemplary impulse response of a filter 20 according to the invention is shown in FIG.5. An impulse response in FIG. 5 is the impulse response of a filter with a cosine window, which is defined by the relation: y(i) (n) = cos(o)(n)*i) * (A-B*cos (2ni/W (n) ) +C*cos ( 4ni/W (n) ) (2) where: ω(η) = 2n*FN(n)/fp W(n), FN(n), fp - are defined above Window type A B C Hann (Hanning) jo.5 10.5
Window type Hamming 0.53836 0.46164 Blackman 10.42 0.5 ,0.08 The operations performed by each of the filters 20 have been shown in FIG.6. The task of the filter bank 2 is to enable the determination of an audio signal's frequency spectrum in the range of frequencies from the lowest audible by human (e.g. C2=16.35 Hz) to ½ fp - sampling rate ( e.g. e7=21096 Hz at 44100 samples per second) . Before each filter begins its operation, parameters of the filter 20 are initiated, the exemplary parameters being the coefficients of particular components of time window function. Then, the current sample PIN of the input signal, having only a real value, is fed to the input of the filter bank 2. Each filter 2, using a recursive algorithm, calculates a new value of components FC(n) and FS (n) based on the previous values of the real component FC(n) and the imaginary component FS (n) , and calculates also values of the sample PIN input to the filter and the sample Ρ0οτ leaving the filter's window and which is stored in an internal shift register. Thanks to the use of a recursive algorithm the number of calculations for each of the filters is constant and does not depend on the filter's window length. The executed operations for a cosine window are defined by the formula: By using trigonometric equations relating to products of trigonometric functions for equations (3) and (4) one obtains a dependence of the components FC(n) and FS (n) on the values of these components for the previous sample of the audio signal and a value of the sample inputted to the filter P™, and the one outputted from the filter P0DT according to the equation shown in FIG.6. In the case of each filter 20 the calculation of the equation for each subsequent sample requires 15 multiplications and 17 additions for Hann or Hamming type windows, or 25 multiplications and 24 additions for a Blackman window. The process of the filter 20 is finished when there are no more audio signal samples at the filter ' s input . [0050] Values of the real component FC(n) and the imaginary component FS (n) of the sample obtained after each subsequent sample of the input signal are forwarded from each filter's 20 output to a system for tracking sound objects 3, and in particular to a spectrum analysing system 31 comprised therein (as shown in FIG.8). Because the spectrum of the
filter bank 2 is calculated after each sample of the input signal, the spectrum analysing system 31 except of the amplitude characteristic can utilize the phase characteristic at the filter bank's 2 output. In particular, in the method and the system according to the invention the change of phase of the current sample of the output signal in relation to the phase for the previous sample is used for precise separation of the frequencies present in the spectrum, what will be described further with reference to FIGS. 7a, 7b, 7c and 7d, and FIG.8. [0051] A spectrum analysing system 31, being a component of the system for tracking objects 3 (as shown in FIG.8) calculates individual components of the signal's spectrum at the filter bank output. To illustrate the operation of this system, an acoustic signal with the following components has been subjected to analysis: There are shown in FIGS. 7a and 7b plots of instantaneous values of quantities obtained at the output of selected group of filters 20 for said signal and values of quantities calculated and analysed by the spectrum analysing system 31. For filters with number n from 266 to 336 with a window having the window width coefficient K=16 there have been represented: the instantaneous value of the real component FC[n], the instantaneous value of the imaginary component FS , which are fed to the input of the spectrum analysis system 31, and the instantaneous value of the spectrum's amplitude FA[n] and the spectrum's phase FF[n], which are calculated by the spectrum analysing system 31. As already mentioned, the spectrum analysing system 31 collects all the possible information necessary to determine the actual frequency of the sound objects present at a given time instant in the signal, including the information about the angular frequency. The correct location of the tone of component frequencies has been shown in FIG.7b, and it is at the intersection of the nominal angular frequency of the filters FQ[n] and the value of the angular frequency at the output of the filters FQ[n], calculated as a derivative of the phase of the spectrum at the output of a particular filter n. Thus, according to the invention, in order to detect a sound object, the spectrum analysing system 31 analyses also the plot of angular frequency F# and FQ[n]. In the case of a signal comprising components which are distant from each other, points which are determined as a
result of analysis of the angular frequency correspond to locations of maxima of the amplitude in FIG. 7a. [0053] Due to some typical phenomena in the signal processing domain basing only on maxima of amplitude of the spectrum is not effective. The presence of a given tone in the input signal affects the value of the amplitude spectrum at adjacent frequencies, leading in consequence to a severely distorted spectrum when the signal comprises two tones close to each other. To illustrate this phenomenon, and to illustrate the functionality of the spectrum analysing system 31 according to the invention, a signal has been subjected also to the analysis, comprising sounds of frequencies: As shown in FIGS. 7c and 7d, in the case of a signal with closely located components, the correct location of a tone determined based on the analysis of angular frequency plots does not correspond to the maximum of amplitude in FIG.7c. Thus, for such a case, thanks to various parameters analysed by the spectrum analysing system 31 it is possible to detect situations which are critical for decomposition of an acoustic signal. In consequence, it is possible to apply specific procedures leading to correct recognition of components, what will be described further with reference to FIG. 8 and FIG.9a, and FIG. b. [0054] The fundamental task of the system for tracking objects 3, a block diagram of which is shown in FIG.8, is to detect at a given time instant all frequency components present in an input signal. As shown in FIG.7b and FIG.7d, the filters adjacent to the input tone have very similar angular frequencies, different from the nominal angular frequencies of those filters. This property is used by another subsystem of the system for tracking objects 3, namely the voting system 32. To prevent incorrect detection of frequency components, the values of the amplitude spectrum FA(n) and angular frequency at the output of filters FQ(n), calculated by the spectrum analysing system 31, are forwarded to the voting system 32 for calculation of their weighted value and detection of its maxima in function of the filter's number (n) . In this way, one obtains a voting system, which takes into account the frequency at the outputs of all the filters 20 adjacent to it in order to determine frequencies present in the input signal for a given frequency at the filter's 2 output. The operation of this, system is shown in
FIGS.9a and 9b. FIG.9a illustrates a relevant case shown in FIGS.7a and 7b, while FIG.9b illustrates a relevant case shown in FIGS.7c and 7d. As it can be seen, the plot of the signal FG(n) (the weighted value calculated by the voting system 32) has distinct peaks in locations corresponding to tones of frequency components present in the input signal. In the case of an input signal comprising components distinctly separated from each other (as shown in FIG.9a) these locations correspond to a maximum of amplitude of the spectrum FA(n) . In the case of a signal comprising components situated too close to each other (as shown in FIG.9b), without the voting system 32 tones reflected in maximum of amplitude of the spectrum would have been detected, which are located in places other than the mentioned peaks in the weighted signal FG(n). [0055] As is shown in FIG.8, the spectrum analysing system 31 and the voting system 32 are connected at their output with a system for associating objects 33. Having at its disposal the list of frequencies detected by the voting system 32 which composes the input signal, and additional parameters, such as amplitude, phase and angular frequency associated to each detected frequency, the system for associating objects 33 combines these parameters in "elements" and next builds sound objects out of them. Preferably, in the system and the method according to the invention, the frequencies (angular frequencies) detected by the voting system 32, and thus "elements", are identified by the filter number n. The system for associating objects 33 is connected to an active objects database 34. The active objects database 34 comprises objects arranged in order depending on the frequency value, wherein the objects have not yet been "terminated". The term "a terminated object" is to be understood as an object such that at a given time instant no element detected by the spectrum analysing system 31 and the voting system 32 can be associated with it. The operation of the system for associating objects 33 has been shown in FIG.10. Subsequent elements of the input signal detected by the voting system 32 are associated with selected active objects in the database 34. To limit the number of required operations, preferably, detected objects of a given frequency are compared only with the corresponding active objects located in a predefined frequency range. At first, the comparison takes into account the angular frequency of an element and an active object. If there is no object sufficiently close to said element (e.g. in the range of distances in frequency corresponding to 0.2 tone) this means that a new object has appeared and it should be added to the active objects 34. If, once associating objects with current elements has been finished, there is no element sufficiently close for an active sound object (e.g. in the range of distances in frequency corresponding to 0.2 tone) this means that no further parameters for the object are detected and it should be terminated. Said terminated
object is taken into account in the association process still for 1 period of its frequency to avoid an accidental termination caused by a temporary interference. During this time it can return to active sound objects in the database 34. After 1 period the object's final point is determined. If the object lasted for a sufficiently long time (e.g. its length was not shorter than the width of the corresponding window W[n]), then this object is transferred to a sound objects database 35. [0056] In the case of associating with each other an active object and an object sufficiently close to, a matching function is further calculated in the system for associating objects 33, which comprises the following weighted values: amplitude matching, phase matching, objects duration time. Such a functionality of the system for associating objects 33 according to the invention is of essential importance in the situation when in a real input signal a component signal from one and the same source has changed frequency. This is because it happens that as a result of frequency changing a number of active objects become closer to each other. Therefore, after calculating the matching function the system for associating objects 33 checks if at a given time instant there is a second object sufficiently close to in the database 34. The system 33 decides which object will be a continuer of the objects which join together. The selection is decided by the result of the matching function comparison. The best matched active object will be continued, and an instruction to terminate will be issued for the remaining ones. Also a resolution improvement system 36 cooperates with the active objects database 34. It tracks the mutual frequency-domain distance of the objects present in the signal. If too close frequencies of active objects are detected the resolution improvement system 36 sends a control signal to start one of the three processes improving the frequency-domain resolution. As mentioned previously, in the case of presence of a few frequencies close to each other, their spectrum overlap. To distinguish them the system has to "listen intently" to the sound. It can achieve this by elongating the window in which the filter samples the signal. In this situation a window adjustment signal 301 is activated, informing the filter bank 2 that in the given range the windows should be elongated. Due to the window elongation the signal dynamics analysis is impeded, therefore if no close objects are detected the resolution improvement system 36 enforces a next shortening of the filter's 20 window. In the solution according to the invention a window with length of 12 to 24 periods of nominal frequency of the filter 20 is assumed. The relation of the frequency-domain resolution with the window's width is shown in FIG.11. The table below illustrates the ability of the system to detect and track at least 4 non-damaged objects subsequently present next to each other, with the minimal distance expressed in
percentage, as a function of the window's width. Window width Detects objects in Tracks objects in th periods) the distance of distance of 12 17.4% 23.2% 16 14.5% 17.4% 20 8.7 14.5% 24 5.9% 11.6% [0057] In another embodiment the system "listens intently" to a sound by modifying the filter bank's spectrum, what is schematically illustrated in FIG.12. The frequency-domain resolution is improved by subtracting from a spectrum at the tracking system's 3 input the expected spectrum of "well localised objects", which are localised in vicinity of new appearing objects. "Well localised objects" are considered as objects the amplitude of which does not vary too quickly (no more than one extreme per window's width) and the frequency of which does not drift too quickly (no more than 10% variation of frequency per window's width) . An attempt to subtract a spectrum of objects varying quicker can lead to the phase inversion at the measurement system input and to a positive feedback resulting in generation of an interfering signal. In practice the resolution improvement system 36 calculates the expected spectrum 303 based on the known instantaneous frequency, amplitude and phase of an object and subtracts them from the real spectrum, causing that the spectrum of adjacent elements will not be interfered so strongly. The spectrum analysing system 31 and the voting system 32 perceive only adjacent elements and a variation of the subtracted object. However, the system for associating objects 33 further takes into account the subtracted parameters while comparing the detected elements with the active objects database 34. Unfortunately, to implement this frequency-domain resolution improvement method a very large number of computations is required and a risk of positive feedback exists. [0058] In a yet another embodiment, the frequency-domain resolution can be improved by subtracting from the input signal an audio signal generated based on well localised (like in the previous embodiment) adjacent objects. Such operation is shown schematically in FIG.13. In practice, this relies on the fact that the resolution improvement system 36 generates an audio signal 302 based on information about frequency, amplitude and phase of the active objects 34, which is forwarded to a differential system 5 at the filter bank's 2 input, as shown schematically in FIG.13. The number of required calculations in an operation of this type is smaller than in the case of the embodiment in FIG.12, however due to an additional delay introduced by the filter bank 2
the risk of system's instability and unintended generation increases. Similarly, also in this case the system for associating objects 33 takes into account the parameters of the subtracted active objects. Due to mechanisms which have been described the method and the system according to the invention provide the frequency-domain resolution of at least ½ semitone (i.e. FN[n+l] / FN = 102,93%) [0059] According to the invention, the information contained in the active objects database 34 is also used by a shape forming system 37. The expected result of the sound signal decomposition according to the invention is to obtain sound objects having the form of sinusoidal waveforms with slowly- varying amplitude envelope and frequency. Therefore, the shape forming system 37 tracks variations of the amplitude envelope and frequency of the active objects in the database 34 and calculates online subsequent characteristic points of amplitude and frequency, which are the local maximum, local minimum and inflection points. Such information allows to unambiguously describe sinusoidal waveforms. The shape forming system 37 forwards these characteristic information in the form of points describing an object online to the active objects database 34. It has been assumed that the distance between points to be determined should be no less than 20 periods of the object's frequency. Distances between points, which are proportional to frequency, are capable to effectively represent dynamics of the objects' variation. Exemplary sound objects have been shown in FIG.14a. This figure illustrates four objects with frequency varying in function of time (sample number) . The same objects have been shown in FIG.14b in the space defined by amplitude and time (sample number). The illustrated points indicate local maxima and minima of the amplitude. The points are connected by a smooth curve, calculated with the use of third order polynomials. Having determined the function of frequency variation and the amplitude envelope it is possible to determine the audio signal. FIG.14c illustrates an audio signal determined based on the shape of the objects defined in FIG.14a and FIG.14b. The object shown in the plots have been described in the form of the table FIG.14d, wherein for each object there are described the parameters of its subsequent characteristic points, including the first point, the last point and the local extrema. Each point has three coordinates, i.e. the position in time expressed by the sample number, the amplitude and the frequency. Such set of points describes unambiguously a slowly-varying sinusoidal waveform [0060] The description of sound objects shown in the table FIG.14d can be written down in the form of a formalized protocol. Standardization of such notation will allow to develop applications using the properties of the sound objects according to the invention. FIG.15 shows an exemplary
format of sound objects notation. 1) Header : The notation starts with a header having as an essential element a header tag comprising a four byte keyword, informing that we deal with the description of sound objects. Next, in two bytes an information about the number of channels (tracks) is specified and two bytes of time unit definition. The header occurs only once at the beginning of a file. 2) Channel: Information about channels (tracks) from this field serves to separate the group of sound objects being in an essential relation, e.g. left or right channel in stereo, vocal track, percussion instruments track, recording from a defined microphone etc. The channel field comprises the channel identifier (number), the number of objects in the channel and the position of the channel from the beginning of an audio signal, measured in defined units. 3) Object: An identifier contained in the first byte decides about the type of the object decides. Identifier "0" denotes a basic unit in the signal record which is the sound object. Value "1" can denote a folder containing a group of objects like, for example, basic tone and its harmonics. Other values can be used to define other elements related to objects. The description of the fundamental sound object includes the number of points. The number of points does not include the first point, which is defined by the object itself. Specifying maximal amplitude in object's parameters allows to control simultaneous amplification of all points of the object. In the case of a folder of objects, this affects the value of amplitude of all the objects contained in the folder. Analogically, specifying information about frequency (applying notation: number of tone*4 of a filter bank = notes*16) allows to simultaneously control the frequency of all the elements related to an object. Furthermore, defining the position of the beginning of an object in relation to a higher level element (e.g. a channel) allows to shift the object in time. 4) Point: Points are used to describe the shape of the sound object in time-frequency-amplitude domain. They have relative value with respect to parameters defined by the sound object. One byte of amplitude defines which part of the maximal amplitude defined by the object the point has. Similarly, tone variation defines by what fraction of tone the frequency has changed. Position of point is defined as relative with respect to the previously defined point in the object. The multilevel structure of recording and relative associations between the fields allow a very flexible operation on sound objects, making them effective tools for designing and modifying audio signals. [0061] Condensed recording of information about sound objects according to the invention, in the format shown in FIG.15,
greatly affects in a positive way the size of registered and transferred files. Taking into account that an audio file can be readily played from this format, we can compare the size of the file shown in FIG.14c, which in .WAV format would contain over 2000 bytes, and in the form of sound objects record "UH0" according to the invention, it would contain 132 bytes. A compression better than 15-fold is not an excellent achievement in this case. In the case of longer audio signals much better results can be achieved. The compression level depends on how much information is contained in the audio signal, i.e. how many and how composed objects can be read from the signal. [0062] Identification of sound objects in an audio signal is not an unambiguous mathematical transformation. The audio signal created as a composition of objects obtained in the result of a decomposition differs from the input signal. The task of the system and the method according to the invention is to minimize this difference. Sources of differences are of two types. Part of them is expected and results from the applied technology, other can result from interference or unexpected properties of input audio signal. To reduce the difference between the audio signal composed of sound objects according to the invention and the input signal a correcting system 4, shown in FIG.l, is used. The system takes parameters of objects from the sound objects database 35 already after terminating the object and performs the operation of modification of selected parameters of objects and points such as to minimize the expected differences or irregularities localised in these parameters. [0063] The first type of correction of sound objects according to the invention, performed by the correcting system 4, is shown in FIG.16. The distortion at the beginning and at the end of the object is caused by the fact that during transient states, when the signal with defined freguency appears or fades, filters with a shorter impulse response react to the change quicker. Therefore, at the beginning the object is bent in the direction of higher frequencies, and at the end it turns towards the lower frequencies. Correction of an object can be based on deforming the object's frequency at the beginning and at the end in the direction defined by the middle section of the obj ect . [0064] A further type of correction according to the invention, performed by the correcting system 4, has been shown in FIG.17. The audio signal samples passing through a filter 20 of the filter bank 2 cause a change at the filter's output, which manifests as a signal shift. This shift has a regular character and is possible to be predicted. Its magnitude depends on the width of the window K of the filter n, the width being in accordance to the invention a function
of frequency. This means that each frequency is shifted by a different value, what affects the sound of the signal perceivably. The magnitude of the shift is ca . 1/2 filter window's width in the area of normal operation of the filter, 1/4 window's width in the initial phase and ca. 3/4 window's width in the case of the objects end. Because for each frequency the magnitude of the shift can be predicted, the task of the correcting system 4 is to properly shift all the points of the object in the opposite direction, so that the dynamics of the representation of the input signal improves. [0065] Yet another type of correction according to the invention, performed by the correcting system 4, is shown in FIG. 18a, FIG. 18B and FIG. 18C. The distortion manifests itself as an object splitting into pieces which are independent objects. This splitting can be caused e.g. by a phase fluctuation in an input signal's component, an interference or mutual influence of closely adjacent objects. The correction of distortions of this type requires the correcting circuit 4 to perform an analysis of the functions of envelope and frequency and to demonstrate that said objects should form an entirety. The correction is simple and is based on combination of the identified objects into one obj ect . [ 0066] A task of the correcting system 4 is also to remove objects having an insignificant influence on the audio signal's sound. According to the invention it was decided, that such objects can be the ones having the maximal amplitude which is lower than 1% of the maximal amplitude present in the whole signal at a given time instant. Change in the signal at the level of 40 dB should not be audible. [ 0067 ] The correcting system performs generally the removal of all irregularities in the shape of sound objects, which operations can be classified as: joining of discontinuous objects, removal of objects' oscillations near the adjacent ones, removal of insignificant objects, as well as the interfering ones, lasting too shortly or audible too weakly. [ 0068 ] To illustrate the results of the use of the method and the system for sound signal decomposition a fragment of stereo audio signal sampled at 44100 samples per second has been tested. The signal is a musical composition including sound of guitar and singing. The plot shown in FIG. 19a illustrating two channel includes ca . 250000 samples ( ca. 5.6 sec.) of the recording. [ 0069] FIG. 19b shows a spectrogram resulting from the operation of the filter bank 2 for the audio signal's left channel (upper plot in FIG.19a). The spectrogram includes the amplitude at the output of 450 filters having frequency from C2 = 16.35 Hz up to e6=10548 Hz. On the left side of the
spectrogram a piano keyboard has been shown as reference points defining the frequency. Furthermore, staffs with bass clef and a staff with treble clef above have been marked. The horizontal axis of the spectrogram corresponds to time instants during a composition, while the darker colour in the spectrogram indicates a higher value of the filtered signal's amplitude . [0070] FIG. 19c shows the result of operation of the voting system 32. Comparing the spectrogram in FIG. 19b with the spectrogram in FIG.19C it can be seen that wide spots representing signal constituent elements have been replaced by distinct lines indicating precise localisation of said constituent elements of the input signal. [0071] FIG.19d shows a cross-section of the spectrogram along the A-A line for the 149008th sample and presents the amplitude in function of frequency. The vertical axis in the middle indicates the real component and the imaginary component and the amplitude of the spectrum. The vertical axis at the right side shows peaks of the voting signal, indicating the temporary localisation of audio signal constituent elements. [0072] FIG. 19e is a cross-section of the spectrogram along the line BB at the frequency of 226.4 Hz. The plot shows the amplitude of the spectrum at the output of the filter 2 with the number n=182. [0073] In FIG. 19f sound objects are shown (without operation of the correcting system 4) . The vertical axis indicates the frequency, while the horizontal axis indicates time expressed by the number of the sample. In the tested fragment of the signal 578 objects have been localised, which are described by 578+995 = 1573 points. To store these objects ca. 9780 bytes are required. The audio signal in FIG. 19a comprising 250000 samples in the left channel requires 500 000 bytes for direct storing, which in the case of using the signal decomposition method and sound objects according to the invention leads to a compression at the level of 49. The use of correcting system 4 further improves the compression level, due to removal of objects having a negligible influence on the signal's sound. [0074] In FIG .19g there are shown amplitudes of selected sound objects, shaped with the use of already determined characteristic points by means of smooth curves created of third order polynomials. In the figure there are shown objects with amplitude higher than 10% of the amplitude of the object with the highest amplitude. [0075] As a result of using the method and the system for signal decomposition according to the invention one obtains
sound objects according to the invention, which can serve for an acoustic signal synthesis. A synthesized test signal in FIG.19a has been shown in FIG. 19h. [0076] The sound objects according to the invention have a number of properties enabling their multiple applications, in particular in processing, analysis and synthesis of sound signals. Sound objects can be acquired with the use of the method for signal decomposition according to the invention as a result of an audio signal decomposition. Sound objects can be also formed analytically, by defining values of parameters shown in FIG.14d. A sound object database can be formed by sounds taken from the surrounding environment or created artificially. Below some advantageous properties of sound objects described by points having three coordinates are listed : 1) Based on parameters describing sound objects it is possible to determine the function of amplitude and frequency variation, and to determine location in respect to other objects, so that an audio signal can be composed of them. 2) One of the parameters which describe sound objects is the time, thanks to which the objects can be shifted, shortened and lengthened in the time domain. 3) A second parameter of sound objects is the frequency, thanks to which the objects can be shifted and modified in the frequency domain. 4) A next parameter of sound objects is the amplitude, thanks to which envelopes of sound objects can be modified. 5) Sound objects can be grouped, by selecting e.g. the ones present in the same time or/and the ones with frequencies being harmonics. 6) Grouped objects can be separated from or appended to an audio signal. This allows to create a new signal from a number of other signals or to split a single signal into a number of independent signals. 7) Grouped objects can be amplified (by increasing their amplitude) or silenced (by decreasing their amplitude) . 8) By modifying proportions of harmonic amplitude included in a group of objects it is possible to modify the timbre of the grouped objects. 9) It is possible to modify the value of all grouped frequencies by increasing or decreasing frequencies of harmonics . 10) It is possible to modify audible emotions contained in sound objects, by modifying the slope (falling or raising) of component frequencies . 11) By presenting an audio signal in the form of objects described by points with three coordinates it is possible to significantly reduce the number of required data bytes without loss of information contained in the signal. [0077] Considering the properties of sound objects, a great deal of applications can be defined for them. The exemplary
ones include: 1) Separation of audio signal sources such as instruments or speakers, based on proper grouping of sound objects present in the signal. 2) Automatic generation of musical notation for individual instruments from an audio signal. 3) Devices for automatic tuning of musical instruments during ongoing musical performance. 4) Forwarding the voice of separated speakers to speech recognition systems. 5) Recognition of emotion contained in separated voices. 6) Identification of separated speakers. 7) Modification of the timbre of recognized instruments. 8) Swapping the instruments (e.g. a guitar playing instead of a piano) ; 9) Modification of a voice of a speaker (raising, lowering, conversion of emotion, intonation) . 10) Swapping of voices of speakers. 11) Synthesis of a voice with the possibility of emotion and intonation control. 12) Smooth joining of speeches. 13) Voice control of devices, even in an environment with interference . 14) Generation of new sounds, "samples", unusual sounds. 15) New musical instruments. 16) Spatial management of sound. 17) Additional possibilities of data compression.
The object of the invention is a method and a system for decomposition of acoustic signal into sound objects having the form of signals with slowly-varying amplitude and frequency, as well as sound objects and their use. The system for decomposition of acoustic signal into sound objects having the form of sinusoidal waveforms with slowly varying amplitude and frequency, comprising an A/D conversion system, connected to a filter bank with central frequencies of filters distributed according to geometrical distribution is characterised in that each filter (20) is adapted to determine the real value FC (n) and the imaginary value FS (n) of the filtered signal, the filter bank (2) is connected to a system for tracking objects (3), wherein the system for tracking objects (3) comprises a spectrum analysing system (31), a voting system (32) adapted to calculate the value of angular frequency weighted by amplitude, a resolution improvement system (36) adapted to adjust the filter's (20) window and/or to determine the difference between the spectrum of the input signal and the spectrum at the filter bank's (2) output, and/or to determine the difference between the input signal and the signal at the filter bank's (2) output, a system for associating objects (33), a shape forming system (37) adapted to determine characteristic points describing slowly-varying sinusoidal waveforms, a database of active objects (34) and a database of sound objects (35). AMENDED CLAIMS received by the International Bureau on 13 July 2015 (13.07.2015) 1. A method for decomposition of acoustic signal into sound objects having the form of sinusoidal wave with slowly- varying amplitude and frequency, comprising a step of determining parameters of short term signal model and a step of determining parameters of long term signal model based on said short term parameters, wherein a step of determining parameters of a short term signal model comprises a conversion of the analogue acoustic signal into a digital input signal PIN, characterised in that in said step of determining parameters of short term signal model the input signal ΡΪΝ is then split into adjacent sub- bands with central frequencies distributed according to logarithmic scale by feeding samples of the acoustic signal to the digital filter bank's input, each digital filter having a window length proportionally to the nominal central frequency - at each filter's (20) output the real value FC (n) and the imaginary value FS (n) of the filtered signal is determined sample by sample, and then based on this - the frequency, the amplitude and the phase of all detected constituent elements of said acoustic signal are determined sample by sample, - an operation improving the frequency-domain resolution of said filtered signal is executed sample by sample and involves at least a step of determining the frequency of all detected constituent elements based on maximum values of the function FG(n) resulting from a mathematical operation reflecting the number of neighbouring filters (20) outputting an angular frequency value substantially similar to an angular frequency value of each consecutive filter (20), and in that in said step of determining parameters of long term signal model: - for each detected element of said acoustic signal an active object in an active objects database (34) is created for its tracking - subsequent detected elements of said acoustic signal are associated sample by sample with at least selected active objects in said active objects database (34) to create a new active object or to append said detected element to an active object, or to close an active object - for each active object in the database (34) values of the envelope of amplitude and values of frequency and their corresponding time instants are determined not less frequently than once per period of duration of a given filter's (20) window W(n) so as to create characteristic points describing slowly-varying sinusoidal waveform of said sound object - at least one selected closed active object is transferred to a database of sound objects (35) to obtain at least one decomposed sound object, defined by a set of characteristic
points with coordinates in time-frequency-amplitude space. 2. The method according to claim 1, characterised in it further comprises a step of correcting selected sound objects which involves a step of correcting of amplitude and/or frequency of selected sound objects as to reduce an expected distortion in said sound objects, the distortion being introduced by said digital filter bank. 3. The method according to claim 1 or 2, characterised in that said operation improving the frequency-domain resolution of said filtered signal further comprises a step of increasing window length of selected filters. 4. The method according to claim 1 or 2, characterised in that said operation improving the frequency-domain resolution of said filtered signal further comprises a step of subtracting an expected spectrum of assuredly located adjacent sound objects from the spectrum at the output of the filters . 5. The method according to claim 1 or 2, characterised in that said operation improving the frequency-domain resolution of said filtered signal further comprises a step of subtracting an audio signal generated based on assuredly located adjacent sound objects from said input signal. 6. A system for decomposition of acoustic signal into sound objects having the form of sinusoidal waveforms with slowly- varying amplitude and frequency, comprising a sub-system for determining parameters of a short term signal model and a sub-system for determining parameters of a long term signal model based on said parameters, wherein said subsystem for determining short term parameters comprises a converter system for conversion of the analogue acoustic signal into a digital input signal PIN characterized in that said subsystem for determining short term parameters further comprises a filter bank (20) with filter central frequencies distributed according to logarithmic distribution, each digital filter having a window length proportionally to the central frequency wherein each filter (20) is adapted to determine a real value FC (n) and an imaginary value FS (n) of said filtered signal, said filter bank (2) being connected to a system for tracking objects (3), wherein said system for tracking objects (3) comprises a spectrum analysing system (31) adapted to detect all constituent elements of the input signal PIN? a voting system (32) adapted to determine the frequency of all detected constituent elements based on maximum values of the function FG(n) resulting from a mathematical operation reflecting the number of neighbouring filters (20) which output an angular frequency value substantially similar to an angular frequency value of each
consecutive filter (20), and in that said subsystem for determining long term parameters comprises a system for associating objects (33), a shape forming system (37) adapted to determine characteristic points describing slowly-varying sinusoidal waveforms, an active objects database (34) and a sound objects database (35) . 7. The system according to claim 6, characterised in that said system for tracking objects (3) is further connected with a correcting system (4) adapted to correct the amplitude and/or the frequency of individual selected sound objects so as to reduce an expected distortion in said sound objects introduced by said digital filter bank and/or adapted to combine discontinuous objects and/or to remove selected sound objects . 8. The system according to claim 7, characterised in that it further comprises a resolution improvement system (36) adapted to increase window length of selected filter and/or to subtract an expected spectrum of assuredly located adjacent sound objects from the spectrum at the output of the filters and/or to subtract an audio signal generated based on assuredly located adjacent sound objects from said input signal 9. A sound object, being a signal with slowly-varying amplitude and frequency, characterised in that it is obtained by the method according to any of claims 1 to 5. 10. A sound object, being a signal with slowly-varying amplitude and frequency, characterised in that it is defined by characteristic points having three coordinates in the time-amplitude-frequency space, wherein each characteristic point is distant from the next one in the time domain by a value proportional to the duration of a filter's (20) window W(n) assigned to the object's frequency. 11. Use of a sound object according to claims 9 to 10 for synthesis of acoustic signal. 12. Use of a sound object for synthesis of acoustic signal according to claim 11, characterised in that said sound object is previously shifted, shortened or lengthened in the time domain and/or is shifted or modified in the frequency domain and/or the envelope of the sound object is modified beforehand . 13. Use of a sound object for synthesis of audio signal according to claim 11 or 12, characterised in that said sound objects are previously grouped with respect to a time instant of their occurrence or with respect to the harmonic contents.
14. Use of a sound object according to claims 9 to _1_0 for automatic generation of musical notation of sounds from at least one source.