[0001]    Many innovations beyond two-channel stereo have failed because of cost, impracticability (e.g. number of loudspeakers), and last but not least a requirement for backwards compatibility. While 5.1 surround multi-channel audio systems are being adopted widely by consumers, also this system is compromised in terms of number of loudspeakers and with a backwards compatibility restriction (the front left and right loudspeakers are located at the same angles as in two-channel stereo, i.e. +/−30°, resulting in a narrow frontal virtual sound stage). 
         [0002]    It is a fact that by far most audio content is available in the two-channel stereo format. For audio systems enhancing the sound experience beyond stereo, it is thus crucial that stereo audio content can be played back, desirably with an improved experience compared to the legacy systems. 
         [0003]    It has long been realized that the use of more front loudspeakers improves the virtual sound stage also for listeners not exactly located in the sweet spot. There has been the aim of playing back stereo signals over more than two loudspeakers for improved results. Especially, there has been a lot of attention on playing back stereo signals with an additional center loudspeaker. However, the improvement of these techniques over conventional stereo playback has not been clear enough that they would have been widely used. The main limitations of these techniques are that they only consider localization and not explicitly other aspects such as ambience and listener envelopment. Further, the localization theory behind these techniques is based a one-virtual-source-scenario, limiting their performance when a number of sources are present at different directions simultaneously. 
         [0004]    These weaknesses are overcome by the techniques proposed in this description by using a perceptually motivated spatial decomposition of stereo audio signals. Given this decomposition, audio signals can be rendered for an increased number of loudspeakers, loudspeaker line arrays, and wavefield synthesis systems. 
         [0005]    The proposed techniques are not limited for conversion of (two channel) stereo signals to audio signals with more channels. But generally, a signal with L channels can be converted to a signal with M channels. The signals can either be stereo or multi-channel audio signals aimed for playback, or they can be raw microphone signals or linear combinations of microphone signals. It is also shown how the technique is applied to microphone signals (a.g. Ambisonics B-format) and matrixed surround downmix signals for reproducing these over various loudspeaker setups. 
         [0006]    When we refer to a stereo or multi-channel audio signal with a number of channels, we mean the same as when we refer to a number of (mono) audio signals. 
       SUMMARY OF THE INVENTION 
       [0007]    According to the main embodiment applying to multiple audio signals, it is proposed to generate multiple output audio signals (y 1 , . . . , y M ) from multiple input audio signals (x 1 , . . . , x L ), in which the number of output is equal or higher than the number of input signals, this method comprising the steps of:
       by means of linear combinations of the input subbands X 1 (i), . . . , X L (i), computing one or more independent sound subbands representing signal components which are independent between the input subbands,   by means of linear combinations of the input subbands X 1 (i), . . . , X L (i), computing one or more localized direct sound subbands representing signal components which are contained in more than one of the input subbands and direction factors representing the ratios with which these signal components are contained in two or more input subbands,   generating the output subband signals, Y 1 (i) . . . Y M (i), where each output subband signal is a linear combination of the independent sound subbands and the localized direct sound subbands   converting the output subband signals, Y 1 (i) . . . Y M (i), to time domain audio signals, y 1  . . . y M .       
 
         [0012]    The index i is the index of the subband considered. According to a first embodiment, this method can be used with only one subband per audio channel, even if more subbands per channel give a better acoustic result. 
         [0013]    The proposed scheme is based on the following reasoning. A number of input audio signals x 1 , . . . , x L  are decomposed into signal components representing sound which is independent between the audio channels and signal components which represent sound which is correlated between the audio channels. This is motivated by the different perceptual effect these two types of signal components have. The independent signal components represent information on source width, listener envelopment, and ambience and the correlated (dependent) signal components represent the localization of auditory events or acoustically the direct sound. To each correlated signal component there is associated directional information which can be represented by the ratios with which this sound is contained in a number of audio input signals. Given this decomposition, a number of audio output signals can be generated with the aim of reproducing a specific auditory spatial image when played back over loudspeakers (or headphones). The correlated signal components are rendered to the output signals (y 1 , . . . , y M ) such that it is perceived by a listener from a desired direction. The independent signal components are rendered to the output signals (loudspeakers) such that it mimics non-direct sound and its desired perceptual effect. This functionality, described on a high level, is taking the spatial information from the input audio signals and transforming this spatial information to spatial information in the output channels with desired properties. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0014]    The invention will be better understood thanks to the attached drawings in which: 
           [0015]      FIG. 1  shows a standard stereo loudspeaker setup, 
           [0016]      FIG. 2  shows the location of the perceived auditory events for different level differences for two coherent loudspeaker signals, the level and time difference between a pair of coherent loudspeaker signals determining the location of the auditory event which appears between the two loudspeakers, 
           [0017]      FIG. 3  ( a ) shows early reflections emitted from the side loudspeakers having the effect of widening of the auditory event. 
           [0018]      FIG. 3  ( b ) shows late reflections emitted from the side loudspeakers relating more to the environment as listener envelopment, 
           [0019]      FIG. 4  shows a way to mix a stereo signal mimicking direct sound and lateral reflections, 
           [0020]      FIG. 5  shows time-frequency tiles representing the decomposition of the signal into subband as a function of time, 
           [0021]      FIG. 6  shows the direction direction factor A and the normalized power of S and AS, 
           [0022]      FIG. 7  shows the least squares estimate weights w 1  and w 2  and the post scaling factor for the computation of the estimate of s, 
           [0023]      FIG. 8  shows the least squares estimate weights w 3  and w 4  and the post scaling factor for the computation of the estimate of N 1 , 
           [0024]      FIG. 9  shows the least squares estimate weights w 5  and w 6  and the post scaling factor for the computation of the estimate of N 2 , 
           [0025]      FIG. 10  shows the estimated s, A, n 1  and n 2 , 
           [0026]      FIG. 11  shows the ±30° virtual sound stage (a) converted to a virtual sound stage with the width of the aperture of a loudspeaker array (b) 
           [0027]      FIG. 12  shows loudspeaker pair selection l and factors a 1  and a 2  as a function of the stereo signal level difference, 
           [0028]      FIG. 13  shows an emission of plane waves through a plurality of loudspeakers, 
           [0029]      FIG. 14  shows the ±30° virtual sound stage (a) converted to a virtual sound stage with the width of the aperture of a loudspeaker array with increased listener envelopment by emitting independent sound from the side loudspeakers (b), 
           [0030]      FIG. 15  shows the eight signals, generated for a setup as in  FIG. 14(   b ), 
           [0031]      FIG. 16  shows each signal corresponding to the front sound stage defined as a virtual source. The independent lateral sound is emitted as plane waves (virtual sources in the far field) 
           [0032]      FIG. 17  shows a quadraphonic sound system (a) extended for use with more loudspeakers (b). 
       
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     Spatial Hearing and Stereo Loudspeaker Playback 
       [0033]    The proposed scheme is motivated an described for the important case of two input channels (stereo audio input) and M audio output channels (M≧2). Later, it is described how to apply the same reasoning as derived at the example of stereo input signals to the more general case of L input channels. 
         [0034]    The most commonly used consumer playback system for spatial audio is the stereo loudspeaker setup as shown in  FIG. 1 . Two loudspeakers are placed in front on the left and right sides of the listener. Usually, these loudspeakers are placed on a circle at angles −30° and +30°. The width of the auditory spatial image that is perceived when listening to such a stereo playback system is limited approximately to the area between and behind the two loudspeakers. 
         [0035]    The perceived auditory spatial image, in natural listening and when listening to reproduced sound, largely depends on the binaural localization cues, i.e. the interaural time difference (ITD), interaural level difference (ILD), and interaural coherence (IC). Furthermore, it has been shown that the perception of elevation is related to monaural cues. 
         [0036]    The ability to produce an auditory spatial image mimicking a sound stage with stereo loudspeaker playback is made possible by the perceptual phenomenon of summing localization, i.e. an auditory event can be made appear at any angle between a loudspeaker pair in front of a listener by controlling the level and/or time difference between the signals given to the loudspeakers. It was Blumlein in the 1930&#39;s who recognized the power of this principle and filed his now-famous patent on stereophony. Summing localization is based on the fact that ITD and ILD cues evoked at the ears crudely approximate the dominating cues that would appear if a physical source were located at the direction of the auditory event which appears between the loudspeakers. 
         [0037]      FIG. 2  illustrates the location of the perceived auditory events for different level differences for two coherent loudspeaker signals. When the left and right loudspeaker signals are coherent, have the same level, and no delay difference, an auditory event appears in the center between the two loudspeakers as illustrated by Region  1  in  FIG. 2 . By increasing the level on one side, e.g. right, the auditory event moves to that side as illustrated by Region  2  in  FIG. 2 . In the extreme case, when only the signal on the left is active, the auditory event appears at the left loudspeaker position as is illustrated by Region  3  in  FIG. 2 . The position of the auditory event can be similarly controlled by varying the delay between the loudspeaker signals. The described principle of controlling the location of an auditory event between a loudspeaker pair is also applicable when the loudspeaker pair is not in the front of the listener. However, some restrictions apply for loudspeakers to the sides of a listener. 
         [0038]    As illustrated in  FIG. 2 , summing localization can be used to mimic a scenario where different instruments are located at different directions on a virtual sound stage, i.e. in the region between the two loudspeakers. In the following, it is described how other attributes than localization can be controlled. 
         [0039]    Important in concert hall acoustics is the consideration of reflections arriving at the listener from the sides, i.e. lateral reflections. It has been shown that early lateral reflections have the effect of widening the auditory event. The effect of early reflections with delays smaller than about 80 ms is approximately constant and thus a physical measure, denoted lateral fraction, has been defined considering early reflections in this range. The lateral fraction is the ratio of the lateral sound energy to the total sound energy that arrived within the first 80 ms after the arrival of the direct sound and measures the width of the auditory event. 
         [0040]    An experimental setup for emulating early lateral reflections is illustrated in  FIG. 3(   a ). The direct sound is emitted from the center loudspeaker while independent early reflections are emitted from the left and right loudspeakers. The width of the auditory event increases as the relative strength of the early lateral reflections is increased. 
         [0041]    More than 80 ms after the arrival of the direct sound, lateral reflections tend to contribute more to the perception of the environment than to the auditory event itself. This is manifested in a sense of “envelopment” or “spaciousness of the environment”, frequently denoted listener envelopment. A similar measure as the lateral fraction for early reflections is also applicable to late reflections for measuring the degree of listener envelopment. This measure is denoted late lateral energy fraction. 
         [0042]    Late lateral reflections can be emulated with a setup as shown in  FIG. 3(   b ). The direct sound is emitted from the center loudspeaker while independent late reflections are emitted from the left and right loudspeakers. The sense of listener envelopment increases as the relative strength of the late lateral reflections is increased, while the width of the auditory event is expected to be hardly affected. 
         [0043]    Stereo signals are recorded or mixed such that for each source the signal goes coherently into the left and right signal channel with specific directional cues (level difference, time difference) and reflected/reverberated independent signals go into the channels determining auditory event width and listener envelopment cues. It is out of the scope of this description to further discuss mixing and recording techniques. 
       Spatial Decomposition of Stereo Signals 
       [0044]    As opposed to using a direct sound from a real source, as was illustrated in  FIG. 3 , one can use direct sound corresponding to a virtual source generated with summing localization. The shaded areas indicate the perceived auditory events. That is, experiments as are shown in  FIG. 3  can be carried out with only two loudspeakers. This is illustrated in  FIG. 4 , where the signal s mimics the direct sound from a direction determined by the factor a. The independent signals, n 1  and n 2 , correspond to the lateral reflections. The described scenario is a natural decomposition for stereo signals with one auditory event, 
         [0000]        x   1 ( n )= s ( n )+ n   1 ( n )  x   2 ( n )= as ( n )+ n   2 ( n )  (1) 
         [0000]    capturing the localization and width of the auditory event and listener envelopment. 
         [0045]    In order to get a decomposition which is not only effective in a one auditory event scenario, but non-stationary scenarios with multiple concurrently active sources, the described decomposition is carried out independently in a number of frequency bands and adaptively in time, 
         [0000]        X   1 ( i,k )= S ( i,k )+ N   1 ( i,k )  X   2 ( i,k )= A ( i,k ) S ( i,k )+ N   2 ( i,k )  (2) 
         [0000]    where i is the subband index and k is the subband time index. This is illustrated in  FIG. 5 , i.e. in each time-frequency tile with indices i and k, the signals S, N 1 , N 2 , and direction factor A are estimated independently. For brevity of notation, the subband and time indices are often ignored in the following. We are using a subband decomposition with perceptually motivated subband bandwidths, i.e. the bandwidth of a subband is chosen to be equal to one critical band. S, N 1 , N 2 , and direction factor A are estimated approximately every 20 ms in each subband. 
         [0046]    Note that more generally one could also consider a time difference of the direct sound in equation (2). That is, one would not only use an direction factor A, but also a direction delay which would be defined as the delay with which S is contained in X 1  and X 2 . In the following description we do not consider such a delay, but it is understood that the analysis can easily be extended to consider such a delay. 
         [0047]    Given the stereo subband signals, X 1  and X 2 , the goal is to compute estimates of S, N 1 , N 2 , and A. A short-time estimate of the power of X 1  is denoted P X     1   (i,k)=E{X 1   2 (i,k)}. For the other signals, the same convention is used, i.e. P X     2   , P s  and P N =P N     1   =P N2  are the corresponding short-time power estimates. The power of N 1  and N 2  is assumed to be the same, i.e. it is assumed that the amount of lateral independent sound is the same for left and right. 
         [0048]    Note that other assumptions than P N =P N     1   =P N2  may be used. For example A 2 P N     1   =P N2    
       Estimating P S , A, and P N    
       [0049]    Given the subband representation of the stereo signal, the power (P X     1   , P X     2   ) and the normalized cross-correlation are computed. The normalized cross-correlation between left and right is: 
         [0000]    
       
         
           
             
               
                 
                   
                     Φ 
                      
                     
                       ( 
                       
                         i 
                         , 
                         k 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       E 
                        
                       
                         { 
                         
                           
                             
                               X 
                               1 
                             
                              
                             
                               ( 
                               
                                 i 
                                 , 
                                 k 
                               
                               ) 
                             
                           
                            
                           
                             
                               X 
                               2 
                             
                              
                             
                               ( 
                               
                                 i 
                                 , 
                                 k 
                               
                               ) 
                             
                           
                         
                         } 
                       
                     
                     
                       
                         E 
                          
                         
                           { 
                           
                             
                               X 
                               1 
                               2 
                             
                              
                             
                               ( 
                               
                                 i 
                                 , 
                                 k 
                               
                               ) 
                             
                           
                           } 
                         
                          
                         E 
                          
                         
                           { 
                           
                             
                               X 
                               2 
                               2 
                             
                              
                             
                               ( 
                               
                                 i 
                                 , 
                                 k 
                               
                               ) 
                             
                           
                           } 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0050]    A, P S , and P N  are computed as a function of the estimated Px 1 , Px 2  and Φ. Three equations relating the known and unknown variables are: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       Px 
                       1 
                     
                     = 
                     
                       
                         P 
                         S 
                       
                       + 
                       
                         P 
                         N 
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       Px 
                       2 
                     
                     = 
                     
                       
                         
                           A 
                           2 
                         
                          
                         
                           P 
                           S 
                         
                       
                       + 
                       
                         P 
                         N 
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     Φ 
                     = 
                     
                       aS 
                       
                         
                           
                             Px 
                             1 
                           
                         
                          
                         
                           Px 
                           2 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   4 
                   ) 
                 
               
             
           
         
       
     
         [0051]    These equations solved for A, P S , and P N , yield 
         [0000]    
       
         
           
             
               
                 
                   
                     A 
                     = 
                     
                       B 
                       
                         2 
                          
                         C 
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       P 
                       S 
                     
                     = 
                     
                       
                         2 
                          
                         
                           C 
                           2 
                         
                       
                       B 
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       P 
                       N 
                     
                     = 
                     
                       
                         X 
                         1 
                       
                       - 
                       
                         
                           2 
                            
                           
                             C 
                             2 
                           
                         
                         B 
                       
                     
                   
                    
                   
                     
 
                   
                    
                   with 
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
             
               
                 
                   
                     B 
                     = 
                     
                       
                         Px 
                         2 
                       
                       - 
                       
                         Px 
                         1 
                       
                       + 
                       
                         
                           
                             
                               ( 
                               
                                 
                                   Px 
                                   1 
                                 
                                 - 
                                 
                                   Px 
                                   2 
                                 
                               
                               ) 
                             
                             2 
                           
                           + 
                           
                             4 
                              
                             
                               Px 
                               1 
                             
                              
                             
                               Px 
                               2 
                             
                              
                             
                               Φ 
                               2 
                             
                           
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     C 
                     = 
                     
                       Φ 
                        
                       
                         
                           
                             Px 
                             1 
                           
                            
                           
                             Px 
                             2 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
       Least Squares Estimation of S, N 1  and N 2    
       [0052]    Next, the least squares estimates of S, N 1  and N 2  are computed as a function of A, P S , and P N . For each i and k, the signal S is estimated as 
         [0000]        Ŝ=ω   1   X   1 +ω 2   X   2 =ω 1 ( S+N   1 )+ω 2 ( AS+N   2 )  (7) 
         [0000]    where ω 1  and ω 2  are real-valued weights. The estimation error is 
         [0000]        E= (1−ω 1 −ω 2   A ) S−ω   1   N   1 −ω 2   N   2   (8) 
         [0053]    The weights ω 1  and ω 2  are optimal in a least mean square sense when the error E is orthogonal to X 1  and X 2 , i.e. 
         [0000]      E{EX 1 }=0 E{EX 2 }=0  (9) 
         [0000]    yielding two equations, 
         [0000]      (1−ω 1 −ω 2   A ) P   s −ω 1   P   N =0, 
         [0000]        A (1−ω 1 −ω 2   A ) P   s −ω 2   P   N =0  (10) 
         [0000]    from which the weights are computed, 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ω 
                       1 
                     
                     = 
                     
                       
                         
                           P 
                           S 
                         
                          
                         
                           P 
                           N 
                         
                       
                       
                         
                           
                             ( 
                             
                               
                                 A 
                                 2 
                               
                               + 
                               1 
                             
                             ) 
                           
                            
                           
                             P 
                             S 
                           
                            
                           
                             P 
                             N 
                           
                         
                         + 
                         
                           P 
                           N 
                           2 
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       ω 
                       2 
                     
                     = 
                     
                       
                         
                           AP 
                           S 
                         
                          
                         
                           P 
                           N 
                         
                       
                       
                         
                           
                             ( 
                             
                               
                                 A 
                                 2 
                               
                               + 
                               1 
                             
                             ) 
                           
                            
                           
                             P 
                             S 
                           
                            
                           
                             P 
                             N 
                           
                         
                         + 
                         
                           P 
                           N 
                           2 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   11 
                   ) 
                 
               
             
           
         
       
     
         [0054]    Similarly, N 1  and N 2 , are estimated. The estimate of N 1  is 
         [0000]        {circumflex over (N)}   1 =ω 3   X   1 +ω 4   X   2 =ω 3 ( S+N   1 )+ω 4 ( AS+N   2 )  (12) 
         [0055]    The estimation error is 
         [0000]        E =(ω 3 −ω 4   A ) S −(1−ω 3 ) N   1 −ω 2   N 2  (13) 
         [0056]    Again, the weights are computed such that the estimation error is orthogonal to X 1  and X 2  resulting in 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ω 
                       3 
                     
                     = 
                     
                       
                         
                           
                             A 
                             2 
                           
                            
                           
                             P 
                             S 
                           
                            
                           
                             P 
                             N 
                           
                         
                         + 
                         
                           P 
                           N 
                           2 
                         
                       
                       
                         
                           
                             ( 
                             
                               
                                 A 
                                 2 
                               
                               + 
                               1 
                             
                             ) 
                           
                            
                           
                             P 
                             S 
                           
                            
                           
                             P 
                             N 
                           
                         
                         + 
                         
                           P 
                           N 
                           2 
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       ω 
                       4 
                     
                     = 
                     
                       
                         
                           - 
                           
                             AP 
                             S 
                           
                         
                          
                         
                           P 
                           N 
                         
                       
                       
                         
                           
                             ( 
                             
                               
                                 A 
                                 2 
                               
                               + 
                               1 
                             
                             ) 
                           
                            
                           
                             P 
                             S 
                           
                            
                           
                             P 
                             N 
                           
                         
                         + 
                         
                           P 
                           N 
                           2 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   14 
                   ) 
                 
               
             
           
         
       
     
         [0057]    The weights for computing the least squares estimate of N 2  are 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         N 
                         ^ 
                       
                       2 
                     
                     = 
                     
                       
                         
                           
                             ω 
                             5 
                           
                            
                           
                             X 
                             1 
                           
                         
                         + 
                         
                           
                             ω 
                             6 
                           
                            
                           
                             X 
                             2 
                           
                         
                       
                       = 
                       
                         
                           
                             ω 
                             5 
                           
                            
                           
                             ( 
                             
                               S 
                               + 
                               
                                 N 
                                 1 
                               
                             
                             ) 
                           
                         
                         + 
                         
                           
                             ω 
                             6 
                           
                            
                           
                             ( 
                             
                               AS 
                               + 
                               
                                 N 
                                 2 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   are 
                 
               
               
                 
                   ( 
                   15 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       ω 
                       5 
                     
                     = 
                     
                       
                         
                           - 
                           
                             AP 
                             S 
                           
                         
                          
                         
                           P 
                           N 
                         
                       
                       
                         
                           
                             ( 
                             
                               
                                 A 
                                 2 
                               
                               + 
                               1 
                             
                             ) 
                           
                            
                           
                             P 
                             S 
                           
                            
                           
                             P 
                             N 
                           
                         
                         + 
                         
                           P 
                           N 
                           2 
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       ω 
                       6 
                     
                     = 
                     
                       
                         
                           
                             P 
                             S 
                           
                            
                           
                             P 
                             N 
                           
                         
                         + 
                         
                           P 
                           N 
                           2 
                         
                       
                       
                         
                           
                             ( 
                             
                               
                                 A 
                                 2 
                               
                               + 
                               1 
                             
                             ) 
                           
                            
                           
                             P 
                             S 
                           
                            
                           
                             P 
                             N 
                           
                         
                         + 
                         
                           P 
                           N 
                           2 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   16 
                   ) 
                 
               
             
           
         
       
     
       Post-Scaling 
       [0058]    Given the least squares estimates, these are (optionally) post-scaled such that the power of the estimates Ŝ, {circumflex over (N)} 1 , {circumflex over (N)} 2  equals to P S  and P N =P N1 =P N2 . The power of Ŝ is 
         [0000]        P   Ŝ =(ω 1   +aω   2 ) 2   P   s +(ω 1   2 +ω 2   2 ) P   N   (17) 
         [0059]    Thus, for obtaining an estimate of S with power P S , Ŝ is scaled 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       S 
                       ^ 
                     
                     ′ 
                   
                   = 
                   
                     
                       
                         
                           P 
                           N 
                         
                       
                       
                         
                           
                             
                               
                                 ( 
                                 
                                   
                                     ω 
                                     1 
                                   
                                   + 
                                   
                                     a 
                                      
                                     
                                         
                                     
                                      
                                     
                                       ω 
                                       2 
                                     
                                   
                                 
                                 ) 
                               
                               2 
                             
                              
                             
                               P 
                               S 
                             
                           
                           + 
                           
                             
                               ( 
                               
                                 
                                   ω 
                                   1 
                                   2 
                                 
                                 + 
                                 
                                   ω 
                                   2 
                                   2 
                                 
                               
                               ) 
                             
                              
                             
                               P 
                               N 
                             
                           
                         
                       
                     
                      
                     
                       S 
                       ^ 
                     
                   
                 
               
               
                 
                   ( 
                   18 
                   ) 
                 
               
             
           
         
       
     
         [0060]    With similar reasoning, {circumflex over (N)} 1  and {circumflex over (N)} 2  are scaled, i.e. 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         N 
                         ^ 
                       
                       1 
                       ′ 
                     
                     = 
                     
                       
                         
                           
                             P 
                             N 
                           
                         
                         
                           
                             
                               
                                 
                                   ( 
                                   
                                     
                                       ω 
                                       3 
                                     
                                     + 
                                     
                                       a 
                                        
                                       
                                           
                                       
                                        
                                       
                                         ω 
                                         4 
                                       
                                     
                                   
                                   ) 
                                 
                                 2 
                               
                                
                               
                                 P 
                                 S 
                               
                             
                             + 
                             
                               
                                 ( 
                                 
                                   
                                     ω 
                                     3 
                                     2 
                                   
                                   + 
                                   
                                     ω 
                                     4 
                                     2 
                                   
                                 
                                 ) 
                               
                                
                               
                                 P 
                                 N 
                               
                             
                           
                         
                       
                        
                       
                         
                           N 
                           ^ 
                         
                         1 
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       
                         N 
                         ^ 
                       
                       2 
                       ′ 
                     
                     = 
                     
                       
                         
                           
                             P 
                             N 
                           
                         
                         
                           
                             
                               
                                 
                                   ( 
                                   
                                     
                                       ω 
                                       5 
                                     
                                     + 
                                     
                                       a 
                                        
                                       
                                           
                                       
                                        
                                       
                                         ω 
                                         6 
                                       
                                     
                                   
                                   ) 
                                 
                                 2 
                               
                                
                               
                                 P 
                                 S 
                               
                             
                             + 
                             
                               
                                 ( 
                                 
                                   
                                     ω 
                                     5 
                                     2 
                                   
                                   + 
                                   
                                     ω 
                                     6 
                                     2 
                                   
                                 
                                 ) 
                               
                                
                               
                                 P 
                                 N 
                               
                             
                           
                         
                       
                        
                       
                         
                           N 
                           ^ 
                         
                         2 
                       
                     
                   
                 
               
               
                 
                   ( 
                   19 
                   ) 
                 
               
             
           
         
       
     
       NUMERICAL EXAMPLES 
       [0061]    The direction factor A and the normalized power of S and AS are shown as a function of the stereo signal level difference and Φ in  FIG. 6 . 
         [0062]    The weights ω 1  and ω 2  for computing the least squares estimate of S are shown in the top two panels of  FIG. 7  as a function of the stereo signal level difference and Φ. The post-scaling factor for Ŝ (18) is shown in the bottom panel. 
         [0063]    The weights ω 3  and ω 2  for computing the least squares estimate of N 1  and the corresponding post-scaling factor (19) are shown in  FIG. 7  as a function of the stereo signal level difference and Φ. 
         [0064]    The weights ω 5  and ω 6  for computing the least squares estimate of N 2  and the corresponding post-scaling factor (19) are shown in  FIG. 7  as a function of the stereo signal level difference and Φ. 
         [0065]    An example for the spatial decomposition of a stereo rock music clips with a singer in the center is shown in  FIG. 10 . The estimates of s, A, n 1  and n 2  are shown. The signals are shown in the time-domain and A is shown for every time-frequency tile. The estimated direct sound s is relatively strong compared to the independent lateral sound n 1  and n 2  since the singer in the center is dominant. 
       Playing Back the Decomposed Stereo Signals Over Different Playback Setups 
       [0066]    Given the spatial decomposition of the stereo signal, i.e. the subband signals for the estimated localized direct sound Ŝ′, the direction factor A, and the lateral independent sound {circumflex over (N)} 1 ′ and {circumflex over (N)} 2 ′, one can define rules on how to emit the signal components corresponding to Ŝ′, {circumflex over (N)} 1 ′ and {circumflex over (N)} 2 ′, from different playback setups. 
       Multiple Loudspeakers in Front of the Listener 
       [0067]      FIG. 11  illustrates the scenario that is addressed. The virtual sound stage of width Φ 0 =30°, shown in Part (a) of the figure, is scaled to a virtual sound stage of width Φ 0 ′ which is reproduced with multiple loudspeakers, shown in Part (b) of the figure. 
         [0068]    The estimated independent lateral sound, {circumflex over (N)}′ 1  and {circumflex over (N)}′ 2 , is emitted from the loudspeakers on the sides, e.g. loudspeakers  1  and  6  in  FIG. 11(   b ). That is, because the more the lateral sound is emitted from the side the more it is effective in terms enveloping the listener into the sound. Given the estimated direction factor A, the angle Φ of the auditory event relative to the ±Φ 0  virtual sound stage is estimated, using the “stereophonic law of sines” (or other laws relating A to the perceived angle), 
         [0000]    
       
         
           
             
               
                 
                   φ 
                   = 
                   
                     
                       sin 
                       
                         - 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         
                           
                             A 
                             - 
                             1 
                           
                           
                             A 
                             + 
                             1 
                           
                         
                          
                         sin 
                          
                         
                             
                         
                          
                         
                           φ 
                           0 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   20 
                   ) 
                 
               
             
           
         
       
     
         [0069]    This angle is linearly scaled to compute the angle relative to the widened sound stage, 
         [0000]    
       
         
           
             
               
                 
                   
                     φ 
                     ′ 
                   
                   = 
                   
                     
                       
                         φ 
                         0 
                         ′ 
                       
                       
                         φ 
                         0 
                       
                     
                      
                     φ 
                   
                 
               
               
                 
                   ( 
                   21 
                   ) 
                 
               
             
           
         
       
     
         [0070]    The loudspeaker pair enclosing Φ′ is selected. In the example illustrated in  FIG. 11(   b ) this pair has indices  4  and  5 . The angles relevant for amplitude panning between this loudspeaker pair, γ 0  and γ 1 , are defined as shown in the figure. If the selected loudspeaker pair has indices l and l+1 then the signals given to these loudspeakers are 
         [0000]      a 1 √{square root over (1+A 2 S)} 
         [0000]      a 2 √{square root over (1+A 2 S)}  (22) 
         [0000]    where the amplitude panning factors a 1  and a 2  are computed with the stereophonic law of sines (or another amplitude panning law) and normalized such that a 1   2 +a 2   2 =1, 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       a 
                       1 
                     
                     = 
                     
                       1 
                       
                         
                           1 
                           + 
                           
                             C 
                             2 
                           
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       a 
                       2 
                     
                     = 
                     
                       C 
                       
                         
                           1 
                           + 
                           
                             C 
                             2 
                           
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   with 
                 
               
               
                 
                   ( 
                   23 
                   ) 
                 
               
             
             
               
                 
                   C 
                   = 
                   
                     
                       sin 
                        
                       
                         ( 
                         
                           
                             γ 
                             0 
                           
                           + 
                           γ 
                         
                         ) 
                       
                     
                     
                       sin 
                        
                       
                         ( 
                         
                           
                             γ 
                             0 
                           
                           - 
                           γ 
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   ( 
                   24 
                   ) 
                 
               
             
           
         
       
     
         [0071]    The factors in √{square root over ( 1 +A 2 )} in (22) are such that the total power of these signals is equal to the total power of the coherent components, S and AS, in the stereo signal. Alternatively, one can use amplitude panning laws which give signal to more than two loudspeakers simultaneously. 
         [0072]      FIG. 12  shows an example for the selection of loudspeaker pairs, l and l+1, and the amplitude panning factors a 1  and a 2  for Φ′ 0 =Φ 0 =30° for M=8 loudspeakers at angles {−30°, −20°, −12°, −4°, 4°, 12°, 20°, 30°}. 
         [0073]    Given the above reasoning, each time-frequency tile of the output signal channels, i and k, is computed as 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       Y 
                       m 
                     
                     = 
                     
                       
                         
                           δ 
                            
                           
                             ( 
                             
                               m 
                               - 
                               1 
                             
                             ) 
                           
                         
                          
                         
                           
                             N 
                             ^ 
                           
                           1 
                           ′ 
                         
                       
                       + 
                       
                         
                           δ 
                            
                           
                             ( 
                             
                               m 
                               - 
                               M 
                             
                             ) 
                           
                         
                          
                         
                           
                             N 
                             ^ 
                           
                           2 
                           ′ 
                         
                       
                       + 
                       
                         
                           ( 
                           
                             
                               
                                 δ 
                                  
                                 
                                   ( 
                                   
                                     m 
                                     - 
                                     l 
                                   
                                   ) 
                                 
                               
                                
                               
                                 a 
                                 1 
                               
                             
                             + 
                             
                               
                                 δ 
                                  
                                 
                                   ( 
                                   
                                     m 
                                     - 
                                     l 
                                     - 
                                     1 
                                   
                                   ) 
                                 
                               
                                
                               
                                 a 
                                 2 
                               
                             
                           
                           ) 
                         
                          
                         
                           
                             1 
                             + 
                             
                               A 
                               2 
                             
                           
                         
                          
                         
                           
                             S 
                             ^ 
                           
                           ′ 
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                       
                   
                    
                   where 
                 
               
               
                 
                   ( 
                   25 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       δ 
                        
                       
                         ( 
                         m 
                         ) 
                       
                     
                     = 
                     
                       { 
                       
                         
                           
                             1 
                           
                           
                             
                               
                                 for 
                                  
                                 
                                     
                                 
                                  
                                 m 
                               
                               = 
                               0 
                             
                           
                         
                         
                           
                             0 
                           
                           
                             otherwise 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   26 
                   ) 
                 
               
             
           
         
       
     
         [0000]    and m is the output channel index 1≦m≦M. The subband signals of the output channels are converted back to the time domain and form the output channels y 1  to y M . In the following, this last step is not always again explicitly mentioned. 
         [0074]    A limitation of the described scheme is that when the listener is at one side, e.g. close to loudspeaker  1 , the lateral independent sound will reach him with much more intensity than the lateral sound from the other side. This problem can be circumvented by emitting the lateral independent sound from all loudspeakers with the aim of generating two lateral plane waves. This is illustrated in  FIG. 13 . The lateral independent sound is given to all loudspeakers with delays mimicking a plane wave with a certain direction, 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       Y 
                       m 
                     
                      
                     
                       ( 
                       
                         i 
                         , 
                         k 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       
                         
                           
                             N 
                             ^ 
                           
                           1 
                           ′ 
                         
                          
                         
                           ( 
                           
                             i 
                             , 
                             
                               k 
                               - 
                               
                                 
                                   ( 
                                   
                                     m 
                                     - 
                                     1 
                                   
                                   ) 
                                 
                                  
                                 d 
                               
                             
                           
                           ) 
                         
                       
                       
                         M 
                       
                     
                     + 
                     
                       
                         
                           
                             N 
                             ^ 
                           
                           2 
                           ′ 
                         
                          
                         
                           ( 
                           
                             i 
                             , 
                             
                               k 
                               - 
                               
                                 
                                   ( 
                                   
                                     M 
                                     - 
                                     m 
                                   
                                   ) 
                                 
                                  
                                 d 
                               
                             
                           
                           ) 
                         
                       
                       
                         M 
                       
                     
                     + 
                     
                       
                         ( 
                         
                           
                             
                               δ 
                                
                               
                                 ( 
                                 
                                   m 
                                   - 
                                   l 
                                 
                                 ) 
                               
                             
                              
                             
                               a 
                               1 
                             
                           
                           + 
                           
                             
                               δ 
                                
                               
                                 ( 
                                 
                                   m 
                                   - 
                                   l 
                                   - 
                                   1 
                                 
                                 ) 
                               
                             
                              
                             
                               a 
                               2 
                             
                           
                         
                         ) 
                       
                        
                       
                         
                           1 
                           + 
                           
                             A 
                             2 
                           
                         
                       
                        
                       
                         
                           S 
                           ^ 
                         
                         ′ 
                       
                     
                   
                 
               
               
                 
                   ( 
                   27 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where d is the delay, 
         [0000]    
       
         
           
             
               
                 
                   d 
                   = 
                   
                     
                       
                         sf 
                         s 
                       
                        
                       sin 
                        
                       
                           
                       
                        
                       α 
                     
                     v 
                   
                 
               
               
                 
                   ( 
                   28 
                   ) 
                 
               
             
           
         
       
     
         [0000]    s is the distance between the equally spaced loudspeakers, v is the speed of sound, f s  is the subband sampling frequency, and ±α are the directions of propagation of the two plane waves. In our system, the subband sampling frequency is not high enough such that d can be expressed as an integer. Thus, we are first converting {circumflex over (N)}′ 1  and {circumflex over (N)}′ 2  to the time-domain and then we add its various delayed versions to the output channels. 
       Multiple Front Loudspeakers Plus Side Loudspeakers 
       [0075]    The previously described playback scenario aims at widening the virtual sound stage and at making the perceived sound stage independent of the location of the listener. 
         [0076]    Optionally one can play back the independent lateral sound, {circumflex over (N)}′ 1  and {circumflex over (N)}′ 2  with separate two loudspeakers located more to the sides of the listener, as illustrated in  FIG. 14 . The ±30° virtual sound stage (a) is converted to a virtual sound stage with the width of the aperture of a loudspeaker array (b). Additionally, the lateral independent sound is played from the sides with separate loudspeakers for a stronger listener envelopment. It is expected that this results in a stronger impression of listener envelopment. In this case, the output signals are also computed by (25), where the signals with index  1  and M are the loudspeakers on the side. The loudspeaker pair selection, l and l+1, is in this case such that Ŝ′ is never given to the signals with index  1  and M since the whole width of the virtual stage is projected to only the front loudspeakers 2≦m≦M −1. 
         [0077]      FIG. 15  shows an example for the eight signals generated for the setup shown in  FIG. 14  for the same music clip for which the spatial decomposition was shown in  FIG. 10 . Note that the dominant singer in the center is amplitude panned between the center two loudspeaker signals, y 4  and y 5 . 
       Conventional 5.1 Surround Loudspeaker Setup 
       [0078]    One possibility to convert a stereo signal to a 5.1 surround compatible multi-channel audio signal is to use a setup as shown in  FIG. 14(   b ) with three front loudspeakers and two rear loudspeakers arranged as specified in the 5.1 standard. In this case, the rear loudspeakers emit the independent lateral sound, while the front loudspeakers are used to reproduce the virtual sound stage. Informal listening indicates that when playing back audio signals as described listener envelopment is more pronounced compared to stereo playback. 
         [0079]    Another possibility to convert a stereo signal to a 5.1 surround compatible signal is to use a setup as shown in  FIG. 11  where the loudspeakers are rearranged to match a 5.1 configuration. In this case, the ±30° virtual stage is extended to a ±110° virtual stage surrounding the listener. 
       Wavefield Synthesis Playback System 
       [0080]    First, signals y 1 , y 2 , . . . y M  are generated similar as for a setup as is illustrated in  FIG. 14(   b ). Then, for each signal, y 1 , y 2 , . . . y M , a virtual source is defined in the wavefield synthesis system. The lateral independent sound, y 1  and y M , is emitted as plane waves or sources in the far field as is illustrated in  FIG. 16  for M=8. For each other signal, a virtual source is defined with a location as desired. In the example shown in  FIG. 16 , the distance is varied for the different sources and some of the sources are defined to be in the front of the sound emitting array, i.e. the virtual sound stage can be defined with an individual distance for each defined direction. 
         [0000]    Generalized Scheme for 2-to-M Conversion 
         [0081]    Generally speaking, the loudspeaker signals for any of the described schemes can be formulated as: 
         [0000]      Y=MN  (29) 
         [0000]    where N is a vector containing the signals {circumflex over (N)}′ 1 , {circumflex over (N)}′ 2 , and Ŝ′. The vector Y contains all the loudspeaker signals. The matrix M has elements such that the loudspeaker signals in vector Y will be the same as computed by (25) or (27). Alternatively, different matrices M may be implemented using filtering and/or different amplitude panning laws (e.g. panning of Ŝ′ using more than two loudspeakers). For wavefield synthesis systems, the vector Y may contain all loudspeaker signals of the system (usually &gt;M). In this case, the matrix M also contains delays, all-pass filters, and filters in general to implement emission of the wavefield corresponding to the virtual sources associated to {circumflex over (N)}′ 1 , {circumflex over (N)} 2  and Ŝ′. In the claims, a relation like (29) having delays, all-pass filters, and/or filters in general as matrix elements of M is denoted a linear combination of the elements in N. 
       Modifying the Decomposed Audio Signals 
     Controlling the Width of the Sound Stage 
       [0082]    By modifying the estimated direction factors, e.g. A(i,k), one can control the width of the virtual sound stage. By linear scaling of the direction factors with a factor larger than one, the instruments being part of the sound stage are moved more to the side. The opposite can be achieved by scaling with a factor smaller than one. Alternatively, one can modify the amplitude panning law ( 20 ) for computing the angle of the localized direct sound. 
       Modifying the Ratio Between Localized Direct Sound and the Independent Sound 
       [0083]    For controlling the amount of ambience one can scale the independent lateral sound signals {circumflex over (N)}′ 1  and {circumflex over (N)}′ 2  for getting more or less ambience. Similarly, the localized direct sound can be modified in strength by means of scaling the S′ signals. 
       Modifying Stereo Signals 
       [0084]    One can also use the proposed decomposition for modifying stereo signals without increasing the number of channels. The aim here is solely to modify either the width of the virtual sound stage or the ratio between localized direct sound and the independent sound. The subbands for the stereo output are in this case 
         [0000]        Y   1   =v   1   {circumflex over (N)}′   1   +v   2   Ŝ′ Y   2   =v   1   {circumflex over (N)}′   2   +v   2   v   3   AŜ′   (30) 
         [0000]    where the factors v 1  and v 2  are used to control the ratio between independent sound and localized sound. For v 3  ≠1 also the width of the sound stage is modified (whereas in this case v 2  is modified to compensate the level change in the localized sound for v 3  ≠1).
 
Generalization to More than Two Input Channels
 
         [0085]    Formulated in words, the generation of {circumflex over (N)}′ 1 , {circumflex over (N)}′ 2  and Ŝ′ for the two-input-channel case is as follows (this was the aim of the least squares estimation). The lateral independent sound {circumflex over (N)}′ 1  is computed by removing from X 1  the signal component that is also contained in X 2 . Similarly, {circumflex over (N)}′ 2  is computed by removing from X 1  the signal component that is also contained in X 1 . The localized direct sound Ŝ′ is computed such that it contains the signal component present in both, X 1  and X 2 , and A is the computed magnitude ratio with which S′ is contained in X 1  and X 2 . A represents the direction of the localized direct sound. 
         [0086]    As an example, now a scheme with four input channels is described. Suppose a quadraphonic system with loudspeaker signals x 1  to x 4 , as illustrated in  FIG. 17(   a ), is supposed to be extended with more playback channels, as illustrated in  FIG. 17(   b ). Similar as in the two-input-channel case, independent sound channels are computed. In this case these are four (or if desired less) signals {circumflex over (N)}′ 1 , {circumflex over (N)}′ 2 , {circumflex over (N)}′ 3 , and {circumflex over (N)}′ 4 . These signals are computed in the same spirit as described above for the two-input-channel case. That is, the independent sound {circumflex over (N)}′ 1  is computed by removing from X 1  the signal components that are either also contained in X 2  or X 4  (the signals of the adjacent quadraphony loudspeakers). Similarly, {circumflex over (N)}′ 2 , {circumflex over (N)}′ 3 , and {circumflex over (N)}′ 4  are computed. Localized direct sound is computed for each channel pair of adjacent loudspeakers, i.e. Ŝ′ 12 , Ŝ′ 23 , Ŝ′ 34 , and Ŝ′ 41 . The localized direct sound Ŝ′ 12  is computed such that it contains the signal component present in both, X 1  and X 2 , and A 12  is the computed magnitude ratio with which Ŝ′ 12  is contained in X 1  and X 2 . A 12  represents the direction of the localized direct sound. With similar reasoning, Ŝ′ 23 , Ŝ′ 34 , Ŝ′ 41 , A 23 , A 34  and A 41  are computed. For playback over the system with twelve channels, shown in  FIG. 17(   b ), {circumflex over (N)}′ 1 , {circumflex over (N)}′ 2 , {circumflex over (N)}′ 3 , and {circumflex over (N)}′ 4  are emitted from the loudspeakers with signals y 1 , y 4 , y 7  and y 12 . To the front loudspeakers, y 1  to y 4 , a similar algorithm is applied as for the two-input-channel case for emitting Ŝ′ 12 , i.e. amplitude panning of Ŝ′ 12  over the loudspeaker pair most close to the direction defined by A 12 . Similarly, Ŝ′ 23 , Ŝ′ 34 , Ŝ′ 41 , are emitted from the loudspeaker arrays directed to the three other sides as a function of A 23 , A 34  and A 41 . Alternatively, as in the two-input-channel case, the independent sound channels may be emitted as plane waves. Also playback over wavefield synthesis systems with loudspeaker arrays around the listener is possible by defining for each loudspeaker in  FIG. 17(   b ) a virtual source, similar in spirit of using wavefield synthesis for the two-input-channel case. Again, this scheme can be generalized, similar to (29), where in this case the vector N contains the subband signals of all computed independent and localized sound channels. 
         [0087]    With similar reasoning, a 5.1 multi-channel surround audio system can be extended for playback with more than five main loudspeakers. However, the center channel needs special care, since often content is produced where amplitude panning between left front and right front is applied (without center). Sometimes amplitude panning is also applied between front left and center, and front right and center, or simultaneously between all three channels. This is different compared to the previously described quadraphony example, where we have used a signal model assuming that there are common signal components only between adjacent loudspeaker pairs. Either one takes this into consideration to compute the localized direct sound accordingly, or, a simpler solution is to downmix the front three channels to two channels and applying afterward the system described for quadraphony. 
         [0088]    A simpler solution for extending the scheme with two input channels for more input channels, is to apply the scheme for two input channels heuristically between certain channels pairs and then combining the resulting decompositions to compute, in the quadraphonic case for example, {circumflex over (N)}′ 1 , {circumflex over (N)}′ 2 , {circumflex over (N)}′ 3 , {circumflex over (N)} 4 , Ŝ′ 12 , Ŝ′ 23 , Ŝ′ 34 , Ŝ′ 41 , A 12 , A 23 , A 34  and A 41 . Playback of these is done as described for the quadraphonic case. 
       Computation of Loudspeaker Signals for Ambisonics 
       [0089]    The Ambisonic system is a surround audio system featuring signals which are independent of the specific playback setup. A first order Ambisonic system features the following signals which are defined relative to a specific point P in space: 
         [0000]      W=S 
         [0000]      X=S cos Ψcos Φ 
         [0000]      Y=S sin Ψcos Φ 
         [0000]      Z=S sin Φ 
         [0000]    where W=S is the (omnidirectional) sound pressure signal in P. The signals X, Y and Z are the signals obtained from dipoles in P, i.e. these signals are proportional to the particle velocity in Cartesian coordinate directions x, y and z (where the origin is in point P). The angles Ψ and Φ denote the azimuth and elevation angles, respectively (spherical polar coordinates). The so-called “B-Format” signal additionally features a factor of √{square root over (2)} for W X, Y and Z. 
         [0090]    To generate M signals, for playback over an M-channel three dimensional loudspeaker system, signals are computed representing sound arriving from the eight directions x, −x, y, −y, z, −z. This is done by combining W X, Y and Z to get directional (e.g. cardioid) responses, e.g. 
         [0000]    
       
      
       x 
       1 
       =W+X x 
       3 
       =W+Y x 
       5 
       =W+Z  
      
     
         [0000]        x   2   =W−X x   4   =W−Y x   6   =W−Z   (31) 
         [0091]    Given these signals, similar reasoning as described for the quadraphonic system above is used to compute eight independent sound subband signals (or less if desired) {circumflex over (N)}′ c  (1≦c≦8). For example, the independent sound {circumflex over (N)} 40   1  is computed by removing from X 1  the signal components that are either also contained in the spatially adjacent channels X 3 , X 4 , X 5  or X 6 . Additionally, between adjacent pairs or triples of the input signals localized direct sound and direction factors representing its direction are computed. Given this decomposition, the sound is emitted over the loudspeakers, similarly as described in the previous example of quadraphony, or in general (29). 
         [0092]    For a two dimensional Ambisonics system, 
         [0000]      W=S 
         [0000]      X=S cos Ψ 
         [0000]      Y=S sin Ψ  (33) 
         [0000]    resulting in four input signals, x 1  to x 4 , the processing is similar to the described quadraphonic system. 
       Decoding of Matrixed Surround 
       [0093]    A matrix surround encoder mixes a multi-channel audio signal (for example 5.1 surround signal) down to a stereo signal. This format of representing multi-channel audio signals is denoted “matrixed surround”. For example, the channels of a 5.1 surround signals may be downmixed by a matrix encoder in the following way (for simplicity we are ignoring the low frequency effects channel): 
         [0000]    
       
         
           
             
               
                 x 
                 1 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 l 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               + 
               
                 
                   1 
                   
                     2 
                   
                 
                  
                 
                   c 
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
               
               + 
               
                 j 
                  
                 
                     
                 
                  
                 
                   1 
                   
                     2 
                   
                 
                  
                 
                   
                     l 
                     s 
                   
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
               
               + 
               
                 j 
                  
                 
                     
                 
                  
                 
                   1 
                   
                     6 
                   
                 
                  
                 
                   
                     r 
                     s 
                   
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
               
             
           
         
       
       
         
           
             
               
                 x 
                 2 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 r 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               + 
               
                 
                   1 
                   
                     2 
                   
                 
                  
                 
                   c 
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
               
               - 
               
                 j 
                  
                 
                     
                 
                  
                 
                   1 
                   
                     2 
                   
                 
                  
                 
                   
                     r 
                     s 
                   
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
               
               - 
               
                 j 
                  
                 
                     
                 
                  
                 
                   1 
                   
                     16 
                   
                 
                  
                 
                   
                     l 
                     s 
                   
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
               
             
           
         
       
     
         [0000]    where l, r, c, l s , and r s  denote the front left, front right, center, rear left, and rear right channels respectively. The j denotes a 90 degree phase shift, and −j is a −90 degree phase shift. Other matrix encoders may use variations of the described downmix. 
         [0094]    Similar as previously described for the 2-to-M channel conversion, one may apply the spatial decomposition to the matrix surround downmix signal. Thus for each subband at each time independent sound subbands, localized sound subbands, and direction factors are computed. Linear combinations of the independent sound subbands and localized sound subbands are emitted from each loudspeaker of the surround system that is to emit the matrix decoded surround signal. 
         [0095]    Note that the normalized correlation is likely to also take negative values, due to the out-of-phase components in the matrixed surround downmix signal. If this is the case, the corresponding direction factors will be negative, indicating that the sound originated from a rear channel in the original multi-channel audio signal (before matrix downmix). 
         [0096]    This way of decoding matrixed surround is very appealing, since it has low complexity and at the same time a rich ambience is reproduced by the estimated independent sound subbands. There is no need for generating artificial ambience, which is very computationally complex. 
       Implementation Details 
       [0097]    For computing the subband signals, a Discrete (Fast) Fourier Transform (DFT) can be used. For reducing the number of bands, motivated by complexity reduction and better audio quality, the DFT bands can be combined such that each combined band has a frequency resolution motivated by the frequency resolution of the human auditory system. The described processing is then carried out for each combined subband. Alternatively, Quadrature Mirror Filter (QMF) banks or any other non-cascaded or cascaded filterbanks can be used. 
         [0098]    Two critical signal types are transients and stationary/tonal signals. For effectively addressing both, a filterbank may be used with an adaptive time-frequency resolution. Transients would be detected and the time resolution of the filterbank (or alternatively only of the processing) would be increased to effectively process the transients. Stationary/tonal signal components would also be detected and the time resolution of the filterbank and/or processing would be decreased for these types of signals. As a criterion for detecting stationary/tonal signal components one may use a “tonality measure”. 
         [0099]    Our implementation of the algorithm uses a Fast Fourier Transform (FFT). For 44.1 kHz sampling rate we use FFT sizes between 256 and 1024. Our combined subbands have a bandwidth which is approximately two times the critical bandwidth of the human auditory system. This results in using about 20 combined subbands for 44.1 kHz sampling rate. 
       Application Examples 
     Television Sets 
       [0100]    For playing back the audio of stereo-based audiovisual TV content, a center channel can be generated for getting the benefit of a “stabilized center” (e.g. movie dialog appears in the center of the screen for listeners at all locations). Alternatively, stereo audio can be converted to 5.1 surround if desired. 
       Stereo to Multi-Channel Conversion Box 
       [0101]    A conversion device would convert audio content to a format suitable for playback over more than two loudspeakers. For example, this box could be used with a stereo music player and connect to a 5.1 loudspeaker set. The user could have various options: stereo+center channel, 5.1 surround with front virtual stage and ambience, 5.1 surround with a ±110° virtual sound stage surrounding the listener, or all loudspeakers arranged in the front for a better/wider front virtual stage. 
         [0102]    Such a conversion box could feature a stereo analog line-in audio input and/or a digital SP-DIF audio input. The output would either be multi-channel line-out or alternatively digital audio out, e.g. SP-DIF. 
         [0000]    Devices and Appliances with Advanced Playback Capabilities 
         [0103]    Such devices and appliances would support advanced playback in terms of playing back stereo or multi-channel surround audio content with more loudspeakers than conventionally. Also, they could support conversion of stereo content to multi-channel surround content. 
       Multi-Channel Loudspeaker Sets 
       [0104]    A multi-channel loudspeaker set is envisioned with the capability of converting its audio input signal to a signal for each loudspeaker it features. 
       Automotive Audio 
       [0105]    Automotive audio is a challenging topic. Due to the listeners&#39; positions and due to the obstacles (seats, bodies of various listeners) and limitations for loudspeaker placement it is difficult to play back stereo or multi-channel audio signals such that they reproduce a good virtual sound stage. The proposed algorithm can be used for computing signals for loudspeakers placed at specific positions such that the virtual sound stage is improved for the listener that are not in the sweet spot. 
       Additional Field of Use 
       [0106]    A perceptually motivated spatial decomposition for stereo and multi-channel audio signals was described. In a number of subbands and as a function of time, lateral independent sound and localized sound and its specific angle (or level difference) are estimated. Given an assumed signal model, the least squares estimates of these signals are computed. 
         [0107]    Furthermore, it was described how the decomposed stereo signals can be played back over multiple loudspeakers, loudspeaker arrays, and wavefield synthesis systems. Also it was described how the proposed spatial decomposition is applied for “decoding” the Ambisonics signal format for multi-channel loudspeaker playback. Also it was outlined how the described principles are applied for microphone signals, ambisonics B-format signals, and matrixed surround signals.
A perceptually motivated spatial decomposition for two-channel stereo audio signals, capturing the information about the virtual sound stage, is proposed. The spatial decomposition allows to re-synthesize audio signals for playback over other sound systems than two-channel stereo. With the use of more front loudspeakers, the width of the virtual sound stage can be increased beyond +/-30deg and the sweet spot region is extended. Optionally, lateral independent sound components can be played back separately over loudspeakers on the two sides of a listener to increase listener envelopment. It is also explained how the spatial decomposition can be used with surround sound and wavefield synthesis based audio system. According to the main embodiment of the invention applying to multiple audio signals, it is proposed to generate multiple output audio signals (y 1... Y M) from multiple input audio signals (x 1,..., x L), in which the number of output is equal or higher than the number of input signals, this method comprising the steps of: --by means of linear combinations of the input subbands X 1 (i),..., X L (i), computing one or more independent sound subbands representing signal components which are independent between the input subbands, --by means of linear combinations of the input subbands X 1 (i),..., X L (i), computing one or more localized direct sound subbands representing signal components which are contained in more than one of the input subbands and direction factors representing the ratios with which these signal components are contained in two or more input subbands, --generating the output subband signals, Y 1 (i)... Y M (i), where each output subband signal is a linear combination of the independent sound subbands and the localized direct sound subbands--converting the output subband signals, Y 1 (i)... Y M (i), to time domain audio signals, y 1... Y M.