Patent Publication Number: US-4149031-A

Title: Multichannel matrix logic and encoding systems

Description:
This application is a division of application Ser. No. 701,228 filed June 30, 1976 (now U.S. Pat. No. 4,085,291), and is a continuation-in-part of my application Ser. No. 468,238 filed May 9, 1974 (now U.S. Pat. No. 3,985,978), Ser. No. 578,078 filed May 16, 1975 (now U.S. Pat. No. 3,970,788), Ser. No. 288,873 filed Sept. 13, 1972 (now U.S. Pat. No. 3,906,156), and Ser. No. 187,065 filed Oct. 6, 1971 (now U.S. Pat. No. 3,856,992), which are hereby incorporated herein by reference. Various subject matter described herein is disclosed in Patent Office Disclosure Documents 035254, 035255, and 035256 filed Sept. 10, 1974. 
    
    
     The present invention relates to apparatus and methods for preparation, recording or transmission, and production of multidirectional audio signals in a manner to provide an enhanced perception of those directional qualities and realism in comparison to conventional stereo practice. More particularly, the present invention is directed to multichannel encoding and decoding of multidirectional audio signals, and to encoding variations and logic systems in connection therewith. 
     Multidirectional audio recording and reproduction systems are often designated as quadrasonic systems, because four loudspeakers in a square array centered on a listener are often used. In such quadrasonic systems, the reproduction loudspeaker locations are conventionally designated LF, RF, LB, and RB, for the left front, right front, left back and right back loudspeaker directions, respectively, with respect to their bearing angles referred to that listener. In one method of preparing audio signals for reproduction through such loudspeakers, four directional microphones are used to generate directional source signals which may be designated S LF , S RF , S LB  and S RB  for source signals from the left front, right front, left back, and right back reference directions. Linear combinations of these signals may be encoded into four transmission signals for broadcast and/or recording purposes. The encoded transmission signals may be broadcast or recorded on independent channels or tracks, or may use subcarrier modulation techniques for broadcasting or disc recording. It is the function of the encoding apparatus to prepare such linear combinations as may be particularly suited to the specific recording or transmission means to be used, and many specific linear combinations have been devised for the encoding of the source signals in the transmission signals. 
     It is frequently desired that only one or more of the four transmission signals be usable for a more limited, but balanced portrayal of the available directional information detected by the source microphones, namely, that the use of some of these transmission signal channels be compatible with existing audio systems using fewer channels. For example, in FM broadcast applications it is desired that one such transmission channel be suitable for monaural reproduction. In this connection, this monaural-compatible transmission signal, which may be designated &#34;M&#34; should desirably consist of the linear sum combination of the source signals, as follows: 
     
         M= S.sub.LF + S.sub.LB + S.sub.RB                          ( 1) 
    
     similarly, it is desired that one further channel serve as a left-right compatible stereo-difference channel. In this connection, the stereo difference channel, which may be designated &#34;Y&#34; may consist of the linear, left-right difference combination of the directional source signals, as follows: 
     
         Y= S.sub.LF - S.sub.RF + S.sub.LB - S.sub.RB               ( 2) 
    
     for modulation of a 38-kHz subcarrier in the usual manner for conventional stereo broadcast which is suitable for compatible reception by a conventional stereo receiver. 
     The remaining transmission channels may be regarded as supplementary to the mono- and stereo-difference channels, and the necessary independent linear combinations for these may be arbitrarily selected, or else devised to suit further compatibility requirements. For example, one such transmission signal may be S LB  with null contribution from the other source signals, while the other may be S RB  with null contributions from the other source signals. After study, it is seen that a suitable decoder, forming linear combinations of M, Y, S LB , S RB  to form each of the source signals separately, is possible. On the other hand, if the four source signals be sequentially sampled, then the encoding and subcarrier modulation may be done as one step. The M and Y transmission channels then remain as above, and the supplementary channels appear as a front-back difference channel (designated &#34;X&#34;): 
     
         x= s.sub.lf + s.sub.rf -(s.sub.lb + s.sub.rb),             (3) 
    
     as a modulation on a quadrature-phase 38-kHz subcarrier, and a diagonal-difference channel (designated &#34;U&#34;): 
     
         u= s.sub.lf -(s.sub.rf + s.sub.lb)+s.sub.rb                ( 4) 
    
     appearing as modulation on a 76-kHz subcarrier. This method of forming linear combinations, or matrixing, as it is called, has the further compatibility advantage in that, if the U channel not be transmitted or received, i.e., be silent, then the decoding means already devised for M, X, Y, and U still provides a satisfactory decode for four-loudspeaker presentation, in that the signal intended for a given loudspeaker exhibits a suppression of the appearance of the other source signals by about 9.5 decibels, a suppression nearly as satisfactory as the practical values, e.g., 20 decibels or better, often obtained if U be not omitted. The resulting matrix is known as a 4-3-4 matrix, to distinguish it from the U-included, or 4-4-4 case. 
     Other 4-4-4 matrixing methods have been employed for disc recording. In the CD-4 system, for example, the stereo-compatible recording channels are shown below as L and R: 
     
         l= s.sub.lf + s.sub.lb                                     (6) 
    
     
         l&#39;= s.sub.lf - s.sub.lb                                    ( 7) 
    
     
         r= s.sub.rf + s.sub.rb                                     ( 8) 
    
     
         r&#39;= s.sub.rf - s.sub.rb                                    ( 9) 
    
     while L&#39; and R&#39; are supplementary channels appearing as modulations of ultrasonic subcarriers on the same groovewalls of the disc. 
     It is characteristic of these conventional systems that they produce stereo-compatible signals that contain only left-right information, and omit any further information allowing for the differentiation of front from back, if only the first two transmission channels be used. Thus, the four-loudspeaker decode from two channels is directionally defective, i.e., an adequate 4-2-4 matrix is not contained in such systems. Such 4-2-4 matrices, of varying degrees of adequacy, however, have been devised specifically for applications in which the transmission or recording of more than two channels would be prohibitively difficult. The well-known QS and SQ matrices are examples of such 4-2-4 systems. 
     In the previously referred to patents and applications of the present inventor, there are described multichannel matrices, here designated by the generic name UMX, that include specific 4-2-4, 4-3-4, and 4-4-4 cases merely by deletion of certain of the transmission channels. These provide for the encoding of directional information into transmission channels which may be in accordance with the following embodiments of specific encode equations: ##EQU1## in which (jθ) designates a leading phase-shift of all frequency components of the audio signal in the amount designated by the phase angle θ. This phase shift is in excess of any frequency-dependent reference phase applied in common to all signals in forming T.sub.Σ. These relative phase shifts θ i  may be made equal to the bearing angles describing the intended bearing-angle locations for each of the source signals S i , and these bearing angles may be measured counterclockwise from the center right (CR or 0°) bearing, the midbearing between the right front (RF) and the right back (RB) directions. 
     The decoding of these signals to provide signals for presentation via loudspeakers, for example, at specific bearing angles φ (as shown in FIG. 2) may be carried out in accordance with the following equation: 
     
         P.sub.φ = T.sub.Σ + T.sub.Δ exp(Jφ)+T.sub.T exp(-jφ)+T.sub.Q  exp( 2jφ)                       (14) 
    
     The resulting overall encode-decode signals obey the equations: ##EQU2## 
     The characteristic kernel of the above equation: ##EQU3## fully characterizes the directional properties of the overll matrix. It is a single-valued, 360°-repetitive function of the single variable α=φ-θ. These directional properties are shown in the polar plot of the magnitude of this kernel as labelled &#34;QMX&#34; in FIG. 3 of the drawings. 
     If the fourth transmission signal T Q  be made null, the overall encode-decode signals obey the equation: ##EQU4## with the characteristic kernel ##STR1## shown plotted in its magnitude as the polar plot labelled &#34;TXM&#34; in FIG. 3 of the drawings. This kernel is a single-valued, 360°-repetitive function of the single variable α=φ-θ. Also, if both transmission channel T Q  and transmission channel T T  be made null, the overall encode-decode signals obey ##EQU5## with the characteristic kernel ##STR2## also shown plotted in its magnitude as the polar plot labelled &#34;BMX&#34; in FIG. 3. This kernel is also a single-valued, 360°-repetitive function of the single variable α=φ-θ. 
     These embodiments of kernel functions have a magnitude which is symmetric about the null difference angle α=φ-θ, namely, about α=0°, whereas the phase of the kernel functions is antisymmetric about α=0°, positive phases for α&gt;0° being matched by negative phases for α&gt;0°. These relations are called axial symmetry. The illustrated TMX kernel also obeys this phase rule albeit vacuously, the phases being zero except for the values of a beyond ±120° where 180° phase is obtained. In all cases, it is noted that a relative maximum is obtained at α=0° and a relatively small value at α=180°. 
     These kernel functions also exhibit rotational symmetry in that, regarded as functions of φ, the orientation of the plot depends upon the value of θ. This is illustrated in FIG. 4 for the QMX kernel. In this illustration of FIG. 4, the four radial lines may be regarded as indicating loudspeaker placements, and the intercept of these lines with the polar curve may be regarded as indicating, by the distance from the center to the intercept, the magnitude of the signal supplied to the loudspeaker from a given source. The dashed line indicates the axis of symmetry of the kernel, and this axis points to the source location for a single source. Eight such source locations are shown, and the rigid body rotation of the kernel to follow the souce location (rotational symmetry) is evident. 
     Because the optimal shape by psycho-acoustic criteria for the kernel function is not precisely known, it may be useful to provide for alternate decoding means by introducing variable coefficients according to the equation: 
     
         P.sub.φ =aT.sub.93 +bT.sub.Δ exp(jφ )+cT.sub.T exp(- jφ)+dT.sub.Q exp( 2 jφ)                           (21) 
    
     so that an overall encode-decode kernel function would be proportional to: 
     
         a+b[exp j(φ-θ)]+c[exp j(θ-φ )]+d[exp2j(φ-θ )](22) 
    
     Examples of the effects of various adjustments of the coefficients a, b, c, and d are discussed in previously referred to patent applications. FIG. 5 schematically shows the introduction of such coefficients. Further discussion is undertaken in the present application. 
     In the previously referred to applications, linear combinations of the T 93  and T 66  transmission signals are disclosed which may be used as the left and right channels for disc recording. These signals, which are designated T L  for the left channel signal and T R  for the right channel signal, may be formed as follows: ##EQU6## The T L  and T R  signals correspond to BMX (two transmission channel) decodes for center left (φ=180°) and center right (φ=0°) and thus form stereo-compatible channels, while being also decodable to yield quadrasonic presentation obeying the BMX kernel. The remaining channels T T  and T Q  then serve the role of supplementary channels as in the examples of FM broadcasting and disc recording discussed above. There remains the improvement, however, in that the stereo-compatible channels are decodable for quadrasonic presentation. 
     When these UMX signals are to be recorded on a disc, the channels T L  and T R  may be recorded in the usual manner together with high-frequency carriers superposed in each of the groovewall modulations, the carriers being phase modulated by frequency emphasized versions of ##EQU7## wherein C L  is the left groovewall carrier channel and C R  is the right groovewall carrier channel in one embodiment. With due regard for the sum-difference combinations, T L  and T R  may be decoded to provide signals obeying the BMX kernel plotted in FIG. 6, whereas the decoding of C L  and C R  provides signals obeying the CMX kernel ##EQU8## shown plotted also in FIG. 6. Half the sum of the two comprises the QMX kernel of FIG. 6. 
     Since the BMX information and the CMX information are transmitted by different means, it is not always possible to ensure their combination in exactly equal proportion, so that the &#34;a&#34; and &#34;b&#34; coefficients of Equation (21) may not be exactly equal to the &#34;c&#34; and &#34;d&#34; coefficients, resulting in a matrix imbalance. As an illustration, the effect of such imbalances which do not disturb axial or rotational symmetry has been calculated and plotted in FIG. 7. In the plots at the left of FIG. 7, the effects of reducing the CMX contribution by varying decibel amounts are shown in the plots, with &#34;backlobes&#34; omitted in the lower portion for the sake of clarity. The dashed curve shows the effect in the BMX kernel along caused by a 3-dB imbalance. 
     The effects of imbalance together with phase shift are also shown in the plots at the right of FIG. 7, the imbalance being related to the phase as in the case of a first-order filter, 3-dB being plotted together with a 45° phase shift, for example. The psychoacoustic consequences of the imbalance shown here or discussed in the copending application remain largely unknown. It is known, however, that more precise localizations obtain with QMX rather than BMX kernels, and that a 10-dB reduction in the CMX contribution has rather little psychoacoustic effect, as does the elimination of TQ altogether. 
     In U.S. Pat. No. 3,906,156 it is observed that the CMX channels may be limited in bandwidth without significant psychoacoustic degradation. If this be done by means of first-order, low-pass filters without phase compensation, then the kernel will progressively change with frequency as in the right-most plot of FIG. 7. In Equation (21), the effect is that the coefficients &#34;c&#34; and &#34;d&#34; become frequency-dependent in magnitude and phase. However, if a matching phase characteristic, but without amplitude variation, be inserted for the coefficients &#34;a&#34; and &#34;b&#34;, then the kernel will progressively change with frequency as in the left-most plot of FIG. 7. 
     Evidently, the smoothly changing kernel, without loss of axial or rotational symmetry, helps minimize the psychoacoustic consequences of the bandwidth reduction for the CMX contribution. However, the full explanation of the psychoacoustic acceptability must also rest upon the properties of human hearing being different at different frequencies, with regard to the perception of directional effects. 
     That the perception of directional effects depends upon frequency is widely believed. For example, at low-frequencies, it is believed that the perception depends upon phase differences in the acoustic signals present at the two ears, while at high frequencies the dependence is mainly upon amplitude differences. However that may be, no precise, reliable mathematical formulation is available. The understanding of these phenomena is further complicated by the fact that the more presence of a listener&#39;s head in the acoustic sound field produces a distortion of the field, so that the relation of the acoustic signals at the ears to the original acoustic-field signals at those same points is an extremely complicated relation. At high frequencies (i.e., wavelengths comparable to, or smaller than typical head dimensions) exceedingly complicated phase and amplitude relations, for example, are the rule. The complexity of these relations may make it credible that the ordinary observation of acute directional perception of single point-source images at high frequencies does not necessarily extend to phantom-source imaging--an imaging possibility that is not prescribed in nature. 
     Even at moderate frequencies, phantom-source imaging is poorly understood. In two-speaker stereo practice, for example, a centered phantom image is usually obtained by providing the two loudspeakers with identical signals. If, however, a phase shift of 180° be introduced between the two signals, most listeners find the image to be disturbed--either unstable in location or else appearing to be greatly spread out--and the localization experience to be productive of a certain degree of auditory distress. With a 90° phase shift, on the other hand, many listeners find no change in image quality, in comparison to 0°, while others require the phase shift to be reduced to 45°, or less, before they report no change. Even at 0° of phase, many listeners find the phantom image to exhibit instability and broadening in comparison to single point-source images. The phase shifts at the listener&#39;s ears, and the amplitude relations as well, bear no simple relation to the loudspeaker-signal phases, of course. 
     In quadrasonic imaging, the introduction of two additional loudspeakers provides the opportunity for more elaborate structuring of the signals presented, via the intervening acoustic space, to the two ears. For example, in seeking to establish a center-front phantom image, the equal excitation of the two front loudspeakers, previously known to be quite effective from two speaker stereo practice, may be supplemented by suitably phased, low amplitude excitation of the back loudspeakers. For example, if a TMX kernel be used, the two back loudspeakers may be supplied with equal signals, as are the two front loudspeakers, but the back signals may be lower in level by about 15.8-dB and show a 180° phasing relative to the front. Experiments have shown that the presence of these back signals stabilizes the center-front phantom image, to the extent of reducing localization disagreements among various listeners by about 50%. 
     Similar experiments have been performed in which the back signals are phased 135° relative to one another at a level of 7.7-dB below the front level, but with 45° phasing relative to the front (BMX kernel) or with 135° phasing relative to the front (QMX kernel). The stabilization of center-front phantom-image localization is similar in these two cases to that observed in the TMX kernel case cited above. Evidently, the phasing of the back signals in these three cases serves to &#34;repel&#34; the image from the back location to stabilize imaging at the front location. Repulsion of the image away from the lower-level loudspeaker in 180° phasing of stereo-pair loudspeaker presentations is a well-known phenomenon, particularly at frequencies such that the acoustic wavelengths are not small compared to head dimensions. 
     The phase and amplitude differences are slight at extremely low frequencies between acoustic signals at the two ears, as caused by simple sources. Directional perception is also not quite so acute as at the middle frequencies, although human hearing may be more sensitive to these phase differences that do remain. The importance and phasing of kernel backlobes for phantom imaging may well be different at such low frequencies, e.g., below about 100 Hz, than appear to be the case at middle frequencies, say 100 Hz to 2000 Hz. For these reasons, it would be useful to provide decode apparatus in which a set of coefficients &#34;a&#34;, &#34;b&#34;, &#34;c&#34;, and &#34;d&#34; of Equation (21) may be provided or selected for signal frequency values below some low frequency such as about 100 Hz which is different from the corresponding set of coefficients used above that frequency. 
     In some signal transmission and recording systems, it is desired to avoid the technological complications and costs of providing more than two transmission channels. In such cases, it is then desired to find a means of enhancing the psychoacoustic quality of the images that may be obtained from the presentation signals of the 4-2-4 matrices that may be employed. This enhancing means has been called &#34;logic&#34;, a means of detecting signal conditions corresponding to a single image and either adjusting the matrix coefficients or adjusting the gains in the presentation channels to reduce the crosstalk into undesired channels. The provision of improvements and advances in the area of &#34;logic&#34; or image location enhancement in multidirectional sound systems would be desirable in respect of providing the performance enhancement for multidirectional sound systems appropriate to such techniques. 
     As indicated, the previously referred to patent applications described generally an encoding matrix system known as the UMX system, and include description of various specific matrix embodiments of the system. However, improvements and variations in such matrix systems for particular purposes would be of interest, such as matrix variations relating to particular aspects of stereo compatibility in tape and disk recording. 
     Accordingly, it is an object of the present invention to provide improvements in multidirectional signal matrixing systems. It is a further object to provide for image location enhancement in multidirectional sound reproducing systems. Another object is the provision of matrix embodiments relating to particular aspects of stereo compatibility in tape and disk recording. 
    
    
     These and other objects will be apparent from the following detailed description and the accompanying drawings of which: 
     FIG. 1 is a general schematic illustration of one type of quadrasonic sound system; 
     FIG. 2 is a schematic view illustrating certain angular relations useful in matrix encoding and decoding; 
     FIG. 3 is a series of three polar plots of the magnitude of the characteristic presentation kernels of certain embodiments of four, three and two transmission channel matrix systems; 
     FIG. 4 is an illustration of the rotational symmetry of the specific four transmission channel presentation kernel function of FIG. 3; 
     FIG. 5 is a block diagram of an embodiment of a 4-4-4 matrixing system; 
     FIG. 6 is a series of three polar plots of the magnitude of several characteristic presentation kernel embodiments; 
     FIG. 7 is a polar plot illustrating the effects of matrix imbalance in certain matrix embodiments; 
     FIG. 8 is a schematic illustration of magnetic tape tracks for multichannel recording; 
     FIG. 9 is a polar plot of the magnitude of a matrix embodiment of a &#34;right front&#34; signal including a matrix-supplement signal; 
     FIG. 10 is schematic diagram of signal synthesis circuitry useful in image enhancement; 
     FIG. 11 is a schematic diagram of circuitry for generating an image enhancement control signal; 
     FIG. 12 is an overall schematic diagram of an embodiment of a matrix logic circuit; 
     FIG. 13 is a schematic diagram of a processing block relating to the circuit of FIG. 12; 
     FIG. 14 is a schematic diagram of a processing block relating to the circuit of FIG. 13; 
     FIG. 15 is a schematic diagram of a processing block relating to the circuit of FIG. 12; 
     FIGS. 16 and 17 are schematic diagrams of envelope generation circuitry adapted for use in connection with the circuit of FIG. 12, and 
     FIGS. 18 and 21 are schematic diagrams of circuitry for processing embodiments of stereo compatible multidirectional signals. 
    
    
     The present invention concerns multidirectional sound systems, including encoding and decoding methods and apparatus, in which multidirectional source signals are transmitted, in matrixed form, in a plurality of channels including two basic transmission channels. In such systems, the multidirectional source signals are generally encoded in the two basic channels with predetermined amplitude and phase relations indicative of source directions, and the encoded channels are adapted to be decoded for formation from the signals in said two channels of more than two loudspeaker presentation signals having source signals corresponding to every source direction appearing in a plurality of such presentation signals but in differing relative amplitude. Advantageously, the two basic transmission channels are monaural and stereo compatible channels, and accordingly may be utilized directly as &#34;left&#34; and &#34;right&#34; presentation signals in accordance with conventional stereo practice and through the use of conventional stereophonic equipment. In this connection, the term &#34;transmitted&#34; as used herein will be understood to include the various forms of recording or signal-storage, such as phonograph records and tapes, and a principal application of the various aspects of the present invention relates directly to the decoding of signals transmitted from the groovewalls of a phonograph recording disk. Such recording disks may also be monaural and stereo-compatible, such that they may be utilized directly with conventional monaural record players or stereophonic playback equipment. 
     As indicated previously, aspects of the present disclosure are directed to presentation image location enhancement, and this aspect of the present disclosure is particularly applicable to multidirectional sound systems in which the source signals are transmitted to the decoder in at least two original transmission channels. In accordance with this aspect of the present invention, at least one synthetic supplementary channel is generated from the original encoded transmission channels, with the source signals in the synthetic supplementary channels having predetermined amplitude and phase relations differently indicative of source directions relative to the amplitude and phase relations of the multidirectional source signals encoded in the original transmission channels. 
     The amplitude and phase relations of the source signals in the synthetic channels are determined by detection of single-image conditions from the signals of the original transmission channels, and determining synthetic re-encode amplitude and phase coefficients differently indicative of source directions in respect of that image location. 
     The signals of the original channels and the synthetic channels are then decoded to form multidirectional presentation signals having image location which is enhanced with respect to the image location which would be provided by decoding the original channels without the synthetic channels. Generally, this decoding of the original transmission channels and the synthetic channels may be carried out by decoding the original channels to provide a first set of more than two loudspeaker presenation signals for presentation at more than two listening space bearing angles and having source signals corresponding to every source direction appearing in a plurality of such presenation signals but in different relative amplitude, and decoding the synthetic channel or channels to provide a second set of a plurality of presentation signals for presentation at listening space bearing angles corresponding to the listening space bearing angle presentation directions of said first set of presentation signals. The two sets of presentation signals are combined by adding to each signal of the first set of presentation signals the signal of said second set of presentation signals for the corresponding listening space bearing angle to provide a plurality of output signals. The output signals are for the presentation directions corresponding to the presentation directions of the first set of presentation signals, but with a sharpened directionality pattern with respect to that of the set of presentation signals obtained from the original signals. The extent to which the synthetic supplementary channels are included in the decode will generally be determined by the degree to which the signals of the original transmission channels are interpreted to represent single image (saliency) conditions. 
     This synthetic approach to image enhancement differs from previous approaches in providing means for generting one or more synthetic channels from the original two channel quad input and quad presentation speaker matrix (4-2-4 matrix), and the insertion of these synthetic channels as supplementary channels into the decode matrix of the correspondingly increased transmission channel matrix, such as a 4-4-4 system where two synthetic supplementary channels are generated. 
     This aspect of the present discosure will be subsequently described in detail in respect of a two original channel matrix, with the understanding that the principles disclosed may be applied to the use of three or more channels. Since there exists no linear means of obtaining mutually-independent 4-channel signals from the original two channels of a two channel matrix, nonlinear means are provided for the detection of single-image conditions, and for determining synthetic re-encode coefficients for that image in a form suitable for representation in the synthetic supplementary channels. The extend to which these synthetic channels are allowed to contribute to the decode is determined by the degree to which the signal conditions in the two original channels may be interpreted as representing single-image conditions, for example, by adjusting the a, b, c, d coefficients of Equation (21) where c and d are the respective coefficients of two synthetically generated supplementary channels which may be reduced or nulled for non-salient conditions and made equal to the a and b coefficients of the two original channels for fully salient conditions. The speed with which the synthetic contribution is provided and the speed with which it is removed are called attack and release speeds corresponding to characteristic attack and release times. These three conditions, the saliency of single-image conditions, attack time, and release time, may be set to match characteristics of human hearing related to the masking of one image by another in psychoacoustics. Means for providing these functions are to be further described with respect to a particular embodiment of a two transmission channel system utilizing mono and stereo compatible BMX channels T L  and T R  in accordance with the UMX matrixing system. 
     The synthetic channel synthesis involves the generation of phase-shifted audio signals in which the amount of phase shift is based on an estimate of the phase shift used in encoding the salient sources in the original channels. For the case of a two transmission channel system in accordance with the UMX system, the original channels T L  and T R  obtained from the recording medium pickup device, being linear combinations of the T.sub.Σ and T.sub.Δ channels as indicated by Equations (23) and (24), may be readily rematrixed to provide the T.sub.Σ and T.sub.Δ channels ##EQU9## by a simple matrix circuit. For the purposes of the following discussion, the resulting T.sub.Σ and T.sub.Δ signals which also may be utilized directly (e.g. as received from an FM broadcast). and representable, as follows: 
     
         T.sub.Σ = Σ.sub.i S.sub.i (t)                  (28) 
    
     
         T.sub.Δ = Σ.sub.i [S.sub.i (t) cos θ.sub.i -S.sub.i (t) sin θ.sub.i ]                                       (29) 
    
     in which S i  (t) is the leading-phase, 90° phase-shifted version (Hilbert transform) of the ith audio signal (S i  (t), and the pair of terms shown in Equation (29) represents a phase shifted version of S i  (t) that lags that audio signal by the phase angle θ i . A synthetic supplementary audio signal channel T ET  may be generated having encoding coefficients differently indicative (from those of the original two channels), as follows: 
     
         T.sub. ET = Σ.sub.i [S.sub.i (t) cos θ.sub.E +S.sub.i (t) sin θ.sub.E ]                                           (30) 
    
     the pair of terms shown in Equation (30 ) forms a phase shifted version of S i  (t) with leading phase angle θ E , which is a phase-angle estimate. The phase-angle estimate θ E  may be formed by phase-lock-loop detection. It will be appreciated that signal T ET  may be regarded as a synthetically generated T T  signal in accordance with the UMX system for the estimated phase angle θ E  of the phase shift used in encoding the salient sources in the original channels. 
     A phase-lock loop detection system for determining θ E  comprises one oscillator that is fixed-tuned and one that is voltage-controlled in frequency. The audio signals T.sub.Σ and T.sub.Δ are used to impress single-sideband (SSB) modulation upon the oscillations, with a quadrature phase relation for the carriers, except for the further difference in phase represented by the difference between the estimated phase angle θ E  and the phase angle to which it is to correspond. Phase detection produces a signal proportional to the sine of this phase difference, and this sine signal may be used as an error signal to control adjustment of the frequency of the voltage-controlled oscillator to bring the error signal to a null. Since phase is the time integral of frequency, this phase-lock is phase-rate controlled. The nominal frequency ω o  for the oscillators may be chosen by convenience but should be more than twice the highest frequency in the audio signals. In this connection, it should be noted that the synthetic channels may have a limited frequency range as described and claimed in parent application Ser. No. 578,078 for the auxiliary channels, and such limited frequency range auxiliary supplementary channels are particularly suited for use with decoders adapted for use with originally transmitted, limited frequency range auxiliary channels. 
     The SSB signal formed from the original channel signal T.sub.Σ is: 
     
         σ.sub.E (t)=Σ.sub.i [S.sub.i (t) cos (ω.sub.o t+θ.sub.E)-S.sub.i (t) sin (ω.sub.o t+θ.sub.E)](31) 
    
     using the voltage-controlled oscillator. The SSB signal δ(t) formed from the original channel signal T.sub.Δ is representable as follows: 
     
         δ(t)=Σ.sub.j [S.sub.j (t) cos θ.sub.j sin ω.sub.o t-S.sub.j (t) sin θ.sub.j sin ω.sub.o t+S.sub.j (t) cos θ.sub.j cos ω.sub.o t+S.sub.j (t) sin θ.sub.j cos ω.sub.o t]=Σ.sub.j [S.sub.j (t) sin (ω.sub.o t+θ.sub.j)+S.sub.j (t) cos (ω.sub.o t+θ.sub.j)](32) 
    
     Upon phase detection, the low-pass-filtered version of the product of these two signals, σ E  (t)δ(t)| LP  may be provided, as follows: ##EQU10## Since the phase-lock loop forms a long-term average and since this average vanishes for S i  S j , for all i and j, the effective error signal σ E  (t)δ(t)| LP ,AV (representing the low-pass filtered, long-term averaged signal) may be represented by: ##EQU11## For independent sources, this error signal may be represented by: 
     
         σ.sub.E (t)δ(t)| .sub.LP,AV =Σ.sub.i S.sub.i.sup.2 |.sub.AV sin (θ.sub.i -θ.sub.E)(35) 
    
     the signal σ E  (t)δ(t)| LP ,AV will be dominated by the salient source, and is to be minimized by the phase-lock-loop control of θ E . The low-pass filtering may have a cutoff frequency of less than 1/2 ω o , but not less than the highest source signal frequency to be transmitted. The long-term averaging time (which may be provided by a low-pass filter) may be selected to be not shorter than the attack time of the control circuitry. 
     A synthetic channel T ET  may be derived from the formation of the following signals: ##EQU12## The sum of signals (36) and (37) is proportional to the synthetic version of a T T  channel in accordance with the UMX matrixing system, namely, that shown in Equation (30). 
     A second synthetic supplementary channel T EQ , which is a synthetic version of T Q  channel in accordance with the UMX matrixing system, may similarly be formed, as follows: 
     
         T.sub.EQ =Σ.sub.i [S.sub.i (t) cos (θ.sub.E +θ.sub.i)-S.sub.i (t) sin (θ.sub.E +θ.sub.i)](38) 
    
     In this connection, a T EQ  synthetic supplementary signal may be obtained from the SSB signal δ E  (t): 
     
         δ.sub.E (t)=Σ.sub.i [S.sub.i (t) cos θ.sub.i cos ω.sub.o t+θ.sub.E)-S.sub.i (t) sin θ.sub.i cos (ω.sub.o t+θ.sub.E)-S.sub.i (t) cos θ.sub.i sin (ω.sub.o t+θ.sub.E)-S.sub.i (t) sin θ.sub.i sin (ω.sub.o t+θ.sub.E)]=Σ.sub.i [S.sub.i (t) cos (ω.sub.o t+θ.sub.E +θ.sub.i)-S.sub.i (t) sin (ω.sub.o t+θ.sub.E +θ.sub.i)]           (39) 
    
     The T EQ  signal is formed from the SSB signal δ E  (t) as follows: 
     
         T.sub.EQ =2δ.sub.E (t) cos ω.sub.o t           (40) 
    
     A signal generator 100 for synthesizing the synthetic supplementary channel T ET  is shown in FIG. 10 of the drawings. The input signals to the signal generator 10 are the original transmission channel signals T.sub.Δ and T.sub.Σ as may be obtained, for example, directly from FM broadcast. The T.sub.Σ and T.sub.Δ input signals may also be obtained by rematrixing from the left and right, UMX transmission channel signals T L  and T R  or other locus-encoded linear combinations of the T.sub.Σ and T.sub.Δ signals, recorded on the respective tape or disk or disk recording tracks of a tape or disk recording as indicated hereinabove. 
     In the signal generator 10, the T.sub.Σ and T.sub.Δ original channel signals are transmitted through wideband 90° phase splitters 102, 104, respectively, which supply the phase related signals as shown in the drawing. The outputs of each phase splitter 102, 104 provide the modulating signals for double sideband modulation. In the illustrated embodiment, the T.sub.Δ audio signal is subjected to double sideband modulation at the fixed carrier frequency ω o  of the fixed frequency oscillator 106, the output of which is transmitted through phase splitter 108 which supplies the phase related signals as shown in the drawing. The T.sub.Σ signal is subjected to double sideband modulation at the variable frequency of the voltage controlled oscillator 110, the output of which is similarly transmitted through a phase splitter 112, to provide a quadrature phase relation for the carriers. 
     The combination of the double sideband modulated T.sub.Δ signals from the respective multiplier circuits 114, 116 at summer 118 produces the single sideband signal of T.sub.Δ which is supplied to automatic gain control circuit 120. Similarly, the combination of the double sideband modulated T.sub.Σ signals from the respective multiplier circuits 122, 124 at summer 126 produces the single sideband signal of T.sub.Σ which is supplied to automatic gain control circuit 128. In this manner, the single sideband signals are provided, with a quadrature phase relation for the carriers, except for the further difference in phase represented by the difference between the estimated phase angle θ E  and the encoding phase angle to which it is to correspond. Phase detection at detector 130 produces a signal proportional to the sine of this difference, and this error signal is in turn applied as a control signal to the voltage controlled oscillator 110 to adjust its frequency to bring the error signal to a null. 
     The signal from multiplier 122 is supplied as an input signal to multiplier 132 which is also supplied with the 0° phase signal of phase splitter 108 as the other input signal. The signal from multiplier 124 is supplied as an input signal to another multiplier 134, which similarly is also supplied with the 0° phase signal of phase splitter 108 as the other input signal. The output signals of multipliers 132, 134 are each subjected to low pass filtering at low pass filters 136, 138, and the filtered signals are combined by summer 140 to provide upon multiplication by two by multiplier 141 the synthetic supplementary channel T ET . The automatic gain-control circuits 120, 128 seek to maintain constant-level signals in well-known manners, and are included to seek to make the phase-lock loop operate nearly equally well with low-level salient signals as well as high-level ones. The second synthetic supplementary channel T EQ  may be generated similarly by circuitry 101 also shown in FIG. 10. In this circuitry, the phase-splitter outputs for the signal T.sub.Δ are supplied to a second single sideband modulator to generate δ E  (t) in a manner similar to the previously described generation of δ(t), but using the quadrature-pair oscillations from the voltage-controlled oscillator 110 instead of the fixed-frequency oscillator 106. However, the fixed frequency oscillator (with amplitude multiplied by 2 by amplifier 103) provides the oscillation used for demodulation. Accordingly, as shown in the drawing, the 0° signal from phase splitter 102 and the 0° signal from the phase splitter 112 are provided as input signals to multiplier 105, while the respective 90° signals from these phase splitters 102, 112 are provided as input signals to multiplier 107. The output signals of the multipliers 105, 107 are respectively added and subtracted at summer 109. The summer output and the amplified 0° signal from amplifier 103 form the inputs to multiplier 111, the output of which is passed through low pass filter 113 to provide the T EQ  signal. 
     Many variations of the indicated phase lock-loop system are possible. The roles of the fixed-frequency and voltage-controlled oscillators may be interchanged; the two components of σ E  (t) may be combined with reversed polarity and multiplied by cos ω o  t to save one balanced modulator in forming T ET , etc. 
     The synthetic channels T ET  and T EQ  will be continuously formed by the signal generator 10 regardless of whether the complex of sources contains a single salient source. On the other hand, the utility of these synthetic channels is questionable in the absence of such a salient source, and in such absence the synthetic channels generally ought best not be used in the decode. Consequently, it is desirable to make a test for saliency and to suppress the use of T ET  and T EQ  if the test be not satisfied. 
     In respect of the illustrated embodiment employing the T.sub.Σ and T.sub.Δ input signals, this saliency test may be cast in terms of the signals T.sub.Δ and T ET  *, where the signal T ET  * is the complex conjugate of signal T ET . In phasor notation, the sum of these two signals is given by: 
     
         T.sub.ET * + T.sub.Δ = Σ.sub.i S.sub.i e.sup.-jθ.sbsp.E + S.sub.i e.sup.-jθ.sbsp.i  = e.sup.-jθ.sbsp.E Σ.sub.i S.sub.i [1+ e.sup.i(θ.sbsp.E.sup.-θ.sbsp.i.sup.) ] (41) 
    
     Similarly, in phase notation, their difference is given by: 
     
         T.sub.ET * - T.sub.Δ =e.sup.-jθ.sbsp.E Σ.sub.i S.sub.i [1- e.sup.j(θ.sbsp.E.sup.-θ.sbsp.i.sup.) ]    (42) 
    
     Except for the phase-factor common to Equations (41) and (42), these describe BMX decodes for the bearing angles θ E  and θ E  + 180°, and the difference in the levels of these two signals can provide the saliency indication. The signal T ET  * may be generated as the low-pass version of σ E  cos ω o  t, and may be obtained by transmitting the input of the automatic gain control block 128 to the input of a multiplier whose other input is supplied from the 0° output of phase splitter 120 and then transmitting the output of the multiplier through a low-pass filter. 
     A control signal generator 150 for generating a control signal C based on this saliency test is shown in FIG 11. The signals T ET  * and T.sub.Δ are provided as essentially constant-level signals through the action of respective automatic gain control circuits 152, 154 indicated. The sum signal indicated by equation (41), and the difference signal indicated by Equation (42) are then formed by summer 156 and (negative) summer 158, respectively, as indicated in the drawing. The differences in the levels of the sum signal from summer 156, and the difference signal from summer 158, are determined by comparing the outputs of level detectors 160, 162 through subtraction of one of the level signals from the other by (negative) summer 164. The saliency condition may be selected, for example, at some predetermined ratio of the level of the sum signal to the level of the difference signal, which is intended to indicate a single image condition. When this condition is met, there would be at least this ratio between the strength of the decodable image at θ E  from the original transmission channels, and the 180° opposed location. In any event, the differences in the levels of these signals may be used to form the control signal C. 
     The level detectors 160, 162 shown may comprise a rectifier circuit and a filter. The filter may be arranged to exhibit two characteristic response times; an attack time and a release time, and these may be adjusted to match well with the characteristics of human hearing. Conventionally, the attack times may be in the order of 0.1 ms., and release times on the order of 100 ms. to one second. Thermoluminescent and electroluminescent devices may be combined with a photo diode for similar effects. Also, the level detectors may be provided wth logarithmic characteristics, so that level differences may correspond to amplitude ratios as indicated hereinabove. 
     The control signal C may be used to actuate gain-controlled amplifiers, and in this regard, the threshold device 166 of FIG. 11 may be set so that rather little change of the value C obtains if the level difference is less than a predetermined threshold corresponding to level differences that may be preassigned arbitrarily, but for which an assignment of at least about 6dB may be preferable. Above this threshold, large variations of signal C correspond to level variations. Except for thee limitations, the variations of signal C provided by the threshold circuit 166 may be monotonic with variations in level difference. 
     The gain-controlled amplifiers actuated by signal C may exhibit gain variations from null transmission to full transmission response to the whole range of the variation of C, in the case of the amplifiers that transmit signals T ET  and T EQ  to the matrix decode system. However, when full transmission of these signals obtains, there will be a 3dB loudness increases for a simple decode of the two original transmission channels plus the two synthetic channels T ET  and T EQ . To compensate for this, the transmission of each of the T ET , T EQ , T.sub.Δ, and T.sub.Σ signals may be reduced proportionately by up to 3dB corresponding to the full transmission of signals T ET  and T EQ , to compensate for the contribution of these channels to the decoded presentation signals. This compensation may be provided in a second set of four gain-controlled amplifiers to provide the 3dB variation, or the two functions may be combined in a single set of four appropriate amplifiers. The T.sub.Σ, T.sub.Δ (or T L  andT R )T ET  and T EQ  signals may be decoded in accordance with the disclosure of U.S. Pat. No. 3,906,156 with the T ET  and T EQ  signals being substituted for the T T  and T Q  signals, respectively. 
     In this connection, it is desirable to be able to form phase shifted versions of T ET  and T EQ  (to avoid the use of further phase shift circuits) for these signals in the decode circuitry from the signals already formed by the signal generator 100 of FIG. 10. For example, ##STR3## may be formed from signals available in the phase-lock loop circuit of FIG. 10, and their sum is the Hilbert transform of T ET . Similarly, the Hilbert transform of T EQ  is proportional to the signal δ E  (t) sin ω o  | LP . With the availability of these phase-shifted versions, including phase-shifted versions of T.sub.Σ and T.sub.Δ, the remaining mixing functions of decoding need involve no further phase shifting, and may be carried out by the appropriate mixing of the indicated phase-shifted signals. 
     It should be pointed out, however, that the functions of modulation and demodulation, as used in the phase-lock loop, involve frequency-dependent phase shifts arising mostly from the use of the necessry band-pass (not shown) and low-pass filters. These may be made to be essentially linear with frequency, to constitute simple delay. Matching delays may be inserted in the transmission path of signals T.sub.Σ and T.sub.Δ, and their phase shifted version before supplying these to the decode matrix. 
     The same phase-lock loop techniques may be used to generate synthetic supplementary channels for the 4-2-4 matrices known as QS (Sansui Electric Co.) and SQ (CBS Laboratories). For the QS matrix, it is known that, upon the insertion of a 90° phase shift between left and right channels, a BMX-type decode matrix may produce very good results. Thus, channels with properties very similar to T.sub.Σ and T.sub.Δ may be derived after such phase insertion has been made, and processed as indicated hereinabove. 
     The encode equation for the two transmission channel SQ matrix is: ##EQU13## from which the sum and difference (L T  - R T ) combinations are: ##EQU14## in which for Equation (45) the number &#34;0.7&#34; is used where the square root of 1/2 is intended, so that all numerical coefficients in Equation (46) have magnitude unity. The two channels of Equation (46) may serve as the signal inputs for the phase-lock loop current 10 of FIG. 10. 
     The phase-lock loop will then form a phase estimate θ E  that is an estimate of the phase difference between Σ T  and Δ T , namely, an estimate of θ in the following equation: ##EQU15## Then, the circuit 10 will form channels by subtracting these phase differences from the phase of Σ T , but adding them to the phase of Δ T , as follows, assuming one-at-a-time saliency for each of these sources: ##EQU16## 
     Upon forming sum and difference between these two signals T ET  and Q ET , the following two channels are derived; ##STR4## These channels may be named L TE  and R Te . 
     Under saliency conditions, there is thus obtained an augmented SQ matrix, as follows: ##EQU17## 
     This matrix equation is completely soluble for the source quantities S LB , S LF , S RF , and S RB . The classical test for this assertion is the calculation of the determinant of the 4× 4 matrix of Equation (50). This determinant does not vanish, verifying the assertion. 
     A saliency test for the SQ matrix will be somewhat different than that proposed for BMX or QS, because of inherent asymmetries in the SQ matrix. Such saliency tests may be devised for the SQ matrix, and have been implemented in various versions of SQ &#34;logic&#34; circuitry that have been devised; such SQ-saliency test methods and apparatus may be utilized to provide a control signal governing the contribution of synthetic channels L TE  and R TE  to the decode matrix. 
     While this aspect of the present invention has been particularly described with respect to systems employing two original encoded transmission channels, it will be appreciated that synthetic supplementary channels may be provided for systems employing more than two original transmission channels. For example, in a system employing three original encoded transmission channels, such as may be exemplified by T L , T R  (or T.sub.Σ, T.sub.Δ) and T T  signals of the UMX matrixing system, a fourth, synthetic supplementary channel, such as signal T EQ , may be generated from two or more of the original channels such as by utilizing the T.sub.Σ and T.sub.Δ channels as discussed in connection with the apparatus of FIGS. 10 and 11 (without utilization of the T ET  signal). However, it will be appreciated that the third original channel may be advantageously utilized with the other two original channels in the saliency test of single image locations, because the addition of the directional information of the third channel permits increased image resolution. 
     Another aspect of the present disclosure is directed to a logic system utilizing envelope generation in respect of at least two encoded original transmission channels for the provision of logic functions of source signal direction, even in the multiple source case where multiple sound sources at different bearing angles are encoded in the original transmission channels. The influence of the directional variation of the envelope-generated logic direction functions is limited to a frequency range below the audio band, and the envelope generated logic functions are utilized in the formation of a separate modulating signal for each presentaton signal decoded from the original transmission signals. 
     The application of the envelope generation logic system to the basic two-channel matrix (BMX) of the universal matrix system (UMX) may hereinafter be referred to as a &#34;BML&#34; system, and this aspect of the present disclosure will now be more particularly described in respect of an embodiment of such a BML system. 
     The theory of such envelope generation rests upon the fact that if a time-dependent function, designated x(t) for purposes of the following discussion, is a bandpass function, so that its spectrum is of negligible contact near zero frequencies, and is also of negligible content for frequencies higher than 2f c , where f c  is a central frequency in the bandpass, then the function x(t) may be represented by the following: 
     
         x(t)= x.sub.c (t) cos ω.sub.c t+ x.sub.s (t) sin ω.sub.c t (51) 
    
     where ω c  = 2πf c  and x c  (t) and x s  (t) are low-pass time functions with frequency content below f c . An envelope function r(t) of the function x(t) is then defined to be the square root of the sum of the squares of the xc(t) and x s  (t) functions as follows: 
     
         r(t)=√x.sub.c.sup.2 (t)+ x.sub.s.sup.2 (t)          52) 
    
     Furthermore, a 90° phase shifter circuit (also known as a Hilbert-transform filter) is realizable for bandpass signals. Transmission of a bandpass signal x(t) through such a phase shifter circuit produces a time dependent signal, which may be described y(t), as follows: 
     
         y(t)=-x.sub. c (t) sin ω.sub. c t+ x.sub.s (t) cos ω.sub.c t (53) 
    
     so that 
     
         x.sub.s (t)= x(t) sin ω.sub.c t+ y(t) cos ωt   (54a) 
    
     
         x.sub.c (t)= x(t) cos ω.sub.c t- y(t) sin ω.sub.c t (54b) 
    
     From Equations (54) it may be seen that the fundamental definition, Equation (52), implies that the envelope function r(t) may be obtained as the square root of the sum of the squares of the bandpass function x(t) and its Hilbert transform signal Y(t), as follows: 
     
         x(t)= √x.sup.2 (t)+ y.sup.2 (t)                     (55) 
    
     Application of such envelope formation as utilized in multidirectional matrix logic systems will now be more particularly described in respect of signal processing for two channel matrix systems such as the two channel BMX matrix system, together with circuitry for accomplishing such processing. 
     For each audio signal [e.g., designated x(t) in accordance with the previous discussion] of a BMX matrix, it will be appreciated that there are also obtainable the 90° phase shifted versions of each of these signals (e.g., y(t) the leading phase version and -y(t) for the lagging phase version). Thus, if x 93  (t) be the sum signal T 93   (i.e., T L  plus T R ) and x.sub.Δ (t) be the difference signal T.sub.Δ  (i.e., T R  -T L ) then for the example of a single source encoded in the BMX transmission channels, the following relationships may be established, and the indicated signals are all directly available or obtainable: 
     
         x.sub.Δ (t)= x(t) cos θ - y(t) sin θ     (56a) 
    
     
         y.sub.Δ (t)= y(t) cos θ+x(t) sin θ       (56b) 
    
     
         x.sub.Σ (t)= x(t)                                    (56c) 
    
     
         y.sub.Σ (t)= y(t)                                    (56d) 
    
     The indicated Equations (56) are given for a single source waveform x(t) at a bearing angle of θ, measured counterclockwise from the right direction. 
     As indicated by the previous discussion, envelope r(t) of a function may be defined from 
     
         r.sup.2 (t)= x.sup.2 (t)+ y.sup.2 (t)                      (57) 
    
     where y(t) is the leading phase, 90° phase shifted version of the x(t) signal. 
     In a similar manner, the envelope function r.sub.Σ (t) may be defined for the sum signal x.sub.Σ (t), and the envelope function r.sub.Δ (t) may be defined for the difference signal x.sub.Δ (t), as follows: 
     
         r.sub.Σ.sup.2 (t)= x.sub.Σ.sup.2 (t)+ y.sub.Σ.sup.2 (t) (58a) 
    
     
         r.sub.Δ.sup.2 (t)= x.sub.Δ.sup.2 (t)+ y.sub.Δ.sup.2 (t) (58b) 
    
     Then, level-envelope signals u(t) and v(t) may be defined from Equations (58) for the sum signal x.sub.Σ (t) and the difference signal, as envelope-normalized versions of these signals and their phase shifted counterparts, respectively, as follows: 
     
         u.sub.Σ (t)= x.sub.Σ (t)/r.sub.Σ (t)     (59a) 
    
     
         v.sub.Σ (t)= y.sub.Σ (t)r.sub.Σ (t)      (59b) 
    
     
         u.sub.Δ (t)= x.sub.Δ (t)/r.sub.Δ (t)     (59c) 
    
     
         v.sub.Δ (t)= y.sub.Δ (t)/r.sub.Δ (t)     (59d) 
    
     In the single source case, which is an important signal circumstances for logic systems, the envelope signal r(t) for the audio signal x(t), the envelope signal r.sub.Σ (t) for the sum signal x.sub.Σ (t), and the envelope signal r.sub.Δ (t) for the difference signal are all equal, as may be readily verified. Thus, for a single source condition, r(t) = r.sub.Δ (t). 
     The sine and cosine functions of the source bearing angle θ may then be formed as functions of the level-envelope signals of Equations (59), as follows: 
     
         cos θ = u.sub.Σ (t)u.sub.Δ (t)+ v.sub.Σ (t)v.sub.Δ (t)                                      (60a) 
    
     
         sin θ = u.sub.Σ (t)v.sub.Δ (t)- v.sub.Σ (t)u.sub.Δ (t)                                      (60b) 
    
     The signal functions of Equations (60) provide for the calculation of logic values (designated cos θ L  and sin θ L ) of cos θ and sin θ in the multiple source signal case, for which case these logic values will be time varying, and should be subjected to low-pass filtering, resulting in the low-pass filtered logic values c(θ L ) and s(θ L ) as more slowly varying quantities: 
     
         c(θ.sub.L)== cos θ.sub.L|.sbsb.LP     (61a) 
    
     
         s(θ.sub.L)= sin θ.sub.L|.sbsb.LP      (61b) 
    
     The square root of the sum of the squares of the sine and cosine logic values defines a logic filter gain function a; 
     
         a.sup.2 = s.sup.2 (θ.sub.L)+ c.sup.2 (θ.sub.L) (62) 
    
     such that the value of a may range from zero to one: 
     
         0≦ a≦ 1                                      (63) 
    
     The exact value of the logic gain function, a, will depend upon the speed of fluctuation of the cos θ L  and sin θ L  signals. 
     The squares complement of the logic-filter gain function, a, may be denoted as a bias signal b: 
     
         b= √1- a.sup.2                                      (64) 
    
     The bias signal b will also vary from one to zero: 
     
         0≦b≦1                                        (65) 
    
     The bias signal, b, and the logic values S(θ L ), C(θ L ), are used in the formation of the individual presentation signals for the different presentation directions. In this connection, a modulating function &#34;M&#34; is provided for each presentation signal, having a presentation bearing angle φ (measured counterclockwise from the right presentation direction), the modulating function being related to the sine and cosine values of the difference between the bearing angle of that presentation signal and the logic value of the encoded signal(s). Further in this regard, for the kth speaker at a bearing angle φ k , there is formed a signal C k  for this speaker, as follows: 
     
         C.sub.k = c(θ.sub.L) cos φ.sub.k + s(θ.sub.L) sin φ.sub.k                                               (66) 
    
     which is the filtered version of cos (φ k  -θ L ). The modulating function M for this kth speaker may be formed upon addition to the bias signal b, of this C k  signal: 
     
         M.sub.k = b+ C.sub.k                                       (67) 
    
     The modulating speaker M k  functions to modulate the signal for that speaker which is supplied from the matrix decode of the original transmission channels. This signal processing, including the provision of a modulating signal M k , is carried out for each presentation signal for each speaker of the presentation array. The modulating signal M k  for each presentation direction is then used to modulate the decode signal for that source direction to produce an overall presentation function having directional sharpness which is enhanced through operations of the logic source encode values determined through the application of the appropriate envelope functions. Further in this regard, if there be but a single source signal encoded in the original transmission channels the signal for a given speaker will be substantially the same as the decode signal provided by decoding of the QMX matrix channels (e.g., signals T.sub.Σ, T.sub.Δ, T T  and T Q ), because for the single source case, b= 0, a=1, and C k  = cos (φ k  - θ). 
     For the case of multiple sources, the logic value of the encoding direction, θ L , is &#34;captured&#34; by the strongest source, and this value changes as one or another source captures it in the course of time. Most of the time, it may be expected that θ L  will be only slowly time varying with characteristic frequencies well below the audio band. Under extremely complex signal conditions, however, fluctuations at an audio-band rate could characterize sin θ L  and cos θ L , causing waveform distortion. Then, the indicated low-pass filtering of sin θ L  and cosθ L  to produce S(θ  L ) and C(θ  L ) acts to progressively disconnect the BML system, by bringing the value of the logic gain function &#34;a&#34; near 0, and the value of the bias signal &#34;b&#34; near 1. 
     This filter presents a problem common to all gain-riding or compander-type circuits. The filter introduces a delay, so that the gain-riding function lags behind the signal conditions causing its actuation. Thus, linear filters could be used only if advance sensing of the signals may be provided. Otherwise, nonlinear filters characterized by differing attack and decay times should be used. Sometimes these should be multiple-stage filters employing a spectrum of decay times. Attack times of 1 ms and decay times grouped in the neighborhood of 100 ms are typical. The choice of these times are necessarily made to suit the specific application, and are based on acceptable performance to the ear, balanced (for decay times) against the waveform distortion introduced at the lowest audio frequency (e.g., 20 Hz). The &#34;pumping&#34; of more-or-less continuous low-level sounds on the part of sporadic high-level sounds is an effect that, in particular, should be kept of minimal aural consequence. 
     Schematically illustrated in FIGS. 12 through 17 is circuitry embodying various aspects of the BML system. FIG. 12 is an overall diagram of the various process blocks of the illustrated embodiment of an envelope-leveler logic system 200. FIGS. 13 and 15 show further detail of two of the processing blocks of FIG. 12, and FIG. 14, in turn, shown further detail of a processing block of FIG. 13. 
     The envelope-leveler logic system 200 of FIG. 12 comprises a logic signal generator 202 which operates on input signals T.sub.Σ (corresponding to signal xΣ(t) of the previous discussion), T.sub.Σ (corresponding to signal y.sub.Σ (t) of the previous discussion), T.sub.Δ (corresponding to signal x.sub.Δ (t) of the previous discussion) and T.sub.Δ  (corresponding to signal y.sub.Δ (t) of the previous discussion). The signals T.sub.Σ and T.sub.Δ, as indicated hereinabove, may be provided from any of their linear combinations, such as signals T L  and T R , and the 90° leading phase shifted versions of these signals, T.sub.Σ and T.sub.Δ, respectively, may be provided by Hilbert transform circuit means (not shown) such as 90° phase shifters. The logic signal generator 202 provides the logic signals C(θ L ) and S(θ L ), which provide the input signals to the bias signal generator 204. The bias signal generator 204 in turn generates the bias signal &#34;b,&#34; which is the complement of the filter gain signal &#34;a&#34; in accordance with the previous discussion. The bias signal generator 204 is illustrated in more detail in FIG. 15. 
     The logic signals S(θ L ) and C(θ L ) from this logic signal generator 202 also function as input signals, together with the bias signal b, for the speaker signal modulator 206. In the interest of representational clarity, only one speaker signal modulator 206 for the &#34;kth&#34; speaker is shown in FIG. 12; however, it will be understood that the logic circuit will include a separate modulator 206 for each output speaker signal, and bypass switching of this plurality of modulators may be mechanically ganged as indicated in the drawing. Each of the speaker signal modulators 206 is provided with the matrix decode signal for the respective speaker served by the modular. The appropriate matrix decode signals may be provided by a decoder such as described in U.S. Pat. No. 3,906,156 which decodes the T.sub.Σ and T.sub.Δ (or T L  and T R ) signals to provide presentation signals in accordance with the UMX matrixing system. Having generally described the processing blocks of the illustrated BML system, the components and operations of the system will now be discussed in more detail. Turning to FIGS. 13 and 14, relating specifically to the logic signal generator 202, it will be seen that the input signals T.sub.Σ and its 90° phase shifted version T.sub.Σ are transmitted to envelope leveler 208, while the signals T.sub.Δ and its 90° phase-shifted version T.sub.Δ are transmitted to another envelope leveler 210. The envelope levelers 208, 210 are each of the type illustrated in FIG. 14, in which input signals (designated x(t) and y(t) in accordance with the nomenclature of the previous discussion of such signals, and which may be either the T.sub.Σ, T.sub.Σ  or T.sub.Δ, T.sub.Δ  signal pair) are respectively transmitted to multipliers 212, 214. The output signals from multipliers 212, 214 are the respective u(t) and v(t) function signals. The other input signal to each of the multipliers 212, 214 is the output signal of a feedback loop comprising multipliers 216, 218 (connected as squaring circuits for the u(t) and v(t) signals), summer 220 (to provide the sum of the squares of the u(t) and v(t) signals), and differential amplifier 222 (to provide a level-set control for the sum of the squares signal). The appropriate four combinations of the V.sub.Σ, VΔ, uΣ and uΔ output signals of envelope levelers 208, 210 are multiplied by multiplier circuits 224, 226, 228, 230 to provide the respective signals u.sub.Σ v.sub.Δ, u.sub.Σ u.sub.Δ, v.sub.Σ v.sub.Δ and v.sub.Σ u.sub.Δ. The u.sub.Σ u.sub.Δ and v.sub.Σ v.sub.Δ signals are combined by summer 232 to produce a cosθ signal, in accordance with Equation (60a) and the v.sub.Σ u.sub.Δ signal is subtracted from the u.sub.Σ v.sub.Δ signal by (negative) summer 234 to produce a sin θ signal in accordance with Equation (60b). The respective cosθ and sinθ signals are transmitted through low pass filters 236, 238, which have a passband below the audio band (e.g., below about 20 Hz) to provide the filtered logic signals C(θ L ) and S(θ L ) in accordance with Equation (61). 
     As shown by FIG. 12, these logic signals are transmitted to the bias signal generator 204, which is shown in more detail in FIG. 15, where the logic signals C(θ L ) and S(θ L ) are squared by multipliers 240, 242, respectively. The phase inverted, squared signals are combined at summer 244, transmitted to differential amplifier 246 having a feedback loop including multiplier (as squarer) 248, the output of which is inverted at summer 250 to provide the bias signal &#34;b&#34; in accordance with Equation (64). 
     As shown in FIG. 12, the bias signal b and the logic signals are transmitted to each of the speaker signal modulators, including the modulator 206 for the kth speaker shown in the drawings. In the modulators, the logic signals are transmitted to sin-cos dividers 252, 254 which are set to sin and cos values of the bearing angle of the particular speaker served by each modulator (i.e., angle θ k  for the kth speaker served by modulator 206). The signal cos θ k  c(θ L ) formed by sin-cos divider 252, the signal sin θ k  S(θ L ) formed by sin-cos divider 254, and the bias signal b are combined by summer 256 to produce the modulating signal M k  in accordance with Equation (67). The modulating signal M k  in turn forms one input to multiplier 258, the other input being the matrix decode signal for the speaker. As indicated by FIG. 12, the illustrated BML system may be readily disconnected, if desired, by a gauged bypass switch 260. 
     In addition to the novelty of the overall design of the BML system, the novelty of the envelope leveler shown in FIG. 14 in respect of its implicit generation of the envelope without the use of low pass filters (except for the filters used to obtain the 90° phase shift) should also be noted. An additional embodiment of such an envelope generator is shown in FIG. 16, and such envelope generators could have other applications such as applications in compander service. For example, such an envelope leveler could be used in the compander of Burwen described at pages 906-911 of the Acoustic Engineering Society Journal for December, 1971. 
     As indicated, an alternate embodiment of an envelope leveler 208, 210 to that of FIG. 14 is shown in FIGS. 16 and 17. The envelope leveler 250 of FIGS. 16 and 17 operates on input signals designated x(t) and y(t), y(t) being the Hilbert transform signal of the x(t) signal in accordance with the previous nomenclature. These input signals are supplied to envelope generator 252, which is shown in more detail in FIG. 17, in order to produce the envelope signal r(t). In the envelope generator, the x(t) and y(t) signals are squared by respective multiplier circuits 254, 256 and the outputs are combined by summer 258. The sum of the squares output is supplied to a differential amplifier 260 with a feedback loop through multiplier 262 in order to produce the signal r(t), as indicated in FIG. 16. 
     In the envelope leveler 250, the x(t) and y(t) signals are also supplied to respective differential amplifiers 264,266 which are each included in a feedback loop through respective multipliers 268,270. The envelope signal r(t) is transmitted to each of the multipliers 268,270 as a second input signal, and the (inverted) output of each of the multiplier 268,270 is transmitted to the respective differential amplifier 264,266 such that the amplifier outputs are, respectively, the desired u(t) and v(t) signals. The usual method of deriving r(t) is by means of full wave rectification of the input signal x(t) followed by low pass filtering. However, such filtering is difficult for the proper removal of ripple for signals of frequencies predominantly near the lower end of the pass-band. This filtering also introduces a delay complicating the compander application--i.e., requiring the use of a nonlinear filter. It is believed that the filtering requirements for the envelope signal r(t) of Equation (55) such as provided by the envelope leveler of FIG. 16, would not be nearly so severe in compander applications. 
     It will be appreciated that in FIGS. 14 through 17, high quality multipliers and/or squarers (indicated by the symbol   are used and shown as elements in high-gain (gain G) feedback loops. It will also be appreciated that the bandwidths of these elements are necessarily signal-level dependent, so that, even when such elements are providing a maximun signal-attenuation effect, care should be taken that the bandwidth is sufficient for loop stability. For equipment optimization, this requirement will imply unusually great bandwidths for the slight-attenuation conditions. While certain previous aspects of the present invention have been described in detail in respect of particular matrix signals such as T.sub.Σ and T.sub.Δ signals, and linear combinations of these signals producing particular T L  and T R  signals, other combinations may also be utilized. For example, in addition to the T.sub.Σ, T.sub.Δ and T L , T R  linear combinations shown in Equations (23) and (24) [and the previously described combination of T T  and T Q  signals shown by Eauations (25) and (26)], there are other combinations serving particular compatibility requirements. For example, in one system for making quadrasonic tape recordings, the number of tracks are doubled, in comparison to stereo practice, to provide tracks bearing the further information required. In particular, each of the original stereo tracks 300,302 is subdivided into two tracks as shown in FIG. 8. The left track 300 is subdivided into left front and left back tracks 304,306, and the right track 302 is similarly subdivided into right front and right back tracks 308,310. The subdivided tracks then provide a kind of compatibility for playback on tape machines provided with only two-track-stereo pickup heads, in that the left track pickup will provide the front-back sum of the left track information, and the right track pickup will similarly provide the front-back sum of the right track information, to provide stereo-compatible signals. 
     As discussed in connection with Equations (6), (7), (8), and (9), however, this means of obtaining stereo-compatible signals omits all information relating to differentiation of frong from back with the exception of four-speaker decoding that presents such differentiation from the two signals along, as may be provided in the manner of 4-2-4 matrices. A simple way to avoid this defect is to let the left-transmit and right-transmit signals of the 4-2-4 matrix appear as common signals in the two left channels and in the two right channels, respectively, the remainder in each being a matrix-supplement channel: 
     
         &#34;LF&#34; = T.sub.L + M.sub.L,                                  (68) 
    
     
         &#34;lb&#34; = t.sub.l - m.sub.l,                                  (69) 
    
     
         &#34;rf&#34; = t.sub.r + m.sub.r,                                  (70) 
    
     
         &#34;rb&#34; = t.sub.r - m.sub.r,                                  (71) 
    
     in which the matrix-supplement channels M L  and M R  are independent linear combinations, independent of T L  and T R , of the quadrasonic information necessary to provide a full discrete decode of that information. 
     A further condition may be placed on the matrix supplement channels, however, namely that &#34;LF&#34; should represent a predominantly left-front signal, and similarly for the other track signals. 
     A signal component &#34;predominantly&#34; appearing in one of four directional signals should have at least a 6 db higher level in that signal than the same signal component in any of the other three directional component signals if an impression approaching &#34;discreteness&#34; is to be the result. 
     For example, in the SQ-type 4-2-4 matrix ##EQU18## the matrix supplement channels may be ##EQU19## The tape channels then are  and the requisite predominance may be observed for the first row of Equation (74) as an example, by noting the numerical values 2.414, 0.414, and 0.707. Although the coefficients for the supplementary channels differ from those of U.S. Pat. No. 3,761,628 directed to the SQ matrix, the entire set of equations, Equations (72) and (73), are soluable for the four quantities S LB , S LF , S RF , and S RB , since the determinant of the coefficients does not vanish, as may be easily verified. 
     When four QMX channels of the UMX system are to be used, then suitable equations may be as follows: ##EQU20## Then &#34;RF&#34;, according to Equation (70), for example, is ##EQU21## The first trigonometric factor provides nulls at LF and RB, and the second at CB (center back), while the maximum is near RF, as may be seen in the plot of magnitudes in FIG. 9. Reflection about the CR (center right) axis provides the plot for &#34;RB&#34;, while reflection about the CF axis produces the plot for &#34;LF&#34;, etc. The requisite predominance is evident from the plot. 
     There are other specific linear combinations that are of interest, relating to stereo compatibility of the BMX matrix. For example, Equations (23) and (24) show linear combinations T L  and T R  of the T.sub.Σ and T.sub.Δ channels, which are each null in respect of signals encoded from a direction corresponding to the encoding azimuth direction of the other of the T L  and T R  signals, which have encoding azimuths of 0° and 180°, respectively, (i.e., at center right and center left when the bearing angle θ is measured counterclockwise from the center right direction). However, many practitioners of quadrasonic disk recording are habituated to regarding left front (LF) and right front (RF) as the encoding points corresponding to stereo left and stereo right (SL and SR), and for a variety of reasons including left-right phase difference considerations, it may be desirable to provide stereo compatible encoding azimuths other than 0° and 180°. In general, if the encoding azimuths for SL and SR are to be moved forward by a given angular increment δ, so that SL= 180°-δ, and SR= δ, with full separation to be obtained between these encoding points, then appropriate linear combinations T SL  and T SR  of T.sub.Σ and T.sub.Δ, which provide for stereo compatible encoding azimuths of δ and 180°-δ, may be represented, as follows: ##EQU22## These functions become respectively null at δ° and 180°-δ; such that a source signal to be encoded from the δ° direction will not appear in the T SL  signal and a source signal to be encoded from the 180°-δ direction will not appear in the T SR  signal, thereby providing signals with stereo compatible encoding azimuths of δ° and 180°-δ, respectively. The mono-compatible signal T S Σ  for the T SL  and T SR  signals is then 
     
         T.sub.SΣ = cos δT.sub.Σ =T.sub.SL + T.sub.SR (83) 
    
     while the difference signal T S Δ is as follows: 
     
         T.sub.SΔ = T.sub.SR = T.sub.SL = j sin δT.sub.Σ + T.sub.Δ                                             (84) 
    
     plainly T S Σ is omnidirectional while the illustrated example of T S Δ shows the greater magnitude for the back half azimuths (θ between 0° and -180°). 
     The total energy of the two channels is 
     
         |T.sub.SL |2+ |T.sub.SR |2= (1-sin θ sin δ)|S.sub.θ |2   (85) 
    
     and the front-back loudness ratio is calculated in the table below, for stereo playback 
     
                       Table 1.                                                    
______________________________________                                    
δ    CB/CF         RB/RF                                            
______________________________________                                    
10° 1.52 dB       0.569 dB                                         
20° 3.10 dB       1.20  dB                                         
22.5°                                                              
           3.50 dB       2.41  dB                                         
30° 4.77 dB       3.21  dB                                         
45° 7.66 dB       4.77  dB                                         
______________________________________                                    
 The left-right separation ratio is
 ##EQU23##
 which for θ = 135° is shown in the table below:
 
    
     
                       Table 2.                                                    
______________________________________                                    
                              CF (PWM)                                    
δ                                                                   
       [T.sub.SL /T.sub. SR ]                                             
                     Phase    Phase                                       
______________________________________                                    
10°                                                                
        9.40 dB      80°                                           
                              56.9°                                
20°                                                                
       11.81 dB      70°                                           
                              42.5°                                
22.5°                                                              
       12.59 dB      67.5°                                         
                              38.7°                                
30°                                                                
       15.67 dB      60°                                           
                              26.9°                                
45°                                                                
       infinite      45°                                           
                              0.00                                        
______________________________________                                    
 
    
     The last column of the above is the left-right phase difference for equal signals encoded at LF and RF, the pairwise mix representation of CF. The formula of this column is ##EQU24## 
     If 5 dB can be regarded as a largest allowable loudness variation, but 3 dB be regarded as more tolerable, then, from Table 1, a suitable range for δ would be in the range of 20° to 30°, providing very desirable separation ratios and phases as may be seen in Table 2. The value δ = 22.5° may be taken as being of special interest. 
     A further adjustment of phase may be made by replacing T SL  by T SL  exp (-jσ) and T SR  by T SR  exp (jσ). It may be shown, however, that for σ ≠ 0°, this introduces a loudness variation into T S Σ while not reducing the stereo loudness variation and not altering the left-right separation ratios. Thus, while not of great interest, this variation may find use in some applications. 
     The decode equations may be formed from ##EQU25## according to Equation (14), omitting T T  and T Q , with the result ##EQU26## to obtain the BMX kernel, Equation (20). Insertion of T SR  and T SL  into a decoder specifically designed for decoding T R  and T L  signals generally will result in a presentation error that is small if δ be small, a backward shift on the right, and a forward shift on the left. 
     Of course, T SL  and T SR  may be used with supplementary channels such as the T T  and T Q  channels to obtain the QMX kernel and the enhanced directionality accruing therefrom. Similarly, sum-difference combinations analogous to Equations (75)-(78) may be devised which may be useful in applications such as tape recording. These are ##EQU27## Then a right front signal &#34;RF&#34; and a right back &#34;RB&#34; signal would be as follows:  with the corresponding &#34;left&#34; sum-differences showing mirror symmetry with these about the CF axis. It is seen that for δ = 22.5°, the two factors in Equation (95) are rotated by only 7.5° in comparison to the two factors of Equation (80) so that very similar predominances may be expected in the corresponding plots. 
     The front-to-back loudness variation exhibited in Equation (85) may be compensated, if desired, in various ways. In this connection, the loudness variation may be diminished by diminishing the contribution of δ at the front and back locations. One way of doing this is to replace δ by d cosθ in Equation (81) and by -d cosθ in Equation (82). However, this introduces the theoretical possibility of a substantial θ dependence not expressible by linear combinations of unity and exp(-jθ), a substantial nondecodable θ dependence in that the decode may not be arranged to produce a dependence solely upon the difference between source and presentation azimuths. In this case, in order to avoid such complications, it is appropriate to provide that d will be small, (e.g., about 0.35 or less), so that the error contribution will be small (e.g., 10 db or more lower than the θ dependent signal component). Then the approximation 
     
         exp(-jd cos θ)≃λ - jd cos θ (97) 
    
     will have the further advantage of simplifying the error term, a term which is expressible in terms of exp(-jθ) and exp(jθ), so that this term would be decodable if a third channel were to be available. In this approach, the polarity pattern of cosθ has much the same effect as the change in polarity for δ between Equations (81) and (82), and λ is to be a parameter to permit T SL  to be null in both its real and imaginary parts. 
     With these modifications, we have the λ formulation, ##EQU28## and T SL  is null for θ i  = θ such that 
     
         d= tan θ                                             (98) 
    
     and 
     
         λ = cos θ 
    
     These two equations are compatible if ##EQU29## in which K would appear in the alternative K formulation ##EQU30## The null for T SR , allowing for a reversal in polarity of cosθ obtains at the same value for sinθ but because of the reversal in polarity of the cosine, the null angle is 180°- tan -1  d, for T SR , but tan -1  d for T SL . Thus, a value of d corresponding to δ = 20° is d= 0.36, and for this value, K= 1.06. The total energy corresponding to Equation (85) does, however, show a loudness variation of the same energy magnitude and pattern as the mono variations, corresponding to Equation (83), namely 1 + d 2  cos 2  θ, with a maximum energy ratio (1+ d 2 ) to 1. Both front and back are diminished in loudness by about 0.5 dB for d= 0.36, relative to the sides. 
     Decode follows from the fact that the three equations ##EQU31## may be solved for u, cos θ, and sin θ in terms of T SL , T SR , -jT ST , the determinant being unity. From these solutions, the combinations 
     
         T.sub.Σ = u                                          (105) 
    
     
         T.sub.Δ = cos θ - j sin θ                (106) 
    
     
         T.sub.T = cos θ + j sin θ                      (107) 
    
     may be assembled for decode in a usual TMX fashion. Then, if T ST  be made zero, an approximate BMX decode enjoying benefits thereof is the result. The above equations are, of course, compatible with the use of a usual UMX fourth channel, T Q , and other linear combinations of these channels may also be provided. 
     The solutions corresponding to Equations (105), (106), and (107) are ##EQU32## as exact expressions for signals suitable for TMX decode, or QMX decode, if T Q  be also available. If T ST  be not available, then Equation (68) becomes simply T.sub.Σ plus the following error terms: ##EQU33## and Equation (109) becomes simply T.sub.Δ plus the following error terms: ##EQU34## as nondecodable terms causing an error in the BMX decode. 
     The utility of the d formulation plainly depends upon d being small (K near unity), since the compromise with the fully decodable δ formulation is for trading the nondecodable terms against the larger loudness variation of the δ formulation. The σ phase shifts may also be introduced in the d formulation with further loudness variation in mono only. 
     Stereo compatible transmission channels having encoding loci at other than 0° and 180° may, of course, be used in systems employing synthetic supplementary channels or logic systems, such as envelope logic systems, as described hereinabove. The stereo compatible transmission channel signals may be recorded on conventional recording media such as magnetic tape and recording disks by means of known apparatus and procedures. For example, the T SL  and T SR  signals may be recorded as the left and right channels of a recording disk by means of a 45--45 cutter, and the recording may be played back by conventional monaural and stereo recording equipment, as well as suitable matrix decode equipment. Furthermore, additional channels such as T T  or T ST  and T Q  may be recorded with the T SL  and T SR  signals as angularly modulated signals in a multiple frequency manner as generally described in application Ser. No. 468,238 referred to hereinabove. 
     Apparatus to provide or decode the mono and stereo-compatible signals may comprise means for producing and combining the various signal components defined by Equations (80), (81) and (100)-(107), respectively, and conventional circuit elements having the appropriate functional parameters may be used for this purpose. Illustrated in FIG. 18 is an embodiment of apparatus 300 for providing signals T SL  and T SR  as defined by Equations (81) and (82). In the encoding apparatus 300, source signals are transmitted to a basic UMX matrix circuit 304, such as described in U.S. Pat. No. 3,906,156 to provide signals T.sub.Σ and T.sub.Δ. The T.sub.Δ signal is transmitted through polarity splitter 306 to provide +T.sub.Δ and -T.sub.Δ versions of this signal. Phase modifier 304 provides phase modified versions of the T.sub.Σ signal which are respectively leading, and lagging the T.sub.Σ signal in phase by δ° (in respect of the T.sub.Δ signal from splitter 306). The lagging T.sub.Σ signal and the -T.sub.Δ are combined in summer 308 to provide the T SL  signal. The +T.sub.Δ signal is combined with the leading T.sub.Σ signal at summer 310 to provide the T SR  signal. The T SL  and T SR  signals are shown with a (δ) subscript to indicate the stereo compatible signals of Equations (81) and (82). 
     Illustrated in the FIG. 19 is apparatus 350 for decoding this δ version of the stereo compatible signals T SL  and T SR . These signals are transmitted to matrix means 352 which provides signals T SR  exp(-jδ) and T SL  exp(jδ), together with the oriinal signals which, in turn, are fed to means 354 which applies the factor 1/cos δ to the signals, appropriate combinations which are combined at summers 358 and 356 to provide UMX sgnals T.sub.Σ and T.sub.Δ, which are fed to a basic UMX decoder 360 to produce speaker presentation signals. 
     Illustrated in FIG. 20, which is analogous to FIG. 7 of U.S. Pat. No. 3,906,156, is circuitry 400 for producing the &#34;K&#34; version of signals T SL  and T SR  as defined by Equations (100) and (101). In this apparatus, each source signal is transmitted to a polarity splitter 402, the outputs of which are in turn transmitted to sin-cos pots 404 adjusted to the respective bearing angle of the source. For each source signal, the -cos and the +cos signals are appropriately modified as indicated by means 406 for introducing the &#34;K&#34; factor, and means 408 for introducing the &#34;d&#34; factor, and the appropriate signals are combined at summers 410. The proper phase relationship is provided by reference phase and +90° phase shifters 412, 416. The appropriate combination of signals are combined at summer 414 to respectively provide the T SL  and T SR  signals, which are shown having a (K) subscript to indicate the &#34;K&#34; version of signals defined by Equations (100) and (101). 
     Illustrated in FIG. 21 is apparatus 450 for decoding stereo compatible signals of the &#34;K&#34; formulation in accordance with Equations 102-107. In this connection, signals T SL , T SR  and -jT ST  are transmitted to a matrix circuit means 452 which solves the matrix of Equations 102-104, the determinant of the matrix being unity, to produce the T.sub.Σ, T.sub.Δ and T T , as set forth in Equations 105-107. As indicated hereinabove, the T ST  signal is a stereo compatible third channel, including signal components adapted to compensate for error terms which would otherwise be introduced in the basic UMX decode. While means for producing T ST  third channel signals are not illustrated, it will be appreciated by one skilled in the art that the defined components of this signal may be combined by appropriate apparatus in view of the teachings of this disclosure and the previously referred to patents and applications. The signals T.sub.Σ, T.sub.Δ and T T  provided by the matrix circuit 452 are transmitted to a basic UMX decoder 454 which provides directional loudspeaker signals for presentation to a listener. 
     It will be appreciated that various sign changes may be made in the definitions of certain of the signals and that signals have generally been defined herein in respect of a unit magnitude which may be varied as desired. Moreover, while aspects of the present invention have been described with respect to various specific embodiments, it will be apparent that numerous variations, modifications and adaptations may be made by one skilled in the art based on the present disclosure, and such modifications, variations and adaptations are intended to be included within the spirit and scope of the present invention as set forth in the following claims. 
     Various of the features of the invention are set forth in the following claims.