Brevet FR2836571A1 — publié le 20030829

La présente invention concerne un procédé et un dispositif de pilotage d'un ensemble de restitution d'un champ acoustique. Le son est un phénomène acoustique ondulatoire qui évolue dans le temps et dans l'espace. Les techniques existantes agissent principalement sur l'aspect temporel des sons, le traitement de l'aspect spatial étant très incomplet. En effet, les systèmes de restitution de haute qualité existants imposent de fait une configuration spatiale prédéterminée de l'ensemble de restitution. Par exemple, les systèmes dits multicanaux adressent des signaux différents et prédéterminés à plusieurs haut-parleurs dont la répartition est fixée et connue. De même, les systèmes dits ambisonic , qui considèrent la direction de provenance des sons qui parviennent à un auditeur, requièrent un ensemble de restitution dont la configuration doit respecter certaines règles de positionnement. Dans ces systèmes, l'environnement sonore est assimilé à une distribution angulaire de sources sonores autour d'un point, correspondant à la position d'écoute. Les signaux correspondent à une décomposition de cette distribution sur une base de fonctions de directivité appelées harmoniques sphériques. Dans l'état actuel de développement de ces systèmes, une restitution de bonne qualité n'est possible qu'avec une répartition sphérique des hautparleurs et une répartition angulaire sensiblement régulière. Ainsi, lorsque les techniques existantes sont mises en oeuvre avec un ensemble de restitution dont la répartition spatiale est quelconque, la qualité de restitution est grandement détériorée, notamment en raison de distorsions angulaires. Des développements techniques récents permettent de considérer une modélisation dans le temps et dans les trois dimensions de l'espace d'un champ acoustique plutôt que la distribution angulaire de l'environnement sonore. Notamment, la thèse de doctorat Représentation de champs acoustiques, application à la transmission et à la reproduction de scènes so- nores complexes dans un contexte multimédia Université Paris VI, Jérôme Daniel, du 11 juillet 2000, définit des fonctions décrivant les caractéristiques ondulatoires d'un champ acoustique et permettant une décomposition sur une base de fonctions de l'espace et du temps qui décrit complètement un champ acoustique tridimensionnel. Cependant, dans ce document, les solutions théoriques s'inspirent des systèmes dits Ambisonic et une restitution de haute qualité ne peut être obtenue que pour les 5 répartitions sphériques régulières existantes. Aucun élément ne permet d'assurer une restitution de haute qualité à partir d'une configuration spatiale quelconque de l'ensemble de restitution. Il apparaît donc qu'aucun système de l'art antérieur ne permet d'effectuer une restitution de qualité à partir d'une configuration spatiale quelconque de l'ensemble de restitution. Le but de l'invention est de remédier à ce problème en fournissant un procédé et un dispositif de détermination de signaux de pilotage d'un ensemble de restitution d'un champ acoustique dont la configuration spatiale est quelconque. L'invention a pour objet un procédé de pilotage d'un ensemble de restitution d'un champ acoustique pour obtenir un champ acoustique restitué de caractéristiques spécifiques sensiblement indépendantes des caractéristiques intrinsèques de restitution dudit ensemble, ledit ensemble de restitution comportant une pluralité d'éléments de restitution, caractérisé en ce qu'il comporte au moins : - une étape d'établissement d'un nombre fini de coefficients représentatifs de la distribution dans le temps et dans les trois dimensions de l'espace dudit champ acoustique à restituer ; - une étape de détermination de filtres de reconstruction représentatifs dudit ensemble de restitution, comportant une sous-étape de prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution ; - une étape de détermination d'au moins un signal de pilotage desdits éléments dudit ensemble de restitution, ledit au moins un signal étant obtenu par l'application, auxdits coefficients, desdits filtres de reconstruction ; et - une étape de délivrance dudit au moins un signal de pilotage, en vue d'une application auxdits éléments de restitution afin d'engendrer ledit champ acoustique restitué par ledit ensemble de restitution. Suivant d'autres caractéristiques : - ladite étape d'établissement d'un nombre fini de coefficients représentatifs de la distribution dudit champ acoustique à restituer comporte : - une étape consistant à fournir un signal d'entrée comportant des informations temporelles et spatiales d'un environnement sonore ; et - une étape de mise en forme dudit signal d'entrée par décomposition desdites informations sur une base de fonctions spatio-temporelles, cette étape de mise en forme permettant de délivrer une représentation dudit champ acoustique à restituer correspondant audit environnement sonore sous la forme d'une combinaison linéaire desdites fonctions ; - ladite étape d'établissement d'un nombre fini de coefficients représentatifs de la distribution dudit champ acoustique à restituer comporte : - une étape consistant à fournir un signal d'entrée comportant un nombre fini de coefficients représentatifs dudit champ acoustique à restituer sous la forme d'une combinaison linéaire de fonctions spatio-temporelles ; - lesdites fonctions spatio-temporelles sont des fonctions dites de Fourier-Bessel et/ou des combinaisons linéaires de ces fonctions ; - ladite sous-étape de prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution est réalisée au moins à partir de paramètres représentatifs, pour chaque élément, des trois coordonnées de sa position par rapport au centre placé dans la zone d'écoute, et/ou de sa réponse spatiotemporelle ; - ladite sous-étape de prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution est réalisée en outre à partir : - de paramètres décrivant, sous forme de coefficients de pondération, une fenêtre spatiale qui spécifie la répartition dans l'espace de contraintes de reconstruction du champ acoustique ; et - d'un paramètre décrivant un ordre de fonctionnement limitant le nombre de coefficients à prendre en compte lors de ladite étape de détermination de filtres de reconstruction ; - ladite sous-étape de prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution est réalisée en outre à partir : - de paramètres constituant une liste de fonctions spatiotemporelles dont la reconstruction est imposée ; et - d'un paramètre décrivant un ordre de fonctionnement limitant le nombre de coefficients à prendre en compte lors de ladite étape de détermination de filtres de reconstruction ; - ladite étape de prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution est réalisée en outre au moins à partir d'un des paramètres choisis dans le groupe constitué : - de paramètres représentatifs d'au moins une des trois coordonnées de la position de chaque ou certains des éléments, par rapport au centre placé dans la zone d'écoute ; - de paramètres représentatifs des réponses spatio-temporelles de chaque ou certains des éléments ; - d'un paramètre décrivant un ordre de fonctionnement limitant le nombre de coefficients à prendre en compte lors de ladite étape de détermination de filtres de reconstruction ; - de paramètres constituant une liste de fonctions spatiotemporelles dont la reconstruction est imposée ; - de paramètres représentatifs des gabarits desdits éléments de restitution ; - d'un paramètre représentatif de la capacité d'adaptation locale souhaitée à l'irrégularité spatiale de la configuration dudit ensemble de restitution ; - d'un paramètre définissant le modèle de rayonnement desdits éléments de restitution ; - de paramètres représentatifs de la réponse en fréquence desdits éléments de restitution ; - d'un paramètre représentatif d'une fenêtre spatiale ; - de paramètres représentatifs d'une fenêtre spatiale sous forme de coefficients de pondération ; et - d'un paramètre représentatif du rayon d'une fenêtre spatiale lorsque celle-ci est une boule ; - le procédé comporte une étape de calibrage permettant de délivrer tout ou partie des paramètres utilisés dans ladite étape de détermination de filtres de reconstruction ; - ladite étape de calibrage comporte, pour au moins l'un des éléments de restitution : - une sous-étape d'acquisition de signaux représentatifs du rayonnement dudit au moins un élément dans le lieu d'écoute ; et - une sous-étape de détermination de paramètres spatiaux et/ou acoustiques dudit au moins un élément ; - ladite étape de calibrage comporte : - une sous-étape d'émission d'un signal spécifique vers ledit au moins un élément dudit ensemble de restitution, ladite sous-étape d'acquisition correspondant à l'acquisition de l'onde sonore émise en réponse par ledit au moins un élément ; et - une sous-étape de transformation desdits signaux acquis en un nombre fini de coefficients représentatifs de l'onde sonore émise, afin de permettre la réalisation de ladite sous-étape de détermination de paramètres spatiaux et/ou acoustiques ; - ladite sous-étape d'acquisition correspond à une sous-étape de réception d'un nombre de coefficients représentatifs du champ acoustique généré par ledit au moins un élément sous la forme d'une combinaison linéaire de fonctions spatio-temporelles, lesquels coefficients sont directement utilisés lors de ladite sous-étape de détermination de paramètres spatiaux et/ou acoustiques dudit au moins un élément ; - ladite sous-étape de calibrage comporte en outre une sous-étape de détermination de la position dans au moins l'une des trois dimensions de l'espace dudit au moins un élément dudit ensemble de restitution ; - ladite étape de calibrage comporte en outre une sous-étape de détermination de la réponse spatio-temporelle dudit au moins un élément dudit ensemble de restitution ; - ladite étape de calibrage comporte en outre une sous-étape de détermination de la réponse en fréquence dudit au moins un élément dudit ensemble de restitution ; - le procédé comporte une étape de simulation de tout ou partie des paramètres nécessaires à la réalisation de ladite étape de détermination de filtres de reconstruction ; - ladite étape de simulation comporte : - une sous-étape de détermination des paramètres manquants parmi les paramètres utilisés lors de ladite étape de détermination de filtres de reconstruction ; -une pluralité de sous-étapes de calcul permettant de déterminer la ou les valeurs du ou des paramètres manquants tels que définis précédemment en fonction des paramètres reçus, de la fréquence, et de valeurs par défaut prédéterminées ; - ladite étape de simulation comporte une sous-étape de détermination d'une liste d'éléments de l'ensemble de restitution actifs en fonction de la fréquence, et lesdites sous-étapes de calcul sont réalisées pour les seuls éléments de ladite liste ; - ladite étape de simulation comporte une sous-étape de calcul d'un paramètre représentatif de l'ordre de fonctionnement limitant le nombre de coefficients à prendre en compte lors de ladite étape de détermination de filtres de reconstruction à partir au moins de la position dans l'espace de tout ou partie des éléments de l'ensemble de restitution ; - ladite étape de simulation comporte une étape de détermination de paramètres représentatifs d'une fenêtre spatiale sous forme de coefficients de pondération à partir d'un paramètre représentatif de la fenêtre spatiale dans le repère sphérique et/ou d'un paramètre représentatif du rayon de ladite fenêtre spatiale lorsque celle-ci est une boule ; - ladite étape de simulation comporte une sous-étape de détermination d'une liste de fonctions spatio-temporelles dont la reconstruction est imposée à partir de la position de tout ou partie des éléments de l'ensemble de restitution ; - le procédé comporte une étape de saisie permettant de déterminer tout ou partie des paramètres utilisés lors de ladite étape de détermination de filtres de reconstruction ; - ladite étape de détermination de filtres de reconstruction comprend : - une pluralité de sous-étapes de calcul réalisées pour un nombre fini de fréquences de fonctionnement et permettant de délivrer une matrice de pondération du champ acoustique, une matrice représentative du rayonnement de l'ensemble de restitution, et une matrice représentative des fonctions spatiotemporelles dont la reconstruction est imposée ; et - une sous-étape de calcul d'une matrice de décodage, réalisée pour un nombre fini de fréquences de fonctionnement, à partir de la matrice de pondération du champ acoustique, de la matrice représentative du rayonnement de l'ensemble de restitution, de la matrice représentative des fonctions spatiotemporelles dont la reconstruction est imposée, et d'un paramètre représentatif de la capacité d'adaptation locale souhaitée à l'irrégularité spatiale de l'ensemble de restitution, représentative des filtres de reconstruction ; - ladite sous-étape de calcul permettant de délivrer une matrice représentative du rayonnement de l'ensemble de restitution, est réalisée à partir de paramètres représentatifs pour chaque élément : - des trois coordonnées de sa position par rapport au centre pla- cé dans la zone d'écoute ; et/ou - de sa réponse spatio-temporelle ; et - ladite sous-étape de calcul permettant de délivrer une matrice représentative du rayonnement de l'ensemble de restitution est réalisée en outre à partir de paramètres représentatifs pour chaque élément de sa réponse en fréquence. L'invention a également pour objet un programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé, lorsque ledit programme est exécuté sur un ordinateur. L'invention a également pour objet un support amovible du type comportant au moins un processeur de traitement et un élément de mémoire non volatile, caractérisé en ce que ladite mémoire comprend un programme comprenant des instructions pour l'exécution des étapes du procédé, lorsque ledit processeur exécute ledit programme. L'invention a également pour objet un dispositif de pilotage d'un ensemble de restitution d'un champ acoustique, comportant une pluralité d'éléments de restitution, caractérisé en ce qu'il comporte au moins : - des moyens de détermination de filtres de reconstruction représentatifs dudit ensemble de restitution adaptés pour permettre une prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution ; et - des moyens de détermination d'au moins un signal de pilotage desdits éléments dudit ensemble de restitution, ledit au moins un signal étant obtenu par application desdits filtres de reconstruction à un nombre fini de coeffi- cients représentatifs de la distribution dans le temps et dans les trois dimensions de l'espace dudit champ acoustique à restituer. Selon d'autres caractéristiques de l'invention : - le dispositif est associé à des moyens de mise en forme d'un signal d'entrée comportant des informations temporelles et spatiales d'un environnement sonore à restituer, adaptés pour décomposer lesdites informations sur une base de fonctions spatio-temporelles afin de délivrer un signal comportant ledit nombre fini de coefficients représentatifs de la distribution dans le temps et dans les trois dimensions de l'espace dudit champ acoustique à restituer, correspondant audit environnement sonore, sous la forme d'une combinaison linéaire desdites fonctions spatio-temporelles ; - lesdites fonctions spatio-temporelles sont des fonctions dites de Fourier-Bessel et/ou des combinaisons linéaires de ces fonctions ; - lesdits moyens de détermination de filtres de reconstruction reçoivent en entrée au moins l'un des paramètres parmi les paramètres suivants : - des paramètres représentatifs d'au moins une des trois coordonnées de la position de chaque ou certains des éléments, par rapport au centre placé dans la zone d'écoute ; - des paramètres représentatifs des réponses spatio-temporelles de chaque ou certains des éléments ; - un paramètre décrivant un ordre de fonctionnement limitant le nombre de coefficients à prendre en compte dans les moyens de détermination de filtres de reconstruction ; - des paramètres représentatifs des gabarits desdits éléments de restitution ; - un paramètre représentatif de la capacité d'adaptation locale souhaitée à l'irrégularité spatiale de la configuration dudit ensemble de restitution ; - un paramètre définissant le modèle de rayonnement desdits éléments de restitution ; - des paramètres représentatifs de la réponse en fréquence desdits éléments de restitution ; - un paramètre représentatif d'une fenêtre spatiale ; - des paramètres représentatifs d'une fenêtre spatiale sous forme de coefficients de pondération ; - de paramètres représentatifs du rayon d'une fenêtre spatiale lorsque celle-ci est une boule ; et - des paramètres représentatifs d'une liste de fonctions spatiotemporelles dont la reconstruction est imposée ; - chacun desdits paramètres reçus par lesdits moyens de détermination de filtres de reconstruction est véhiculé par l'un des signaux parmi le groupe de signaux suivants : - un signal de définition comportant des informations représentatives des caractéristiques spatiales de l'ensemble de restitution ; - un signal supplémentaire comportant des informations représentatives des caractéristiques acoustiques associées aux éléments de l'ensemble de restitution ; et - un signal d'optimisation comportant des informations relatives à une stratégie d'optimisation, afin de délivrer, à l'aide des paramètres contenus dans ces signaux, un signal représentatif desdits filtres de reconstruction représentatifs dudit ensemble de restitution ; - le dispositif est associé à des moyens de détermination de tout ou partie des paramètres reçus par lesdits moyens de détermination de filtres de reconstruction, lesdits moyens comportant au moins l'un des éléments suivants : - des moyens de simulation ; - des moyens de calibrage ; - des moyens de saisie de paramètres ; - lesdits moyens de détermination de filtres de reconstruction sont adaptés pour déterminer un ensemble de filtres représentatifs de la position dans l'espace des éléments de l'ensemble de restitution ; et - lesdits moyens de détermination de filtres de reconstruction sont adaptés pour déterminer un ensemble de filtres représentatifs de l'effet de salle induit par la zone d'écoute. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés, sur lesquels : - la Fig. 1 est une représentation d'un repère sphérique ; - la Fig. 2 est un schéma d'un système de restitution selon l'invention ; - la Fig. 3 est un schéma synoptique du procédé de l'invention ; - ta Fig. 4 est un schéma détaillant les moyens de calibrage ; - la Fig. 5 est un schéma détaillant l'étape de calibrage ; - la Fig. 6 est un schéma de l'étape de simulation ; - la Fig. 7 est un schéma des moyens de détermination de filtres de reconstruction ; - la Fig. 8 est un schéma de l'étape de détermination de filtres de reconstruction ; - la Fig. 9 est un mode de réalisation de l'étape de mise en forme du signal d'entrée ; et - la Fig. 10 est un mode de réalisation de l'étape de détermination de signaux de pilotage. Sur la figure 1, on a représenté un repère sphérique classique, de manière à préciser le système de coordonnées auquel il est fait référence dans le texte. Ce repère est un repère orthonormal, d'origine 0 et comportant trois axes (OX), (OY) et (OZ). Dans ce repère, une position notée x est décrite au moyen de ses coordonnées sphériques (r, 0, 0), où r désigne la distance par rapport à l'origine 0, SI'orientation dans le plan vertical et hui'orientation dans le plan horizontal. Dans un tel repère, un champ acoustique est connu si l'on définit en tout point à chaque instant t la pression acoustique notée p (r, #,#,t), dont la transformée de Fourier temporelle est notée P (r, où f désigne la fréquence. La figure 2 est une représentation d'un système de restitution selon l'invention. Ce système comporte un décodeur 1 pilotant un ensemble de restitution 2 qui comporte une pluralité d'éléments 31 à 3N, tels que des haut-parleurs, des enceintes acoustiques ou toute autre source sonore, agencés de manière quelconque dans un lieu d'écoute 4. On place arbitrairement, dans le lieu d'écoute 4, l'origine 0 du repère que l'on appelle centre 5 de l'ensemble de restitution L'ensemble des caractéristiques spatiales, acoustiques et électrodynamiques est considéré comme les caractéristiques intrinsèques de restitution. Le système comporte également des moyens 6 de mise en forme d'un signal d'entrée SI et des moyens 7 de génération de paramètres comportant des moyens 8 de simulation, des moyens 9 de calibrage et des moyens 10 de saisie de paramètres. Le décodeur 1 comporte des moyens 11 de détermination de signaux de pilotage et des moyens 12 de détermination de filtres de reconstruction. Le décodeur 1 reçoit en entrée un signal SIFB comportant des informations représentatives du champ acoustique tridimensionnel à restituer, un signal de définition SL comportant des informations représentatives des caractéristiques spatiales de l'ensemble de restitution 2, un signal supplémentaire RP comportant des informations représentatives des caractéristiques acoustiques associées aux éléments 31 à 3N et un signal d'optimisation OS comportant des informations relatives à une stratégie d'optimisation. Le décodeur émet à l'attention de chacun des éléments 31 à 3N de l'ensemble de restitution 2, un signal sci à SCN de pilotage spécifique. Sur la figure 3, on a représenté schématiquement les étapes principales du procédé mis en oeuvre dans un système selon l'invention tel que décrit en référence à la figure 2. Le procédé comporte une étape 20 de saisie de paramètres d'optimisation, une étape 30 de calibrage permettant de mesurer certaines caractéristiques de l'ensemble de restitution 2 et une étape 40 de simulation. Lors de l'étape 20 de saisie de paramètres mise en oeuvre par les moyens 10 d'interface, certains paramètres du fonctionnement du système peuvent être définis manuellement par un opérateur ou être délivrés par un dispositif adapté. Lors de l'étape 30 de calibrage, décrite plus en détails en référence aux figures 4 et 5, les moyens de calibrage 9 sont reliés tour à tour avec chacun des éléments 31 à 3N de l'ensemble de restitution 2 afin de mesurer des paramètres associés à ces éléments. L'étape 40 de simulation, mise en oeuvre par les moyens 8, permet de simuler les signaux de paramètres nécessaires au fonctionnement du système qui ne sont ni saisis lors de l'étape 20 ni mesurés lors de l'étape 30. Les moyens 7 de génération de paramètres délivrent alors en sortie le signal de définition SL, le signal supplémentaire RP et le signal d'optimisation OS. Ainsi, les étapes 20,30 et 40 permettent de déterminer l'ensemble des paramètres nécessaires à la mise en oeuvre de l'étape 50. Suite à ces étapes, le procédé comporte une étape 50 de détermination de filtres de reconstruction mise en oeuvre par les moyens 12 du décodeur 1 et permettant de délivrer un signal FD représentatif des filtres de reconstruction. Cette étape 50 de détermination de filtres de reconstruction permet de prendre en compte les caractéristiques au moins spatiales de l'ensemble de restitution 2 définies lors des étapes 20 de saisie, 30 de calibrage ou 40 de simulation. L'étape 50 permet également de prendre en compte les caractéristiques acoustiques associées aux éléments 31 à 3N de l'ensemble de restitution 2 et les informations relatives à une stratégie d'optimisation. Les filtres de reconstruction obtenus à l'issue de l'étape 50 sont par la suite mémorisés dans le décodeur 1 de sorte que les étapes 20,30, 40 et 50 ne sont répétées qu'en cas de modification de l'ensemble de restitution 2 ou des stratégies d'optimisation. En fonctionnement, le signal SI comportant des informations temporelles et spatiales d'un environnement sonore à restituer, est fourni aux moyens 6 de mise en forme, par exemple par une acquisition directe ou par lecture d'un enregistrement ou par synthèse à l'aide de logiciels informatiques. Ce signal SI est mis en forme lors d'une étape 60 de mise en forme. A l'issue de cette étape, les moyens 6 délivrent au décodeur 1 un signal SIFB comportant un nombre fini de coefficients représentatifs, sur une base de fonctions spatio-temporelles, de la distribution dans le temps et dans les trois dimensions de l'espace, d'un champ acoustique à restituer correspondant à l'environnement sonore à restituer. En variante, le signal SIFB est fourni par des moyens extérieurs, par exemple un micro-ordinateur comportant des moyens de synthèse. L'invention se base sur l'utilisation d'une famille de fonctions spatiotemporelles permettant de décrire les caractéristiques de tout champ acoustique. Dans le mode de réalisation décrit, ces fonctions sont des fonctions dites de Fourier-Bessel sphériques de première espèce appelées par la suite fonctions de Fourier-Bessel. Dans une zone vide de sources sonores et vide d'obstacles, les fonctions de Fourier-Bessel sont solutions de l'équation des ondes et constituent une base qui engendre tous les champs acoustiques produits par des sources sonores situées à l'extérieur de cette zone. Tout champ acoustique tridimensionnel s'exprime donc par une combinaison linéaire des fonctions de Fourier-Bessel, selon l'expression de la transformée de Fourier-Bessel inverse qui s'exprime : Dans cette équation, les termes Pl,m(f) sont, par définition, les coeffi- 2f cients de Fourier-Bessel du champ p (/, t), k= -, c est la célérité du son dans c l'air (340 ms-1), jz (kr) est la fonction de Bessel sphérique de première espèce d'ordre 1 définie par il (X) =J--M/2 (x) où Jv (x) est la fonction de Bessel de pre- 2x mière espèce d'ordre v, et y/m () est l'harmonique sphérique réelle d'ordre 1 et de terme m, avec m allant de-1 à 1, définie par : - cos 0 yl (R (cosS) pour m=O 21r - = ' (cos sin (M ) pour m Dans cette équation, les Pt (x) sont les fonctions de Legendre associées définies par : avec PI (x) les polynômes de Legendre, définis par : Les coefficients de Fourier-Bessel s'expriment aussi dans le domaine temporel par les coefficients PZ, m (t) correspondant à la transformée de Fourier temporelle inverse des coefficients Pl, (f). En variante, le procédé de l'invention utilise des bases de fonctions s'exprimant comme des combinaisons linéaires, éventuellement infinies, de fonctions de Fourier-Bessel. Lors de l'étape 60 de mise en forme, réalisée par les moyens 6, le signal d'entrée SI est décomposé en coefficients de Fourier-Bessel Pl, m (t) de manière à établir les coefficients formant le signal SIFB- La décomposition en coefficients de Fourier-Bessel est conduite jusqu'à un ordre limite L défini préalablement à cette étape 60 de mise en forme lors de l'étape 20 de saisie. A l'issue de l'étape 60, le signal SIFB délivré par les moyens 6 de mise en forme est introduit dans les moyens 11 de détermination des signaux de pilotage. Ces moyens 11 reçoivent également le signal FD représentatif des filtres de reconstruction définis en prenant en compte notamment la configuration spatiale de l'ensemble de restitution 2. Les coefficients du signal SIFB, délivrés à l'issue de l'étape 60, sont utilisés par les moyens 11 lors d'une étape 70 de détermination des signaux de pilotage SC1 à scN des éléments de l'ensemble de restitution 2 à partir de l'application des filtres de reconstruction déterminés lors de l'étape 50 à ces coefficients. Les signaux sci à SCN sont alors délivrés afin d'être appliqués aux éléments 31 à 3N de l'ensemble de restitution 2 qui restituent le champ acoustique dont les caractéristiques sont sensiblement indépendantes des caractéristiques intrinsèques de restitution de l'ensemble de restitution 2. Grâce au procédé de l'invention, les signaux SC1 à sCNde pilotage sont adaptés pour permettre une restitution optimale du champ acoustique qui exploite au mieux les caractéristiques spatiales et/ou acoustiques de l'ensemble de restitution 2, notamment l'effet de salle, et qui intègre la stratégie d'optimisation choisie. Ainsi, du fait de la quasi-indépendance entre les caractéristiques intrinsèques de restitution de l'ensemble de restitution 2 et du champ acoustique restitué, il est possible de rendre ce dernier sensiblement identique au champ acoustique correspondant à l'environnement sonore représenté par les informations temporelles et spatiales reçues en entrée. On va maintenant décrire plus en détail les étapes principales du procédé de l'invention. Lors de l'étape 20 de saisie de paramètres un opérateur ou un système de mémoire adapté peut spécifier tout ou partie des paramètres de calcul et notamment : - ; n, représentatif de la position de l'élément 3n par rapport au centre d'écoute 5 ; ; n s'exprime dans le repère sphérique au moyen des coordonnées , Le signal de définition SL véhicule les paramètres Xn, le signal supplémentaire RP, les paramètres Hn (f) et Nl,m,n(f) et le signal d'optimisation OS, les paramètres Gn (/),/. (/), { (4, )} (/), Z (/), /), (/), (/) et 7 (/). Les moyens 10 d'interface mettant en oeuvre cette étape 20 sont des moyens de type classiques tels qu'un micro-ordinateur ou tout autre moyen approprié. On va maintenant décrire plus en détail l'étape 30 de calibrage et les moyens 9 qui la mettent en oeuvre. Sur la figure 4 on a représenté le détail des moyens 9 de calibrage. Ils comportent un module 91 de décomposition, un module 92 de détermination de réponse impulsionnelle et un module 93 de détermination de paramètres de calibrage. Les moyens 9 de calibrage sont adaptés pour être connectés à un dispositif d'acquisition sonore 100 tel qu'un microphone ou tout autre dispositif adapté, et pour être connectés tour à tour à chaque élément 3n de l'ensemble de restitution 2 afin de prélever des informations sur cet élément. Sur la figure 5, on a représenté le détail d'un mode de réalisation de l'étape 30 de calibrage mise en oeuvre par les moyens 9 de calibrage et permettant de mesurer des caractéristiques de l'ensemble de restitution 2. Lors d'une sous-étape 32, les moyens 9 de calibrage émettent un signai spécifique un (t) tel qu'une séquence pseudo-aléatoire MLS (Maximum Length Sequence) à l'attention d'un élément 3n. Le dispositif d'acquisition 100 reçoit, lors d'une sous-étape 34, l'onde sonore émise par l'élément 3n en réponse à la réception du signal un (t) et transmet des signaux CI, (t) représentatifs de l'onde reçue au module de décomposition 91. Lors d'une sous-étape 36, le module de décomposition 91 décompose les signaux captés par le dispositif d'acquisition 100 en un nombre fini de coefficients de Fourier-Bessel (t). Par exemple, le dispositif 100 délivre des informations de pression p (t) et de vitesse v (t) au centre 5 de l'ensemble de restitution. Dans ce cas, les coef- ficients qo, o (t) à ql, 1 (t) représentatifs du champ acoustique se déduisent des signaux co, o (t) à cl, l (t) selon les relations suivantes : qo, o (t) =-74Li-r co, o (t) avec co, o (t) =p (t) ql l (t) =fc+cl l (t) avec q, (t) = vy (t) ? t. o (=-/x'c) avec c,, ot) = vz (t) lx-4r-e ql l (t) =-pc ; cl l (t) avec cl l (t) = vx (t) r Dans ces équations, vx (t), vy (t) et vz (t) désignent les composantes du vecteur vitesse v (t) dans le repère orthonormal considéré et p désigne la masse volumique de l'air. Lorsque ces coefficients sont définis par le module 91, ils sont adressés au module 92 de détermination de réponse. Lors d'une sous-étape 38, le module 92 de détermination de réponse détermine les réponses impulsionnelles hpl m (t) qui relient les coefficients de Fourier-Bessel ql, m (t) et le signal émis un (t). La réponse impulsionnelle délivrée par le module 92 de détermination de réponse est adressée au module 93 de détermination de paramètres. Lors d'une sous-étape 39, le module 93 déduit des informations sur des éléments de l'ensemble de restitution. Dans le mode de réalisation décrit, le module 93 de détermination de paramètres détermine la distance rn entre l'élément 3n et le centre 5 à partir de sa réponse hpo, o (t) et de la mesure du temps mis par le son pour se propager de l'élément 3n au dispositif d'acquisition 100, grâce à des méthodes d'estimation de retard sur la réponse hpo, o (t). Dans le mode de réalisation décrit, le dispositif d'acquisition 100 est en mesure d'encoder sans ambiguïté l'orientation d'une source dans l'espace. Ainsi, il apparaît pour chaque instant t des relations trigonométriques entre les 3 réponses hpl,-i (t), hpl, o (t) et hpl, l (t) impliquant les coordonnées an et ,. Le module 93 détermine les valeurs hpi,. i, hpl, o et Api, 1 correspondant aux valeurs prises par les réponses hpl,-i (t), hpi, o (t) et Api, (t) à un instant t choisi arbitrairement tel que par exemple l'instant pour lequel hpo, o (t) atteint son maximum. Par la suite, le module 93 estime des coordonnées an et à partir des valeurs hpl,-i, hpi, o et hpl, l au moyen des relations trigonométriques suivantes : - pourio > 0 : =arctan '.-' '. o ( J hp, ~, + hp 2 - pouro hp1,-1 - pour hp1,1 > 0 : #n=-arctan hp1,1 hp1,-1 - pour hp1,1 0 : #n=#-arctan hp1,1 Ces relations admettent les cas particuliers suivants : - pour hp1,0=0 et hp1,1=0 : #n=#/2, - pour hp1,1=0 et hp1,-1=0 et hp1,0=0 : #n et #n sont indéfinis - pour hp1,1=0 et hp1,-1#0 et hp1,0=0 : #n=#/2 et #n=-signe(hp1,-1)#/2 - pour hp1,1=0 et hp1,-1#0 et hp1,0#0 : #n=-signe(hp1,-1)#/2 Avantageusement, les coordonnées On et , sont estimées sur plusieurs instants. La détermination finale des coordonnées On et #n est obtenue au moyen de techniques de moyennage entre les différentes estimations. En variante, les coordonnées fin et #n sont estimées à partir d'autres réponses parmi les hpl. (t) disponibles ou sont estimées dans le domaine fréquentiel à partir des réponses HPl,m(f). Ainsi définis, les paramètres ru #n et sont transmis au décodeur 1 par le signal de définition SL. Dans le mode de réalisation décrit, le module 93 délivre également la fonction de transfert Hn (f) de chaque élément 3n, à partir des réponses hplm (t) issues du module 92 de détermination de réponse. Une solution consiste à construire la réponse hp'o, o (t) correspondant à la sélection de la partie de la réponse hpo, o (t) qui comporte un signal non nul et dénué des réflexions introduites par le lieu d'écoute 4. La réponse en fréquence Hn (j) est déduite par transformée de Fourier de la réponse hp o o (t) préalablement fenêtrée. La fenêtre peut être choisie parmi les fenêtres classiques de lissage, comme par exemple rectangulaire, Hamming, Hanning, et Blackman. Les paramètres Hn(f) ainsi définis sont transmis au décodeur 1 par le signal supplémentaire RP. Dans le mode de réalisation décrit, le module 93 délivre également la réponse spatio-temporelle NZ, m, n (j) de chaque élément 3n de l'ensemble de restitution 2, déduite en appliquant un ajustement de gain et un alignement temporel des réponses impulsionnelles hplm (t) à partir de la mesure de la distance rn de l'élément 3n de la manière suivante : La réponse spatio-temporelle contient une grande quantité d'informations caractérisant l'élément 3n, notamment sa position et sa réponse en fréquence. Elle est représentative également de la directivité de l'élément 3n, de sa non ponctualité, ainsi que de l'effet de salle résultant du rayonnement de l'élément 3n dans le lieu d'écoute 4. Le module 93 applique un fenêtrage temporel à la réponse pour ajuster la durée de prise en compte de l'effet de salle. La réponse spatiotemporelle exprimée dans le domaine fréquentiel N/, m, n (j) est obtenue par transformée de Fourier de la préponse #l,m,n(t). La réponse spatio-temporelle Nlmn (f) est ensuite fenêtrée fréquentiellement afin d'ajuster la bande de fréquences sur laquelle l'effet de salle est pris en compte. Le module 93 délivre alors les paramètres Nlm, n (/) ainsi mis en forme qui sont fournis au décodeur 1 par le signal supplémentaire RP. Les sous-étapes 32 à 39 sont répétées pour tous les éléments 31 à 3N de l'ensemble de restitution 2. En variante, les moyens 9 de calibrage sont adaptés pour recevoir d'autre types d'informations référentes à l'élément 3n. Par exemple, ces informations sont introduites sous la forme d'un nombre fini de coefficients de Fourier- Bessel représentatifs du champ acoustique produit par l'élément 3n dans le lieu d'écoute 4. De tels coefficients peuvent notamment être délivrés par des moyens de simulation acoustique mettant en oeuvre une modélisation géométrique du lieu d'écoute 4 pour déterminer la position des sources images induites par les réflexions dues à la position de l'élément 3n et à la géométrie du lieu d'écoute 4. Les moyens de simulation acoustique reçoivent en entrée le signal un (t) émis par le module 92 et délivrent, à l'aide du signal Cl, (t), les coefficients de Fourier-Bessel déterminés par superposition du champ acoustique émis par l'élément 3n et des champs acoustiques émis par les sources images lorsque l'élément 3n reçoit le signal un (t). Dans ce cas le module 91 de décomposition effectue uniquement une transmission du signal clm (t) au module 92. En variante, les moyens 9 de calibrage comportent d'autres moyens d'acquisition d'informations référentes aux éléments 31 à 3N, tels que des moyens de mesure de position au laser, des moyens de traitement du signal mettant en oeuvre des techniques de formation de voies ou tout autre moyen approprié. Les moyens 9 mettant en oeuvre l'étape 30 de calibrage sont constitués par exemple d'une carte électronique ou d'un programme d'ordinateur ou de tout autre moyen approprié. On va maintenant décrire le détail de l'étape 40 de simulation de paramètres et les moyens 8 qui la mettent en oeuvre. Cette étape est réalisée pour chaque fréquence fade fonctionnement. Les modes de réalisation décrits nécessitent de connaître pour chaque élément 3n sa position complète décrite par les paramètres rn, #n et #n et/ou sa réponse spatio-temporelle décrite par les paramètres Nl,m,n(f). Dans un premier mode de réalisation, décrit en référence à la figure 6, les paramètres qui ne sont ni saisis, par un opérateur ou par des moyens extérieurs, ni mesurés, sont simulés. L'étape 40 débute par une sous-étape 41 de détermination des paramètres manquants dans les signaux RP, SL et OS reçus. Lors d'une sous-étape 42, le paramètre H, (A représentatif de la réponse des éléments de l'ensemble de restitution 2 prend la valeur par défaut 1. Lors d'une sous-étape 43, le paramètre Gn (f) représentatif des gabarits des éléments de l'ensemble de restitution 2 est déterminé par seuillage sur le paramètre H, (f) dans le cas où celui-ci est mesuré, défini par l'utilisateur, ou fourni par des moyens extérieurs, sinon, Gn (f) prend la valeur par défaut 1. L'étape 40 comporte ensuite une sous-étape 44 de détermination des éléments actifs à la fréquence f considérée. Lors de cette sous-étape, une liste {n*} (f) d'éléments de l'ensemble de restitution actifs à la fréquence f est déterminée, ces éléments étant ceux dont le gabarit Gn(f) est non nul pour cette fréquence. La liste {*} (/) comprend Nf élé- ments et elle est transmise au décodeur 1 par le signal d'optimisation OS. Elle est utilisée pour sélectionner les paramètres correspondant aux éléments actifs à chaque fréquencefparmi l'ensemble des paramètres. Les paramètres d'indice n* correspondent au nème élément actif à la fréquence/ Lors d'une sous-étape 45, le paramètre L (/) représentatif de l'ordre de fonctionnement du module de détermination des filtres à la fréquence courante/, est déterminé de la manière suivante : - les moyens de simulation 8 calculent le plus petit angle amm formé par une paire d'éléments de l'ensemble de restitution au moyen d'une relation trigonométrique, telle que par exemple : parmi l'ensemble des couples (nl*, n2*) tels que nl* n2* ; - les moyens 9 de simulation déterminent l'ordre maximum L (j) qui est le plus grand entier respectant la relation Lors d'une sous-étape 46, le paramètre RM(f) définissant le modèle de rayonnement des éléments constituant l'ensemble de restitution, est déterminé automatiquement en prenant par défaut le modèle de rayonnement sphérique. Lors d'une sous-étape 47, le paramètre Wz (f) qui décrit la fenêtre spatiale représentative de la répartition dans l'espace de contraintes de reconstruction du champ acoustique sous forme de pondération de coefficients de FourierBessel est déterminé de la manière suivante : - si le paramètre W(r,f) représentatif de la fenêtre spatiale dans le repère sphérique est fourni ou saisi, Wl (/) est déduit de sa valeur en appliquant l'expression : - si le paramètre R (f), qui représente un rayon lorsque la fenêtre spa- tiale est une boule de rayon R (f), est fourni par des moyens extérieurs ou saisi, Wz (f) est déduit de sa valeur en appliquant l'expression : sinon, Wz (A est déduit de Lffl, en appliquant l'expression : +, (/) = 8 (/ ? () + (kR) --7, (kR) ()) avec =- ; W, (f) = 8 irIR3 "12 (kR) +j, 2 kR 2 ; zf ; - En variante, si la fenêtre spatiale n'est pas spécifiée, les moyens de simulation 8 attribuent au paramètre Wiff, une valeur par défaut, par exemple une fenêtre de Hamming de taille 2L (f) + 1, évaluée en 1. Le paramètre Wl (f) est déterminé pour les valeurs de l allant de 0 à L (f). Lors d'une sous-étape 48, le paramètre { (l", mk)} (f) est déduit à partir des paramètres L et et de la manière suivante : Dans un premier temps, les moyens 9 calculent les coefficients où (#n *, #n*) est la direction de l'élément de restitution 3n*. Dans un second temps, les moyens 9 calculent les coefficients Dans un troisième temps, les moyens 8 calculent, à l'aide d'un paramètre supplémentaire a, la liste des paramètres {(lk, mk)}(f), que l'on appelle Cet qui est initialement vide. Pour chaque valeur de l'ordre l, en commençant à 0, les moyens 8 réalisent les sous-étapes suivantes : - recherche de Gl= max (Gl,m) ; -détermination de la liste C, des coefficients (1, m) tels que Gl,m (en dB) soit compris entre Gl - # (en dB) et Gl (en dB). Si la sommé du nombre de termes dans C et du nombre de termes dans Ci est supérieure ou égale au nombre Nf d'éléments de restitution actifs à la fréquence f, la liste C est complète, sinon, on ajoute CI à C et on recommence la recherche de Gl pour l+1. Dans le cas où tous les éléments 31* à 3Nj sont dans un plan horizontal et où la liste des { (l, ? K/c)} (f) n'est ni saisie, ni fournie, les moyens 8 de simulation effectuent un traitement simplifié : La liste des coefficients {(lk, mk)}(f) prend la forme : où Ll est choisi de sorte que le nombre d'éléments dans cette liste soit inférieur au nombre Nf d'éléments 3n* actifs à la fréquence f Ll peut prendre pour valeur la partie entière de (Nf-1)/2, mais il est préférable de prendre pour, une valeur plus faible. Lors d'une sous-étape 49, le paramètre u (f), qui représente à la fréquence courante fIa capacité d'adaptation locale souhaitée, variant entre 0 et 1, est déterminé automatiquement en prenant par exemple la valeur par défaut 0,7. Ainsi, les moyens 9 de simulation permettent, lors de l'étape 40, de compléter les signaux SL, RP et OS de manière à délivrer aux moyens 12 de détermination de filtres de reconstruction l'ensemble des paramètres nécessaires à leur mise en oeuvre. En fonction des paramètres saisis ou mesurés, certaines des sous- étapes de simulation décrites ne sont pas réalisées. L'étape 40 de simulation constituée de l'ensemble des sous-étapes 41 à 49, est répétée pour toutes les fréquences considérées. En variante, chaque sous-étape est réalisée pour toutes les fréquences avant de passer à la sous- étape suivante. Dans un autre mode de réalisation, tous les paramètres intervenants sont fournis au décodeur 1 et l'étape 40 ne comporte alors que la sous étape 41 de réception et vérification des signaux SL, RP et OS et la sous étape 44 de détermination des éléments actifs à la fréquence f considérée. Les moyens 8 de simulation mettant en oeuvre l'étape 40 sont par exemple des programmes d'ordinateurs ou des cartes électroniques dédiées pour une telle application ou tout autre moyen approprié. On va maintenant décrire plus en détail l'étape 50 de détermination de filtres de reconstruction et les moyens 12 qui la mettent en oeuvre. Sur la figure 7, on a représenté les moyens 12 de détermination de filtres de reconstruction qui comportent un module 82 de détermination de matrices de transfert à partir des paramètres des signaux SL, RP et OS ainsi que des moyens 84 de détermination d'une matrice de décodage D*. Les moyens 12 comportent également un module 86 de stockage de la réponse des filtres de reconstruction et un module 88 de paramétrage de filtres de reconstruction. Sur la figure 8, on a représenté le détail de l'étape 50 de détermination de filtres de reconstruction. L'étape 50 est répétée pour chaque fréquence de fonctionnement et comporte une pluralité de sous-étapes de détermination de matrices représentatives des paramètres définis préalablement. L'étape 50 de détermination de filtres de reconstruction comprend une sous-étape 51 de détermination d'une matrice W de pondération du champ acoustique à partir des signaux l(f) et Wl(f). W est une matrice diagonale de taille (L(f)+1)2 contenant les coefficients de pondération WI (l) et dans laquelle chaque coefficient Who se trouve 21+ 1 fois à la suite sur la diagonale. La matrice Wa donc la forme suivante : De même, l'étape 50 comporte une sous-étape 52 de détermination d'une matrice M représentative du rayonnement de l'ensemble de restitution à partir des paramètres Nl,m,n*(f), RM(f), Hn*(f), #n* et L (f). M est une matrice de taille (L (n+1) 2 sur Nf, constituée d'éléments M/, m*, ies indices/, yM désignant) a iigne +/+M et n* désignant la colonne n. La matrice M a donc la forme suivante : Les éléments M,* s'obtiennent en fonction du modèle de rayonnement RM (f) : - si RM (j) définit un modèle de rayonnement en ondes planes, - si RM(f) définit un modèle de rayonnement en ondes sphériques, - si RM (j) définit un modèle utilisant les mesures effectuées des réponses spatio-temporelles, avec recours au modèle d'ondes planes pour les mesures manquantes, alors Ml,m,n* = Nl,m,n*(f) pour les indices l,m,n* fournis et la fréquence courante f. Le reste des Ml,m,n* est déterminé selon la relation : - si RM (j) définit un modèle utilisant les mesures effectuées des réponses spatio-temporelles, avec recours au modèle d'ondes sphériques pour les mesures manquantes, alors M* = NV/, m. M* (/) pour les indices l,m,n* fournis et la fréquence courante f. Le reste des M* est déterminé selon la relation : Dans ces expressions #l(rn*,f) est défini par l'expression : La matrice M ainsi définie est représentative du rayonnement de l'ensemble de restitution. En particulier, M est représentative de la configuration spatiale de l'ensemble de restitution. Lorsque le procédé utilise les coefficients Nl,m,n(f), la matrice M est représentative des réponses spatio-temporelles des éléments 31 à 3N et donc notamment de l'effet de salle induit par le lieu d'écoute 4. L'étape 50 comporte également une sous-étape 53 de détermination d'une matrice F représentative des fonctions de Fourier-Bessel dont on exige une reconstruction parfaite. Cette matrice est déterminée à partir du paramètre L (f), ainsi que des paramètres { (lk, mk) (f) de la manière suivante. A partir de la liste { ()} (/), en appelant K le nombre d'éléments (Ik, mk) de la liste { mu)} (0, la matrice F constituée est de taille K sur (lu+1) 2. Chaque ligne k de la matrice F contient un 1 sur la colonne ++mk, et des 0 ailleurs. Par exemple, pour une configuration de l'ensemble de restitution de type dit 5. 1 , dont la liste { (4, )} (/) peut prendre la forme { (0, 0), (1,-1), (1, 1)}, la matrice F s'écrit : Lorsque le paramètre 0 est nul, le décodeur 1 ne reproduit que les fonctions de Fourier-Bessel énumérées par les paramètres { (lk, rnk) l (f), les autres étant ignorées. Lorsque pif) est réglé à 1, le décodeur reproduit parfaitement les fonctions de Fourier-Bessel désignées par { mu)} 0 mais reproduit en plus partiellement de nombreuses autres fonctions de Fourier-Bessel parmi celles dispo- nibles jusqu'à l'ordre L (t) de sorte que le champ reconstruit est globalement plus proche de celui décrit en entrée. Cette reconstruction partielle permet au décodeur 1 de s'adapter à des configurations de restitution très irrégulières dans leur répartition angulaire. Les sous-étapes 51 à 53 mises en oeuvre par le module 82 peuvent être exécutées séquentiellement ou simultanément. L'étape 50 de détermination de filtres de reconstruction comprend ensuite une sous-étape 54 de prise en compte de l'ensemble des paramètres déterminés précédemment, mise en oeuvre par le module 84 afin de délivrer une matrice D* de décodage représentative des filtres de reconstruction. Cette matrice D* est délivrée à partir des matrices M, F, Wet du pa- ramètre, pif) selon l'expression suivante : D* = pA MT W+A MT FT (F MA MT FTrl F (I (L+l) 2 - pM A MT W) avec A = ( (I-&num;) IN +, uMT WM)-' 1 T où MT désigne la matrice transposée conjuguée de M. Les éléments D*,,,,, de la matrice D* sont organisés de la manière suivante : La matrice D* est donc représentative de la configuration de l'ensemble de restitution, des caractéristiques acoustiques associées aux éléments 31 à 3N et des stratégies d'optimisation. Dans le cas où le procédé utilise les coefficients Nl,m,n(f), la matrice D* est représentative notamment de l'effet de salle induit par le lieu d'écoute 4. Par la suite, lors d'une sous-étape 55, le module 86 de stockage de la réponse des filtres de reconstruction à la fréquence courante f complète pour la fréquence f la matrice D (j) représentative de la réponse en fréquence des filtres de reconstruction, en recevant en entrée la matrice D*. Les éléments de la matrice D* sont stockés dans la matrice D (f), en inversant le procédé de détermination de la liste {n*} (f) décrit précédemment en référence à la figure 6. Plus précisément, chaque élément D*nlm de le matrice D* est stocké dans l'élément Dn*, (/) de la matrice DW. Les éléments de D(f) non déterminés à l'issue de cette sous-étape sont fixés à 0. Une telle utilisation de la liste {n*} (j) permet la prise en compte de gabarits hétérogènes des éléments de restitution 31 à 3N. Les éléments D,,, m (/) de la matrice D (l) sont organisés de la manière suivante : L'ensemble des sous-étapes 51 à 55 est répété pour toutes les fré- quences f considérées et les résultats sont stockés dans le module de stockage 86. A l'issue de ce traitement, la matrice D (f) représentative des réponses en fréquence de l'ensemble des filtres de reconstruction est adressée au module 88 de paramétrage de filtres de reconstruction. Lors d'une sous-étape 58, le module 88 de paramétrage de filtres de reconstruction fournit alors le signal FD représentatif des filtres de reconstruction, en recevant en entrée la matrice D (f). Chaque élément Dn, I, m (f) de la matrice Do est un filtre de reconstruction qui est décrit dans le signal FD au moyen de paramètres qui peuvent prendre différentes formes. Par exemple, les paramètres du signal FD associés à chaque filtre Dn, (f) peuvent prendre les formes suivantes : - une réponse en fréquence, dont les paramètres sont directement les valeurs de Dn, ffl pour certaines fréquencesf ; - une réponse impulsionnelle finie, dont les paramètres dn,l,m(t) sont calculés par transformée de Fourier temporelle inverse de D,,,,,, (f). Chaque réponse impulsionnelle dn, (t) est échantillonnée puis tronquée à une longueur propre à chaque réponse ; ou - des coefficients d'un filtre récursif à réponse impulsionnelle infinie calculé à partir des Dn, (f) avec des méthodes d'adaptation classiques. Ainsi, les moyens 12 de détermination de filtres de reconstruction déli- vrent à l'issue de l'étape 50 un signal FD aux moyens 11 de détermination de signaux de pilotage. Dans ce mode de réalisation, ce signal FD est représentatif des paramètres suivants : - configuration spatiale des éléments de l'ensemble de restitution ; - caractéristiques acoustiques associées aux éléments de l'ensemble de restitution, notamment les réponses en fréquence et les réponses spatiotemporelles représentatives, entre autre, de l'effet de salle induit par le lieu d'écoute 4 ; - stratégies d'optimisation, notamment les fonctions spatio-temporelles dont on impose la reconstruction, la répartition dans l'espace de contraintes de reconstruction du champ acoustique et la capacité d'adaptation locale souhaitée à l'irrégularité spatiale de la configuration de l'ensemble de restitution 2. Les moyens 12 de détermination de filtres de reconstruction peuvent être réalisés sous la forme d'un logiciel dédié à cette fonction ou encore être intégrés dans une carte électronique ou tout autre moyen approprié. On va maintenant décrire plus en détail l'étape 60 de mise en forme du signal d'entrée. Lorsque le système est mis en oeuvre, il reçoit le signal d'entrée SI qui comporte des informations temporelles et spatiales d'un environnement sonore à restituer. Ces informations peuvent être de plusieurs natures, notamment : - un environnement sonore codé selon une distribution angulaire tel que par exemple le format communément appelé format B ; - une description d'un environnement sonore au moyen d'informations de position de sources virtuelles qui composent l'environnement sonore et des signaux émis par ces sources ; - un environnement sonore codé en multicanal, c'est-à-dire au moyen de signaux destinés à alimenter des haut-parleurs dont la répartition angulaire est fixée et connue et qui inclut notamment les techniques dites 7.1 , 5.1 , quadriphoniques, stéréophoniques, et monophoniques. - un environnement sonore donné par son champ acoustique sous la forme de coefficients de Fourier-Bessel. Ainsi que cela a été dit en référence à la figure 3, lors de l'étape 60, les moyens 6 de mise en forme reçoivent le signal d'entrée SI et le décomposent en coefficients de Fourier-Bessel représentatifs d'un champ acoustique correspondant à l'environnement sonore décrit par le signal SI. Ces coefficients de Fourier-Bessel sont délivrés au décodeur 1 par le signal SIFB- En fonction de la nature du signal d'entrée SI, l'étape 60 de mise en forme varie. En référence à la figure 9 on va maintenant décrire la décomposition en coefficients de Fourier-Bessel dans le cas où l'environnement sonore est codé dans le signal SI sous la forme de la description d'une scène sonore au moyen d'informations de position des sources virtuelles qui la composent et des signaux émis par ces sources. Une matrice E permet d'attribuer à chaque source virtuelle s un mo- dèle de rayonnement, par exemple en onde sphérique. E est une matrice de taille (L+1) 2 sur S, où S est le nombre de sources présentes dans la scène et L est l'ordre auquel est conduite la décomposition. La position d'une source s est désignée par ses coordonnées sphériques rs, On introduit également le vecteur Y qui contient les transformées de Fourier temporelles Ys(f) des signaux ys(t) émis par les sources. Y s'écrit : Les coefficients de Fourier-Bessel Pl,m(f) sont placés dans un vecteur P de taille (L+1/, où les 2l+1 termes d'ordre 1 sont placés les uns à la suite des autres par ordre 1 croissant. Le coefficient Pl, m (j) est ainsi l'élément d'indice 12+1+m du vecteur P qui s'écrit : Ainsi que cela est représenté en référence à la figure 9, l'obtention des coefficients Pli (/) de Fourier-Bessel, constituant le signal SIFB, correspond à un filtrage de chaque signal Ys(f) au moyen du filtre Elms (/) puis par sommation des résultats. Les coefficients Plm (f) s'expriment donc de la façon suivante : L'implantation des filtres Elms (f) peut être réalisée selon des méthodes classiques de filtrage, telles que par exemple : - le filtrage dans le domaine fréquentiel ; - le filtrage à l'aide d'un filtre à réponse impulsionnelle finie ; ou - le filtrage à l'aide d'un filtre à réponse impulsionnelle infinie. Il s'agit de la méthode la plus directe qui consiste à déduire de l'expression El,m,s(f) un filtre récursif, par exemple à l'aide d'une transformée bilinéaire. Dans le cas où le signal SI correspond à la représentation d'un environnement sonore selon un format multicanal, les moyens 6 de mise en forme effectuent les opérations décrites ci-après. Une matrice S permet d'attribuer à chaque canal c une source de rayonnement, par exemple en onde plane de direction de provenance (bye, c) correspondant à la direction de l'élément de restitution associé au canal c dans le format multicanal considéré. S est une matrice de taille (L+ 1) 2 sur C, où C est le nombre de canaux. Les éléments SZ, m, c de la matrice S s'écrivent : On définit également le vecteur Y qui contient les signaux Ye (t) correspondant à chaque canal. Y s'écrit : Les coefficients de Fourier-Bessel pl. m (t) regroupés comme précédemment dans le vecteur P sont obtenus par la relation : Chaque coefficient de Fourier-Bessel Pl, (t) constituant le signal SIFB est obtenu par combinaison linéaire des signaux y, (t) : Dans le cas où le signal SI correspond à la description angulaire d'un environnement sonore selon le format-B, les quatre signaux W (t), X (t), Y (t) et Z (t) de ce format se décomposent en appliquant de simples gains : Enfin, dans le cas où le signal SI correspond à une description du champ acoustique sous la forme de coefficients de Fourier-Bessel, l'étape 60 consiste en une simple transmission de signal. Ainsi, à l'issue de l'étape 60 de mise en forme, les moyens 6 délivrent, à l'attention des moyens 11 de détermination de signaux de pilotage, un signal SIFB correspondant à la décomposition du champ acoustique à restituer en un nombre fini de coefficients de Fourier-Bessel. Les moyens 6 peuvent être réalisés sous la forme de logiciels informatiques dédiés ou encore être réalisés sous la forme d'une carte informatique dédiée ou tout autre moyen approprié. On va maintenant décrire plus en détail l'étape 70 de détermination de signaux de pilotage. Les moyens 11 de détermination de signaux de pilotage reçoivent en entrée le signal SIFB correspondant aux coefficients de Fourier-Bessel représentatifs du champ acoustique à restituer et le signal FD représentatif des filtres de reconstruction provenant des moyens 12. Ainsi que cela a été dit précédemment, le signal FD intègre des paramètres caractéristiques de l'ensemble de restitution 2. A partir de ces informations, lors de l'étape 70, les moyens 11 déter- minent les signaux sci (t) à sev (t) délivrés à l'attention des éléments 31 à 3N. Ces signaux sont obtenus par l'application au signal SIFB des filtres de reconstruction, de réponse en fréquence D,,,,, (n, et transmis dans le signal FD. Les filtres de reconstruction sont appliqués de la façon suivante : avec PZ, m (j) les coefficients de Fourier-Bessel constituant le signal SIFB et Vn (j) défini par : où SC (/) est la transformée de Fourier temporelle de son (t). Selon la forme des paramètres du signal FD, chaque filtrage des Pl, m (j) par Dn,l,m(f) peut être réalisé selon des méthodes classiques de filtrage, telles que par exemple : - le signal FD fournit directement les réponses en fréquence Dn,k,m(f), et le filtrage est effectué dans le domaine fréquentiel, par exemple, à l'aide des techniques usuelles de convolution par blocs ; - le signal FD fournit les réponses impulsionnelles finies dn,l,m(t), et le filtrage est effectué dans le domaine temporel par convolution ; ou - le signal FD fournit les coefficients de filtres récursifs à réponses impulsionnelles infinies, et le filtrage est effectué dans le domaine temporel au moyen des relations de récurrence. Sur la figure 10, on a représenté le cas du filtre à réponse impulsionnelle finie. On définit Tn, le nombre d'échantillons propre à chaque réponse dn, l, m (t), ce qui conduit à l'expression de convolution suivante : L'étape 70 se termine par un ajustement des gains et l'application de retards afin d'aligner temporellement les fronts d'onde des éléments 31 à 3N de l'ensemble de restitution 2 par rapport à l'élément le plus éloigné. Les signaux sci (t) à sev (t) destinés à alimenter les éléments 31 à 3N se déduisent des signaux vi (t) à VM. (t) selon l'expression : Chaque élément 31 à 3N reçoit donc un signal de pilotage spécifique SC1 à SCN et émet un champ acoustique qui contribue à la reconstruction optimale du champ acoustique à restituer. Le pilotage simultané de l'ensemble des éléments 31 à 3N permet une reconstruction optimale du champ acoustique à restituer. Par ailleurs, le système décrit peut également fonctionner dans des modes simplifiés. Par exemple, dans un premier mode de réalisation simplifié, lors de l'étape 50, le module 12 de détermination de filtres reçoit uniquement les paramètres suivants : - , représentatif de la position de l'élément 3n de l'ensemble de restitution 2 ; - Wl, décrivant, directement sous forme de pondération des coefficients de Fourier-Bessel, une fenêtre spatiale représentative de la répartition dans l'espace de contraintes de reconstruction du champ acoustique ; et - L, imposant l'ordre limite de fonctionnement des moyens 12 de détermination de filtres de reconstruction. Dans ce mode simplifié, ces paramètres sont indépendants de la fréquence et les éléments 31 à 3N de l'ensemble de restitution sont actifs et supposés idéaux pour toutes les fréquences. Les sous-étapes de l'étape 50 ne sont donc réalisées qu'une seule fois. Lors de la sous-étape 52, la matrice M est construite à partir d'un modèle de rayonnement en ondes planes. Les éléments Ml, m, n de la matrice M se simplifient en : Dans ce mode simplifié, /= 1 et la liste { (} (/) ne contient aucun terme. Lors de la sous-étape 54, le module 84 détermine alors directement la matrice D selon l'expression simplifiée : Le stockage de la réponse des filtres de reconstruction n'est plus nécessaire, et la sous-étape 55 n'est pas réalisée. De même, les filtres décrits dans la matrice D étant de simples gains, la sous-étape 58 n'est pas non plus réalisée et le module 84 fournit directement le signal FD. Lors de l'étape 70, la détermination des signaux de commande s'effectue dans le domaine temporel et correspond à de simples combinaisons linéaires des coefficients Pl. m (t) suivies d'un alignement temporel selon l'expression : f. max (i%)- 1 c) == vj -----'-avec ) = p () D, 1=0 m=-l Le module 11 fournit alors les signaux de commande sci (t) à scat) destinés à l'ensemble de restitution. Dans un autre mode de réalisation simplifié, lors de l'étape 50, le module 12 de détermination de filtres reçoit en entrée les paramètres suivants : - xi, représentatif de la position de l'élément 3n de l'ensemble de restitution 2 ; - t (lkmk) l, constituant la liste des fonctions spatio-temporelles dont la reconstruction est imposée ; et - L, imposant l'ordre de fonctionnement des moyens 12 de détermination de filtres de reconstruction. Dans ce mode simplifié les paramètres sont indépendants de la fréquence et les éléments 31 à 3N de l'ensemble de restitution sont actifs et supposés idéaux pour toutes les fréquences. Les sous-étapes de l'étape 50 ne sont donc réalisées qu'une seule fois. Lors de la sous-étape 52, la matrice M est construite à partir d'un modèle de rayonnement en ondes planes. Les éléments Ml, m, n de la matrice M se simplifient en : La sous-étape 53 de détermination de la matrice F reste inchangée. Dans ce mode simplifié p = 0 et lors de la sous-étape 54, le module 84 détermine directement la matrice D selon l'expression simplifiée : Le stockage de la réponse des filtres de reconstruction n'est plus nécessaire, et la sous-étape 55 n'est pas réalisée. Les filtres décrits dans la matrice D étant de simples gains, la sous-étape 58 n'est pas non plus réalisée. C'est donc le module 84 qui fournit directement le signal FD. Lors de l'étape 70, la détermination des signaux de commande s'effectue dans le domaine temporel et correspond à de simples combinaisons linéaires des coefficients plm (t) suivies d'un alignement temporel selon l'expression : sc. - (t max (rn) - rn) avec - rn Vn - 1=0 mol /=0m=- ; Le module 11 fournit alors les signaux de commande sci (t) à sc (t) des- tinés à l'ensemble de restitution. Il apparaît que selon l'invention, les signaux de pilotage sci à SCN sont adaptés pour exploiter au mieux les caractéristiques spatiales de l'ensemble de restitution 2, les caractéristiques acoustiques associées aux éléments 31 à 3N et les stratégies d'optimisation de manière à reconstruire un champ acoustique de haute qualité. Il apparaît donc que le procédé mis en oeuvre permet notamment d'obtenir une restitution optimum d'un champ acoustique tridimensionnel quelle que soit la configuration spatiale de l'ensemble de restitution 2. L'invention n'est pas limitée aux modes de réalisation décrits. Notamment, le procédé de l'invention peut être mis en oeuvre par des calculateurs numériques tels qu'un ou plusieurs processeurs informatiques ou des processeurs de signaux numériques (DSP). Il peut également être mis en oeuvre à partir d'une plate-forme géné- rale telle qu'un ordinateur personnel. Il est également possible de concevoir une carte électronique destinée à être insérée dans un autre élément et adaptée pour mémoriser et exécuter le procédé de l'invention. Par exemple, une telle carte électronique s'intègre dans un ordinateur. Dans d'autres modes de réalisation, tout ou partie des paramètres nécessaires à l'exécution de l'étape de détermination de filtres de reconstruction est extrait de mémoires préenregistrées ou est délivré par un autre appareil dédié à cette fonction. REVENDICATIONS 1. Procédé de pilotage d'un ensemble de restitution (2) d'un champ acoustique pour obtenir un champ acoustique restitué de caractéristiques spécifiques sensiblement indépendantes des caractéristiques intrinsèques de restitution dudit ensemble (2), ledit ensemble de restitution (2) comportant une pluralité d'éléments de restitution (31 à 3N), caractérisé en ce qu'il comporte au moins : - une étape d'établissement d'un nombre fini de coefficients représentatifs de la distribution dans le temps et dans les trois dimensions de l'espace dudit champ acoustique à restituer ; - une étape (50) de détermination de filtres de reconstruction représentatifs dudit ensemble de restitution (2), comportant une sous-étape (54) de prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution (2) ; - une étape (70) de détermination d'au moins un signal de pilotage (sci à SCN) desdits éléments (31 à 3N) dudit ensemble de restitution (2), ledit au moins un signal étant obtenu par l'application, auxdits coefficients, desdits filtres de reconstruction ; et - une étape de délivrance dudit au moins un signal de pilotage (si à SCN), en vue d'une application auxdits éléments de restitution (31 à 3N) afin d'engendrer ledit champ acoustique restitué par ledit ensemble de restitution (2). 2. Procédé selon la revendication 1, caractérisé en ce que ladite étape d'établissement d'un nombre fini de coefficients représentatifs de la distribution dudit champ acoustique à restituer comporte : - une étape consistant à fournir un signal d'entrée (SI) comportant des informations temporelles et spatiales d'un environnement sonore ; et - une étape de mise en forme (60) dudit signal d'entrée (SI) par décomposition desdites informations sur une base de fonctions spatio-temporelles, cette étape de mise en forme (60) permettant de délivrer une représentation dudit champ acoustique à restituer correspondant audit environnement sonore sous la forme d'une combinaison linéaire desdites fonctions. 3. Procédé selon la revendication 1, caractérisé en ce que ladite étape d'établissement d'un nombre fini de coefficients représentatifs de la distribution dudit champ acoustique à restituer comporte : - une étape consistant à fournir un signal d'entrée (SIFB) comportant un nombre fini de coefficients représentatifs dudit champ acoustique à restituer sous la forme d'une combinaison linéaire de fonctions spatio-temporelles. 4. Procédé selon l'une quelconque des revendications 2 ou 3, caractérisé en ce que lesdites fonctions spatio-temporelles sont des fonctions dites de Fourier-Bessel et/ou des combinaisons linéaires de ces fonctions. 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que ladite sous-étape (54) de prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution (2) est réalisée au moins à partir de paramètres représentatifs, pour chaque élément (3n), des trois coordonnées de sa position (Xn) par rapport au centre (5) placé dans la zone d'écoute (4), et/ou de sa réponse spatio-temporelle (NI, m, n (f)). 6. Procédé selon la revendication 5, caractérisé en ce que ladite sous- étape (54) de prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution (2) est réalisée en outre à partir : - de paramètres (W1 (f)) décrivant, sous forme de coefficients de pondération, une fenêtre spatiale qui spécifie la répartition dans l'espace de contraintes de reconstruction du champ acoustique ; et - d'un paramètre (L (f)) décrivant un ordre de fonctionnement limitant le nombre de coefficients à prendre en compte lors de ladite étape (50) de détermination de filtres de reconstruction. 7. Procédé selon l'une quelconque des revendications 5 ou 6, caractérisé en ce que ladite sous-étape (54) de prise en compte de caractéristiques dudit ensemble de restitution (2) est réalisée en outre à partir : - de paramètres ( { (ik, mk)} (f)) constituant une liste de fonctions spatiotemporelles dont la reconstruction est imposée ; et - d'un paramètre (L (f)) décrivant un ordre de fonctionnement limitant le nombre de coefficients à prendre en compte lors de ladite étape (50) de détermination de filtres de reconstruction. 8. Procédé selon l'une quelconque des revendications 5 à 7, caractérisé en ce que ladite étape (54) de prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution (2) est réalisée en outre au moins à partir d'un des paramètres choisis dans le groupe constitué : - de paramètres (Xn) représentatifs d'au moins une des trois coordonnées de la position de chaque ou certains des éléments (31 à 3N), par rapport au centre (5) placé dans la zone d'écoute (4) ; - de paramètres (nul, m, n (fil) représentatifs des réponses spatiotemporelles de chaque ou certains des éléments (31 à 3N) ; - d'un paramètre (L (f)) décrivant un ordre de fonctionnement limitant le nombre de coefficients à prendre en compte lors de ladite étape (50) de détermination de filtres de reconstruction ; - de paramètres ( { ( 9. Procédé selon l'une quelconque des revendications 5 à 8, caractérisé en ce qu'il comporte une étape de calibrage (30) permettant de délivrer tout ou partie des paramètres utilisés dans ladite étape (50) de détermination de filtres de reconstruction. 10. Procédé selon la revendication 9, caractérisé en ce que ladite étape de calibrage (30) comporte, pour au moins l'un des éléments de restitution (3n) : - une sous-étape d'acquisition (34) de signaux représentatifs du rayonnement dudit au moins un élément (3n) dans le lieu d'écoute (4) ; et - une sous-étape (39) de détermination de paramètres spatiaux et/ou acoustiques dudit au moins un élément (3n). 11. Procédé selon la revendication 10, caractérisé en ce que ladite étape de calibrage (30) comporte : - une sous-étape d'émission (32) d'un signal spécifique (un (t)) vers ledit au moins un élément (3n) dudit ensemble de restitution (2), ladite sous-étape (34) d'acquisition correspondant à l'acquisition de l'onde sonore émise en réponse par ledit au moins un élément (3n) ; et - une sous-étape de transformation (36) desdits signaux acquis en un nombre fini de coefficients représentatifs de l'onde sonore émise, afin de permettre la réalisation de ladite sous-étape (39) de détermination de paramètres spatiaux et/ou acoustiques. 12. Procédé selon la revendication 10, caractérisé en ce que ladite sous-étape (34) d'acquisition correspond à une sous-étape de réception d'un nombre de coefficients représentatifs du champ acoustique généré par ledit au moins un élément (3n) sous la forme d'une combinaison linéaire de fonctions spatio-temporelles, lesquels coefficients sont directement utilisés lors de ladite sous- étape (39) de détermination de paramètres spatiaux et/ou acoustiques dudit au moins un élément (3n). 13. Procédé selon l'une quelconque des revendications 9 à 12, caractérisé en ce que ladite sous-étape de calibrage (30) comporte en outre une sous- étape de détermination de la position dans au moins l'une des trois dimensions de l'espace dudit au moins un élément (3n) dudit ensemble de restitution (2). 14. Procédé selon l'une quelconque des revendications 9 à 13, caractérisé en ce que ladite étape de calibrage (30) comporte en outre une sous-étape de détermination (38) de la réponse spatio-temporelle (Nl, m, n (fl) dudit au moins un élément (3n) dudit ensemble de restitution. 15. Procédé selon l'une quelconque des revendications 9 à 14, caractérisé en ce que ladite étape de calibrage (30) comporte en outre une sous-étape de détermination de la réponse en fréquence (Hn (f)) dudit au moins un élément (3n) dudit ensemble de restitution (2). 16. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte une étape de simulation (40) de tout ou partie des paramètres nécessaires à la réalisation de ladite étape (50) de détermination de filtres de reconstruction. 17. Procédé selon la revendication 16, caractérisé en ce que ladite étape de simulation (40) comporte : - une sous-étape (41) de détermination des paramètres manquants parmi les paramètres utilisés lors de ladite étape (50) de détermination de filtres de reconstruction ; - une pluralité de sous-étapes de calcul (42,43, 44,45, 46,47, 48,49) permettant de déterminer la ou les valeurs du ou des paramètres manquants tels que définis précédemment en fonction des paramètres reçus, de la fréquence, et de valeurs par défaut prédéterminées. 18. Procédé selon la revendication 17, caractérisé en ce que ladite étape de simulation (40) comporte une sous-étape (44) de détermination d'une liste ( {n*} (f)) d'éléments de l'ensemble de restitution actifs en fonction de la fréquence, et en ce que lesdites sous-étapes de calcul sont réalisées pour les seuls éléments de ladite liste. 19. Procédé selon l'une quelconque des revendications 17 ou 18, caractérisé en ce que ladite étape de simulation (40) comporte une sous-étape (45) de calcul d'un paramètre (L (f)) représentatif de l'ordre de fonctionnement limitant le nombre de coefficients à prendre en compte lors de ladite étape (50) de détermination de filtres de reconstruction à partir au moins de la position dans l'espace de tout ou partie des éléments (3n) de l'ensemble de restitution. 20. Procédé selon l'une quelconque des revendications 17 à 19, caractérisé en ce que ladite étape de simulation comporte une étape de détermination (47) de paramètres (W1 (f)) représentatifs d'une fenêtre spatiale sous forme de coefficients de pondération à partir d'un paramètre (W (r, f)) représentatif de la fenêtre spatiale dans le repère sphérique et/ou d'un paramètre (R (f)) représentatif du rayon de ladite fenêtre spatiale lorsque celle-ci est une boule. 21. Procédé selon l'une quelconque des revendications 17 à 20, caractérisé en ce que ladite étape de simulation (40) comporte une sous-étape de détermination (43) d'une liste ( {t) 22. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte une étape de saisie (20) permettant de déterminer tout ou partie des paramètres utilisés lors de ladite étape (50) de détermination de filtres de reconstruction. 23. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que ladite étape (50) de détermination de filtres de reconstruction comprend : - une pluralité de sous-étapes de calcul (51,52, 53) réalisées pour un nombre fini de fréquences de fonctionnement et permettant de délivrer une matrice (W) de pondération du champ acoustique, une matrice (M) représentative du rayonnement de l'ensemble de restitution (2), et une matrice (F) représentative des fonctions spatio-temporelles dont la reconstruction est imposée ; et - une sous-étape (54) de calcul d'une matrice (D*) de décodage, réalisée pour un nombre fini de fréquences de fonctionnement, à partir de la matrice (W) de pondération du champ acoustique, de la matrice (M) représentative du rayonnement de l'ensemble de restitution (2), de la matrice (F) représentative des fonctions spatio-temporelles dont la reconstruction est imposée, et d'un paramètre (u (f)) représentatif de la capacité d'adaptation locale souhaitée à l'irrégularité spatiale de l'ensemble de restitution, représentative des filtres de reconstruction. 24. Procédé selon la revendication 23, caractérisé en ce que ladite sous-étape de calcul (52) permettant de délivrer une matrice (M) représentative du rayonnement de l'ensemble de restitution (2), est réalisée à partir de paramètres représentatifs pour chaque élément (3n) : - des trois coordonnées de sa position (Xn) par rapport au centre (5) placé dans la zone d'écoute (4) ; et/ou - de sa réponse spatio-temporelle (Ni, m, n (f)). 25. Procédé selon la revendication 24, caractérisé en ce que ladite sous-étape de calcul (52) permettant de délivrer une matrice (M) représentative du rayonnement de l'ensemble de restitution (2) est réalisée en outre à partir de paramètres représentatifs pour chaque élément (3n) de sa réponse en fréquence (Hn (f)). 26. Programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 25, lorsque ledit programme est exécuté sur un ordinateur. 27. Support amovible du type comportant au moins un processeur de traitement et un élément de mémoire non volatile, caractérisé en ce que ladite mémoire comprend un programme comprenant des instructions pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 25, lorsque ledit processeur exécute ledit programme. 28. Dispositif de pilotage d'un ensemble de restitution (2) d'un champ acoustique, comportant une pluralité d'éléments de restitution (31 à 3N), caractérisé en ce qu'il comporte au moins : - des moyens (12) de détermination de filtres de reconstruction représentatifs dudit ensemble de restitution (2) adaptés pour permettre une prise en compte de caractéristiques au moins spatiales dudit ensemble de restitution (2) ; et - des moyens (11) de détermination d'au moins un signal de pilotage (SC1 à SCN) desdits éléments (31 à 3N) dudit ensemble de restitution (2), ledit au moins un signal étant obtenu par application desdits filtres de reconstruction à un nombre fini de coefficients représentatifs de la distribution dans le temps et dans les trois dimensions de l'espace dudit champ acoustique à restituer. 29. Dispositif selon la revendication 28, caractérisé en ce qu'il est associé à des moyens (6) de mise en forme d'un signal d'entrée (SI) comportant des informations temporelles et spatiales d'un environnement sonore à restituer, adaptés pour décomposer lesdites informations sur une base de fonctions spatiotemporelles afin de délivrer un signal (SIFB) comportant ledit nombre fini de coefficients représentatifs de la distribution dans le temps et dans les trois dimensions de l'espace dudit champ acoustique à restituer, correspondant audit environnement sonore, sous la forme d'une combinaison linéaire desdites fonctions spatiotemporelles. 30. Dispositif selon la revendication 29, caractérisé en ce que lesdites fonctions spatio-temporelles sont des fonctions dites de Fourier-Bessel et/ou des combinaisons linéaires de ces fonctions. 31. Dispositif selon l'une quelconque des revendications 28 à 30, caractérisé en ce que lesdits moyens (12) de détermination de filtres de reconstruction reçoivent en entrée au moins l'un des paramètres parmi les paramètres suivants : - des paramètres (Xn) représentatifs d'au moins une des trois coordonnées de la position de chaque ou certains des éléments (31 à 3N), par rapport au centre (5) placé dans la zone d'écoute (4) ; - des paramètres (Ni, m, n (f)) représentatifs des réponses spatiotemporelles de chaque ou certains des éléments (31 à 3N) ; - un paramètre (L (f)) décrivant un ordre de fonctionnement limitant le nombre de coefficients à prendre en compte dans les moyens (12) de détermination de filtres de reconstruction ; - des paramètres (Gn (f)) représentatifs des gabarits desdits éléments de restitution (31 à 3N) ; - un paramètre (u (f)) représentatif de la capacité d'adaptation locale souhaitée à l'irrégularité spatiale de la configuration dudit ensemble de restitution (2) ; - un paramètre (RM (f)) définissant le modèle de rayonnement desdits éléments de restitution (31 à 3N) ; - des paramètres (Hn (f)) représentatifs de la réponse en fréquence desdits éléments de restitution (31 à 3N) ; - un paramètre (W (r, f)) représentatif d'une fenêtre spatiale ; - des paramètres (WI (f)) représentatifs d'une fenêtre spatiale sous forme de coefficients de pondération ; - un paramètre (R (f)) représentatif du rayon d'une fenêtre spatiale lorsque celle-ci est une boule ; et - des paramètres ( { (ik, mk)} (f)) constituant une liste de fonctions spatiotemporelles dont la reconstruction est imposée. 32. Dispositif selon l'une quelconque des revendications 28 à 31, caractérisé en ce que chacun desdits paramètres reçus par lesdits moyens (12) de détermination de filtres de reconstruction est véhiculé par l'un des signaux parmi le groupe de signaux suivants : - un signal de définition (SL) comportant des informations représentatives des caractéristiques spatiales de l'ensemble de restitution (2) ; - un signal supplémentaire (RP) comportant des informations représentatives des caractéristiques acoustiques associées aux éléments (31 à 3N) de l'ensemble de restitution (2) ; et - un signal d'optimisation (OS) comportant des informations relatives à une stratégie d'optimisation, afin de délivrer, à l'aide des paramètres contenus dans ces signaux, un signal (FD) représentatif desdits filtres de reconstruction représentatifs dudit ensemble de restitution (2). 33. Dispositif selon la revendication 32, caractérisé en ce qu'il est associé à des moyens (7) de détermination de tout ou partie des paramètres reçus par lesdits moyens (12) de détermination de filtres de reconstruction, lesdits moyens (7) comportant au moins l'un des éléments suivants : - des moyens (8) de simulation ; - des moyens (9) de calibrage ; - des moyens (10) de saisie de paramètres. 34. Dispositif selon l'une quelconque des revendications 28 à 33, caractérisé en ce que lesdits moyens (12) de détermination de filtres de reconstruction sont adaptés pour déterminer un ensemble de filtres représentatifs de la position dans l'espace des éléments (31 à 3N) de l'ensemble de restitution (2). 35. Dispositif selon l'une quelconque des revendications 28 à 34, caractérisé en ce que lesdits moyens (12) de détermination de filtres de reconstruction sont adaptés pour déterminer un ensemble de filtres représentatifs de l'effet de salle induit par la zone d'écoute (4).