WO2009103624A1 - Procede et dispositif pour ameliorer la bande passante d'un systeme physique - Google Patents

Procede et dispositif pour ameliorer la bande passante d'un systeme physique Download PDF

Info

Publication number
WO2009103624A1
WO2009103624A1 PCT/EP2009/051424 EP2009051424W WO2009103624A1 WO 2009103624 A1 WO2009103624 A1 WO 2009103624A1 EP 2009051424 W EP2009051424 W EP 2009051424W WO 2009103624 A1 WO2009103624 A1 WO 2009103624A1
Authority
WO
WIPO (PCT)
Prior art keywords
value
filter
impulse response
sample
samples
Prior art date
Application number
PCT/EP2009/051424
Other languages
English (en)
Inventor
Michel Ayraud
Nathalie Pascal
Original Assignee
E2V Semiconductors
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by E2V Semiconductors filed Critical E2V Semiconductors
Priority to EP09712207A priority Critical patent/EP2243218A1/fr
Priority to US12/918,337 priority patent/US20100332577A1/en
Publication of WO2009103624A1 publication Critical patent/WO2009103624A1/fr

Links

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H03H17/0294—Variable filters; Programmable filters
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H03H17/06—Non-recursive filters

Definitions

  • the invention relates to a method and an electronic circuit for improving the bandwidth of a device which, for physical reasons, has a bandwidth less than that which one would like.
  • any electronic system made from physical elements has a limited bandwidth due to the non-ideal characteristics of the system. This is the case in particular for all electronic acquisition systems of physical quantities (sensors) but also for all systems for processing electronic signals or for transmitting electronic or optical signals. It follows from this limitation of non-ideal systems that any physical system behaves, with respect to an input quantity, as a filter which is in general a low-pass filter, but which could be more sophisticated than a simple low-pass filter. This is reflected in practice by the fact that the output of the system fails to follow the too fast variations of an input quantity of the system. When we speak here of variation, it can be mainly a temporal variation but also a spatial variation.
  • the input variation may be a time variation of light and the output variation is a time variation of voltage that represents this light variation; we then realize that, just because the sensor is a sensor made from physical elements (photodiodes, circuits for collecting electrical charges, amplifiers, transmission circuits, etc.), the electronic signal can not follow instantly. a very sudden variation of light at the entrance; the sensor response includes a low-pass time filter function.
  • the sensor can not generate an electron signal spatially varying as abruptly as the pattern.
  • the senor produces a certain filtering function, but this time in the spatial domain, therefore, as a function of a spatial variable.
  • the bandwidth in the spatial domain can be as important as in the time domain and the correction principles according to the present invention apply in both cases.
  • the invention proposes a much simpler solution for increasing the bandwidth of a physical system.
  • a finite impulse response filter is used which is calculated in the following manner, from the behavior (observed or known) of the physical system: the impulse response of the physical system is determined according to a temporal or spatial variable; the impulse response is the response to a normalized input pulse of duration tending to zero and vertical rising and falling edges; sample-by-sample is calculated an impulse response of similar shape but compressed according to the scale of the variable in a ratio corresponding to a desired bandwidth increase and expanded in amplitude in the same ratio, and the coefficients of a filter at finite impulse response able to provide at its output a succession of samples of the compressed response when applying a succession of corresponding samples of the response of the physical system to its input.
  • This finite impulse response filter is incorporated in the physical system, output (preferably) or input or inside the system, to improve the bandwidth in a ratio n.
  • an improved bandwidth signal is collected at the output of the corrected system in a factor n.
  • the choice of the interval [0, T] is such that it encompasses the most significant part of the impulse response a (t), for example the part of the response for which the amplitude of the response is greater than at least 5% of the maximum amplitude of the response.
  • Samples in the interval [T, nT] may very well be arbitrarily set to zero, even if the values of a (t) are not zero in this interval because these values are small. This may result in a poorer result but a result that is entirely acceptable.
  • the coefficients of the filter are preferably coefficients Ci where i is the index varying from 0 to N, the value of the coefficient Ci being defined by the following iteration:
  • Ci (bi - ai) C 0 - ai - i Ci .. .. - a ⁇ j.Cj .. .. - a - ⁇ Ci - ⁇ ) / a 0 ,
  • the value of the curve a (t) at time 0 is zero or close to zero, it is replaced by an arbitrary value at 0 which is not zero, preferably less than or equal to a ⁇ , for example a- ⁇ / 2.
  • the invention relates to an electronic system comprising a set of physical elements whose nature induces a bandwidth limitation of the system, and an electronic filter.
  • compensation system for giving the system a bandwidth n times greater than the system without this filter, n being a number greater than 1, in which the system without the filter has a pulse response a (t) in function a variable t, characterized in that the compensation filter is a finite impulse response filter with N + 1 coefficients Ci where i is an index varying from 0 to N, the value of the coefficient Ci being defined by the following iteration :
  • Ci (bi - ai.C 0 - a ⁇ i Ci .. .. - a ⁇ j.Cj .. .. - ai Ci-i) / a 0 , where ai is a sample value of index i varying from 1 at
  • a 0 is zero, we choose for a 0 other non-zero value, preferably between 0 and ai.
  • FIG. 1 symbolically represents a physical system receiving an input quantity e and providing an output quantity F (e);
  • FIG. 2 represents the impulse response F [x (t)] of the physical system
  • FIG. 3 represents the principle of frequency response compensation by an inverse transfer function in the frequency domain
  • FIG. 4 represents the principle of establishing a desired impulse response from the actual impulse response of the system
  • FIG. 5 represents the principle according to the invention of correction of the bandwidth by a finite impulse response filter associated with the physical system to be corrected
  • FIG. 6 represents the general structure of a finite impulse response filter
  • FIG. 7 represents the detail of the elaboration of the desired impulse response b (t) from the natural response a (t) of the system
  • FIG. 8 represents the frequency behavior of an uncorrected physical system in the presence of an input variable e (t) in slots of increasing frequency
  • FIG. 9 represents the behavior of the corrected physical system, under the same conditions as in FIG. 8.
  • FIG. 1 represents in a symbolic form (a simple rectangle) any physical system SP which, when an input input quantity e is applied to it, provides an output quantity F (e) at the output.
  • the input quantity e is an electronic signal and the output quantity F (e) is another electronic signal.
  • the input quantity e is a luminance and the output quantity F (e) is an electronic signal, a luminance sensor being present in the system.
  • the system In many applications, it would be desirable for the system to behave in an ideal manner, that is, the magnitude F (e) varies according to a determined function of e both when e varies slowly and when e varies rapidly.
  • the light sensor produce an electronic signal of amplitude proportional to the luminance with the same coefficient of proportionality when the luminance varies slowly and when it varies rapidly.
  • the system is a physical system and has limitations specific to any physical system, in particular limitations of ability to respond instantaneously to an instantaneous variation of the input quantity e. It introduces function deformations F (e) when e varies too rapidly.
  • rapid variation we must understand that we can speak of speed in time or in space. Indeed, we must consider the problems of limitation of physical systems both in terms of variation over time and variations in space as will be detailed below.
  • the variable considered for the impulse response is time t but the reasoning would be the same if the variable was spatial; this is all the more true that the spatial variable can be reduced to a temporal variable, for example when the spatially varying signals are read sequentially.
  • temporal variable t is generic and can be directly or indirectly transposed to a spatial variable.
  • a compensation filter having the transfer function H (f) in the frequency space, thus a filtering curve like that of FIG.
  • FIG. 4 shows at 4a the input pulse x (t) (as in FIG. 2a), at 4b the impulse response a (t) of the initial physical system (as in FIG. 2b), and at 4c an answer impulse b (t) desired, similar in shape to that of Figure 4b but compressed according to the time scale and expanded according to the scale of the amplitudes.
  • the width at half height of the response b (t) is that of the curve a (t) but divided by a factor n.
  • the maximum height of the response b (t) is that of the curve a (t) but multiplied by n.
  • the impulse response of the physical system is determined, an improved impulse response is deduced which is a response similar to the initial response, amplified in amplitude and compressed according to the time scale, and the characteristics of the impulse response are calculated.
  • a compensation filter which, when it receives at its input the initial impulse response produces at its output the improved impulse response.
  • the filter is a finite impulse response filter whose coefficients are easy to calculate once the desired output response for a given input signal form is known.
  • FIG. 5 represents the corrected SPC physical system whose impulse response is the impulse response of FIG. 4c.
  • This system comprises the initial physical system SP and a finite impulse response filter FIR whose function is to convert the waveform a (t) of FIG. 4b into the waveform b (t) of FIG. 4c, and at the same time to improve the bandwidth of any waveform F (e) out of the physical system SP.
  • a finite impulse response filter it is a circuit that processes an analog signal sampled at a sampling frequency F, and cyclically adds a series of weighted successive samples. by individual coefficients chosen according to the filtering function to be obtained.
  • the number of samples in the series can be of ten or several tens, or even several hundreds depending on the complexity of the filtering function to be performed.
  • the sampling of the signal to be processed by the filter can be done in the filter or upstream of the filter depending on whether the physical system delivers unsampled or sampled signals.
  • Figure 6 shows a finite impulse response filter.
  • the signal to be filtered is a signal u (t), applied in the form of successive samples at the frequency F and the output signal of the filter is a signal u '(t) output in the form of samples.
  • the filter is represented as a set of delay circuits represented by the boxes z "1 and a set of multiplier circuits making it possible to apply to multiply the signal by a weighting coefficient d specific to this circuit
  • the letter z is the variable conventionally used for the representation mathematical sampling of sampled systems and the function z "1 represents a unit delay applied to a sample, therefore a delay 1 / F if the sampling frequency is F.
  • the filter further comprises an adder ADD.
  • the output signal u '(t) of the filter is a sampled signal which is the weighted sum of the series of N + 1 last successive samples received.
  • u '(t) U 0 -C 0 + ui. Ci + U 2 -C 2 + .... U N -C N if there are N + 1 samples by counting the sample U 0 .
  • the successive samples U 0 , u- ⁇ , U 2 , etc. until u N respectively represent the current sample u (t 0 ) at a time t 0 and the previous samples u (t o -dt), u (t o -2dt), etc.
  • the coefficients of the FIR filter are calculated according to the invention from the knowledge of the impulse response a (t) of the physical system SP. This impulse response is measurable or computable from the known constitution of the physical system.
  • the calculation of the filter coefficients is done by determining which coefficients are to be used in the filter for a sampled input signal having the shape of the impulse response a (t) of Figure 4b so that the filter with the sampled form of the answer b (t) of Figure 4c.
  • n integer and greater than 1, the bandwidth of the physical system SP.
  • n can be 2 or 3.
  • the duration T may be the duration during which the impulse response is at least equal to 5% of its maximum value.
  • N is an arbitrarily chosen integer, which is all the higher since we want to represent more finely the impulse response b (t).
  • the improved response b (t) is in dashed lines: it is a curve established by smoothing from the N + 1 samples bj calculated according to the formula above. From the samples ai and bj the coefficients of the finite impulse response filter are calculated which makes it possible to obtain at its output the response b (t) when a signal whose time variation is that of the response a is applied to its input. (t).
  • b n a n .C 0 + a n -i .Ci + a o .C n
  • Ci (bi - a -
  • C 0 ) / a 0 (2.a 2 - a -
  • Ci (na n i - ai, C 0 - ai.-i .Ci - ai .Ci-i) / a 0
  • CN (na n N - aN-Co - aN-i .Ci - ai.C ⁇ ⁇ ⁇ -i) / ao
  • this filter receives the output F [e (t)] of the physical system whose bandwidth is to be increased, and it provides an output F '[e (t)] with increased bandwidth.
  • the FIR compensation filter is placed in the physical system rather than at the output thereof, that is to say that the filter can be placed upstream of the physical elements that tend to reduce bandwidth.
  • n is a non-integer value, for example 2.5 or 3.5 (or even any number that this has in practice little interest).
  • Samples values a n i are no longer available when neither is nor an integer. In this case, it will be possible to use as sample value a n i an interpolated value between the two real samples closest to the non-existent sample of rank
  • FIG. 8 represents a simulation performed for a physical system SP whose transfer function is assimilated to a first-order low-pass filter.
  • An input rectangular slot waveform (solid line) at a frequency that is greater than a filter cut-off frequency and which is variable in time (the duration of the slots decreases) is applied as a magnitude.
  • e (t) at the input of the physical system SP.
  • the output of the system is a magnitude F [e (t)] which is represented by a dashed curve.
  • FIG. 9 represents the output F '[e (t)] of the system corrected by the FIR filter calculated in the manner explained above, the same slots e (t) being applied to the input. We see that this exit follows much better slots and only deteriorates for the shortest slots.
  • the coefficients calculated are as follows:
  • This filter is only given as an illustrative example showing the considerable improvement that can be achieved in the bandwidth.
  • the physical system will not be as simple as a first-order filter, the number of coefficients N of the filter will be rather 32 or 64 and the n-factor of bandwidth increase will be rather between 2 and 4.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Filters That Use Time-Delay Elements (AREA)
  • Complex Calculations (AREA)

Abstract

L'invention concerne l'amélioration de la bande passante des systèmes physiques. On utilise un filtre à réponse impulsionnelle finie qui est calculé de la manière suivante, à partir du comportement (observé ou connu) du système physique (SP) : on détermine la réponse impulsionnelle a(t) du système physique selon une variable temporelle ou spatiale; on calcule échantillon par échantillon une réponse impulsionnelle b(t) de forme semblable mais comprimée selon l'échelle de la variable t dans un rapport n et dilatée en amplitude dans le même rapport, et on calcule les coefficients d'un filtre à réponse impulsionnelle finie apte à fournir à sa sortie le signal b(t) lorsqu'on applique le signal a(t) à son entrée. Ce filtre à réponse impulsionnelle finie est incorporé au système physique, de préférence à la sortie, pour en améliorer la bande passante dans le rapport n.

Description

PROCEDE ET DISPOSITIF POUR AMELIORER LA BANDE PASSANTE D'UN SYSTEME PHYSIQUE
L'invention concerne un procédé et un circuit électronique pour améliorer la bande passante d'un dispositif qui, pour des raisons physiques, possède une bande passante inférieure à celle que l'on souhaiterait.
On sait que tout système électronique réalisé à partir d'éléments physiques possède une bande passante limitée due aux caractéristiques non idéales du système. C'est le cas en particulier de tous les systèmes d'acquisition électronique de grandeurs physiques (capteurs) mais aussi de tous les systèmes de traitement de signaux électroniques ou de transmission de signaux électroniques ou optiques. II résulte de cette limite des systèmes non idéaux que tout système physique se comporte, vis-à-vis d'une grandeur d'entrée, comme un filtre qui est en général un filtre passe-bas, mais qui pourrait être plus sophistiqué qu'un simple filtre passe-bas. Cela se traduit en pratique par le fait que la sortie du système ne parvient pas à suivre les variations trop rapides d'une grandeur d'entrée du système. Quand on parle ici de variation, ce peut être surtout une variation temporelle mais aussi une variation spatiale.
Par exemple, pour un capteur d'image, la variation d'entrée peut être une variation temporelle de lumière et la variation de sortie est une variation temporelle de tension électrique qui représente cette variation de lumière ; on s'aperçoit alors que, du seul fait que le capteur est un capteur réalisé à partir d'éléments physiques (photodiodes, circuits de recueil de charges électriques, amplificateurs, circuits de transmission, etc), le signal électronique ne peut pas suivre instantanément une variation très brusque de lumière à l'entrée ; la réponse du capteur inclut une fonction de filtrage temporel passe-bas. De même, si on considère un capteur matriciel comportant de nombreux pixels très rapprochés, éclairé par un motif d'image spatiale très fin, on s'aperçoit que le capteur ne peut pas engendrer un signal électronique variant spatialement aussi brusquement que le motif d'image qui l'éclairé ; là encore le capteur produit une certaine fonction de filtrage, mais cette fois dans le domaine spatial donc en fonction d'une variable spatiale. La bande passante dans le domaine spatial peut être aussi importante que dans le domaine temporel et les principes de correction selon la présente invention s'appliquent dans les deux cas.
Pour analyser le comportement du capteur aussi bien en réponse temporelle qu'en réponse spatiale, on peut essayer de déterminer sa fonction de transfert temporelle ou spatiale. Plusieurs sortes de fonction de transfert, utilisant des variables de temps ou de fréquence, ou encore la variable de Laplace p des systèmes analogiques ou la variable z des systèmes échantillonnés, sont classiquement utilisées pour représenter des fonctions de transfert de systèmes physiques.
Lorsqu'on sait qu'un système physique présente une fonction de transfert limitée en fréquence, on sait qu'on peut améliorer la réponse globale du système en incluant dans celui-ci un filtre inverse qui tend à compenser la fonction de filtrage naturelle du système. Ainsi, si le système agit comme filtre temporel passe-bas, ce qui est presque toujours le cas, on peut lui incorporer un filtre passe-haut qui tend à compenser l'effet de filtrage passe-bas qu'on a préalablement identifié.
Si on connaît la fonction de filtrage passe-bas du système physique, on essaye dans ce cas de réaliser un filtre analogique ou digital qui possède une fonction de filtrage à peu-près inverse. Les fonctions de filtrage sophistiquées qui en résultent en général sont très difficiles à réaliser de manière analogique ; elles ne peuvent guère être exécutées autrement que par calcul et en utilisant de puissants processeurs. On peut par exemple déterminer la fonction de transfert du système physique dans l'espace fréquentiel par une transformation de Fourier ; cette transformation établit une courbe de réponse dans le domaine fréquentiel. Puis on calcule la fonction de filtrage inverse ; dans le domaine fréquentiel la fonction inverse est tout simplement l'inverse de la fonction de transfert initiale du système physique. Lors de l'utilisation du système, on applique une transformation de Fourier au signal de sortie du système physique pour définir des composantes dans le domaine fréquentiel, on lui applique la fonction de filtrage inverse dans le domaine fréquentiel, et on repasse dans le domaine temporel par une transformation inverse de Fourier pour obtenir un signal électronique temporel filtré, ayant une bande passante meilleure que celle du système physique non corrigé. Ce type de procédé requiert des puissances de calcul très élevées pour effectuer la transformée de Fourier et la transformée inverse à chaque instant.
L'invention propose une solution beaucoup plus simple pour augmenter la bande passante d'un système physique. Dans cette solution, on utilise un filtre à réponse impulsionnelle finie qui est calculé de la manière suivante, à partir du comportement (observé ou connu) du système physique : on détermine la réponse impulsionnelle du système physique selon une variable temporelle ou spatiale ; la réponse impulsionnelle est la réponse à une impulsion d'entrée normalisée de durée tendant vers zéro et à fronts de montée et descente verticaux ; on calcule échantillon par échantillon une réponse impulsionnelle de forme semblable mais comprimée selon l'échelle de la variable dans un rapport correspondant à une augmentation de bande passante désirée et dilatée en amplitude dans le même rapport, et on calcule les coefficients d'un filtre à réponse impulsionnelle finie apte à fournir à sa sortie une succession d'échantillons de la réponse comprimée lorsqu'on applique à son entrée une succession d'échantillons correspondants de la réponse du système physique. Ce filtre à réponse impulsionnelle finie est incorporé au système physique, en sortie (de préférence) ou en entrée ou à l'intérieur du système, pour en améliorer la bande passante dans un rapport n.
Ainsi, on ne recherche pas un filtre qui tente de compenser intégralement la fonction de transfert du système physique mais on recherche un filtre qui donne au système physique corrigé une réponse impulsionnelle de même forme générale que sa réponse naturelle mais comprimée selon l'échelle de la variable temporelle ou spatiale et dilatée en amplitude.
Plus précisément, le procédé selon l'invention est un procédé de réalisation d'un système comprenant un ensemble d'éléments physiques limitant la bande passante du système, caractérisé en ce qu'on détermine la réponse impulsionnelle a(t) du système en fonction d'une variable t au moins entre une valeur 0 et une valeur T, on détermine au moins N échantillons successifs ai d'indice i variant de 1 à N (N entier supérieur à 1 ) dont la succession représente la forme de la réponse impulsionnelle a(t), répartis avec un pas T/N depuis la valeur t=T/N jusqu'à la valeur t=T, ainsi qu'une valeur initiale d'échantillon arbitraire a0 et des valeurs d'échantillons, éventuellement nuls, représentant approximativement cette courbe entre l'instant T et un instant n.T où n est un coefficient représentant un facteur d'augmentation de bande passante désirée, on détermine une courbe b(t) = n.a(n.t) qui est une réplique approximative de la courbe a(t), dilatée en amplitude dans le rapport n et comprimée selon l'échelle de la variable t dans le même rapport n, on établit une valeur initiale d'échantillon b0 = n.aO et on prend N échantillons bj dont la succession représente cette courbe b(t), répartis avec un pas T/N depuis la valeur t=T/N jusqu'à la valeur t=T, on réalise un filtre à réponse impulsionnelle finie à N+1 coefficients, ayant une entrée et une sortie et apte à fournir sur sa sortie les valeurs successives bj de b0 à bN des échantillons bj lorsque les échantillons successifs ai de a0 à aN sont appliqués à son entrée, et on incorpore ce filtre au système.
Lorsqu'on utilise ensuite le système, on recueille en sortie du système corrigé un signal de bande passante améliorée dans un facteur n. Le choix de l'intervalle [0,T] est tel qu'il englobe la partie la plus significative de la réponse impulsionnelle a(t), par exemple toute la partie de la réponse pour laquelle l'amplitude de la réponse est supérieure à au moins 5% du maximum de l'amplitude de la réponse. Les échantillons compris dans l'intervalle [T, nT] peuvent très bien être pris arbitrairement à zéro, même si les valeurs de a(t) ne sont pas nulles dans cet intervalle, car ces valeurs sont faibles. Il peut en résulter un moins bon résultat mais un résultat tout-à-fait acceptable.
Les coefficients du filtre sont de préférence des coefficients Ci où i est l'indice variant de 0 à N, la valeur du coefficient Ci étant définie par l'itération suivante :
C0 = bo/ao ou une valeur de préférence égale à n ou proche de n, et
Ci = (bi - ai. C0 - ai-i Ci .. .. - aμj.Cj .. ..- a-ι Ci--ι )/a0,
Si la valeur de la courbe a(t) au temps 0 est nulle ou proche de zéro, elle est remplacée par une valeur arbitraire a0 non nulle, de préférence inférieure ou égale à a^, par exemple a-ι/2.
Outre le procédé de réalisation, l'invention concerne un système électronique comportant un ensemble d'éléments physiques dont la nature induit une limitation de bande passante du système, et un filtre électronique de compensation permettant de donner au système une bande passante n fois plus grande que celle qu'aurait le système sans ce filtre, n étant un nombre supérieur à 1 , dans lequel le système dépourvu du filtre possède une réponse impulsionnelle a(t) en fonction d'une variable t, caractérisé en ce que le filtre de compensation est un filtre à réponse impulsionnelle finie à N+1 coefficients Ci où i est un indice variant de 0 à N, la valeur du coefficient Ci étant définie par l'itération suivante :
Co = bo/ao et
Ci = (bi - ai. C0 - aμi Ci .. .. - aμj.Cj .. ..- ai Ci-i )/a0, dans lequel ai est une valeur d'échantillon d'indice i variant de 1 à
N prise parmi N échantillons successifs de la réponse impulsionnelle a(t), répartis avec un pas T/N entre une valeur t= T/N et une valeur t=T, a0 étant remplacé par une valeur arbitraire non nulle si a(t) est nulle pour t=0, et dans lequel bi est une valeur d'échantillon égale à n fois la valeur d'un échantillon an i = a(n.i.TVN) de la réponse impulsionnelle a(t).
Là encore, si a0 est nul, on choisit pour a0 une autre valeur non nulle, de préférence comprise entre 0 et a-i.
Le facteur n est de préférence un nombre entier. Il est de préférence égal à 2, 3 ou 4. II peut également être non entier et dans ce cas l'indice n.i peut ne pas être entier de sorte qu'il n'existe pas véritablement d'échantillon an i. On attribue alors à l'échantillon fictif an i qui sert pour le calcul de l'échantillon bi = n.an i une valeur interpolée entre deux échantillons réels encadrant cet échantillon fictif.
D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description détaillée qui suit et qui est faite en référence aux dessins annexés dans lesquels :
- la figure 1 représente symboliquement un système physique recevant une grandeur d'entrée e et fournissant une grandeur de sortie F(e) ;
- la figure 2 représente la réponse impulsionnelle F[x(t)] du système physique ;
- la figure 3 représente le principe d'une compensation de réponse en fréquence par une fonction de transfert inverse dans le domaine des fréquences ; - la figure 4 représente le principe d'établissement d'une réponse impulsionnelle désirée à partir de la réponse impulsionnelle réelle du système ;
- la figure 5 représente le principe selon l'invention de correction de la bande passante par un filtre à réponse impulsionnelle finie associé au système physique à corriger ;
- la figure 6 représente la structure générale d'un filtre à réponse impulsionnelle finie ;
- la figure 7 représente le détail de l'élaboration de la réponse impulsionnelle désirée b(t) à partir de la réponse naturelle a(t) du système ;
- la figure 8 représente le comportement en fréquence d'un système physique non corrigé en présence d'une grandeur d'entrée e(t) en créneaux de fréquence croissante ;
- la figure 9 représente le comportement du système physique corrigé, dans les mêmes conditions qu'à la figure 8.
La figure 1 représente sous une forme symbolique (un simple rectangle) un système physique quelconque SP qui, lorsqu'on lui applique une grandeur d'entrée e en entrée fournit une grandeur de sortie F(e) en sortie. Par exemple, la grandeur d'entrée e est un signal électronique et la grandeur de sortie F(e) est un autre signal électronique. Ou bien, la grandeur d'entrée e est une luminance et la grandeur de sortie F(e) est un signal électronique, un capteur de luminance étant présent dans le système.
Dans beaucoup d'applications, on souhaiterait que le système se comporte de manière idéale, c'est-à-dire que la grandeur F(e) varie selon une fonction déterminée de e aussi bien lorsque e varie lentement que lorsque e varie rapidement. Par exemple on veut que le capteur de lumière produise un signal électronique d'amplitude proportionnelle à la luminance avec le même coefficient de proportionnalité lorsque la luminance varie lentement et lorsqu'elle varie rapidement. Mais le système est un système physique et a des limitations propres à tout système physique, en particulier des limitations d'aptitude à répondre instantanément à une variation instantanée de la grandeur d'entrée e. Il introduit des déformations de fonction F(e) lorsque e varie trop rapidement. Lorsqu'on parle ici de variation rapide, il faut comprendre qu'on peut parler de rapidité dans le temps ou dans l'espace. En effet, on doit considérer les problèmes de limitation des systèmes physiques aussi bien en termes de variation dans le temps que de variations dans l'espace comme on le détaillera plus loin.
La plupart du temps les systèmes physiques se comportent comme des filtres passe-bas, c'est-à-dire qu'ils ont une bande passante temporelle ou spatiale limitée vers le haut. Leur courbe de réponse en fréquence temporelle ou spatiale est plate depuis les fréquences basses puis chute fortement au-delà d'une certaine fréquence, alors qu'on souhaiterait idéalement qu'elle reste plate pour des fréquences bien supérieures. Mais ils peuvent avoir aussi une bande passante limitée à la fois vers le haut et vers le bas, et même avoir une courbe de réponse en fréquence qui n'est plate nulle part. Une manière de représenter le comportement de filtre du système physique consiste à déterminer sa courbe de réponse impulsionnelle, c'est-à- dire d'appliquer une impulsion normalisée de Dirac à l'entrée du système et d'observer la déformation de cette impulsion à la sortie du système. Quelle que soit la qualité du système, il y a toujours une déformation. La courbe d'amplitude de réponse en fonction du temps (si le temps est la variable concernée) ou d'une variable spatiale est une représentation du comportement général du système physique.
La figure 2 représente en 2a l'impulsion d'entrée normalisée (amplitude unitaire, durée tendant vers zéro, fronts de montée et de descente verticaux), représentée comme une fonction x(t), et en 2b la réponse impulsionnelle a(t) du système ; la réponse a(t) est la courbe F[x(t)] produite par le système lorsqu'on lui applique le signal d'entrée e=x(t). La variable considérée pour la réponse impulsionnelle est le temps t mais le raisonnement serait le même si la variable était spatiale ; ceci est d'autant plus vrai que la variable spatiale peut être ramenée à une variable temporelle par exemple lorsque les signaux variant spatialement sont lus séquentiellement. On considérera donc ici et dans toute la suite que le terme variable temporelle t est générique et est transposable directement ou indirectement à une variable spatiale. Sur la figure 3 on a représenté cette même réponse impulsionnelle dans l'espace des fréquences, où l'impulsion d'entrée x(f) a un spectre normalisé uniforme en amplitude sur toute la bande de fréquences positives et négatives (figure 3a) ; le système physique est par exemple un filtre passe-bas (figure 3b) de fonction G(f) en fonction de la fréquence, plate pour les fréquences basses, chutant pour les fréquences plus hautes ; on sait alors tracer une courbe inverse de celle de ce filtre, figure 3c, correspondant à une fonction de transfert inverse H(f) = 1 /G(f). En plaçant dans le système physique, en amont, en aval, ou au cœur du système, un filtre de compensation ayant la fonction de transfert H(f) dans l'espace des fréquences, donc une courbe de filtrage comme celle de la figure 3c, le système global aura pour fonction de transfert le produit des fonctions de la transfert du système physique G(f) et de la fonction de transfert du filtre H(f). Ce produit (figure 3d) est égal à 1 et on devrait retrouver l'impulsion d'entrée normalisée en sortie du système ainsi compensé.
La réalité n'est pas si simple car il est très difficile de réaliser un tel filtre inverse de fonction de transfert H(f) = 1 /G(f) connaissant seulement la forme de la courbe de réponse G(f). On peut faire alors un filtrage par le calcul en faisant la transformée de Fourier à chaque instant du signal électronique, en l'inversant, et en faisant la transformée de Fourier inverse, mais cela ne peut être fait qu'avec des processeurs très puissants et très rapides.
Selon l'invention on va utiliser la connaissance de la courbe de réponse impulsionnelle temporelle (ou spatiale mais, comme on l'a dit, le principe est le même) du système physique pour calculer un filtre à réponse impulsionnelle finie qui va faire correspondre à cette réponse impulsionnelle une autre réponse impulsionnelle améliorée ; la forme de la réponse améliorée correspond à une fonction de filtrage semblable à celle du système physique initial mais de bande passante plus large. La forme de réponse impulsionnelle améliorée est obtenue en dilatant l'amplitude de la réponse initiale a(t) mais en comprimant d'autant l'échelle des temps.
La figure 4 représente en 4a l'impulsion d'entrée x(t) (comme à la figure 2a), en 4b la réponse impulsionnelle a(t) du système physique initial (comme à la figure 2b), et en 4c une réponse impulsionnelle b(t) désirée, semblable en forme à celle de la figure 4b mais comprimée selon l'échelle des temps et dilatée selon l'échelle des amplitudes. La largeur à mi-hauteur de la réponse b(t) est celle de la courbe a(t) mais divisée par un facteur n. La hauteur maximale de la réponse b(t) est celle de la courbe a(t) mais multipliée par n.
Selon l'invention on détermine donc la réponse impulsionnelle du système physique, on en déduit une réponse impulsionnelle améliorée qui est une réponse semblable à la réponse initiale, dilatée en amplitude et comprimée selon l'échelle des temps, et on calcule les caractéristiques d'un filtre de compensation qui, lorsqu'il reçoit à son entrée la réponse impulsionnelle initiale produit à sa sortie la réponse impulsionnelle améliorée. Le filtre est un filtre à réponse impulsionnelle finie dont les coefficients sont faciles à calculer dès lors qu'on connaît la réponse de sortie souhaitée pour une forme de signal d'entrée donnée. La figure 5 représente le système physique SPC corrigé ayant pour réponse impulsionnelle la réponse impulsionnelle de la figure 4c. Ce système comporte le système physique initial SP et un filtre à réponse impulsionnelle finie FIR dont la fonction est de convertir la forme d'onde a(t) de la figure 4b en la forme d'onde b(t) de la figure 4c, et du même coup d'améliorer la bande passante de toute forme d'onde F(e) sortant du système physique SP.
On rappelle ce qu'est un filtre à réponse impulsionnelle finie : c'est un circuit qui traite un signal analogique échantillonné à une fréquence d'échantillonnage F, et qui fait de manière cyclique l'addition d'une série d'échantillons successifs pondérés par des coefficients individuels choisis en fonction de la fonction de filtrage à obtenir. Le nombre d'échantillons dans la série peut être d'une dizaine ou plusieurs dizaines, voire même plusieurs centaines selon la complexité de la fonction de filtrage à réaliser. L'échantillonnage du signal à traiter par le filtre peut être fait dans le filtre ou en amont du filtre selon que le système physique délivre des signaux non échantillonnés ou échantillonnés.
La figure 6 représente un filtre à réponse impulsionnelle finie. Le signal à filtrer est un signal u(t), appliqué sous forme d'échantillons successifs à la fréquence F et le signal de sortie du filtre est un signal u'(t) sortant sous forme d'échantillons. Le filtre est représenté sous forme d'un ensemble de circuits à retard représentés par les boîtes z"1 et d'un ensemble de circuits multiplicateurs permettant d'appliquer de multiplier le signal par un coefficient de pondération d propre à ce circuit. La lettre z est la variable servant classiquement à la représentation mathématique des systèmes échantillonnés et la fonction z"1 représente un retard unitaire appliqué à un échantillon, donc un retard 1/F si la fréquence d'échantillonnage est F.
Le filtre comporte encore un additionneur ADD. Comme on le comprend aisément en regardant la figure 6, le signal de sortie u'(t) du filtre est un signal échantillonné qui est la somme pondérée de la série des N+1 derniers échantillons successifs reçus. u'(t) = U0-C0+ u-i. Ci + U2-C2 +.... UN-CN s'il y a N+1 échantillons en comptant l'échantillon U0. Les échantillons successifs U0, u-ι, U2, etc jusqu'à uN représentent respectivement l'échantillon courant u(t0) à un instant t0 et les échantillons précédents u(to-dt), u(to-2dt), etc. jusqu'à u(to-N.dt). Les coefficients du filtre FIR sont calculés selon l'invention à partir de la connaissance de la réponse impulsionnelle a(t) du système physique SP. Cette réponse impulsionnelle est mesurable ou calculable à partir de la constitution connue du système physique. Le calcul des coefficients du filtre est fait en déterminant quels coefficients doivent être utilisés dans le filtre pour qu'un signal d'entrée échantillonné ayant la forme de la réponse impulsionnelle a(t) de la figure 4b sorte du filtre avec la forme échantillonnée de la réponse b(t) de la figure 4c.
Le calcul va être exposé d'abord dans un exemple simple où on souhaite multiplier par un nombre n, entier et supérieur à 1 , la bande passante du système physique SP. Typiquement, n peut être égal à 2 ou 3. Les figures qui suivent maintenant se réfèrent à un exemple dans lequel n est égal à 2.
On considère une durée T sur laquelle la réponse impulsionnelle a(t) est significative ; elle s'étend théoriquement à l'infini, mais en pratique elle tend toujours vers zéro et on ne s'intéresse pas à la zone où elle devient très faible. A titre d'exemple, la durée T peut être la durée pendant laquelle la réponse impulsionnelle est au moins égale à 5% de sa valeur maximale.
On échantillonne cette courbe par pas de T/N, N étant un nombre entier choisi arbitrairement, qui est d'autant plus élevé qu'on veut représenter plus finement la réponse impulsionnelle b(t). N peut être compris par exemple entre 10 et 100, et un choix de N=15 ou 31 ou 63 est particulièrement indiqué, pour que le nombre N+1 soit une puissance de deux.
On note a0,
Figure imgf000013_0001
, a2, ....ai, ... aN les valeurs de la réponse impulsionnelle à des instants respectifs 0, T/N, 2T/N, etc. séparés par des intervalles de temps égaux à T/N, à partir du début (temps zéro) de la réponse (qui est l'instant d'application de l'impulsion d'entrée normalisée qui donne lieu à cette réponse). L'indice i représente le numéro de l'échantillon entre 0 et N. Les valeurs de ces échantillons doivent être précisément mesurées ou calculées car ce sont ces échantillons qui constituent la base de la représentation de la réponse du système physique et aussi de la réponse améliorée comme on va le voir. L'échantillon a0 d'indice i=0 est traité différemment comme on le verra et n'est pas forcément égal à la valeur a(t) pour t=0, en particulier dans le cas fréquent où cette valeur est nulle, de sorte que la succession d'échantillons de la courbe a(t) entre les instants t=0 et t=T peut être considérée comme composée d'une valeur arbitraire d'échantillon initial a0 et de N valeurs d'échantillons a\ dont la succession représente la courbe a(t). On détermine aussi des échantillons an i de la courbe a(t) entre le temps T et le temps n.T. Mais les échantillons au-delà du temps T peuvent aussi être mis arbitrairement à zéro car la queue de la réponse impulsionnelle, au delà du temps T n'est pas significative. Ces échantillons au-delà du temps T seront utilisés pour déterminer la réponse impulsionnelle améliorée, même si leur valeur est nulle.
A partir de la courbe a(t) on construit une courbe b(t) qui est la réponse impulsionnelle améliorée désirée.
Pour que la courbe b(t) soit une réponse impulsionnelle semblable à la courbe a(t) mais dilatée en amplitude dans le rapport n et comprimée en temps dans le même rapport n, il faut que la valeur b(t) à un instant t soit n fois la valeur de la courbe a(t) à l'instant n.t.
Par conséquent, à partir d'un échantillon an i qui est la valeur de la courbe a(t) à l'instant t=n.i.T/N, on établit un échantillon bj qui est la valeur de la courbe b(t) à l'instant i.T/N en prenant pour valeur bj la valeur: bi = n.an i La figure 7 illustre cette constitution de la courbe b(t) à partir d'échantillons de la courbe a(t) dans le cas où n=2 et N=1 6. La réponse initiale a(t), connue initialement, est en traits pleins ; elle comporte un échantillon initial arbitraire a0 à l'instant t=0 et N échantillons entre TVN et T et d'autres échantillons entre T et nT, ces derniers pouvant être mis à zéro éventuellement ; la réponse améliorée b(t) est en traits pointillés : c'est une courbe établie par lissage à partir des N+1 échantillons bj calculés selon la formule ci-dessus. A partir des échantillons ai et bj on calcule les coefficients du filtre à réponse impulsionnelle finie qui permet d'obtenir à sa sortie la réponse b(t) lorsqu'on applique à son entrée un signal dont la variation temporelle est celle de la réponse a(t). Lorsqu'on utilisera un tel filtre en lui appliquant non plus un signal a(t) mais un signal F[e(t)] provenant de la sortie du système physique, ce filtre délivrera un signal F'[e(t)] de bande passante améliorée dans le rapport n.
Le calcul des coefficients C0, C-i , ...C1, ... CN du filtre est un calcul itératif reposant sur le fait que la sortie b(t) du filtre est fonction de la série d'échantillons successifs aj appliqués à l'entrée : b(t) = ao.Co+ a-ι.C-1 + a2.C2 +.... aN.CN
Avant l'instant 0, les échantillons d'entrée sont tous nuls. Il en résulte les égalités suivantes : bo = ao Co + O + 0 + 0 + 0
Figure imgf000014_0001
bj = ai. C0 + aj.-i .Ci + +ao.Cj + 0 +0
bn = an.C0 + an-i .Ci + ao.Cn
et par conséquent :
C0 = bo/ao cette valeur est en principe égale à n du fait de la dilatation de la courbe dans le rapport n. C'est parce qu'on ne pourrait pas calculer bo/ao si a0 était nul ou trop proche de zéro qu'on donne alors à l'échantillon a0 une valeur arbitraire ou qu'on choisit tout simplement pour le premier coefficient du filtre la valeur C0 = n (ou à la rigueur une valeur proche de n) sans avoir même besoin de choisir une valeur d'échantillon a0.
Ci = (bi - a-|.C0)/a0 = (2.a2 - a-|.C0)/a0 dans le cas où n=2 d'après la formule bi = n.an i
Ci = (n.an i - ai. C0 - ai.-i .Ci - ai .Ci-i)/a0
CN = (n.an N - aN-Co - aN-i .Ci - ai.Cι\ι-i)/ao
On rappelle que pour ce calcul les échantillons ai au-delà du temps T peuvent être arbitrairement pris à zéro si on ne souhaite pas utiliser les échantillons réels de la courbe a(t) qui sont de toutes façons faibles puisque la partie significative de la réponse se situe entre les instants 0 et T.
Il peut arriver que la valeur de la courbe a(t) à l'instant t=0 soit nulle ou très proche de zéro. Dans ce cas, on forcera la valeur d'échantillon a0 à une valeur différente de zéro, par exemple à la valeur aV2 qui est une valeur intermédiaire entre 0 et ai. Et on prendra encore b0 = n fois la valeur forcée de a0 pour conserver le rapport bo/ao égal à n.
On aboutit ainsi à un filtre de compensation à réponse impulsionnelle finie à N+1 coefficients qui augmente dans un facteur n la bande passante du système physique initial. En fonctionnement, ce filtre reçoit la sortie F[e(t)] du système physique dont on veut augmenter la bande passante, et il fournit une sortie F'[e(t)] à bande passante augmentée. On pourrait envisager aussi que le filtre de compensation FIR soit placé dans le système physique plutôt qu'à la sortie de celui-ci, c'est-à-dire qu'on peut éventuellement placer le filtre en amont des éléments physiques qui tendent à réduire la bande passante.
Dans ce qui précède n a été considéré comme étant un entier, de préférence égal à 2 ou 3 ou 4. On pourrait envisager aussi que n soit une valeur non entière, par exemple 2,5 ou 3,5 (ou même un nombre quelconque bien que cela ait en pratique peu d'intérêt).
On conserve le même principe d'établissement d'une courbe b(t) correspondant à la courbe a(t) mais dilatée en amplitude d'un facteur n et comprimée dans le temps du même facteur n. La différence réside dans l'obtention des échantillons bj qui vont servir à constituer la courbe b(t) et qui vont servir au calcul des coefficients du filtre.
On ne dispose plus de valeurs d'échantillons an i lorsque n.i n'est pas un nombre entier. Dans ce cas, on pourra utiliser comme valeur d'échantillon an i une valeur interpolée entre les deux échantillons réels les plus proches de l'échantillon inexistant de rang n.i
Par exemple, si n=2,5 et i = 3, on établira un échantillon fictif a7,5 qui sera une interpolation entre l'échantillon a7 et l'échantillon a8, par exemple simplement la moyenne a7,5 = (a7 + a8)/2 des deux échantillons connus adjacents. Par conséquent, dans cet exemple, on trouvera b3 =
2,5(a7 + as)/2 et, dans le calcul des coefficients du filtre, on utilisera les valeurs d'échantillons connus de rang i = 0 à N pour les échantillons ai comme pour les échantillons bj.
La figure 8 représente une simulation effectuée pour un système physique SP dont la fonction de transfert est assimilée à un filtre passe-bas du premier ordre. Une forme d'onde en créneau rectangulaire d'entrée (en trait plein), à une fréquence qui est supérieure à une fréquence de coupure du filtre et qui est variable dans le temps (la durée des créneaux va en diminuant) est appliquée comme grandeur e(t) à l'entrée du système physique SP. La sortie du système est une grandeur F[e(t)] qui est représentée par une courbe tiretée. Du fait du filtrage passe-bas, on voit que la sortie n'arrive pas à suivre fidèlement les fronts de variations en créneau. La figure 9 représente la sortie F'[e(t)] du système corrigé par le filtre FIR calculé de la manière expliquée précédemment, les mêmes créneaux e(t) étant appliqués à l'entrée. On voit que cette sortie suit beaucoup mieux les créneaux et ne se détériore que pour les créneaux les plus brefs.
Le filtre FIR qui sert dans l'exemple de la figure 9 a été calculé pour n=10, avec N=15 (filtre à 16 coefficients). Les coefficients calculés sont les suivants :
C0 = I O C8 = - 0.0000176
C1 = -6.376 C9 = -3.917x10"5
C2 = -1 .422 C10 = -8.741 x10"6
C3 = -0.317 C11 = -1 .950 x10"6 C4 = -0,070 C12 = -4.352x10"7 C5 = -0.0158 C13 = -9.710x10"8
C6 = -0.00353 C14 = -2.1 67x10"8
C7 = -0.000786 C15 = -4.835x10"9
Ce filtre n'est donné qu'à titre d'exemple illustratif montrant bien l'amélioration considérable qu'on peut obtenir sur la bande passante. Dans la pratique le système physique ne sera pas aussi simple qu'un filtre du premier ordre, le nombre de coefficients N du filtre sera plutôt de 32 ou 64 et le facteur n d'augmentation de bande passante sera plutôt compris entre 2 et 4.

Claims

REVENDICATIONS
1. Procédé d'augmentation de la bande passante d'un système physique comprenant un ensemble d'éléments physiques (SP) limitant la bande passante du système, caractérisé en ce qu'on détermine la réponse impulsionnelle a(t) du système physique en fonction d'une variable temporelle ou spatiale t entre des valeurs t=0 et t=T de la variable, on calcule échantillon par échantillon une réponse impulsionnelle de forme semblable mais comprimée selon l'échelle de la variable dans un rapport n correspondant à une augmentation de bande passante désirée et dilatée en amplitude dans le même rapport n, on détermine les coefficients d'un filtre à réponse impulsionnelle finie apte à fournir à sa sortie une succession d'échantillons de la réponse comprimée lorsqu'on applique à son entrée une succession d'échantillons correspondants de la réponse du système physique, et on incorpore ce filtre à réponse impulsionnelle finie au système physique.
2. Procédé selon la revendication 1 , caractérisé en ce qu'on détermine au moins N échantillons successifs ai d'indice i variant de 1 à N (N entier supérieur à 1 ) dont la succession représente la forme de la réponse impulsionnelle a(t), répartis avec un pas T/N depuis la valeur t=T/N jusqu'à la valeur t=T, ainsi qu'une valeur initiale d'échantillon arbitraire a0 et des valeurs d'échantillons, éventuellement nuls, représentant approximativement cette courbe entre l'instant T et un instant n.T où n est un coefficient représentant un facteur d'augmentation de bande passante désirée, on détermine une courbe b(t) = n.a(n.t) qui est une réplique approximative de la courbe a(t), dilatée en amplitude dans le rapport n et comprimée selon l'échelle de la variable t dans le même rapport n, on établit une valeur initiale d'échantillon b0 = n.aO et on prend Néchantillons bj dont la succession représente cette courbe b(t), répartis avec un pas T/N depuis la valeur t=T/N jusqu'à la valeur t=T, et en ce que le filtre à réponse impulsionnelle finie est un filtre à N+1 coefficients, ayant une entrée et une sortie et apte à fournir sur sa sortie les valeurs successives bj de b0 à bN des échantillons bj lorsque les échantillons successifs ai de a0 à aN sont appliqués à son entrée.
3. Procédé selon la revendication 1 , caractérisé en ce que les coefficients du filtre sont des coefficients Ci où i est l'indice variant de 0 à N, la valeur du coefficient Ci étant définie par l'itération suivante : C0 = bo/ao ou une valeur proche de n, et
Ci = (bi - ai. C0 - ai-i Ci .. .. - aμj.Cj .. ..- a-ι Ci--ι )/a0,
4. Procédé selon l'une des revendications 1 et 2, caractérisé en ce que la valeur a0 choisie comme première valeur d'échantillon pour calculer les coefficients est choisie égale à une valeur comprise entre la valeur de a(t) pour t=0 et la valeur du deuxième échantillon a^ dans le cas où a(t) pour t=0 est une valeur nulle ou proche de zéro.
5. Procédé selon l'une des revendications précédentes caractérisé en ce que n est choisi entre 2 et 4.
6. Procédé selon l'une des revendications précédentes, caractérisé en ce que n est un nombre entier.
7. Procédé selon l'une des revendications 1 à 5, caractérisé en ce que n n'est pas entier et les valeurs d'échantillons an i d'indice n.i non entier sont des valeurs interpolées entre des échantillons voisins d'indice entier.
8. Système électronique comportant un ensemble d'éléments physiques (SP) dont la nature induit une limitation de bande passante du système, et un filtre électronique de compensation (FIR) permettant de donner au système une bande passante n fois plus grande que celle qu'aurait le système sans ce filtre, n étant un nombre supérieur à 1 , dans lequel le système dépourvu du filtre possède une réponse impulsionnelle a(t) en fonction d'une variable t, caractérisé en ce que le filtre de compensation est un filtre à réponse impulsionnelle finie à N+1 coefficients Ci où i est un indice variant de 0 à N, la valeur du coefficient C1 étant définie par l'itération suivante : C0 = bo/ao ou une valeur proche de n Ci = (bi - ai. C0 - ai-i Ci .. .. - aμj.Cj .. ..- a-ι Ci--ι )/a0, dans lequel a0 est une valeur arbitraire, ai est une valeur d'échantillon d'indice i prise parmi N échantillons successifs de la réponse impulsionnelle a(t), répartis avec un pas T/N entre une valeur t= T/N et une valeur t=T, et dans lequel bi est une valeur d'échantillon égale à n fois la valeur d'un échantillon an i = a(n.i.TVN) de la réponse impulsionnelle a(t).
9. Système électronique selon la revendication 8, caractérisé en ce que la valeur choisie pour a0 dans l'itération est égale à une valeur comprise entre la valeur de a(t) pour t=0 et la valeur du deuxième échantillon a-\ dans le cas où a(t) pour t=0 est une valeur nulle ou proche de zéro.
10. Système selon la revendication 9, caractérisé en ce que le facteur n est un nombre entier de préférence égal à 2, 3 ou 4.
11. Système selon la revendication 9, caractérisé en ce que n n'est pas entier et les valeurs d'échantillons ani d'indice n.i non entier sont des valeurs interpolées entre des échantillons voisins d'indice entier.
PCT/EP2009/051424 2008-02-19 2009-02-09 Procede et dispositif pour ameliorer la bande passante d'un systeme physique WO2009103624A1 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
EP09712207A EP2243218A1 (fr) 2008-02-19 2009-02-09 Procede et dispositif pour ameliorer la bande passante d'un systeme physique
US12/918,337 US20100332577A1 (en) 2008-02-19 2009-02-09 Method and device for improving the passband of a physical system

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0800885A FR2927741B1 (fr) 2008-02-19 2008-02-19 Procede et dispositif pour ameliorer la bande passante d'un systeme physique
FR0800885 2008-02-19

Publications (1)

Publication Number Publication Date
WO2009103624A1 true WO2009103624A1 (fr) 2009-08-27

Family

ID=39791282

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2009/051424 WO2009103624A1 (fr) 2008-02-19 2009-02-09 Procede et dispositif pour ameliorer la bande passante d'un systeme physique

Country Status (4)

Country Link
US (1) US20100332577A1 (fr)
EP (1) EP2243218A1 (fr)
FR (1) FR2927741B1 (fr)
WO (1) WO2009103624A1 (fr)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5208596A (en) * 1992-04-10 1993-05-04 Rca Thomson Licensing Corporation DAC distortion compensation
US20030043894A1 (en) * 2001-08-24 2003-03-06 Sestok Charles K. Time domain equalizer for DMT modulation
US20060224365A1 (en) * 2005-04-01 2006-10-05 Pickerd John J Oscilloscope having an enhancement filter

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4703447A (en) * 1985-04-05 1987-10-27 The Grass Valley Group, Inc. Mixer controlled variable passband finite impulse response filter
US5168459A (en) * 1991-01-03 1992-12-01 Hewlett-Packard Company Adaptive filter using continuous cross-correlation
US5982305A (en) * 1997-09-17 1999-11-09 Microsoft Corporation Sample rate converter
JP4113758B2 (ja) * 2002-10-29 2008-07-09 株式会社イシダ 重量測定装置、ノイズ除去方法およびディジタルフィルタの設計方法
US7584235B2 (en) * 2004-03-18 2009-09-01 Tektronix, Inc. Variable passband autoregressive moving average filter

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5208596A (en) * 1992-04-10 1993-05-04 Rca Thomson Licensing Corporation DAC distortion compensation
US20030043894A1 (en) * 2001-08-24 2003-03-06 Sestok Charles K. Time domain equalizer for DMT modulation
US20060224365A1 (en) * 2005-04-01 2006-10-05 Pickerd John J Oscilloscope having an enhancement filter

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
AURAND J F ET AL: "Bandwidth enhancement of single-shot data acquisition systems by deconvolution of system response", PRECISION ELECTROMAGNETIC MEASUREMENTS, 1994. DIGEST., 1994 CONFERENCE ON BOULDER, CO, USA 27 JUNE-1 JULY 1994, NEW YORK, NY, USA,IEEE, 27 June 1994 (1994-06-27), pages 282, XP010123952, ISBN: 978-0-7803-1984-4 *
FREDDY BA HLI: "Network Synthesis by Impulse Response for Specified Input and Output in the Time Domain", July 1953 (1953-07-01), Cambridge, MA, USA, pages 1 - 65, XP002536082, Retrieved from the Internet <URL:http://dspace.mit.edu/bitstream/handle/1721.1/4810/RLE-TR-261-14267258.pdf?sequence=1> [retrieved on 20090708] *
SEDKI M. RIAD: "The Deconvolution Problem: An Overview", PROCEEDINGS OF THE IEEE, vol. 74, no. 1, January 1986 (1986-01-01), USA, pages 82 - 85, XP002536083 *
STEVEN W. SMITH: "The Scientist and Engineer's Guide to Digital Signal Processing; Chapter 17 - Custom Filters", 1997, XP002499860 *
TAPAN K. SARKAR; DONALD D. WEINER; VIJAY K. JAIN; SOHEIL A. DIANAT: "Impulse Response Determination in the Time Domain - Theory", IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, vol. AP-30, no. 4, July 1982 (1982-07-01), USA, pages 657 - 663, XP002536081 *

Also Published As

Publication number Publication date
FR2927741B1 (fr) 2011-08-05
EP2243218A1 (fr) 2010-10-27
FR2927741A1 (fr) 2009-08-21
US20100332577A1 (en) 2010-12-30

Similar Documents

Publication Publication Date Title
CA2207670A1 (fr) Procede de synthese d&#39;un filtre numerique a reponse impulsionnelle finie et filtre obtenu selon le procede
WO2012019251A2 (fr) Procede d&#39;analyse de signaux fournissant des frequences instantanees et des transformees de fourier glissantes et dispositif d&#39;analyse de signaux
EP0131521B1 (fr) Procédé et installation d&#39;analyse et de restitution de signal à échantillonnage et interpolation
EP2177020B1 (fr) Procédé de correction de sensibilité et capteur d&#39;image matriciel pour la mise en oeuvre de ce procédé
FR2714549A1 (fr) Filtre numérique et convertisseur analogique-numérique à suréchantillonnage employant ce filtre.
WO2009103624A1 (fr) Procede et dispositif pour ameliorer la bande passante d&#39;un systeme physique
EP2887307B1 (fr) Procédé de traitement d&#39;images, notamment issues de systèmes de visualisation de nuit et système associé
EP0688476B1 (fr) Procede et dispositif de filtrage d&#39;un signal temporel numerique, et application a la correction d&#39;echos dans un canal de transmission
EP1067392B1 (fr) Procédé de tests de composants électroniques
EP0970562B1 (fr) Filtre numerique pour retards fractionnaires
EP0732809B1 (fr) Dispositif de filtrage digital
US11069331B2 (en) Noise reduction filter for signal processing
EP0140396B1 (fr) Procédé et dispositif de reconstitution d&#39;un signal analogique
EP2394424B1 (fr) Procédé de correction d&#39;images obtenues par une matrice de points photosensibles
FR2840087A1 (fr) Depollution d&#39;image numerique bruitee et comprimee
EP4113131B1 (fr) Estimation, malgre une fraude, de la puissance consommee sur une phase
EP0014604B1 (fr) Filtre à échantillonnage
EP1346231A1 (fr) Procedes et appareils pour le traitement analogique du signal emis par un detecteur de particules
EP2135220B1 (fr) Procede de correction de bruit spatial de capteur d&#39;image par bornage des luminances
EP0031270A1 (fr) Dispositif de filtrage à faible distorsion de phase et circuit de traitement de signaux de télévision couleur comportant un tel dispositif
EP4020978A1 (fr) Procede d&#39;etalonnage d&#39;une matrice de photodetecteurs, dispositif d&#39;etalonnage et systeme d&#39;imagerie associes
EP0332515A1 (fr) Dispositif à transfert de charges pour le traitement de blocs de données analogiques
BE1020218A3 (fr) Procede pour ameliorer la resolution temporelle des informations fournies par un filtre compose et dispositif correspondant.
FR2885758A1 (fr) Systeme d&#39;acquisition multimode
FR3056858A1 (fr) Dispositif de generation d&#39;un signal porteur continu a partir d&#39;un signal impulsionnel de reference

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09712207

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2009712207

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 12918337

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE