Brevet FR2874274A1 — publié le 20060217

Procédé pour localiser un impact sur une surface et dispositif pour la mise en oeuvre de ce procédé. La présente invention est relative aux procédés 5 pour localiser un impact sur une surface et aux dispositifs pour la mise en oeuvre de ces procédés. Plus particulièrement, l'invention concerne un procédé dans lequel on localise un impact sur une surface appartenant à un objet doté d'au moins NCAPT capteurs acoustiques (l'objet formant interface acoustique peut être fait d'une seule pièce ou de plusieurs éléments, assemblés ou au moins en contact mutuel), NCAPT étant un entier naturel au moins égal à 2, procédé dans lequel: - on fait capter par les capteurs respectivement NCAPT signaux ski (t) à partir d'ondes acoustiques générées dans l'objet par ledit impact, i étant un indice compris entre 1 et NCAPT qui désigne le capteur correspondant, - on calcule, pour chaque point de référence d'indice j, au moins un paramètre de validation représentatif d'au moins une intercorrélation de phases complexes représentatives desdits signaux captés ski(t) et de signaux de référence rji(t), chaque signal de référence rji(t) correspondant au signal qui serait reçu par le capteur i en cas d'impact en un point de référence j parmi NREF points de référence appartenant à ladite surface, NREF étant un entier naturel au moins égal à 1 et j étant un indice compris entre 1 et NREF, - et on localise l'impact en déterminant au moins un point de référence le plus proche du point d'impact, en appliquant au moins un critère de validation au paramètre de validation. Le document WO-A-03/107261 décrit un exemple d'un tel procédé qui donne déjà toute satisfaction. La présente invention a notamment pour but de perfectionner encore ce procédé, notamment pour obtenir un traitement du signal encore plus stable et plus fiable. A cet effet, selon l'invention, un procédé du genre en question est caractérisé en ce que ledit paramètre de validation est représentatif d'au moins une intercorrélation: PRODkj1112 '12p (w) çSki1 (w) goRji1 (w) *vski2 (w) *gRji2 (w) É É É ÇOSki2p (w) *çoRji2p (w) Ou: Ski (w) est la phase complexe de Ski (w) , pour i = il, i2, ..., i2p, - çpRj;. (w) est la phase complexe de Rai (w), pour i = tir i2r 12p, - * désigne l'opérateur conjugué, appliqué à çoSki(w) pour i = 12m et à çRii (w) pour i = i2m 1, m étant un entier 15 compris entre 1 et p; - Ski (w) est la transformée de Fourier de ski (t) , (w) est la transformée de Fourier de rai (t) , - p est un entier naturel non nul inférieur ou égal à NCAPT/2 r i2p sont 2p indices désignant 2p capteurs, compris chacun entre 1 et NCAPTÉ Grâce à ces dispositions, on obtient un procédé de positionnement d'impact qui est particulièrement fiable, notamment parce que l'intercorrélation susmentionnée est indépendante du type d'impact et de la réponse des capteurs. Dans divers modes de réalisation de l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes: - le procédé comporte une étape de calcul de l'intercorrélation PRODkj1l12"12(w), puis une étape de transformée de Fourier inverse de cette intercorrélation pour obtenir une fonction temporelle prodki1112"12p (t) à partir de laquelle on calcule ensuite ledit paramètre de validation; - on normalise l' intercorrélation PRODV-112É12p (w) avant de procéder à la transformée de Fourier inverse; - pour chaque point de référence j, on calcule une 5 fonction de ressemblance Vk](t) choisie parmi: É prOdkj1112..12p (t) , et une combinaison linéaire de plusieurs fonctions prodki1112..12p (t) correspondant à plusieurs sous-ensembles de 2p capteurs parmi NCAPT; - NCAPT est un nombre pair et Vkj(t) est proportionnelle à prodki1112.. 12p (t) , où p = NCAPT /2; - NCAPT = 3 et pour chaque point de référence j, on calcule une fonction de ressemblance Vkj (t) choisie parmi: Vkj (t) = a3. [prodkj12 (t) + prodkj23 (t) ], É et Vkj (t) = b3. [prodki12 (t) + prodki23 (t) + prodkJ13 (t) 1, a3 et b3 étant des constantes; - a3 = et b3 = 1/3; NCAPT = 4 et pour chaque point de référence j, on 20 calcule une fonction de ressemblance Vki(t) choisie parmi: Vki(t) = a4. prodk,1234 (t) , Vkj (t) = b4. [prodk312 (t) + prodki34 (t) ] , Vkj (t) = c4. [prodki12 (t) + prodki23 (t) + prodk,34 (t) + prodkj14 (t) 1, a4, b4 et c4 étant des constantes; - a4 = 1, b4 = et c4 = ; - le paramètre de validation est choisi parmi: MAXIMUMoi = I Vkj (t=0) I, É MAXIMUM1 = Max (IVkJ (t) I), Max (MAXIMUMoj) CONTRASTE1j = (Ei MAXIMUMIi - Max (MAXIMUMu) ) / (NREF -1) ENERGIE = EI (Et [ski (t) 12) on calcule uniquement Vki (t) pour l'instant t = 0 correspondant sensiblement à l'impact, comme étant la partie réelle de PRODki1112"12p (rw) , et on utilise MAXIMUMoi comme paramètre de validation; - on utilise au moins un critère de validation 5 choisi parmi les critères suivants: CONTRASTER > SEUIL1, avec SEUIL1 > 1, É CONTRASTER > SEUIL1 et MAXIMUMoi / MAXIMUMIiBRV"T > SEUIL2, où SEUILS > 1 et MAXIMUMI,BRU'T correspond au paramètre MAXIMUMIi des signaux traités antérieurement et n'ayant pas donné lieu à validation, É MAXIMUMoj > SEUIL3, avec SEUIL3 > 0 et MAXIMUMoj / MAXIMUMoiBRa'T > SEUIL4 avec SEUIL4 > 1, où MAXIMUMoiBR zT correspond au paramètre MAXIMUMoi des signaux traités antérieurement et n'ayant pas donné lieu à validation, É MAXIMUMoi / Moyenne (MAXIMUMoiBR'T) > SEUIL5 avec SEUIL5 > 1, ^ ENERGIE / ENERGIEBRUIT > SEUIL6, où SEUIL6 > 1 et ENERGIEBRVIT correspond au paramètre ENERGIE des signaux traités antérieurement et n'ayant pas donné lieu à validation; - les signaux de référence sont prédéterminés de manière théorique; - les signaux de référence sont utilisés avec ledit objet sans phase d'apprentissage; - les signaux de référence sont préalablement appris sur un dispositif de référence identique audit objet, puis sont utilisés avec ledit objet sans phase d'apprentissage. Par ailleurs, invention a également pour objet un dispositif spécialement adapté pour mettre en oeuvre un procédé tel que défini ci-dessus. D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante d'une de ses formes de réalisation, donnée à titre d'exemple non limitatif, en regard des dessins joints. Sur les dessins: - la figure 1 est une vue schématique en perspective montrant un exemple de dispositif comprenant une interface acoustique adaptée pour mettre en oeuvre un procédé selon une forme de réalisation de l'invention, - et la figure 2 est un schéma bloc du dispositif de la figure 1. Sur les différentes figures, les mêmes références 10 désignent des éléments identiques ou similaires. La figure 1 représente un dispositif 1 destiné à mettre en oeuvre la présente invention, qui comporte par exemple: - une unité centrale 2 de micro-ordinateur, - un écran 3 relié à l'unité centrale 2, - et un objet solide 5 sur lequel l'unité centrale 2 peut repérer la position d'un impact, comme il sera expliqué ci-après. L'objet 5 peut être de tout type (table, tablette, vitre, mur, porte, fenêtre, écran d'ordinateur, panneau d'affichage, borne interactive, jouet, tableau de bord de véhicule, arrière de dossier de siège, sol, pare-chocs de véhicule, etc.) dans lequel on peut faire propager des ondes acoustiques (notamment des ondes de Lamb) en générant des impacts sur sa surface 9, comme il sera expliqué ci-après. Au moins deux capteurs acoustiques 6 sont fixés à l'objet 5 et sont reliés par exemple à des entrées microphoniques 7 de l'unité centrale 2, par l'intermédiaire de câbles 8 ou par tout autre moyen de transmission (radio, infrarouge ou autre), de façon à capter lesdites ondes acoustiques et les transmettre à l'unité centrale 2. On appellera par la suite NCAPT le nombre de capteurs et chaque capteur sera identifié par un indice i compris entre 1 et NCAPT Les capteurs acoustiques 6 peuvent être par exemple des capteurs piézoélectriques, ou autres (par exemple, des capteurs capacitifs, des capteurs magnétostrictifs, des capteurs électromagnétiques, des vélocimètres acoustiques, des capteurs optiques [interféromètres laser, vibromètres laser, ...], etc.). Ils peuvent être adaptés pour mesurer par exemple les amplitudes des déplacements dus à la propagation des ondes acoustiques dans l'objet 5, ou encore la vitesse ou l'accélération de tels déplacements, ou bien encore il peut s'agir d'un capteur de pression mesurant les variations de pression dues à la propagation des ondes acoustiques dans l'objet 5. Pour permettre à l'unité centrale 2 de localiser un impact sur la surface 9, on détermine tout d'abord les signaux reçus par les capteurs i lorsqu'un impact est généré en un certain nombre NREF de points de référence 10 (identifiés chacun par un indice j compris entre 1 et NREF) de la surface 9. Dans l'exemple représenté sur la figure 1, la surface 9 constitue une interface acoustique en forme de clavier, et les zones constituant les points de référence 10 peuvent éventuellement être matérialisés par des marquages délimitant ces zones et les informations qui leur sont associées. A cet effet, on peut procéder tout d'abord à une étape d'apprentissage au cours de laquelle on génère des impacts au niveau des points de référence j de la surface 9. Ces impacts peuvent être générés par exemple en excitant successivement les points de référence j avec un outil quelconque (y compris une partie du corps humain, tel qu'un ongle), avantageusement un outil dont la surface de contact reste constante dans le temps. La force de l'impact est par exemple normale à la surface 9 ou orientée suivant une direction constante. A chaque impact, les réponses impulsionnelles sont 35 détectées par les capteurs 6 et mémorisées par l'unité centrale 2 pour former une banque de signaux dits de référence, notés ri, (t) (signal de référence détecté par le capteur numéroté i pour une excitation au point de référence numéroté j). Ces signaux de référence sont au nombre de NCAPT. NREF En variante, les signaux de référence sont prédéterminés de manière théorique, et le cas échéant utilisés ensuite avec ledit objet 5 sans phase d'apprentissage. Selon une autre variante, les signaux de référence peuvent être préalablement appris sur un dispositif de référence identique audit objet 5, puis sont le cas échéant utilisés avec ledit objet 5 sans phase d'apprentissage. Les signaux de référence (t) peuvent s'exprimer de manière suivante: r;i (t) = e; (t) * (t) * mi (t) (1) Avec: e3(t) fonction temporelle de la force d'excitation, h;i(t) réponse impulsionnelle (fonction de Green) pour une force appliquée au point j et détectée par le capteur i, mi (t) réponse impulsionnelle du capteur i, * symbole représentant l'opérateur de convolution temporelle. En passant dans le domaine fréquenciel, l'équation (1) 30 devient: R;i (w) = E; (w) H;i (w) Mi (w) (2) Avec: Ej(w) : transformée de Fourier de ej(t) Hji (w) : transformée de Fourier de hji (t) Mi (w) : transformée de Fourier de mat) Après l'étape d'apprentissage, le dispositif 1 est utilisé pour localiser n'importe quel impact en un point numéroté k de la surface 9. Les NCAPT signaux nouvellement détectés par les capteurs i, notés ski(t) sont alors comparés avec les signaux de référence ri, (t) précédemment mémorisés, de façon à déterminer si le point d'impact k correspond à un point de référence j connu. A cet effet, on détermine tout d'abord Ski (w) , les NCAPT transformées de Fourier des signaux ski(t), puis on détermine les produits Mkji des phases exponentielles des nouveaux signaux Ski (w) avec les phases des signaux Rai (w) (l'indice i correspond au numéro du capteur et le signe indique le complexe conjugué). Pour simplifier la notation, on utilisera ci-après le préfixe qP pour indiquer la phase exponentielle d'une variable complexe, par exemple: çoSki(w) 20 est la phase exponentielle de Ski (w) et gRji (w)) est la phase exponentielle de Rai (w) . Chaque produit Mkji peut donc s'écrire: Mkji (CO) = Ski (CO) ÇoRji (w) * (3) Si l'on décompose Ski (w) et Rji (w) d'après l'équation (2) : Ski (co) = Ek (w) Hki (co) Mi ((o) 30 et Rai (w) = EJ (w) Hui (w) Mi (w) D'où . MkJi ((D) = ÇOEk (co) eki (w) çPMi (w) çoEj (w) *çPHJi (w) *ei (w) * ( 3' ) On calcule alors un ou plusieurs produits de corrélation PRODkJ1112"12p (w) , respectivement par corrélation d'un nombre pair 2p de signaux MkJi (w) provenant de 2p capteurs il, i2, ... i2p parmi les NCAPT capteurs (p est un entier compris entre 1 et NCAPT/ 2) . PRODkjili2..i2p (w) = Mkjil (w) Mkji2 (CO) *. MkJi2p (w) * (4) On notera que l'opérateur de conjugaison * est appliqué à Mkji (CO) çoSki (CO) avec i i 2m, m étant un entier compris entre 1 et p; Si l'on développe la formule (4), compte tenu du fait que: ÇOEk (w) ÇOEk (w) * çpEJ (w) 49E; (w) * = 1 el (w) el (co) * çoM2 (w) 0M2 (w) * = 1, on obtient: PRODkj1l12..12p (w) = CpHji1 (w) *vHki1 (w) çoHji2 (w) *çpHki2 (w) çOHj2p (w) *çoHk2p (w) On peut noter que type d'excitation ni de la réponse des capteurs, ce qui en fait une variable extrêmement intéressante pour comparer le signal reçu à partir de l'impact au point k avec les signaux de référence de la banque de signaux, afin de déterminer la position du point k. On peut normaliser la ou les valeurs PRODkj1112"12p (w) , par exemple comme suit: PRODkj1112..12p N (CO) = (Npts / Ikj11i2.. 2p) É PRODkj1112.. 2p (a)) (5) et (4') PRODkj1112..12p (CO) ne dépend pas du Où . PRODkj1l12"12p N(w) est la valeur normalisée de PRODkjili2..i2o (w) , Npts = est la durée des signaux temporels ski (t) et rji (t) captés par les capteurs 6 et mémorisés par l'unité centrale 2, c'est-à-dire le nombre de points de chacun de ces signaux après échantillonnage et numérisation, Ikjili2..i2p = 10 I PRODkj1112"'12p (w) I (6) (Ik1112"12p est l'intégrale de I PRODkj1112"12p (w) I sur les fréquences). On notera que la valeur normalisée peut éventuellement être calculée directement, sans déterminer séparément PRODkj1112..'2p (CO) . On revient ensuite dans le domaine temporel en calculant la transformée de Fourier inverse de PRODkjili2..i2pN(w), soit prodkj1112..12p (t) . Pour déterminer si le point k correspond à un des points de référence jo, on utilise des fonctions de ressemblance Vkj (t) qui peuvent, selon le cas, être égales (ou plus généralement proportionnelles) respectivement: - à prodkj1112..12p (t), - ou à une combinaison linéaire (notamment une moyenne) de plusieurs fonctions prodkj1112"12p (t) correspondant à plusieurs sous-ensembles de 2p capteurs. A titre d'exemple, pour NCAPT = 2 capteurs, les fonctions de ressemblance Vki(t) seront égales (ou plus généralement proportionnelles) respectivement à prodkj12 (t) (i1=1, i2=2 et p=1). Pour NCAPT = 3 capteurs, les fonctions de ressemblance Vki(t) peuvent être choisies égales (ou plus généralement 30 proportionnelles) respectivement: - soit à 1/2. [prodkj12 (t) + prodk323 (t) 1 (respectivement i1=1, i2=2 et p=1 pour prodkj12 (t) , i1=2, i2=3 et p=l pour prodkj23 (t)) , - soit à 1/3. [prodkj12 (t) + prodkj23 (t) + prodkj13(t)] (respectivement i1=1, i2=2 et p=l pour prodkj12 (t) , i1=2, i2=3 et p=l pour prodkj23 (t) , i1=1, i2=3 et p=l pour prodkj13 (t)) . Pour NCAPT = 4 capteurs, les fonctions de ressemblance 5 Vki(t) peuvent être choisies égales (ou plus généralement proportionnelles) respectivement: - soit à prodki1234 (t) (i1=1, i2=2, i3=3, i4=4 et p=2), - soit à 1/2. [prodkj12 (t) + prodkj34 (t) ] 10 (respectivement i1=1, i2=2 et p=l pour prodkjl2 (t) , i1=3, i2=4 et p=l pour prodkj34 (t)) , - soit à 1/4. [prodkjl2 (t) + prodkj23 (t) + prodkj34 (t) + prodkjl4 (t) ] (respectivement i1=1, i2=2 et p=l pour prodkjl2 (t) , i1=2, 12=3 et p=l pour prodkj23 (t) , i1=3, i2=4 et p=l pour prodkj34 (t) , i1=1, i2=4 et p=l pour prodk,14 (t)) Pour un nombre de capteurs supérieur à 4, on procède de façon similaire aux cas décrits ci-dessus, comme expliqué précédemment, en calculant les fonctions Vk3(t) comme étant égales ( ou plus généralement proportionnelles) respectivement: - à prodk31112"IZp (t) (notamment si NCAPT = 2p, et i1=1, i2=2, i3=3, Ér 12p-NCAPT) - ou à une combinaison linéaire de plusieurs fonctions prodk31112"12p (t) correspondant à plusieurs sous-ensembles de 2p capteurs parmi NCAPT: notamment, Vkj (t) peut être égale (ou plus généralement proportionnelle) à 1/n fois la somme de n fonctions prodkil112-i2p (t) correspondant à n sous-groupes différents de 2p capteurs parmi NCAPT, chaque capteur i étant de préférence compris au moins une fois dans ces sous-groupes. Dans ces différents exemples, chaque fonction de ressemblance Vkj(t) est comprise entre -1 et 1. Les fonctions de ressemblance étant déterminées, on 35 calcule à partir de ces fonctions un ou plusieurs paramètres de validation. A titre d'exemple, on peut calculer un ou plusieurs des paramètres de validation suivants: - MAXIMUMoj = 1% (t=0) I 5 - MAXIMUMR = Max (I % (t) I) Max (MAXIMUMR) CONTRASTER = (Ej MAXIMUM1j - Max (MAXIMUM1j) ) / (Nt-et -1) ENERGIE = E1 (Et [ski (t) ] 2) Ainsi, on peut valider un impact situé au point k et affirmer qu'il est situé à un point de référence j0, par exemple s'il satisfait au moins l'un des groupes de critères suivants. 1) Groupe 1: L'impact k est considéré situé au point j0 si: CONTRASTE1jo > CONTRASTER pour j différent de jo et CONTRASTE1io > SEUIL1, avec SEUIL.' > 1, par exemple SEUIL1 = 2. 2) Groupe 2: L'impact k est considéré situé au point j0 si 25 CONTRASTE1io,CONTRASTER pour j différent de jo et CONTRASTE1jo > SEUIL1 et MAXIMUMoio / MAXIMUM1joBRUIT > SEUIL2 où SEUIL2 > 1, par exemple SEUIL2 = 2, et MAXIMUM1joBRUIT correspond au paramètre MAXIMUM 1j0 (moyenné ou non) des signaux traités antérieurement et n'ayant pas donné lieu à validation. Ce critère évite de valider du bruit ou des impacts successifs aléatoires sur l'objet 5. 3) Groupe 3: L'impact k est considéré situé au point jO si 10 MAXIMUMojo > MAXIMUMoj pour j différent de jo et MAXIMUMoio > SEUIL3, avec SEUIL3 > 0, par exemple SEUIL3 = 0.5, et MAXIMUMojo / MAXIMUMojoBRUIT > SEUIL4 avec SEUIL4 > 1, par exemple SEUIL4 = 2. MAXIMUMo.oBRUIT correspond au paramètre MAXIMUM 03o o3o (moyenné ou non) des signaux traités antérieurement et n'ayant pas 25 donné lieu à validation. En outre, pour affiner, on peut ajouter le critère suivant: MAXIMUMoio / Moyenne (MAXIMUMojoBRUIT) > SEUIL5 (7) 30 avec SEUIL5 > 1, par exemple SEUIL5 = 4 En plus de ces critères, on peut également ajouter le critère énergétique suivant: ENERGIE / ENERGIEBRUIT > SEUILy (8) où SEUILy > 1, par exemple SEUILy = 2 et ENERGIEBRUIT correspond au paramètre ENERGIE (moyenné ou non) des signaux traités antérieurement et n'ayant pas donné lieu à validation. on pourra donc valider un point point de A titre d'exemple, d'impact k et déterminer qu'il correspond à un référence jO en utilisant l'une des combinaisons suivantes de critères et groupes de critères: (groupe 1), (groupe 2), (groupe 3), (groupe 1) ou (groupe 2), - (groupe 2) ou (groupe 3), (groupe 1) ou (groupe 3), - (groupe 1) ou (groupe 2) ou (groupe 3), (groupe 1) et (7), (groupe 2) et (7), - (groupe 3) et (7), [(groupe 1) ou (groupe 2)] et (7), - [(groupe 2) ou (groupe 3)] et (7), [(groupe 1) ou (groupe 3)] et (7), [(groupe 1) ou (groupe 2) ou (groupe 3)] et (7), (groupe 1) et (8), (groupe 2) et (8), - (groupe 3) et (8), [(groupe 1) ou (groupe 2)] et (8), - [(groupe 2) ou (groupe 3)] et (8), [(groupe 1) ou (groupe 3)] et (8), - [(groupe 1) ou (groupe 2) ou (groupe 3)] et (8), - (groupe 1) et (7) et (8), (groupe 2) et (7) et (8), (groupe 3) et (7) et (8), - [(groupe 1) ou (groupe 2)] et (7) et (8), [(groupe 2) ou (groupe 3) ] et (7) et (8), [(groupe 1) ou (groupe 3)] et (7) et (8), - [(groupe 1) ou (groupe 2) ou (groupe 3)] et (7) et (8). L'unité centrale 2 peut ainsi localiser le point d'impact k sur la surface 9 de l'objet 5. La détermination de ce point d'impact peut éventuellement être la seule information recherchée par l'unité centrale, ou bien elle peut éventuellement permettre à ladite unité centrale 2 d'en déduire une autre information, par exemple une information prédéterminée affectée à un emplacement de la surface 9 (la surface 9 peut ainsi constituer un clavier acoustique). Ladite information affectée à un emplacement de la surface 9 peut être une information prédéterminée affectée à l'avance audit emplacement, ou encore une information déterminée dynamiquement à chaque nouvel impact sur la surface 9, en fonction des impacts précédemment détectés. On notera que, lorsqu'on utilise le paramètre MAXIMUMo3 pour valider les impacts, par exemple lorsqu'on utilise le (groupe 3) comme critère de validation, il est possible de calculer des valeurs approchées de MAXIMUMoi de façon très simple et rapide. En effet, lorsque k est égal à j, les termes de type ej1 (w) *çoHk1 (w) de l'équation (4') sont purement réels. Or sommer simplement la partie réelle du produit PRODki1112"12p (w) revient à calculer la transformée de Fourier inverse de l'équation (4') au temps t=0. Selon cette variante, on calcule donc la partie réelle de chaque produit PRODkiili2 i2p (w) qui donne une valeur approchée de la valeur de prodki1112..12p (t) à t = O. On en déduit ensuite la valeur de Vkj (0) comme expliqué précédemment, ce qui donne la valeur du paramètre MAXIMUMoj, puis on applique le ou les critères de validation voulus, par exemple le groupe 3 susmentionné. Avec cette variante, la charge de calcul est beaucoup moins importante et un suivi en continu est plus facilement possible, même avec un nombre important de points de références. Par ailleurs, il est possible de positionner le point d'impact k sur la surface 9 même lorsqu'il n'est pas sur un des points de référence, par interpolation, comme expliqué dans le document WO-A-03/107261 susmentionné. Par ailleurs les signaux de référence peuvent être théoriquement modélisés et appliqués à des objets réels dont les caractéristiques acoustiques sont identiques à celles de l'objet considéré et utilisé dans la modélisation. Les signaux de référence appris ou théoriquement modélisés peuvent être appliqués sans phase d'apprentissage aucune à des objets ayant des caractéristiques acoustiques identiques à celles de l'objet utilisé pour l'apprentissage des signaux de référence ou considéré et utilisé dans la modélisation. REVENDICATIONS 1. Procédé dans lequel on localise un impact sur une surface (9) appartenant à un objet (5) doté d'au moins NCAPT capteurs acoustiques (6), NCAPT étant un entier naturel 5 au moins égal à 2, procédé dans lequel: on fait capter par les capteurs (6) respectivement NCAPT signaux ski (t) à partir d'ondes acoustiques générées dans l'objet par ledit impact, i étant un indice compris entre 1 et NCAPT qui désigne le capteur (6) correspondant, - on calcule, pour chaque point de référence (10) d'indice j, au moins un paramètre de validation représentatif d'au moins une intercorrélation de phases complexes représentatives desdits signaux captés ski(t) et de signaux de référence r3i(t), chaque signal de référence r3i(t) correspondant au signal qui serait reçu par le capteur i (6) en cas d'impact en un point de référence j (10) parmi NREF points de référence (10) appartenant à ladite surface (9), NREF étant un entier naturel au moins égal à 1 et j étant un indice compris entre 1 et NREF, et on localise l'impact en déterminant au moins un point de référence (10) le plus proche du point d'impact, en appliquant au moins un critère de validation au paramètre de validation, caractérisé en ce que ledit paramètre de validation est représentatif d'au moins une intercorrélation: PRODk,1112. 12p (w) çoSkil (co) ÇoRJi1 (w) *eki2 (CO) *ÇORJ12 (w) . . . ÇOSki2p (CO) *CPRJi2p (co) Oë . - Ski (r)) est la phase complexe de Ski (a)), pour i 11, 12, ..., 12p, - rpRii (ro) est la phase complexe de Rji (cc)), pour i ilr 12r..., i2p, désigne l'opérateur conjugué, appliqué à vSki(w) 2874274 18 pour i = 12m et à çoRJi (w) pour i = i2m_1i m étant un entier compris entre 1 et p; - Ski (w) est la transformée de Fourier de ski (t) , - RJi (w) est la transformée de Fourier de rJi (t) , - p est un entier naturel non nul inférieur ou égal à NCAPT/2r i2r... i2p sont 2p indices désignant 2p capteurs (6), compris chacun entre 1 et NCAPTÉ 2. Procédé selon la revendication 1, comportant une étape de calcul de l'intercorrélation PRODk1112I2p (w) , puis une étape de transformée de Fourier inverse de cette intercorrélation pour obtenir une fonction temporelle prodkJ1112"12p (t) à partir de laquelle on calcule ensuite ledit paramètre de validation. 3. Procédé selon la revendication 2, dans lequel on normalise l'intercorrélation PRODk,1112"12p(w) avant de procéder à la transformée de Fourier inverse. 4. Procédé selon la revendication 2 ou la revendication 3, dans lequel pour chaque point de référence j (10), on calcule une fonction de ressemblance Vk](t) choisie parmi: - prodkJ2112..12p (t) , et une combinaison linéaire de plusieurs fonctions prodkJ1122"12p (t) correspondant à plusieurs sous-25 ensembles de 2p capteurs (6) parmi NCAPT 5. Procédé selon la revendication 4, dans lequel NCAPT est un nombre pair et Vk] (t) est proportionnelle à prodkJ1122"12p (t) , où p = NCAPT /2. 6. Procédé selon la revendication 5, dans lequel NCAPT = 2 et Vkj (t) est proportionnelle à prodkJ12 (t) . 7. Procédé selon la revendication 4, dans lequel NCAPT = 3 et pour chaque point de référence j (10), on calcule une fonction de ressemblance Vk3(t) choisie parmi: - Vkj (t) = a3. [prodkJ12 (t) + prodkJ23 (t) ] , - et Vki(t) = b3. [prodk3l2 (t) + prodki23 (t) + prodki13 (t) ] , a3 et b3 étant des constantes. 8. Procédé selon la revendication 7, dans lequel a3 = et b3 = 1/3. 9. Procédé selon la revendication 4, dans lequel NCAPT = 4 et pour chaque point de référence j (10), on calcule une fonction de ressemblance VkJ(t) choisie parmi: - Vkj (t) = a4. prodki1234 (t) , - Vkj (t) = b4. [prodkil2 (t) + prodki34 (t) ], - Vkj (t) = c4. [prodkj12 (t) + prodk323 (t) + prodkj34 (t) + prodkj14 (t) ] , a4, b4 et c4 étant des constantes. 12. Procédé selon la revendication 9, dans lequel 15 a4 = 1, b4 = et c4 = 14. 13. Procédé selon l'une quelconque des revendications 4 à 10, dans lequel le paramètre de validation est choisi parmi: MAXIMUMoi = I Vk (t=0) I, MAXIMUMIi = Max (I VkJ (th), Max (MAXIMUM0j) CONTRASTEZ j = (E MAXIMUM13 Max (MAXIMUMM,) ) / (NREF -1) ENERGIE = EZ (Et [ski (t) ]2) . 12. Procédé selon la revendication 11, dans lequel on calcule uniquement Vk](t) pour l'instant t = 0 correspondant sensiblement à l'impact, comme étant la partie réelle de PRODki1112"12p (Co) , et on utilise MAXIMUMoj comme paramètre de validation. 13. Procédé selon la revendication 11 ou la revendication 12, dans lequel on utilise au moins un critère de validation choisi parmi les critères suivants: - CONTRASTER] > SEUIL1, avec SEUIL1 > 1, - CONTRASTER > SEUIL1 et MAXIMUMoj / MAXIMUM1sav'T > SEUIL2, où SEUIL2 > 1 et MAXIMUMIiBRVIT correspond au paramètre MAXIMUM1i des signaux traités antérieurement et n'ayant pas donné lieu à validation, - MAXIMUMoi > SEUIL3, avec SEUIL3 > 0 et MAXIMUMoj / MAXIMUMoiBRVIT > SEUIL4 avec SEUIL4 > 1, où MAXIMUMoiBRVIT correspond au paramètre MAXIMUMoi des signaux traités antérieurement et n'ayant pas donné lieu à validation, - MAXIMUMoi / Moyenne (MAXIMUMoiBRVIT) > SEUIL5, avec 10 SEUIL5 > 1, - ENERGIE / ENERGIEBRVIT > SEUIL6, où SEUIL6 > 1 et ENERGIEBRVIT correspond au paramètre ENERGIE des signaux traités antérieurement et n'ayant pas donné lieu à validation. 14. Procédé selon l'une quelconque des revendications précédentes, dans lequel les signaux de référence sont prédéterminés de manière théorique. 15. Procédé selon la revendication 14, dans lequel les signaux de référence sont utilisés avec ledit objet (5) sans phase d'apprentissage. 16. Procédé selon l'une quelconque des revendications 1 à 13, dans lequel les signaux de référence sont préalablement appris sur un dispositif de référence identique audit objet (5), puis sont utilisés avec ledit objet (5) sans phase d'apprentissage. 17. Dispositif spécialement adapté pour mettre en oeuvre un procédé selon l'une quelconque des revendications précédentes, ce dispositif comprenant: - un objet (5) doté d'au moins NCAPT capteurs acoustiques (6), NCAPT étant un entier naturel au moins égal à 2, pour capter respectivement NCAPT signaux ski (t) à partir d'ondes acoustiques générées dans l'objet (5) par un impact sur une surface (9) appartenant audit objet, i étant un indice compris entre 1 et NCAPT qui désigne le capteur (6) correspondant, - des moyens pour calculer, pour chaque point de référence (10) d'indice j, au moins un paramètre de validation représentatif d'au moins une intercorrélation de phases complexes représentatives desdits signaux captés ski (t) et de signaux de référence rji (t) , chaque signal de référence rii(t) correspondant au signal qui serait reçu par le capteur i (6) en cas d'impact en un point de référence j (10) parmi NREF points de référence (10) appartenant à ladite surface (9), NREF étant un entier naturel au moins égal à 1 et j étant un indice compris entre 1 et NREF, et des moyens pour localiser l'impact en déterminant au moins un point de référence (10) le plus proche du point d'impact, en appliquant au moins un critère de validation au paramètre de validation, caractérisé en ce que ledit paramètre de validation est représentatif d'au moins une intercorrélation: PRODkoili2..i2p (w) =çoSkil (w) eji1 (w) *eki2 (w) *çRji2 (w) . . . çoSki2p (w) *eji2p (CO) où . - tpSki (w) est la phase complexe de Ski (w) , pour i 11, i2, ..., i2p, - çRji (w) est la phase complexe de Rji (w) , pour i = i2, ..., i2p, - désigne l'opérateur conjugué, appliqué à q'Ski(w) 25 pour i = i2m et à el,: (co) pour i = i2m 1i m étant un entier compris entre 1 et p; Ski (w) est la transformée de Fourier de ski (t) , - Rji (co) est la transformée de Fourier de rji (t), - p est un entier naturel non nul inférieur ou égal à NCAPT/2, i2, i2p sont 2p indices désignant 2p capteurs (6), compris chacun entre 1 et NCAPTÉ