Brevet FR3060830A1 — publié le 20180622

Traitement en sous-bandes d'un contenu ambisonique réel pour un décodage perfectionné La présente invention se rapporte au domaine du traitement de signal audio ou acoustique, et plus particulièrement au traitement de contenus sonores multicanal réels au format ambiophonique (ou « ambisonique » ci-après). La technique ambisonique consiste à exploiter dans chaque bande de fréquences un sous-ensemble de canaux qui possèdent des caractéristiques de directivité recherchées. A titre d'exemple d'application, on peut citer : - La séparation de sources sonores : o Pour le divertissement (karaoké : suppression de la voix), o Pour la musique (mixage des sources séparées dans un contenu multicanal), o Pour les télécommunications (rehaussement de la voix, débruitage), o Pour la domotique (commande vocale), o Le codage audio multicanal. - Le décodage pour une diffusion multicanal : o Pour le cinéma, o Pour la musique, o Pour la réalité virtuelle. L'ambisonie consiste en une projection du champ acoustique sur une base de fonctions harmoniques sphériques (base illustrée sur la figure 1), pour obtenir une représentation spatialisée de la scène sonore. La fonction Ymcrn (0,0) est l'harmonique sphérique d'ordre m et d'indice nu, dépendant des coordonnées sphériques (O, 0), définie avec la formule suivante : cos né si o- = 1 Yr°1132 (0, = Pmri(cos 0). { sin nO si o- = -1 et n > 1 où Pmn(cos (p) est une fonction pôlaire impliquant le polynome de Legendre : (m-n)! 5 \ dn 1-mn (X) - n (m+n)! , - ' - cos2 x) 2 dxn rm avec cc) = 1 et co = 2 pour n > 1 n et Pm(x) = -a 2m 1m dan (x2 - 1)m ! 3060830 2 Dans la représentation de la figure 1, le premier « vecteur » de la base d'harmoniques sphériques (en haut de la figure 1) correspond à l'ordre m=0, les trois « vecteurs » dans la ligne suivante correspondent à l'ordre m=1 (orientés suivant les trois directions de l'espace), etc. En pratique, un encodage ambisonique réel se fait à partir d'un réseau de capteurs, généralement 5 répartis sur une sphère, qui sont combinés pour synthétiser un contenu ambisonique dont les canaux respectent au mieux les directivités des harmoniques sphériques (comme illustré sur la figure 2). En référence à la figure 2, un microphone MIC comporte une pluralité de capsules piézoélectriques Cl, C2, ... qui reçoivent des ondes sonores selon différentes directions d'arrivée de l'espace. Une unité de traitement UT recevant les signaux issus de ces capsules réalise un encodage ambisonique à l'aide 10 d'une matrice de filtres présentée ci-après, et délivre des signaux ambisoniques (formalisés dans une base d'harmoniques sphériques du type illustré sur la figure 1). On décrit ci-après les principes de base de l'encodage ambisonique. Le formalisme ambisonique, initialement limité à la représentation de fonctions harmoniques sphériques d'ordre 1, a par la suite été étendu aux ordres supérieurs. Le formalisme ambisonique 15 avec un nombre de composantes plus important est communément nommé « Higher Order Ambisonics » (ou « HOA » ci-après). A chaque ordre m correspondent 2m+1 fonctions harmoniques sphériques, comme illustré sur la figure 1. Ainsi, un contenu d'ordre M contient un total de (M+1)2 canaux (4 canaux à l'ordre 1, 9 canaux à l'ordre 2, 16 canaux à l'ordre 3, et ainsi de suite). 20 On entend ci-après par « composantes ambisoniques » le signal ambisonique dans chaque canal ambisonique, en référence aux « composantes vectorielles » dans une base vectorielle qui serait formée par chaque fonction harmonique sphérique. Ainsi par exemple, on peut compter : - une composante ambisonique pour l'ordre m=0, - trois composantes ambisoniques pour l'ordre m=1, 25 - cinq composantes ambisoniques pour l'ordre m=2, - sept composantes ambisoniques pour l'ordre m=3, etc. Les signaux ambisoniques captés pour ces différentes composantes sont alors répartis sur un nombre N de canaux qui se déduit de l'ordre maximum m qu'il est prévu de capter dans la scène sonore. Par exemple, si une scène sonore est captée avec un microphone ambisonique à 20 capsules piézoélectriques, alors l'ordre ambisonique maximum capté est M=3, afin qu'il n'y ait pas plus de 20 3060830 3 canaux N=(M+1)2, le nombre de composantes ambisoniques considérées est 7+5+3+1 = 16 et le nombre N de canaux est N=16, donné par ailleurs par la relation N=(M+1)2, avec M=3. La captation ambisonique x(t) d'ordre M et composée de N sources sonores s, d'incidence (Oi,cpi) se propageant en champ libre peut s'écrire alors mathématiquement sous la forme 5 matricielle suivante : 1 1 x(t) = As(t) =s(t) Yen(01, 01) ON) Où A est une matrice dite « matrice de mélange », de dimensions (M+1)2 x N et dont chaque colonne A, contient les coefficients de mélange de la source i. Physiquement, cette matrice A correspond aux coefficients d'encodage de chaque source i, associés 10 à chaque direction de chaque source i. Pour extraire les sources d'un tel contenu, il faut procéder à l'estimation d'une matrice B dite « matrice de séparation », inverse de la matrice A. Pour obtenir la matrice B, une étape de séparation aveugle de sources peut être mise en oeuvre, par exemple en utilisant un algorithme d'analyse en composantes indépendantes (ou « ACI » ci-après), ou encore un algorithme d'analyse en composantes principales. La matrice B=A-1 permet l'extraction des sources 15 par l'opération suivante : s(t) = Bx(t) Cette étape revient à faire de la formation de voies (ou « beamforming » ci-après), c'est-à-dire à combiner différents canaux ayant des directivités distinctes, afin de créer une nouvelle composante possédant la directivité souhaitée. Un exemple de beamforming pour extraire trois composantes, pour un contenu HOA d'ordre 2, 4 ou 6, est illustré figure 3. Plus l'ordre est élevé, plus le 20 beamforming est directif et le nombre de composantes pouvant être extraites est élevé. En pratique, la génération des signaux ambisoniques x(t)=As(t) passe par une étape intermédiaire de captation microphonique tel qu'illustré sur la figure 2, où les sources s(t) sont captées par les capsules du microphone MIC pour former les signaux pl, p2, p3... On formalise ensuite la matrice d'encodage microphonique E telle que x(t)=E.p(t), pour obtenir les composantes ambisoniques xl, 25 x2, xN (dans N canaux ambisoniques comme illustré sur la figure 4). En référence maintenant à la figure 4, on estime, comme présenté ci-avant, la matrice de décodage B inverse de la matrice A, pour déterminer les signaux de sources sl, s2, s3 : s(t) = Bx(t) 3060830 4 Pour décoder un contenu HOA sur un système de haut-parleurs, la démarche est similaire. On acquière des signaux ambisoniques dans N canaux x1, x2, ..., xN, mais, ici, au lieu de considérer s(t) comme la somme des contributions de sources, on considère s(t) comme la somme des signaux émis par un jeu de haut-parleurs (ce qui permet d'alimenter alors effectivement ces haut-parleurs avec les 5 signaux s1, s2, s3...). On formule donc ici la matrice de décodage B à partir des positions des haut- parleurs d'un système de restitution sonore et on extrait les signaux destinés aux haut-parleurs selon le même procédé que celui utilisé pour la séparation de sources. En réalité, les capteurs utilisés possèdent des limitations physiques qui entrainent une dégradation de l'encodage microphonique, et donc une dégradation de la directivité des composantes 10 ambisoniques. Par exemple, l'encodage des hautes fréquences se dégrade lorsque l'espacement inter-capteurs devient approximativement plus grand qu'une demi-longueur d'onde : ceci est dû au phénomène de repliement spatial. En basses fréquences, les capsules microphoniques tendent à devenir omnidirectionnelles et il devient impossible d'obtenir les directivités recherchées. Plus précisément, les dégradations en basses fréquences sont plus marquées lorsqu'il s'agit de synthétiser 15 des composantes ambisoniques d'ordre élevé. De manière générale, les directivités associées sont plus complexes et donc plus sensibles aux variations de propriétés des capteurs. La figure 5 illustre le degré de corrélation entre un encodage théorique et un encodage réel à partir d'un microphone sphérique à 32 capsules, en fonction de la fréquence et de l'ordre ambisonique. La figure 5 montre que le degré de corrélation le plus élevé est généralement atteint pour des fréquences comprises 20 entre 1 kHz et 10 kHz. Néanmoins, pour les autres gammes de fréquences (excepté pour les ordres ambisoniques 0 et 1), l'extraction de sources ne conduirait pas toujours au même résultat pour un encodage théorique et pour un encodage réel de ces mêmes sources. Plus précisément, pour des fréquences en-dehors de l'intervalle [1 kHz-10 kHz], Les composantes extraites sont potentiellement dégradées. 25 La figure 6 montre la directivité réelle dans le plan horizontal des premières composantes des ordres 0, 1, 2 et 3 en fonction de la fréquence sonore. Il apparaît, sur la figure 6, que les composantes réelles ne sont pas convenablement encodées. En effet, si on prend l'exemple de la composante de l'ordre 0 à la fréquence de 10 kHz, on constate qu'elle n'est pas circulaire, contrairement à la composante théorique et à la même composante calculée aux fréquences entre 300 et 1000Hz. Ainsi, 30 la directivité de cette composante à la fréquence de 10kHz n'est plus respectée, ce qui pourrait induire un rendu spatial dégradé. Par ailleurs, les composantes à l'ordre 1, 2 et 3 ont également des directivités biaisées pour des fréquences plus basses que 10 kHz. 3060830 5 Plus généralement, dès lors que la directivité théorique n'est plus respectée, le beamforming effectué ne permet plus d'extraire convenablement les composantes recherchées. Par exemple, cela se traduit par l'apparition d'interférences pendant la séparation des sources. Cela peut également se traduire par une dégradation du rendu spatial dans des bandes de fréquences concernées par une 5 diffusion multicanal. Plus particulièrement, on constate une perte d'énergie en basses fréquences dans les ordres élevés pendant l'encodage. Cela induit que les sources extraites grâce à des canaux d'ordres élevés peuvent perdre une partie de leur énergie dans les fréquences concernées. L'utilisation du beamforming pour la séparation de sources ou la restitution d'un contenu ambisonique idéal ou d'une captation multicanal est déjà utilisée notamment pour la séparation, ou 10 encore pour le décodage multicanal. Pour la séparation de sources, une inversion de la matrice de mélange estimée par analyse en composantes indépendantes est utilisée pour extraire les sources. Pour le décodage multicanal, la matrice des coefficients ambisoniques relatifs aux haut-parleurs peut être inversée. En revanche, le traitement d'un contenu ambisonique réel, affecté par les limitations physiques du système d'enregistrement, n'est pas abordé dans l'art antérieur. La seule solution 15 proposée actuellement est de limiter la bande-passante totale des sources extraites, ce qui n'est pas satisfaisant. La présente invention vient améliorer cette situation. Elle propose à cet effet un procédé, mis en oeuvre par des moyens informatiques, de traitement d'un contenu ambisonique comportant une pluralité de composantes ambisoniques d'une pluralité 20 d'ordres définissant une succession de canaux ambisoniques dans chacun desquels est représentée une composante ambisonique, le procédé comportant : - un filtrage fréquentiel des composantes ambisoniques dans une pluralité de bandes de fréquences, - une élaboration d'une matrice de décodage ambisonique, 25 - un traitement de la matrice de décodage ambisonique pour extraire, par réduction de dimension de matrice, une pluralité de sous-matrices de décodage ambisonique chacune associée à un ordre ambisonique et à une bande de fréquences choisie pour cet ordre ambisonique, - des applications respectives des sous-matrices de décodage aux composantes 30 ambisoniques dans chaque bande de fréquences choisie, et une reconstruction bandes-à- bandes des résultats desdites applications respectives, pour délivrer une pluralité de signaux décodés, associés chacun à une source sonore. 3060830 6 On entend ici par « source sonore » aussi bien : - une source sonore effectivement identifiée et localisée dans l'espace tridimensionnel (en technique d'extraction de source), auquel cas la matrice de décodage est une matrice de séparation de sources, ou 5 - un haut-parleur parmi plusieurs haut-parleurs, de position bien identifiée dans l'espace, et alimenté en particulier par l'un des signaux décodés précités. L'élaboration de sous-matrices de décodage ambisonique pour chaque bande de fréquences, et pour chaque ordre ambisonique, permet de tirer parti dans chaque bande de fréquences d'un nombre 10 maximum de canaux ambisoniques qui sont réellement valides dans chaque sous-matrice, afin de restituer un signal décodé peu ou pas dégradé. Selon une réalisation, chaque sous-matrice de décodage ambisonique est associée à une bande de fréquences choisie en fonction d'un critère de validité des composantes ambisoniques de l'ordre auquel est associée ladite sous-matrice, dans ladite bande de fréquences choisie. 15 Une telle réalisation permet d'isoler les composantes ambisoniques constituant chaque ordre, afin de les traiter dans la plage de fréquences dans laquelle elles sont valides. Par « valides », on entend un respect de la représentation ambisonique théorique, comme par exemple l'ordre m=4 dans la bande de fréquences 4000 à 6000 Hz dans l'exemple de la figure 5, ou encore l'ordre m=3 dans la bande de fréquences 2000 à 9000 Hz. 20 Ainsi, dans une réalisation, le critère de validité des composantes peut être défini par des conditions de captation desdites composantes ambisoniques, par au moins un microphone ambisonique. Dans cette réalisation par exemple, le procédé peut comprendre en outre : - une réception de données d'au moins un microphone ambisonique utilisé pour capter lesdites composantes ambisoniques ; 25 - une détermination des bandes de fréquences choisies pour construire lesdites sous- matrices, en fonction desdites données de microphone ambisonique. La connaissance des données du microphone ambisonique utilisé pour la captation ambisonique permet d'affiner la détermination des bandes de fréquences choisies pour l'élaboration des sous-matrices. En effet, le traitement ambisonique est fait sur des sous-matrices dont les composantes 30 ambisoniques répondent strictement au critère de validité dans les bandes de fréquences associées. 3060830 7 Toutefois, les données du microphone ambisonique utilisé pour la captation ne sont pas toujours accessibles. En variante, on peut donc prévoir la détermination des bandes de fréquences à l'aide d'un abaque préalablement établi à partir de mesures effectuées sur une pluralité de microphones ambisoniques, afin d'établir des plages de fréquences « moyennes », associées à un ordre 5 ambisonique, dans lesquelles les composantes ambisoniques de chaque ordre ambisonique répondent généralement au critère de validité précité. Ainsi, selon une réalisation, chaque sous-matrice de décodage ambisonique étant associée à un ordre ambisonique et à une bande de fréquences choisie pour cet ordre ambisonique, une bande de fréquences peut être choisie dans une plage de 100Hz à 10kHz pour l'ordre 10 ambisonique m=1, une bande de fréquences peut être choisie dans une plage de 500Hz à 10kHz pour l'ordre ambisonique m=2, une bande de fréquences peut être choisie dans une plage de 2000Hz à 9000Hz pour l'ordre ambisonique m=3, 15 une bande de fréquences peut être choisie dans une plage de 3000Hz à 7000Hz pour l'ordre ambisonique m=4. Ainsi, il s'avère qu'il est possible de définir un critère de validité sur la base de valeurs moyennes des bandes de fréquence pour chaque ordre ambisonique, même si les données du microphone ambisonique utilisé pour la captation de composantes ambisoniques sont inaccessibles. 20 Selon une réalisation particulière, le traitement de la matrice de décodage ambisonique comporte : - une inversion de la matrice élaborée de décodage ambisonique, pour obtenir une matrice de mélange dont : * les lignes correspondent à des canaux ambisoniques respectifs, et * les colonnes correspondant à des sources sonores, 25 - un traitement de la matrice de mélange pour extraire, par réduction de dimension de matrice, une pluralité de sous-matrices de mélange chacune associée à un ordre ambisonique et à une bande de fréquences choisie, et - une inversion des sous-matrices de mélange pour obtenir respectivement lesdites sous-matrices de décodage ambisonique. 3060830 8 On comprend ainsi qu'un filtrage fréquentiel des composantes d'ordre m=4 entre 4000 à 6000 Hz, dans l'exemple de la figure 5, permet de construire une sous-matrice, en particulier de mélange (matrice notée A ci-avant), à N=(m+1)2=25 lignes, en retenant les 25 premiers canaux ambisoniques. Néanmoins, à cet effet, il est préférable que le signal ambisonique soit suffisamment représenté dans 5 cette bande de fréquences 4-6 kHz, comme on le verra plus loin. Par ailleurs, si le signal ambisonique est bien représenté aussi dans les basses fréquences, par exemple entre 100 et 200Hz, on peut construire en outre une sous-matrice pour l'ordre m=1 par exemple, à N=4 lignes. On peut donc obtenir finalement une pluralité de sous-matrices de mélange, chacune associée à un ordre ambisonique m, et comportant chacune un nombre de lignes correspondant à un nombre de canaux 10 ambisoniques valides pour cet ordre m et dans la bande de fréquences à laquelle cette sous-matrice est associée. Dans une réalisation, le traitement du contenu ambisonique est mené pour une séparation de source et ladite matrice de décodage est une matrice de séparation aveugle de sources élaborée à partir des composantes ambisoniques. 15 Par exemple, la matrice de séparation peut être élaborée à partir des composantes ambisoniques filtrées à une bande de fréquences choisie et préférentiellement dans laquelle le nombre de canaux ambisoniques valides selon le critère précité est maximum. Ainsi, les canaux sont retenus pour une exactitude de représentation à un tel ordre ambisonique le plus élevé, mais aussi pour conserver un maximum de canaux représentés correctement dans cette 20 bande de fréquences, à des ordres ambisoniques moins élevés. Dans cette réalisation, on peut simplifier des sous-matrices de mélange avant leur inversion, par réduction d'un nombre de colonnes de chaque sous-matrice, les colonnes restantes des sous-matrices étant choisies de manière à conserver des signaux de plus grandes énergies après application des sous-matrices de décodage. 25 En effet, conserver les signaux de plus grande énergie permet de mieux représenter, et donc de mieux restituer, le champ sonore. En complément ou en variante, on peut choisir de privilégier des signaux extraits les plus décorrélés, ou les plus indépendants suivant un critère d'indépendance choisi. Ainsi, dans cette réalisation, on simplifie des sous-matrices de mélange avant leur inversion, par 30 réduction d'un nombre de colonnes de chaque sous-matrice, les colonnes restantes des sous- 3060830 9 matrices étant choisies de manière à conserver des signaux les moins corrélés après application des sous-matrices de décodage. Par ailleurs, dans un environnement réverbérant, le signal est constitué des champs directs issus de la propagation équivalente « champ libre » de chaque source et de réflexions sur des parois de 5 l'environnement acoustique. Ainsi, dans une réalisation alternative ou complémentaire, on simplifie des sous-matrices de mélange avant leur inversion, par réduction d'un nombre de colonnes de chaque sous-matrice, les colonnes restantes des sous-matrices étant choisies de manière à conserver des signaux correspondant à des champs sonores directs après application des sous-matrices de décodage. 10 Bien entendu, dans une réalisation où le traitement du contenu ambisonique est mené pour une restitution ambisonique sur une pluralité de haut-parleurs, la matrice de décodage précitée peut être une matrice inverse de positions spatiales relatives des haut-parleurs. Dans une réalisation illustrée plus loin en référence à la figure 9, le procédé comprend en particulier, pour un contenu ambisonique décomposé en sous-bandes de fréquences, une application de sous- 15 matrices de décodage, obtenues par : - Pour chaque ordre ambisonique du contenu, une détermination d'une bande de fréquences sur laquelle ledit ordre respecte un critère de validité prédéterminé d'encodage ambisonique, - Sur la base desdites bandes de fréquences, une application d'un banc de filtres au contenu ambisonique pour produire une pluralité de signaux en sous-bandes, de dimensions variables 20 correspondant à des canaux ambisoniques valides dans cette sous-bande, - Une détermination d'une matrice de décodage de taille maximale dans la bande de fréquence de l'ordre ambisonique maximal et d'une matrice de mélange associée, inverse ou pseudo-inverse de ladite matrice de décodage, - Pour chaque autre bande de fréquences, une détermination d'une matrice de mélange de taille 25 réduite, sous-matrice de ladite matrice de mélange, et d'une sous-matrice de séparation, inverse ou pseudo-inverse de ladite sous-matrice de mélange, - Une reconstruction des signaux séparés pleine-bande par application d'un banc de filtre de synthèse aux signaux séparés issus de la multiplication desdits signaux par lesdites matrices. La présente invention vise aussi un programme informatique comportant des instructions pour la 30 mise en oeuvre du procédé lorsque ce programme est exécuté par un processeur. Un exemple 3060830 10 d'ordinogramme de l'algorithme général d'un tel programme est illustré sur la figure 7 commentée ci-après, laquelle est précisée dans les figures 8 et 9. La présente invention vise aussi un dispositif informatique comportant : - une interface d'entrée pour recevoir des signaux de composantes ambisoniques, 5 - une interface de sortie pour délivrer des signaux décodés, associés chacun à une source sonore, - et un programme informatique pour la mise en oeuvre du procédé. Un exemple d'un tel dispositif est illustré sur la figure 10 commentée plus loin. La présente invention propose ainsi d'utiliser la formation de voies à partir d'un encodage 10 ambisonique réel en tirant parti, dans chaque bande de fréquences, de tous les canaux dont la directivité respecte le formalisme ambisonique. Une forme de réalisation présentée ci-avant permet alors de déterminer une ou plusieurs matrices de mélange Ak, correspondant à des sous-matrices obtenues à partir de la matrice théorique A, et formulées chacune dans une bande de fréquences, puis inversée pour donner des matrices de décodage Bk. 15 Ainsi, l'invention offre un traitement générique d'un contenu ambisonique quelconque, et notamment réel, possiblement affecté par des limitations physiques d'un système d'enregistrement, et ce sans aucune contrainte visant à limiter la bande-passante totale des sources extraites. D'autres avantages et caractéristiques de l'invention apparaitront à la lecture de la description détaillée ci-après d'exemples de réalisation de l'invention, et à l'examen des dessins annexés sur 20 lesquels : - la figure 1 illustre une base de fonctions harmoniques sphériques d'ordre 0 (première ligne) à 3 (dernière ligne), avec en gris clair les valeurs positives, et en gris foncé les valeurs négatives, - la figure 2 illustre un système d'encodage ambisonique à partir d'un microphone sphérique, - la figure 3 illustre la formation de voies pour l'extraction de trois composantes, pour différents 25 ordres ambisoniques, - la figure 4 illustre très schématiquement un système de décodage ambisonique à partir de composantes ambisoniques, - la figure 5 illustre la corrélation entre un encodage ambisonique idéal et un encodage réel, 3060830 11 -la figure 6 illustre la directivité dans le plan horizontal, mesurée pour un encodage ambisonique réel (avec de gauche à droite successivement les composantes des ordres 0, 1, 2 et 3), - la figure 7 illustre les principales étapes d'un exemple de procédé au sens de l'invention, - la figure 8 illustre les étapes d'un mode de réalisation particulier du procédé selon l'invention, 5 - la figure 9 est un schéma-bloc d'un algorithme de traitement correspondant au mode de réalisation illustré sur la figure 7, et - la figure 10 illustre schématiquement un dispositif possible pour la mise en oeuvre de l'invention. Le schéma d'ensemble d'un procédé de traitement ambisonique global au sens de l'invention est présenté figure 7. Il s'agit par exemple d'un procédé de décodage ambisonique. On entend par les 10 termes « décodage ambisonique » aussi bien la fourniture de signaux décodés par exemple destinés à alimenter des haut-parleurs respectifs pour une restitution ambiophonique, qu'une fourniture, de façon plus générale, de signaux associés chacun à une source sonore, notamment dans la technique de séparation de sources. A l'étape S1, on dispose d'un contenu ambisonique x(t) comprenant une pluralité de composantes 15 ambisoniques CA, d'ordres successifs m=0, 1, ..., M (avec par exemple M=4) et, issu d'un enregistrement, ou d'une « captation », par au moins un microphone ambisonique MIC. Un microphone ambisonique est un microphone composé d'une pluralité de capsules microphoniques généralement réparties de manière sphérique et de manière la plus régulière possible. Ces capsules jouent le rôle de capteurs de signaux sonores. Les capsules microphoniques sont agencées sur le 20 microphone ambisonique de manière à capter des signaux sonores selon leur directivité dans l'espace. Comme illustré sur la figure 5, l'ensemble des capsules formant un tel microphone ambisonique peut acquérir différentes composantes ambisoniques à des ordres ambisoniques jusqu'à M, mais l'exactitude de la représentation ambisonique pour ces différents ordres n'est pas réellement respectée pour toutes les fréquences du spectre audio entre 0 et 20kHz. Néanmoins, 25 l'invention propose ici d'isoler certaines fréquences du spectre pour lesquelles les composantes ambisoniques, pour des ordres donnés, sont exactes (comme par exemple dans la plage de fréquences entre 4000 et 6000Hz pour l'ordre m=4 sur la figure 5, ou plus largement la plage entre 2000Hz et 9000 Hz pour l'ordre m=3, etc.). Néanmoins, les variations fréquentielles de l'exactitude de représentation ambisonique de chaque 30 ordre de la figure 5 sont obtenues pour un microphone particulier ayant des dimensions et un 3060830 12 nombre donné de capsules. Ainsi, pour un autre microphone, d'autres variations spectrales peuvent être attendues. L'étape S2 vise donc à récupérer les données caractérisant le microphone ambisonique MIC (et éventuellement les conditions de captation du contenu ambisonique c(t), et/ou encore les conditions 5 de réverbération pendant la captation, ou autres). Plus généralement, une donnée caractérisante du microphone ambisonique MIC peut être l'espacement inter-capsules. En effet, l'encodage des hautes fréquences se dégrade lorsque l'espacement inter-capteurs devient plus grand qu'une demi-longueur d'onde. Ceci est dû au phénomène de repliement spatial (ou « aliasing »). A l'inverse, pour un signal basse fréquence, des 10 capsules microphoniques trop rapprochées ne peuvent générer la directivité souhaitée. A l'étape S3, on peut appliquer un banc de filtre d'analyse BFA au contenu ambisonique x(t) afin de sélectionner ensuite, à l'étape S31, des signaux de composantes ambisoniques filtrés dans des plages de fréquences dans lesquelles la représentation ambisonique pour un ordre donné m est la plus exacte (respectant ainsi un « critère de validité » de la représentation ambisonique), et ce en 15 fonction des données du microphone définies ci-dessus. En fonction du type de traitement appliqué au contenu ambisonique x(t), entre un traitement de séparation de sources SAS ou un traitement en vue d'une restitution sur haut-parleurs RES, l'étape S4 vise l'obtention d'une matrice de décodage B, en fonction du type de traitement choisi. Dans le cas d'une restitution ambisonique sur haut-parleurs, la matrice de décodage B est l'inverse d'une 20 matrice A contenant des coefficients propres à des positions spatiales de haut-parleurs utilisés pour la restitution. Dans le cas d'une séparation de sources, la matrice de décodage B est élaborée initialement à l'étape S4 en vue d'un traitement de séparation aveugle de sources à partir des composantes ambisoniques filtrées et sélectionnées. Plus particulièrement, cette matrice de décodage B est élaborée pour la 25 bande de fréquences contenant le plus grand nombre de canaux ambisoniques valides (et le plus grand ordre susceptible d'être obtenu M). La détermination des bandes de fréquences de validité des différents ordres ambisoniques peut être adaptée au microphone ambisonique ayant servi à la captation des composantes ambisoniques à décoder. Pour ce faire, il est possible par exemple de se baser sur les variations fréquentielles de 30 l'exactitude de la représentation ambisonique pour différents ordres m, du type illustré sur la figure 5. 3060830 13 Plus généralement, il peut être déterminé encore une allure « moyenne » des variations fréquentielles de l'exactitude de la représentation ambisonique pour les différents ordres m pour différents modèles de microphones ambisoniques, et se servir de ces allures moyennes si ces données ne sont pas disponibles, au décodage. 5 A l'étape S7, on détermine au moins deux matrices B1, B2, issues d'une réduction matricielle de la matrice de décodage B pour chaque sous-bande de fréquences (dans l'exemple illustré les sous-bandes de fréquences fl et f2). Un exemple de réalisation plus précis de cette réduction matricielle sera décrit plus loin en référence à la figure 8. Puis, à l'étape S8, on effectue le produit de chaque matrice B1 et B2 obtenues à l'étape précédente par les signaux ambisoniques filtrés dans les sous- 10 bandes fl, f2 correspondantes. On obtient ainsi, dans chaque sous-bande k (k=1,2), un ensemble de signaux extraits sk. A l'étape S9, on combine les vecteurs de signaux extraits sl (1 pour k=1) et s2 (2 pour k=2) afin d'obtenir les signaux reconstruits pleine bande (par application par exemple d'un banc de filtre de synthèse). 15 La figure 8 illustre les étapes d'un mode de réalisation particulier du procédé selon l'invention. Plus précisément, la figure 8 présente des étapes du procédé qui peuvent être mises en oeuvre entre les étapes S4 et S7 de la figure 7. A l'étape S4, comme décrit ci-dessus, on obtient la matrice de décodage B définie plus haut. A l'étape S5 on peut effectuer une inversion de cette matrice de décodage B (ou de façon équivalente, une 20 détermination de sa pseudo-inverse) afin d'obtenir la matrice de mélange A correspondante (étape S51). Dans le cas d'une séparation de sources, la matrice de mélange A peut ainsi contenir des coefficients relatifs à des positions respectives de sources sonores à extraire. Dans le cas d'une restitution sur haut-parleurs, la matrice de mélange A peut contenir des coefficients relatifs à la position des haut-parleurs sur lesquels on souhaite restituer les signaux décodés. Plus précisément, 25 les lignes de la matrice de mélange A correspondent aux canaux ambisoniques successifs (définissant successivement les ordres m=0 à m=M, où M est l'ordre ambisonique maximum disponible) et ses colonnes correspondent aux sources ou aux haut-parleurs. A l'étape S6, on peut effectuer une diminution des dimensions de la matrice de mélange A, pour obtenir des sous-matrices Al, A2. Il s'agit d'une réduction matricielle dont le nombre de lignes 30 correspond aux nombres de canaux ambisoniques pour chaque ordre. Typiquement, si les signaux ambisoniques sont bien encodés dans la bande de 100 à 1000Hz, où l'ordre m=1 est bien respecté (au moins pour le microphone ambisonique de la figure 5), il est extrait déjà de la matrice A une 3060830 14 sous-matrice Al à N=4 lignes associée à l'ordre m=1 et à la bande de fréquences 100-1000Hz. Ensuite, si les signaux ambisoniques sont bien représentés dans la bande de 1000 à 10 000Hz, où l'ordre m=2 est bien respecté, il est extrait ensuite de la matrice A une matrice A2 à N=9 lignes et associée à l'ordre m=2 et à la bande de fréquences 1000-10 000Hz, et ainsi de suite. Le nombre de 5 sous-matrices dépend ainsi de l'ordre du contenu ambisonique x(t) dont les composantes sont retenues comme valides à l'étape S31. Chaque sous-matrice correspond alors à une bande de fréquences, et peut ainsi contenir un nombre de lignes correspondant au nombre de canaux valides pour cette bande de fréquences. Plus précisément, comme illustré sur la figure 8, pour chaque sous-bande, on identifie le nombre de canaux valides correspondant. Par exemple, pour une sous-bande 10 f1 choisie pour l'ordre m=1 du contenu ambisonique x(t), on extrait une matrice Al comportant quatre lignes (N1=(m+1)2) correspondant aux quatre canaux ambisoniques à l'ordre 1, et le nombre de « sources » (sources à extraire ou haut-parleurs) en colonnes. Comme illustré sur la figure 8, les quatre lignes retenues pour la construction de la sous-matrice Al sont les coefficients de la matrice initiale globale A : 15 - C11, C12, C13, - C21, C22, C23, - C31, C32, C33, et - C41, C42, C43. Concernant la sous-matrice A2, ces lignes de la matrice globale A peuvent être reprises, ainsi que les 20 suivantes, jusqu'à la ligne : - C91, C92, C93. Pour la matrice de mélange A2, correspondant à l'ordre 2 du contenu ambisonique x(t), et donc à la sous-bande f2, on conserve donc neuf lignes, correspondant aux neuf canaux de l'ordre 2, et le nombre de sources à extraire en colonnes. 25 Chaque sous-matrice de mélange ainsi obtenue est de dimension N x Ntarget, avec Ntarget le nombre de sources issues de la séparation de sources aveugle ou le nombre de haut-parleurs prévus pour une restitution. Dans le cas d'une restitution sur haut-parleurs, le nombre de haut-parleurs est préférentiellement égal ou supérieur au nombre de lignes. Par exemple, pour la matrice de mélange Al de quatre lignes, 30 on peut ne conserver qu'un jeu de quatre colonnes. Dans le cas d'une séparation de sources, le 3060830 15 nombre de colonnes peut être inférieur ou égal au nombre de lignes. Par exemple, pour la matrice de mélange Al de quatre lignes, on peut supprimer des colonnes et garder par exemple des sources dont les signaux sont de plus grandes énergies et/ou celles qui sont les moins corrélés (sources les moins « mélangées » possibles) et/ou les signaux correspondent au champ direct des sources, ou 5 autres. A l'étape S71 on effectue une inversion de chaque sous-matrice de mélange Al, A2 afin d'obtenir respectivement les sous-matrices de décodage Bl, B2 présentées plus haut (étape S7). Le passage par la matrice de mélange A permet en particulier de conserver des niveaux satisfaisants d'énergie des composantes ambisoniques liées à chaque ordre, malgré les réductions matricielles. En d'autres 10 termes, les étapes S5 à S71 permettent « d'affiner » le décodage du contenu ambisonique x(t). La figure 9 est un schéma-bloc d'un algorithme de traitement correspondant au mode de réalisation illustré sur les figures 7 et 8. On a repris les mêmes références d'étapes Si, S2, etc., pour désigner des étapes identiques ou similaires et présentées ci-avant en référence aux figures 7 et 8. On nomme « canaux » les signaux microphoniques ambisoniques et « sources » les signaux à extraire 15 (sources effectivement à extraire ou les signaux d'alimentation des haut-parleurs). A l'étape Si, on dispose d'un contenu ambisonique x(t) d'ordre M, comprenant une pluralité de canaux ambisoniques N enregistrés à traiter. De manière générale, le nombre de canaux ambisoniques enregistrés est égale à N=(M+1)2. A l'étape S2, on dispose des données relatives à la captation ambisonique du contenu x(t) (données relatives au microphone ambisonique MIC utilisé, etc.). 20 Connaissant les limites de validité de l'encodage microphonique, on détermine une bande de fréquences pour chaque ordre ambisonique. Un banc de filtre permettant une reconstruction est appliqué aux N canaux ambisoniques à l'étape S3, pour donner K sous-bandes notées xk. Les sous-bandes sont choisies pour correspondre aux différentes plages de validité de l'encodage microphonique. 25 Dans une réalisation particulière à l'étape S4A illustrée en trait plein, on utilise une matrice B de séparation de sources élaborée en fonction des composantes ambisoniques filtrées en fréquence (flèche de dessus venant sur le rectangle S4A). Plus particulièrement, une méthode de séparation aveugle de sources est appliquée dans la sous-bande contenant le plus de canaux valides, pour obtenir une matrice de séparation B de dimensions Ntarget x N, Ntarget étant le nombre de sources 30 obtenues par la méthode de séparation aveugle dans la sous-bande de fréquences choisie. Les canaux valides sont déterminés à partir d'un critère de validité relatif à chaque ordre du contenu ambisonique x(t) en fonction de chaque bande de fréquences du banc de filtres. Plus généralement, 3060830 16 afin de maximiser la qualité de la séparation de sources, on choisit une bande de fréquences comprenant le plus de composantes ambisoniques valides. On entend par « valides » des composantes dont les critères énergétiques ou la directivité n'ont pas été biaisés lors de la captation ambisonique, comme présenté ci-avant en référence à la figure 5. La validité de chaque ordre dans 5 des bandes de fréquences du domaine audio peut être établie en connaissant les limites du microphone ambisonique utilisé lors de la captation du contenu ambisonique x(t), ou encore à l'aide d'un abaque établi sur la base de mesures effectuées sur une pluralité de microphones ambisoniques, permettant de réaliser une moyenne de la validité de chaque ordre ambisonique dans chaque bande de fréquences. 10 Par exemple, les canaux ambisoniques d'ordre 1 ont tendance à être valides dans une bande de fréquences allant de 100HZ à environ 10kHz. La bande de fréquences dans laquelle les canaux ambisoniques d'ordre 2 peuvent être plus généralement valides peut par exemple aller de 1kHz à 9kHz, etc. Dans une réalisation variante en vue d'une restitution d'une scène sonore sur plusieurs haut-parleurs 15 (plus de deux en général), à l'étape S4B (illustrée par les traits en pointillés sur la figure 9, pour désigner cette variante), la matrice de décodage est construite en fonction de la position des haut-parleurs sur lesquels le contenu doit être restitué. Plus exactement, cette matrice B de décodage correspond à l'inverse d'une matrice de mélange A qui est définie par les positions spatiales respectives des haut-parleurs. 20 En revenant au traitement général (pour une restitution ou pour une séparation de sources), à l'étape S5, la matrice de mélange « théorique » A (pour les deux variantes précitées) est construite par inversion de B. Pour la séparation de sources, la matrice de mélange est composée de N lignes et de Ntarget colonnes, la ième colonne contenant les coefficients harmoniques sphériques, relatifs aux coordonnées (01,cp1) de la source s,. Ci-dessous se trouve un exemple de matrice de mélange A dans 25 le cas d'une séparation de sources pour un contenu ambisonique d'ordre 2 composé de cinq sources : 1 1 1 1 L Ordre 0 Yit(3) Y111(4) Y22(5) Y2'2-(1) 3060830 17 Pour la diffusion sur haut-parleurs, A est composée de N lignes et d'un minimum de N colonnes, la ième colonne contenant les coefficients harmoniques sphériques, relatifs aux coordonnées (Bi, cpi) 5 du haut-parleur i. A l'étape S6, et pour chaque sous-bande k, une sous-matrice de mélange Ak est construite, telle que Ak est une version tronquée de la matrice A, en ne conservant que les Nk lignes correspondant aux canaux effectivement valides dans cette sous-bande k. Pour la séparation de sources, si Nk est inférieur au nombre de sources Ntarget recherchées dans la 10 sous-bande, on ne conserve qu'un jeu de Ntarget,k, colonnes (avec Ntarget,k inférieur ou égal à Nk), choisies suivant des critères énergétiques (par exemple en séparant les sources ayant la plus grande contribution) ou suivant d'autres critères d'intérêt tels que définis précédemment. La matrice Ak a ainsi pour dimensions Nk x Ntarget,k, avec Ntarget,k = min(Nk, Ntarget) par exemple. Ci-dessous se trouve un exemple d'une matrice Ak(4x4) tronquée à l'ordre 1 ambisonique : 15 Composantes conservées (1) Yi11(2) (3) Y 1(4) Y11(5) non valide _Y,1, (1) Y212(5) Pour la restitution sur haut-parleurs, un jeu de Nk haut-parleurs est sélectionné pour la restitution, et Ak a donc pour dimensions Nk x Nk. 20 A l'étape S7, la matrice Ak est inversée pour donner Bk. Lorsque la sous-matrice Ak n'est pas une matrice carrée, une infinité de possibilités existe pour l'inversion. Une pseudo-inversion peut être appliquée, ou encore une inversion en appliquant des contraintes supplémentaires (par exemple choix de la solution donnant le beamforming le plus directif, ou minimisant les lobes secondaires). De manière générale, on entend par « inversion de matrice », aussi bien une inversion classique de 25 matrice, qu'une pseudo-inversion comme présenté ci-avant. A , I- Ordre - Ordre 1 valide 3060830 18 Puis, à l'étape S8, Bk est appliquée à la sous-bande xk pour obtenir les signaux sk tels que sk = Bk. xk Une fois que des sources ont été extraites dans chaque sous-bande, les signaux pleine-bande correspondants sont reconstruits par un filtre de synthèse à partir des signaux de sous-bandes de 5 même direction, à l'étape S9. Ci-dessous, un exemple de mise en oeuvre du procédé selon un mode de réalisation particulier de l'invention est décrit à titre d'exemple. On dispose d'un contenu ambisonique d'ordre 2 (9 canaux) échantillonné à 16kHz, noté x(t) constitué de 3 sources que l'on veut extraire. L'encodage ambisonique aux ordres 0 et 1 est valide entre 200Hz 10 et 8000Hz. L'encodage de l'ordre 2 est valide entre 900Hz et 8000Hz. Un banc de filtre est implémenté, constitué de deux bandes de fréquences, 200Hz-900Hz (jusqu'à l'ordre 1) et 900Hz-8000Hz (utilisation de l'ordre 2) Le banc de filtre est appliqué à x(t), pour former x1(t) et x2(t). x1(t) est constitué de 4 canaux (ambisonie d'ordre 1) et x2(t) contient 9 canaux (ambisonie d'ordre 2). 15 Une matrice de séparation B de dimensions 3x9 est estimée par analyse en composantes indépendantes effectuée dans la sous-bande 900Hz-8000Hz c'est-à-dire x2(t). Une matrice de mélange théorique A, de dimensions 9x3, est déduite par inversion de B, chaque colonne i contenant les coefficients harmoniques sphériques de la source i. Dans le même temps, les matrices Al et A2 sont calculées à partir de A pour extraire les sources dans 20 chaque sous-bande : - Al contient uniquement les coefficients jusqu'à l'ordre 1 pour les trois sources, soit : Al= A (les quatre premières lignes, les trois premières colonnes), - A2 contient les coefficients relatifs aux neufs canaux pour les trois sources, on a donc : A2=A Al et A2 sont inversées pour former les matrices de séparation B1 et B2. 25 Les trois sources sont extraites dans chaque sous-bande d'indices respectifs 1 et 2: s1=B1.xl et s2=B2.x2 3060830 19 Puis, les sources pleine-bande sont reconstituées par application du filtre de synthèse aux signaux en sous-bandes s1 et s2, par exemple par sommation bandes à bandes (si le banc de filtres d'analyse a opéré en bande de base) : s= s1 + s2 5 En référence à la figure 10, la présente invention vise en outre un dispositif DIS pour la mise en oeuvre de l'invention. Ce dispositif DIS peut comporter une interface d'entrée IN pour recevoir des signaux ambisoniques x(t). Le dispositif DIS peut comprendre une mémoire MEM pour stocker des instructions d'un programme informatique au sens de l'invention. Les instructions du programme informatique sont des instructions de traitement des signaux ambisoniques x(t). Elles sont mises en 10 oeuvre par un processeur PROC, afin de délivrer, via une interface de sortie OUT, des signaux décodés s(t). Bien entendu, la présente invention ne se limite pas aux formes de réalisation décrites ci-avant à titre d'exemple ; elle s'étend à d'autres variantes. Typiquement, les plages de fréquences pour lesquelles la représentation ambisonique est valide sont 15 données ci-avant à titre d'exemple et peuvent différer selon la nature du ou des microphones ambisoniques utilisés pour la captation, voire des conditions de captation elles-mêmes. REVENDICATIONS1. Procédé, mis en oeuvre par des moyens informatiques, de traitement d'un contenu ambisonique comportant une pluralité de composantes ambisoniques d'une pluralité d'ordres définissant une succession de canaux ambisoniques dans chacun desquels est représentée une composante ambisonique, le procédé comportant : - un filtrage fréquentiel des composantes ambisoniques dans une pluralité de bandes de fréquences, - une élaboration d'une matrice (B) de décodage ambisonique, - un traitement de la matrice (B) de décodage ambisonique pour extraire, par réduction de dimension de matrice, une pluralité de sous-matrices (B1, B2) de décodage ambisonique chacune associée à un ordre ambisonique et à une bande de fréquences choisie pour cet ordre ambisonique, - des applications respectives des sous-matrices de décodage aux composantes ambisoniques dans chaque bande de fréquences choisie, et une reconstruction bandes-à- bandes des résultats desdites applications respectives, pour délivrer une pluralité de signaux décodés, associés chacun à une source sonore. 2. Procédé selon la revendication 1, dans lequel chaque sous-matrice est associée à une bande de fréquences choisie en fonction d'un critère de validité des composantes ambisoniques de l'ordre auquel est associée ladite sous-matrice, dans ladite bande de fréquence choisie. 3. Procédé selon la revendication 2, dans lequel le critère de validité des composantes est défini par des conditions de captation desdites composantes ambisoniques, par au moins un microphone ambisonique. 4. Procédé selon la revendication 3, comportant : - une réception de données d'au moins un microphone ambisonique utilisé pour capter lesdites composantes ambisoniques ; - une détermination des bandes de fréquences choisies pour construire lesdites sous- matrices (B1, B2), en fonction desdites données de microphone ambisonique. 5. Procédé selon l'une des revendications précédentes, dans lequel, chaque sous-matrice de décodage ambisonique (B1, B2) étant associée à un ordre ambisonique et à une bande de fréquences choisie pour cet ordre ambisonique, 3060830 21 une bande de fréquences est choisie dans une plage de 100Hz à 10kHz pour l'ordre ambisonique m=1, une bande de fréquences est choisie dans une plage de 500Hz à 10kHz est choisie pour l'ordre ambisonique m=2, 5 une bande de fréquences est choisie dans une plage de 2000Hz à 9000Hz est choisie pour l'ordre ambisonique m=3, une bande de fréquences est choisie dans une plage de 3000Hz à 7000Hz est choisie pour l'ordre ambisonique m=4. 10 6. Procédé selon l'une des revendications précédentes, dans lequel le traitement de la matrice (B) de décodage ambisonique comporte : - une inversion de la matrice élaborée (B) de décodage ambisonique, pour obtenir une matrice (A) de mélange dont : 15 * les lignes correspondent à des canaux ambisoniques respectifs, et * les colonnes correspondant à des sources sonores, - un traitement de la matrice de mélange (A) pour extraire, par réduction de dimension de matrice, une pluralité de sous-matrices de mélange (Al, A2) chacune associée à un ordre ambisonique et à une bande de fréquences choisie, et 20 - une inversion des sous-matrices de mélange (Al, A2) pour obtenir respectivement lesdites sous-matrices (B1, B2) de décodage ambisonique. 7. Procédé selon l'une des revendications précédentes, dans lequel le traitement du contenu ambisonique est mené pour une séparation de source et ladite matrice de décodage (B) est 25 une matrice de séparation aveugle de sources élaborée à partir des composantes ambisoniques (S4A). 8. Procédé selon la revendication 7, prise en combinaison avec la revendication 2, dans lequel la matrice de séparation (B) est élaborée à partir des composantes ambisoniques filtrées à 30 une bande de fréquences choisie et dans laquelle le nombre de canaux ambisoniques valides selon ledit critère est maximum. 3060830 22 9. Procédé selon l'une des revendications 7 et 8, prises en combinaison avec la revendication 6, comportant en outre une simplification des sous-matrices de mélange (Al, A2) avant leur inversion, par réduction d'un nombre de colonnes de chaque sous-matrice, les colonnes restantes des sous-matrices étant choisies de manière à conserver des signaux de plus 5 grandes énergies après application des sous-matrices de décodage. 10. Procédé selon l'une des revendications 7 à 9, prises en combinaison avec la revendication 6, comportant en outre une simplification des sous-matrices de mélange (Al, A2) avant leur inversion, par réduction d'un nombre de colonnes de chaque sous-matrice, les colonnes 10 restantes des sous-matrices étant choisies de manière à conserver des signaux les moins corrélés après application des sous-matrices de décodage. 11. Procédé selon l'une des revendications 7 à 10, prises en combinaison avec la revendication 6, comportant en outre une simplification des sous-matrices de mélange (Al, A2) avant leur 15 inversion, par réduction d'un nombre de colonnes de chaque sous-matrice, les colonnes restantes des sous-matrices étant choisies de manière à conserver des signaux correspondant à des champs sonores directs après application des sous-matrices de décodage. 20 12. Procédé selon l'une des revendications 1 à 6, dans lequel le traitement du contenu ambisonique est mené pour une restitution ambisonique sur une pluralité de haut-parleurs et ladite matrice de décodage (B) est une matrice inverse de positions spatiales relatives des haut-parleurs (S4B). 25 13. Procédé selon l'une des revendications précédentes, comportant, pour un contenu ambisonique (x) décomposé en sous-bandes de fréquences (k), une application de sous-matrices de décodage (Bk), obtenues par : - Pour chaque ordre ambisonique du contenu, une détermination d'une bande de fréquences sur laquelle ledit ordre respecte un critère de validité prédéterminé d'encodage 30 ambisonique, - Sur la base desdites bandes de fréquences, une application d'un banc de filtres au contenu ambisonique (x) pour produire une pluralité de signaux en sous-bandes (xk), de dimensions variables correspondant à des canaux ambisoniques valides dans cette sous-bande (k), 3060830 23 - Une détermination d'une matrice de décodage (B) de taille maximale dans la bande de fréquence de l'ordre ambisonique maximal et d'une la matrice de mélange associée (A), inverse ou pseudo-inverse de ladite matrice de décodage (B), - Pour chaque autre bande de fréquences (k), une détermination d'une matrice de mélange 5 (Ak) de taille réduite, sous-matrice de ladite matrice de mélange (A), et d'une sous-matrice de décodage (Bk), inverse ou pseudo-inverse de ladite sous-matrice de mélange (Ak), - Une reconstruction des signaux séparés pleine-bande (s) par application d'un banc de filtre de synthèse aux signaux séparés (sk) issus de la multiplication desdits signaux (xk) par lesdites matrices (Bk). 10 14. Programme informatique caractérisé en ce qu'il comporte des instructions pour la mise en oeuvre du procédé selon l'une des revendications 1 à 13, lorsque ce programme est exécuté par un processeur. 15 15. Dispositif informatique comportant : - une interface d'entrée pour recevoir des signaux de composantes ambisoniques, - une interface de sortie pour délivrer des signaux décodés, associés chacun à une source sonore, - et un circuit de traitement pour la mise en oeuvre du procédé selon l'une des 20 revendications 1 à 13.