DETERMINATION DE LA FREQUENCE FONDAMENTALE D'UN SIGNAL PERIODIQUE INCLUANT DES COMPOSANTES HARMONIQUES DETERMINATION OF THE FUNDAMENTAL FREQUENCY OF A PERIODIC SIGNAL INCLUDING HARMONIC COMPONENTS
La présente invention concerne un procédé permettant de déterminer de façon très précise la fréquence fondamentale d'un signal périodique incluant des composantes harmoniques. L'invention concerne également un produit logiciel permettant la mise en œuvre du procédé, ainsi qu'un compteur d'énergie électrique comportant un tel produit logiciel. The present invention relates to a method for very accurately determining the fundamental frequency of a periodic signal including harmonic components. The invention also relates to a software product for implementing the method, as well as an electrical energy meter comprising such a software product.
Un domaine d'application particulièrement envisagé dans l'invention, bien que non limitatif, est celui des compteurs électriques utilisés chez des particuliers ou en industrie pour comptabiliser l'énergie électrique délivrée par un réseau de distribution d'énergie électrique et effectivement utilisée par un utilisateur. A field of application particularly contemplated in the invention, although not limiting, is that of electric meters used in private homes or in industry to account for the electrical energy delivered by a distribution network of electrical energy and actually used by a customer. user.
Aujourd'hui, les réseaux électriques connus utilisent un courant alternatif sinusoïdal, monophasé ou triphasé, à une fréquence fondamentale déterminée, égale à 50 Hz pour les réseaux européens, et à 60 Hz pour les réseaux américains. Néanmoins, les signaux qui transitent dans un réseau de distribution d'électricité subissent des perturbations, dont certaines, de nature continue, et connues sous le nom de pollution harmonique, peuvent être très gênantes. Les signaux électriques du réseau ne sont ainsi pas des signaux sinusoïdaux purs à la fréquence fondamentale, mais des signaux périodiques comprenant non seulement la fréquence fondamentale, mais aussi des harmoniques, c'est-à-dire des multiples de la fréquence fondamentale, générées par des équipements électroniques ou électriques (appareils domestiques : téléviseurs, ordinateurs, imprimantes, fours à micro-ondes, lampes à décharge..., charges industrielles : variateurs de vitesse, fours à arc, postes à souder...) raccordés réseau. Today, the known electrical networks use a sinusoidal alternating current, single-phase or three-phase, at a determined fundamental frequency, equal to 50 Hz for European networks, and 60 Hz for US networks. Nevertheless, the signals that pass through an electricity distribution network undergo disturbances, some of which, of a continuous nature, and known as harmonic pollution, can be very troublesome. The electrical signals of the network are thus not pure sinusoidal signals at the fundamental frequency, but periodic signals comprising not only the fundamental frequency, but also harmonics, that is to say multiples of the fundamental frequency, generated by electronic or electrical equipment (domestic appliances: televisions, computers, printers, microwave ovens, discharge lamps ..., industrial loads: variable speed drives, arc furnaces, welding machines, etc.) connected to the network.
Pour pouvoir être en mesure de déterminer la présence de pollueurs harmoniques, certains compteurs électriques sont dotés de moyens logiciels permettant de faire une analyse harmonique, typiquement par transformée de Fourier discrète, des signaux réseau en courant et en tension reçus sur chacune de leurs phases, et de déduire une valeur appelée taux de distorsion ou THD (initiales anglo-saxonnes mises pour Total Harmonie Distorsion),
représentant le rapport de la valeur efficace des harmoniques à celle de la grandeur alternative fondamentale (courant ou tension). In order to be able to determine the presence of harmonic polluters, some electric meters are equipped with software means for performing a harmonic analysis, typically by discrete Fourier transform, network signals in current and voltage received on each of their phases, and to deduce a value called distortion rate or THD (Anglo-Saxon initials set for Total Harmony Distortion), representing the ratio of the rms value of the harmonics to that of the fundamental alternating quantity (current or voltage).
Aujourd'hui, il serait souhaitable d'aller plus loin dans l'analyse de façon à pouvoir notamment faire une discrimination entre les équipements pollueurs et les équipements pollués. Cela ne peut se faire qu'en déterminant avec précision les amplitudes, et surtout les phases à toutes les fréquences harmoniques. Pour déterminer les fréquences harmoniques présentes dans un signal, il convient de déterminer au préalable de façon très précise la fréquence fondamentale du signal. Today, it would be desirable to go further in the analysis so as to discriminate between polluting equipment and polluted equipment. This can only be done by accurately determining the amplitudes, and especially the phases at all harmonic frequencies. In order to determine the harmonic frequencies present in a signal, the fundamental frequency of the signal must be very precisely determined beforehand.
A ce jour, certains compteurs peuvent estimer la fréquence fondamentale à partir d'une détection des passages à zéro du signal. Cette technique garantit une précision maximale au dixième d'Hertz près. Ceci est néanmoins insuffisant pour déterminer les amplitudes, et surtout les phases des harmoniques avec la précision nécessaire. To date, some counters can estimate the fundamental frequency from a detection of zero crossings of the signal. This technique guarantees maximum precision to the nearest tenth of Hertz. This is nevertheless insufficient to determine the amplitudes, and especially the phases of the harmonics with the necessary precision.
La présente invention a pour but de proposer un procédé permettant de déterminer de façon beaucoup plus précise la fréquence fondamentale d'un signal périodique ayant un contenu harmonique. The present invention aims to provide a method for determining much more accurately the fundamental frequency of a periodic signal having a harmonic content.
Pour ce faire, l'invention a pour objet un procédé de détermination de la fréquence fondamentale d'un signal périodique incluant des composantes harmoniques, caractérisé en ce qu'il comporte les étapes successives suivantes : To do this, the subject of the invention is a method for determining the fundamental frequency of a periodic signal including harmonic components, characterized in that it comprises the following successive steps:
- Une première étape d'échantillonnage et de pondération dudit signal à une fréquence d'échantillonnage prédéfinie pour délivrer un nombre déterminé N d'échantillons du signal; A first step of sampling and weighting said signal at a predefined sampling frequency to deliver a determined number N of samples of the signal;
- Une deuxième étape d'estimation grossière de la fréquence fondamentale dans laquelle une première valeur estimée de la fréquence fondamentale est calculée à partir desdits échantillons du signal, la première valeur étant exprimée dans une unité de fréquence choisie de sorte que la partie entière de la première valeur estimée comporte au plus trois chiffres ; A second step of rough estimation of the fundamental frequency in which a first estimated value of the fundamental frequency is calculated from said samples of the signal, the first value being expressed in a selected frequency unit so that the entire portion of the first estimated value has not more than three digits;
- Une troisième étape d'estimation au dixième près, dans ladite unité, de la fréquence fondamentale dans laquelle une deuxième valeur estimée au
dixième près, dans ladite unité, de la fréquence fondamentale est déterminée à partir de ladite première valeur, la troisième étape consistant à : o rechercher, parmi un premier ensemble de valeurs de fréquence possibles correspondant à des variations, par pas constant de 0,1, de la première valeur estimée, une première valeur de fréquence pour laquelle l'amplitude correspondante du signal, calculée à partir des coefficients de Fourier discrets de rang 1, est maximale, et o faire correspondre la deuxième valeur estimée au dixième près, dans ladite unité, à ladite première valeur de fréquence. A third step of estimating, to the nearest tenth, in said unit, the fundamental frequency in which a second value estimated at tenth, in said unit, the fundamental frequency is determined from said first value, the third step of: o finding, out of a first set of possible frequency values corresponding to variations, by constant step of 0.1 of the first estimated value, a first frequency value for which the corresponding amplitude of the signal, calculated from the discrete Fourier coefficients of rank 1, is maximum, and o to match the second estimated value to the nearest tenth, in said unit, at said first frequency value.
Ainsi, dans le cas où l'estimation grossière ne permet d'obtenir qu'une valeur décimale de fréquence fondamentale avec au plus un chiffre après la virgule, le procédé selon l'invention offre la possibilité d'affiner le résultat et de déterminer la fréquence fondamentale au dixième près, dans l'unité considérée. Thus, in the case where the rough estimate only makes it possible to obtain a decimal value of fundamental frequency with at most one digit after the decimal point, the method according to the invention offers the possibility of refining the result and determining the fundamental frequency to the nearest tenth, in the unit under consideration.
Dans une implémentation préférée, l'étape de recherche de la première valeur de fréquence de la troisième étape d'estimation au dixième près, dans ladite unité, comportera avantageusement les sous étapes successives suivantes : In a preferred implementation, the step of finding the first frequency value of the third estimation step to the nearest tenth, in said unit, will advantageously comprise the following successive substeps:
1) Initialisation de la valeur courante de la fréquence fondamentale à la première valeur estimée et d'une valeur de fréquence possible dudit premier ensemble à ladite première valeur estimée incrémentée de 0,1 ; 1) initializing the current value of the fundamental frequency to the first estimated value and a possible frequency value of said first set to said first estimated value incremented by 0.1;
2) Calcul d'une première amplitude du signal à la valeur courante et d'une deuxième amplitude du signal à la valeur de fréquence possible, à partir des coefficients de Fourier discrets de rang 1 ; 2) calculating a first amplitude of the signal at the current value and a second amplitude of the signal at the possible frequency value, from the discrete Fourier coefficients of rank 1;
3) Comparaison de la première amplitude et de la deuxième amplitude calculées ; 3) Comparison of the first amplitude and the second amplitude calculated;
4) Si la deuxième amplitude est supérieure à la première amplitude,
4i) remplacement de la valeur courante par la valeur de fréquence possible, puis incrémentation de la valeur de fréquence possible dudit premier ensemble de 0,1; 4) If the second amplitude is greater than the first amplitude, 4i) replacing the current value by the possible frequency value, then incrementing the possible frequency value of said first set of 0.1;
4ii) calcul d'une troisième amplitude du signal à la valeur courante et d'une quatrième amplitude du signal à la valeur de fréquence possible, à partir des coefficients de Fourier discrets de rang 1 ; 4ii) calculating a third amplitude of the signal at the current value and a fourth amplitude of the signal at the possible frequency value, from the discrete Fourier coefficients of rank 1;
4iii) réitération des sous-étapes 4i) à 4ii) tant que la quatrième amplitude est supérieure à la troisième amplitude calculées ; 4iii) reiterating the substeps 4i) to 4ii) as long as the fourth amplitude is greater than the calculated third amplitude;
5) Sinon, 5) Otherwise,
5i) remplacement de la valeur de fréquence possible par la valeur courante décrémentée de 0,1; 5i) replacement of the possible frequency value by the current value decremented by 0.1;
5ii) calcul d'une troisième amplitude du signal à la valeur courante et d'une quatrième amplitude du signal à la valeur de fréquence possible, à partir des coefficients de Fourier discrets de rang 1 ; 5ii) calculating a third amplitude of the signal at the current value and a fourth amplitude of the signal at the possible frequency value, from the discrete Fourier coefficients of rank 1;
5iii) tant que la quatrième amplitude est supérieure à la troisième amplitude calculées, remplacement de la valeur courante par la valeur de fréquence possible, décrémentation de la valeur de fréquence possible dudit premier ensemble de 0,1 et réitération de la sous-étape 5ii); 5iii) as long as the fourth amplitude is greater than the calculated third amplitude, replacing the current value by the possible frequency value, decrementing the possible frequency value of said first set of 0.1 and reiterating substep 5ii) ;
6) Faire correspondre la première valeur de fréquence recherchée à la dernière valeur courante à l'issue des sous étapes 4) ou 5) pour laquelle la comparaison est positive. 6) Match the first desired frequency value to the last current value at the end of sub-steps 4) or 5) for which the comparison is positive.
On optimise ainsi le nombre de calculs nécessaires. This optimizes the number of calculations required.
Dans le cas où l'on souhaite affiner encore le résultat, le procédé selon l'invention comporte avantageusement une quatrième étape d'estimation au centième près, dans ladite unité, de la fréquence fondamentale dans laquelle une troisième valeur estimée au centième près, dans ladite unité, de la fréquence fondamentale est déterminée à partir de ladite deuxième valeur estimée au dixième près, dans ladite unité, la quatrième étape consistant à :
o rechercher, parmi un deuxième ensemble de valeurs de fréquence possibles correspondant à des variations, par pas constant de 0,01, de la deuxième valeur estimée au dixième près, dans ladite unité, une deuxième valeur de fréquence pour laquelle l'amplitude correspondante du signal, calculée à partir des coefficients de Fourier discrets de rang 1, est maximale, et o faire correspondre la troisième valeur estimée au centième près, dans ladite unité, à ladite deuxième valeur de fréquence. In the case where it is desired to further refine the result, the method according to the invention advantageously comprises a fourth step of estimating to the nearest hundredth, in said unit, the fundamental frequency in which a third value estimated to the nearest hundredth, in said unit, the fundamental frequency is determined from said second value estimated to the nearest tenth, in said unit, the fourth step consisting of: o searching, among a second set of possible frequency values corresponding to variations, by constant pitch of 0.01, of the second value estimated to the nearest tenth, in said unit, a second frequency value for which the corresponding amplitude of the signal, calculated from the discrete Fourier Fourier coefficients of rank 1, is maximum, and o match the third value estimated to the nearest hundredth, in said unit, to said second frequency value.
L'étape de recherche de la deuxième valeur de fréquence de la quatrième étape d'estimation au centième près, dans ladite unité, comportera avantageusement des sous étapes similaires à celles implémentées pour la troisième étape d'estimation, ici encore de façon à optimiser le nombre de calculs nécessaires. The step of searching for the second frequency value of the fourth estimation step to the nearest hundredth, in said unit, will advantageously comprise substeps similar to those implemented for the third estimation step, again in order to optimize the number of calculations needed.
Une précision au millième près, dans ladite unité, peut être encore avantageusement obtenue en prévoyant une cinquième étape d'estimation au millième près, dans ladite unité, de la fréquence fondamentale dans laquelle une quatrième valeur estimée au millième près, dans ladite unité, de la fréquence fondamentale est déterminée à partir de ladite troisième valeur estimée au centième près, dans ladite unité, la cinquième étape consistant à : An accuracy to the nearest one thousandth in said unit can be further advantageously obtained by providing a fifth step of estimating to the nearest thousandth, in said unit, the fundamental frequency in which a fourth value estimated to within a thousandth, in said unit, of the fundamental frequency is determined from said third value estimated to the nearest hundredth, in said unit, the fifth step consisting in:
o rechercher, parmi un troisième ensemble de valeurs de fréquence possibles correspondant à des variations, par pas constant de 0,001 de la troisième valeur estimée au centième près, dans ladite unité, une troisième valeur de fréquence pour laquelle l'amplitude correspondante du signal, calculée à partir des coefficients de Fourier discrets de rang 1, est maximale, et o faire correspondre la quatrième valeur estimée au millième près, dans ladite unité, à ladite troisième valeur de fréquence. o searching, among a third set of possible frequency values corresponding to variations, in a constant pitch of 0.001 of the third value estimated to the nearest hundredth, in said unit, a third frequency value for which the corresponding amplitude of the signal, calculated from the discrete Fourier coefficients of rank 1, is maximum, and o match the fourth estimated value to the nearest thousandth, in said unit, said third frequency value.
La présente invention a également pour objet un produit logiciel destiné à être mis en œuvre par un microprocesseur ou un microcontrôleur, et réalisant les étapes du procédé selon l'invention. The present invention also relates to a software product intended to be implemented by a microprocessor or a microcontroller, and carrying out the steps of the method according to the invention.
Une application particulière de l'invention concerne la détermination de
la fréquence fondamentale de la tension de phase et des courants de phase et de neutre d'un réseau de distribution d'énergie électrique comportant au moins une phase et un neutre. A cet effet, un autre objet de l'invention est un compteur d'énergie électrique recevant une tension entre au moins une phase et un neutre ainsi que des courants de phase et de neutre d'un réseau de distribution d'énergie électrique, caractérisé en ce qu'il comporte des moyens de conversion analogiques-numériques de la tension et des courants de phase, un microcontrôleur et des moyens logiciels mettant en œuvre le procédé selon l'invention pour la détermination de la fréquence fondamentale de la tension et des courants de phase et de neutre A particular application of the invention concerns the determination of the fundamental frequency of the phase voltage and the phase and neutral currents of an electrical energy distribution network comprising at least one phase and one neutral. For this purpose, another object of the invention is an electrical energy meter receiving a voltage between at least one phase and a neutral as well as phase and neutral currents of an electrical energy distribution network, characterized in that it comprises means for analog-to-digital conversion of the voltage and phase currents, a microcontroller and software means implementing the method according to the invention for determining the fundamental frequency of the voltage and currents. live and neutral
La présente invention sera mieux comprise au vu de la description suivante faite en référence aux figures annexées, dans lesquelles : The present invention will be better understood from the following description given with reference to the appended figures, in which:
- la figure 1 illustre sous forme de synoptique les différentes étapes du procédé selon l'invention pour la détermination de la fréquence fondamentale au millième près, dans l'unité considérée ; FIG. 1 illustrates in the form of a block diagram the different steps of the method according to the invention for determining the fundamental frequency to the nearest thousandth, in the unit under consideration;
- la figure 2 illustre des résultats de simulation ; FIG. 2 illustrates simulation results;
- la figure 3 illustre sous forme de synoptique différentes sous étapes préférentiellement mises en œuvre pour la détermination de la fréquence fondamentale au dixième près, dans l'unité considérée ; FIG. 3 illustrates in the form of a different block diagram in the steps preferably implemented for determining the fundamental frequency to the nearest tenth, in the unit under consideration;
- la figure 4 illustre sous forme de synoptique différentes sous étapes préférentiellement mises en œuvre pour la détermination de la fréquence fondamentale au centième près, dans l'unité considérée ; FIG. 4 illustrates in the form of a different block diagram, in steps preferentially implemented for determining the fundamental frequency to the nearest hundredth, in the unit under consideration;
- la figure 5 illustre sous forme de synoptique différentes sous étapes préférentiellement mises en œuvre pour la détermination de la fréquence fondamentale au millième près, dans l'unité considérée ; FIG. 5 illustrates in the form of a block diagram various sub-steps preferentially implemented for the determination of the fundamental frequency to the nearest thousandth, in the unit under consideration;
- la figure 6 illustre un exemple d'échantillons d'un signal à contenu harmonique dont on veut déterminer la fréquence fondamentale avec une grande précision ; FIG. 6 illustrates an example of samples of a harmonic content signal whose fundamental frequency is to be determined with great precision;
- la figure 7 illustre schématiquement un compteur d'énergie électrique implémentant le procédé selon l'invention. FIG. 7 schematically illustrates an electric energy meter implementing the method according to the invention.
La présente invention est fondée sur le principe connu que l'on peut obtenir une représentation spectrale discrète de tout signal échantillonné, périodique ou non, à partir d'une transformée discrète de Fourier. Ainsi, le
signal peut être décomposé sous forme d'une somme de signaux purs (sinus et cosinus) pondérés par des coefficients appelés coefficients de Fourier. The present invention is based on the known principle that one can obtain a discrete spectral representation of any sampled signal, periodic or not, from a discrete Fourier transform. So, the signal can be decomposed as a sum of pure signals (sine and cosine) weighted by coefficients called Fourier coefficients.
Pour un signal périodique à une fréquence fondamentale donnée et présentant un contenu harmonique, les deux coefficients de Fourier discrets de rang h peuvent s'exprimer selon les relations suivantes : For a periodic signal at a given fundamental frequency and having a harmonic content, the two discrete Fourier coefficients of rank h can be expressed according to the following relations:
(Relations ( 1))
(Relationships (1))
dans lesquelles : in which :
- ah (F) et bh (F) représentent les coefficients de Fourier discrets de rang h à la fréquence F ; - a h (F) and b h (F) represent the discrete Fourier coefficients of rank h at the frequency F;
- T est la période du signal, inverse de la fréquence F ; - T is the period of the signal, inverse of the frequency F;
- Te est la période d'échantillonnage ; - T e is the sampling period;
- N est le nombre total d'échantillons considérés ; - N is the total number of samples considered;
- NT est le nombre total de périodes du signal sur lequel les échantillons ont été prélevés ; - NT is the total number of periods of the signal on which the samples were taken;
- i est la référence d'un échantillon courant ; i is the reference of a current sample;
- ,5( est un échantillon du signal . -, 5 (is a sample of the signal.
L'amplitude ^( 0) du signal à la fréquence fondamentale F0 peut être exprimée en fonction des coefficients de Fourier discrets de rang 1 calculés pour cette fréquence selon la relation : The amplitude ^ ( 0 ) of the signal at the fundamental frequency F 0 can be expressed as a function of the discrete Fourier coefficients of rank 1 calculated for this frequency according to the relation:
^( 0) = /ai2 ( 0) + ô1 2( 0) (Relation (2)) ^ ( 0 ) = / ai 2 ( 0 ) + ô 1 2 ( 0 ) (Relation (2))
Avec, d'après les relations (1 ) ci-dessus With, according to the relations (1) above
i N-\2T i N- \ 2T
«l( b) =—∑-^S{i)œs(2 F0iTe) " L (b) = -Σ- ^ S {i) œs (2 F 0 iT e )
(Relations (3))
Dans une mise en œuvre préférée, le procédé selon l'invention
consiste, comme cela va à présent être détaillé en référence à la figure 1, à effectuer une première estimation grossière de la fréquence fondamentale du signal, exprimée dans une unité de fréquence choisie de sorte que la partie entière de la première valeur estimée comporte au plus trois chiffres, puis de rechercher successivement une estimation plus précise, typiquement au dixième près, dans l'unité considérée, puis au centième près, dans l'unité considérée, puis au millième près, dans l'unité considérée, en recherchant à chaque fois la valeur de fréquence pour laquelle l'amplitude correspondante du signal, calculée à partir des coefficients de Fourier discrets de rang 1 selon les relations (2) et (3) ci-dessus, est maximale. (Relationships (3)) In a preferred embodiment, the process according to the invention consists, as will now be detailed with reference to FIG. 1, of making a first rough estimate of the fundamental frequency of the signal, expressed in a frequency unit chosen so that the integer part of the first estimated value comprises at most three digits, then search successively for a more precise estimate, typically to the nearest tenth, in the unit under consideration, then to the nearest hundredth, in the unit under consideration, then to the nearest thousandth, in the unit under consideration, looking for each time the frequency value for which the corresponding amplitude of the signal, calculated from the discrete Fourier coefficients of rank 1 according to the relations (2) and (3) above, is maximum.
Dans la suite, on conviendra des notations suivantes : In the following, we will agree the following notations:
- s(t) représente le signal périodique ayant un contenu harmonique, dont on souhaite déterminer la fréquence fondamentale ; - s (t) represents the periodic signal having a harmonic content, whose fundamental frequency is to be determined;
- F i est une estimation de la fréquence fondamentale au dixième io_i - F i is an estimate of the fundamental frequency at the tenth io _i
près, dans l'unité considérée ; near, in the considered unit;
- F o est une estimation de la fréquence fondamentale au centième près, dans l'unité considérée ; - F o is an estimate of the fundamental frequency to the nearest hundredth, in the considered unit;
- F est une estimation de la fréquence fondamentale au millième près, dans l'unité considérée. - F is an estimate of the fundamental frequency to the nearest thousandth, in the considered unit.
Conformément au synoptique simplifié représenté sur la figure 1, une première étape 100 du procédé selon l'invention consiste à échantillonner et à pondérer le signal s(t) à une fréquence d'échantillonnage Fe prédéfinie de manière à délivrer un nombre déterminé N d'échantillons du signal . Le but de la pondération est de borner le signal dans le temps. On utilise de préférence une fenêtre de pondération de type Hanning qui permet d'obtenir la limitation voulue avec peu d'influence sur le signal. According to the simplified block diagram shown in FIG. 1, a first step 100 of the method according to the invention consists of sampling and weighting the signal s (t) at a predefined sampling frequency F e so as to deliver a determined number N d samples of the signal. The purpose of weighting is to limit the signal in time. A Hanning weighting window is preferably used which makes it possible to obtain the desired limitation with little influence on the signal.
Une deuxième étape 200 d'estimation grossière de la fréquence fondamentale est alors réalisée, dans laquelle une première valeur estimée F de la fréquence fondamentale est calculée à partir desdits échantillons du signal. La détection des passages par zéro (« zéro crossing » en terminologie anglo saxonne) est préférentiellement utilisée pour cette deuxième étape.
Conformément à cette méthode connue, un test est effectué sur deux échantillons successifs de signal pour déterminer s'ils sont de signe opposé. Les passages à zéro sur front montant et front descendant sont comptabilisés sur l'horizon d'une seconde. Pour affiner le résultat obtenu, une interpolation linéaire est de préférence réalisée sur la dernière détection d'un passage à zéro. A l'issue de la deuxième étape 200, la première valeur estimée F de la fréquence fondamentale est délivrée. Cette valeur est cependant peu précise (précision maximum au dixième de l'unité considérée) . Il est important de noter, pour la compréhension, que la première valeur F est exprimée dans une unité de fréquence choisie de sorte que la partie entière de cette première valeur comprend au plus trois chiffres. Ainsi, pour un signal de fréquence fondamentale de l'ordre de 50 ou 60 Hertz, l'unité considérée dans la suite du procédé est le Hertz. En revanche, si l'on est en présence de signaux de l'ordre de 50 000 Hertz, l'unité considérée dans la suite sera le KiloHertz. A second step 200 of rough estimation of the fundamental frequency is then performed, in which a first estimated value F of the fundamental frequency is calculated from said samples of the signal. The detection of zero crossings ("zero crossing" in English terminology) is preferably used for this second step. According to this known method, a test is performed on two successive samples of signal to determine if they are of opposite sign. Zero crossings on rising edge and falling edge are counted over the one second horizon. To refine the result obtained, a linear interpolation is preferably performed on the last detection of a zero crossing. At the end of the second step 200, the first estimated value F of the fundamental frequency is delivered. This value is however not very precise (maximum precision to one tenth of the considered unit). It is important to note, for the sake of understanding, that the first value F is expressed in a chosen frequency unit so that the integer part of this first value comprises at most three digits. Thus, for a fundamental frequency signal of the order of 50 or 60 Hertz, the unit considered in the following process is Hertz. On the other hand, if one is in the presence of signals of the order of 50 000 Hertz, the unit considered in the continuation will be KiloHertz.
Une troisième étape 300 d'estimation au dixième près, dans l'unité considérée, de la fréquence fondamentale est alors réalisée. Lors de cette étape 300, une deuxième valeur estimée au dixième près, dans l'unité considérée, de la fréquence fondamentale est déterminée à partir de la première valeur estimée F issue de la première étape 200. Plus précisément, conformément à l'invention, la troisième étape 300 consiste à : A third step 300 estimate to the nearest tenth, in the unit considered, of the fundamental frequency is then performed. During this step 300, a second value estimated to within one-tenth, in the unit under consideration, of the fundamental frequency is determined from the first estimated value F resulting from the first step 200. More precisely, according to the invention, the third step 300 consists of:
- rechercher, parmi un premier ensemble de valeurs de fréquence possibles correspondant à des variations, par pas constant de 0, 1, de la première valeur estimée, une première valeur de fréquence pour laquelle l'amplitude correspondante du signal, calculée à partir des coefficients de Fourier discrets de rang 1, est maximale, et searching, among a first set of possible frequency values corresponding to variations, by constant step of 0, 1, of the first estimated value, a first frequency value for which the corresponding amplitude of the signal, calculated from the coefficients discrete Fourier of rank 1, is maximal, and
- faire correspondre la deuxième valeur estimée au dixième près, dans l'unité considérée à ladite première valeur de fréquence. - Match the second estimated value to the nearest tenth, in the unit considered at said first frequency value.
La troisième étape 300 revient ainsi à rechercher, parmi toutes les valeurs possibles F telles que The third step 300 thus returns to search, among all the possible values F such that
F = F x +p x O,lF = F x + p x O, l
ïo-1 ïo -1
p = 0 ; ± 1 ; ± 2 ...
la valeur pour laquelle on obtient la relation p = 0; ± 1; ± 2 ... the value for which we get the relationship
Max relation (4)
Max relation (4)
A titre d'exemple, si une fréquence grossière de 51,3 Hz est trouvée à l'issue de l'étape 200, l'unité considérée est le Hertz, et le premier ensemble de valeurs de fréquence possibles, dans lequel on recherchera une estimation au dixième de Hertz près, de la fréquence fondamentale, comprendra au plus les dix-neuf valeurs suivantes : 50,1 Hz; 50,2 Hz; 50,3 Hz; 50,4 Hz ; 50,5 Hz; 50,6 Hz ; 50,7 Hz; 50,8 Hz; 50,9 Hz ; 51,0 Hz; 51,1 Hz; 51,2 Hz; 51,3 Hz ; 51,4 Hz; 51,5 Hz ; 51,6 Hz; 51,7 Hz; 51,8 Hz et 51,9 Hz. For example, if a coarse frequency of 51.3 Hz is found at the end of step 200, the considered unit is Hertz, and the first set of possible frequency values, in which a search will be made for estimate to one-tenth of a Hertz, of the fundamental frequency, shall include not more than nineteen values: 50.1 Hz; 50.2 Hz; 50.3 Hz; 50.4 Hz; 50.5 Hz; 50.6 Hz; 50.7 Hz; 50.8 Hz; 50.9 Hz; 51.0 Hz; 51.1 Hz; 51.2 Hz; 51.3 Hz; 51.4 Hz; 51.5 Hz; 51.6 Hz; 51.7 Hz; 51.8 Hz and 51.9 Hz.
A l'issue de l'étape 300, une deuxième valeur F _i estimée au dixième près, dans l'unité considérée, de la fréquence fondamentale est délivrée. At the end of step 300, a second value F i, estimated to the nearest tenth, in the unit in question, of the fundamental frequency is delivered.
Le procédé de l'invention peut s'arrêter à ce niveau dans tous les cas où une précision supérieure n'est pas recherchée. The method of the invention can stop at this level in all cases where a higher accuracy is not sought.
Dans les autres cas, le procédé se poursuit par une quatrième étape 400 d'estimation au centième près, dans l'unité considérée, de la fréquence fondamentale. Cette quatrième étape 400 est similaire à la troisième étape 300, à ceci près que la valeur initiale de la fréquence que l'on cherche à affiner ici correspond à la deuxième valeur estimée ^10-i issue de l'étape 300. Ici, une troisième valeur estimée F 9 est déterminée à partir de la deuxième valeur F i en : In other cases, the method continues with a fourth step 400 of estimate to the nearest hundredth, in the unit considered, of the fundamental frequency. This fourth step 400 is similar to the third step 300, except that the initial value of the frequency which one seeks to refine here corresponds to the second estimated value ^ 10 -i resulting from the step 300. Here, a third estimated value F 9 is determined from the second value F i in:
io_i io _i
- recherchant, parmi un deuxième ensemble de valeurs de fréquence possibles correspondant à des variations, par pas constant de 0,01, de la deuxième valeur estimée ^10-i , une deuxième valeur de fréquence pour laquelle l'amplitude correspondante du signal, calculée à partir des coefficients de Fourier discrets de rang 1, est maximale, et seeking, among a second set of possible frequency values corresponding to variations, by a constant pitch of 0.01, of the second estimated value 10 10 -i, a second frequency value for which the corresponding amplitude of the signal, calculated from the discrete Fourier coefficients of rank 1, is maximal, and
- en faisant correspondre la troisième valeur estimée F 2 à la deuxième valeur de fréquence issue de la recherche. by matching the third estimated value F 2 to the second frequency value resulting from the search.
La troisième étape 400 revient ainsi à rechercher, parmi toutes les
valeurs possibles F telles que The third step 400 thus returns to search, among all the possible values F such that
-2 -2
F F + gx 0,01 q = 0 ; ± 1 ; ± 2 ... F F + gx 0.01 q = 0; ± 1; ± 2 ...
la valeur pour laquelle on obtient la relation the value for which we get the relationship
Max relation (5).Max relation (5).
A titre d'exemple, si une fréquence estimée au dixième de Hz près de For example, if a frequency estimated at one-tenth of a Hz near
51,3 Hz est trouvée à l'issue de l'étape 300, le deuxième ensemble de valeurs de fréquence possibles, dans lequel on recherchera une estimation au centième de Hertz près de la fréquence fondamentale, comprendra au plus les dix-neuf valeurs suivantes : 51,21 Hz; 51,22 Hz; 51,23 Hz; 51,24 Hz ; 51,25 Hz; 51,26 Hz ; 51,27 Hz; 51,28 Hz; 51,29 Hz ; 51,30 Hz; 51,31 Hz; 51,32 Hz; 51,33 Hz ; 51,34 Hz; 51,35 Hz ; 51,36 Hz; 51,37 Hz; 51,38 Hz et 51,39 Hz. 51.3 Hz is found at the end of step 300, the second set of possible frequency values, in which an estimate hundredths of Hertz is sought near the fundamental frequency, will include at most the nineteen following values : 51.21 Hz; 51.22 Hz; 51.23 Hz; 51.24 Hz; 51.25 Hz; 51.26 Hz; 51.27 Hz; 51.28 Hz; 51.29 Hz; 51.30 Hz; 51.31 Hz; 51.32 Hz; 51.33 Hz; 51.34 Hz; 51.35 Hz; 51.36 Hz; 51.37 Hz; 51.38 Hz and 51.39 Hz.
A l'issue de l'étape 400, une troisième valeur F 9 estimée au centième près, dans l'unité considérée, de la fréquence fondamentale est ainsi délivrée. At the end of step 400, a third value F 9 estimated to the nearest hundredth, in the considered unit, of the fundamental frequency is thus delivered.
Si la précision obtenue n'est pas satisfaisante pour l'application envisagée, le procédé peut se poursuivre par une cinquième étape 500 d'estimation au millième près, dans l'unité considérée, de la fréquence fondamentale. Ici encore, la cinquième étape est très similaire aux deux étapes 300 et 400 décrites précédemment, à ceci près que la valeur initiale de la fréquence que l'on cherche à affiner correspond à la troisième valeur estimée F 2 issue de l'étape 400. Par suite, une quatrième valeur estimée If the accuracy obtained is not satisfactory for the intended application, the method can be continued by a fifth step 500 estimate to the nearest thousandth, in the unit considered, of the fundamental frequency. Here again, the fifth step is very similar to the two steps 300 and 400 described above, except that the initial value of the frequency that one seeks to refine corresponds to the third estimated value F 2 resulting from step 400. As a result, a fourth estimated value
F est déterminée à partir de la troisième valeur estimée F 9 en : F is determined from the third estimated value F 9 in:
10 J 10 2 10 D 10 2
- recherchant, parmi un troisième ensemble de valeurs de fréquence possibles correspondant à des variations, par pas constant de 0,001, de la troisième valeur estimée _2 , une troisième valeur de fréquence pour
laquelle l'amplitude correspondante du signal, calculée à partir des coefficients de Fourier discrets de rang 1, est maximale, et searching, among a third set of possible frequency values corresponding to variations, by a constant pitch of 0.001, of the third estimated value _ 2 , a third frequency value for which the corresponding amplitude of the signal, calculated from the discrete Fourier coefficients of rank 1, is maximum, and
- faisant correspondre la quatrième valeur estimée ^10-3 à la troisième valeur de fréquence issue de la recherche. - matching the fourth estimated value ^ 10 -3 to the third frequency value from the search.
La quatrième étape 500 revient ainsi à rechercher, parmi toutes les valeurs possibles F telles que The fourth step 500 thus returns to search, among all the possible values F such that
£,10 J £ 10 J
F F _2 + £x 0,001 FF _ 2 + £ x 0.001
10 l 10 l
£,10~J £, 10 ~ J
k = 0 ; ± 1 ; ± 2 ... k = 0; ± 1; ± 2 ...
la valeur pour laquelle on obtient la relation the value for which we get the relationship
Max relation (6).Max relation (6).
A titre d'exemple, si une fréquence estimée au centième de Hz près de For example, if a frequency estimated at one hundredth of a Hz near
51,32 Hz est trouvée à l'issue de l'étape 400, le troisième ensemble de valeurs de fréquence possibles, dans lequel on recherchera une estimation au millième de Hertz près de la fréquence fondamentale, comprendra au plus les dix-neuf valeurs suivantes : 51,311 Hz ; 51,312 Hz ; 51,313 Hz ; 51,314 Hz ; 51,315 Hz ; 51,316 Hz ; 51,317 Hz ; 51,318 Hz ; 51,319 Hz ; 51,320 Hz; 51,321 Hz; 51,322 Hz; 51,323 Hz ; 51,324 Hz; 51,325 Hz ; 51,326 Hz; 51,327 Hz; 51,328 Hz et 51,329 Hz. 51.32 Hz is found at the end of step 400, the third set of possible frequency values, in which an estimate of one thousandth of Hertz is sought near the fundamental frequency, will comprise at most the nineteen following values : 51.311 Hz; 51.312 Hz; 51.313 Hz; 51.314 Hz; 51.315 Hz; 51.316 Hz; 51.317 Hz; 51.318 Hz; 51.319 Hz; 51.320 Hz; 51.321 Hz; 51.322 Hz; 51.323 Hz; 51.324 Hz; 51.325 Hz; 51.326 Hz; 51.327 Hz; 51.328 Hz and 51.329 Hz.
A l'issue de l'étape 500, une quatrième valeur ^10-3 estimée au millième près, dans l'unité considérée, de la fréquence fondamentale est ainsi délivrée. At the end of step 500, a fourth value ^ 10 -3 estimated to the nearest thousandth, in the unit in question, of the fundamental frequency is thus delivered.
On pourrait penser que le procédé pourrait se poursuivre ainsi. Néanmoins, les simulations effectuées par la Demanderesse ont permis de montrer que le procédé d'estimation présente une limite, typiquement au millième près, dans l'unité considérée. En effet, l'estimation de la fréquence fondamentale dépend des différents paramètres suivants : One might think that the process could continue like this. Nevertheless, the simulations carried out by the Applicant have made it possible to show that the estimation method has a limit, typically to the nearest thousandth, in the unit under consideration. Indeed, the estimation of the fundamental frequency depends on the following different parameters:
- la période T du signal dont on recherche la fréquence
fondamentale ; the period T of the signal whose frequency is sought fundamental ;
- la variation x sur T correspondant à la précision - the variation x on T corresponding to the precision
- la période d'échantillonnage Te - the sampling period T e
- le nombre total NT de périodes du signal sur lequel les échantillons ont été prélevés - the total number N T of periods of the signal on which the samples were taken
- le nombre total N d'échantillons considérés, donné par la relation
- the total number N of samples considered, given by the relation
Des simulations ont été réalisées pour un signal sinusoïdal de fréquence fondamentale égale à 50Hz et pour une période d'échantillonnage Te égale à 256 psecondes, pour mesurer l'influence du nombre total N d'échantillons considérés sur la précision de l'estimation . Simulations were carried out for a sinusoidal signal of fundamental frequency equal to 50 Hz and for a sampling period T e equal to 256 psec, to measure the influence of the total number N of samples considered on the accuracy of the estimate.
Plus précisément, la figure 2 représente la variation de la précision sur la fréquence fondamentale du signal périodique ( proche de 50Hz) identifiée à la fréquence donnant le maximum d'amplitude de la Transformée de Fourier Discrète au rang fondamental, en fonction du nombre d'échantillons constituant la fenêtre d'acquisition dudit signal échantillonné à une période d'échantillonnage Te , la fenêtre de pondération appliquée étant une fenêtre deMore precisely, FIG. 2 represents the variation of the accuracy on the fundamental frequency of the periodic signal (close to 50 Hz) identified at the frequency giving the maximum amplitude of the Discrete Fourier Transform to the fundamental rank, as a function of the number of samples constituting the acquisition window of said sampled signal at a sampling period T e , the weighting window applied being a window of
Hanning . Hanning.
On note la relation : We note the relationship:
Tf = NT x T = Nx Te T f = N T x T = N x T e
avec : with:
Tf , la durée en secondes de la fenêtre d'acquisition Tf, the duration in seconds of the acquisition window
NT , le nombre de périodes du fondamental réel de période T et N , le nombre d'échantillons acquis à la période d'échantillonnage Te N T , the number of periods of the real fundamental period T and N, the number of samples acquired at the sampling period T e
Ainsi, si N est égal à 1250, on obtient : Thus, if N is equal to 1250, we obtain:
Tf = 1250 x 256 ps = 0,32 s qui correspond à un nombre NT de 16 périodes de 50 Hz, Tf = 1250 x 256 ps = 0.32 s which corresponds to a number N T of 16 periods of 50 Hz,
et pour N égal à 5000, on obtient :
= > T = 5000 x 256 ps = 1 ,28 s, soit un nombre NT de 64 périodes de 50 Hz, and for N equal to 5000, we obtain: => T = 5000 x 256 ps = 1, 28 s, which is a number N T of 64 periods of 50 Hz,
Sur le graphique de la figure 2, l'axe des ordonnées est représenté en échelle logarithmique de 1 Hertz à 10"9 Hertz, et l'axe des abscisse est linéaire, de 0 à 6000 échantillons. In the graph of Figure 2, the ordinate axis is shown on a logarithmic scale from 1 Hz at 10 "9 Hertz, and the axis of abscissa is linear, from 0-6000 samples.
Il apparaît clairement sur la figure 2 que, pour obtenir une précision plus fine que 0,1 Hertz, il faut choisir un nombre d'échantillons supérieur à un certain seuil noté NQ I , correspondant ici à 150 échantillons. Pour obtenir une précision plus fine que 0,01 Hertz il faut choisir un nombre d'échantillons supérieur à un second seuil, lui-même supérieur au premier et noté No,oi ' correspondant ici à 400 échantillons. Enfin, pour obtenir une précision encore plus fine que 0,001 Hertz, il faut choisir un nombre d'échantillons supérieur à un troisième seuil supérieur aux deux autres et noté No,ooi ' correspondant ici à 600 échantillons. It appears clearly in FIG. 2 that, to obtain a precision that is finer than 0.1 Hertz, it is necessary to choose a number of samples greater than a certain threshold, denoted NQ I, here corresponding to 150 samples. To obtain an accuracy that is finer than 0.01 Hertz, it is necessary to choose a number of samples greater than a second threshold, itself greater than the first, and denoted No, where here 400 samples correspond. Finally, to obtain a precision even finer than 0.001 Hertz, it is necessary to choose a number of samples greater than a third threshold higher than the other two and noted No, where here corresponds to 600 samples.
II résulte notamment de ce graphique qu'une précision au millième d'Hertz près pour un signal de fréquence fondamental voisine de 50 Hertz peut être obtenue en choisissant un nombre d'échantillons N égal à au moins 600. It follows in particular from this graph that a precision to the nearest thousandth of a Hertz for a fundamental frequency signal close to 50 Hertz can be obtained by choosing a number of samples N equal to at least 600.
Bien entendu, le nombre d'échantillons minimum doit être recalculé pour chaque fréquence fondamentale à estimer et pour chaque fréquence d'échantillonnage prévue, en fonction de la précision souhaitée. Of course, the minimum number of samples must be recalculated for each fundamental frequency to be estimated and for each expected sampling frequency, depending on the desired accuracy.
Dans les étapes du procédé décrites en référence à la figure 1, chacune des étapes 300, 400 et 500 comportent une sous étape lors de laquelle on recherche une valeur plus précise parmi respectivement un premier, deuxième et troisième ensemble d'au plus dix-neuf valeurs possibles de fréquence. In the steps of the method described with reference to FIG. 1, each of the steps 300, 400 and 500 comprises a sub-step during which a more precise value is sought from among a first, second and third set of at most nineteen, respectively. possible values of frequency.
Une mise en œuvre préférée de la sous étape de recherche au sein des étapes 300, 400 et 500 permettant de réduire le nombre de calculs effectués va à présent être décrite en référence aux figures 3 à 5 : A preferred implementation of the sub-step of searching within steps 300, 400 and 500 for reducing the number of calculations performed will now be described with reference to FIGS. 3 to 5:
En référence tout d'abord à la figure 3, la recherche de la première valeur de fréquence ^10-i estimée au dixième près, dans l'unité considérée, de la fréquence fondamentale comporte une première sous-étape
d'initialisation 301 dans laquelle la valeur courante F0 de la fréquence fondamentale est initialisée à la première valeur estimée F , issue de l'étape 200, et une valeur de fréquence possible du premier ensemble, notée dans la suite ^10-i est initialisée à la première valeur estimée F incrémentée de 0, 1. Referring firstly to Figure 3, the search for the first frequency value ^ 10 -i estimated to the nearest tenth, in the unit considered, of the fundamental frequency comprises a first substep boot 301 in which the current value F 0 of the fundamental frequency is set to the first estimated value F, after step 200, and a frequency value as possible to the first set, denoted in the remainder 10 ^ -i is initialized to the first estimated value F incremented by 0, 1.
Pour ces deux valeurs 0 et ^10-i , on calcule alors (sous-étape 302) les deux amplitudes A(F0) et A{F γ ) obtenues à partir des coefficients de For these two values 0 and ^ 10 -i, we then calculate (sub-step 302) the two amplitudes A (F 0 ) and A {F γ) obtained from the coefficients of
Fourier discrets de rang 1, soit, conformément aux relations (2) et (3) précitée Discrete Fourier of rank 1, that is, in accordance with the above-mentioned relations (2) and (3)
Ces deux amplitudes sont ensuite comparées (sous étape 303) . These two amplitudes are then compared (in step 303).
Dans le cas où l'amplitude A{F γ ) est supérieure à l'amplitude A(F0) , cela signifie, conformément à l'invention, que la valeur F i constitue uneIn the case where the amplitude A {F γ) is greater than the amplitude A (F 0 ), this means, according to the invention, that the value F i constitutes a
3 io_i 3 io _i
meilleure estimation de la fréquence fondamentale que la valeur 0 . On remplace donc (sous étape 304) la valeur courante 0 par la valeur de fréquence possible ^ ι , puis on recherche s'il existe, dans le premier ensemble, une autre valeur encore meilleure en incrémentant la valeur de fréquence possible F ι dudit premier ensemble de 0, 1. Ceci peut s'exprimer par les relations : best estimate of the fundamental frequency than the value 0 . So we replace (in step 304) the current value 0 by the possible frequency value ^ ι, then we search if there exists, in the first set, another even better value by incrementing the possible frequency value F ι of said first set of 0, 1. This can be expressed by the relationships:
10 1 10 1
F = F + 0,1 F = F + 0.1
10"1 10"1 10 "1 10 " 1
On calcule alors (sous-étape 305), pour ces deux nouvelles valeurs 0 et^10-i les deux amplitudes A(F0) et A(F^_ ) obtenues à partir des coefficients de Fourier discrets de rang 1, et on procède à une nouvelle comparaison (sous étape 306) des amplitudes obtenues. Les sous étapes 304
et 305 sont réitérées tant que l'amplitude A(F ι ) est supérieure à l'amplitude A(F0) calculée. La recherche s'arrête dès que le résultat de la comparaison (sous étape 306) est négatif. Dans ce cas, la valeur recherchée correspond à la dernière valeur courante à l'issue des sous étapes 304 à 306, pour laquelle la comparaison était positive, c'est-à-dire la valeur avant incrémentation. Cela revient à décrémenter à nouveau la valeur courante de 0,1 (sous étape 307). Then (sub-step 305), for these two new values 0 and ^ 10 -i are computed the two amplitudes A (F 0 ) and A (F ^ _) obtained from the discrete Fourier coefficients of rank 1, and proceeds to a new comparison (in step 306) of the amplitudes obtained. Sub-steps 304 and 305 are reiterated as long as the amplitude A (F ι) is greater than the calculated amplitude A (F 0 ). The search stops as soon as the result of the comparison (under step 306) is negative. In this case, the value sought corresponds to the last current value at the end of substeps 304 to 306, for which the comparison was positive, that is to say the value before incrementation. This amounts to decrementing the current value again by 0.1 (in step 307).
Si au contraire, à l'issue de la sous étape de comparaison 303, l'amplitude A{F γ ) est inférieure à l'amplitude A(F0) , cela signifie, conformément à l'invention, que la valeur F0 constitue une meilleure estimation de la fréquence fondamentale que la valeuri7 i . Il est inutile d'aller rechercher dans ce cas les valeurs du premier ensemble supérieures à ^10-ι■ On conserve donc une valeur courante égale à 0 , mais on remplace If, on the other hand, at the end of the comparison sub-step 303, the amplitude A {F γ) is smaller than the amplitude A (F 0 ), this means, according to the invention, that the value F 0 is a better estimate of the fundamental frequency that valeuri 7 i. It is useless to go in this case to find the values of the first set greater than ^ 10 -ι ■ So we keep a current value equal to 0 , but we replace
(sous étape 308) la valeur de fréquence possible ^10-i par la valeur courante F0 décrémentée de 0,1. Ceci peut s'exprimer par les relations : (In step 308) the possible frequency value 10 10 -i by the current value F 0 decremented by 0.1. This can be expressed by relationships:
Fo = ·^ο F o = · ^ ο
F = F0 - o,\ F = F 0 - o, \
10 1 10 1
On calcule alors (sous-étape 309), pour ces deux nouvelles valeurs 0 et F ! , les deux amplitudes A(F0) et A(F^_ ) obtenues à partir des coefficients de Fourier discrets de rang 1, et on procède à une nouvelle comparaison (sous étape 310) des amplitudes obtenues. We then calculate (sub-step 309), for these two new values 0 and F! , the two amplitudes A (F 0 ) and A (F ^ _) obtained from the discrete Fourier coefficients of rank 1, and a new comparison is made (in step 310) of the amplitudes obtained.
Dans le cas où l'amplitude A{F γ ) est supérieure à l'amplitude A(F0) , cela signifie, conformément à l'invention, que la valeur F i constitue uneIn the case where the amplitude A {F γ) is greater than the amplitude A (F 0 ), this means, according to the invention, that the value F i constitutes a
3 io_i 3 io _i
meilleure estimation de la fréquence fondamentale que la valeur 0 . On remplace donc (sous étape 311) la valeur courante 0 par la valeur de fréquence possible^ ι , puis on recherche s'il existe, dans le premier ensemble, une autre valeur encore meilleure en décrémentant à nouveau la
valeur de fréquence possible F ι dudit premier ensemble de 0,1 - Ceci peut s'exprimer par les relations : best estimate of the fundamental frequency than the value 0 . So we replace (under step 311) the current value 0 by the possible frequency value ^ ι, then we search if there exists, in the first set, another even better value by decrementing again the possible frequency value F ι of said first set of 0.1 - This can be expressed by the relationships:
F0 = F . F 0 = F.
10 1 10 1
F = F - 0,1 F = F - 0.1
10"1 10"1 10 "1 10 " 1
Les sous étapes 309 et 310 sont réitérées tant que l'amplitude A(F^_i ) est supérieure à l'amplitude A(F0) calculées. La recherche s'arrête dès que le résultat de la comparaison (sous étape 310) est négatif. Dans ce cas, la valeur recherchée correspond à la dernière valeur courante à l'issue des sous étapes 311, 309 et 310, pour laquelle la comparaison était positive, c'est-à-dire la valeur avant décrémentation. Cela revient à incrémenter à nouveau la valeur courante de 0,1 (sous étape 312). Sub-steps 309 and 310 are reiterated as long as the amplitude A (F ^ _i) is greater than the calculated amplitude A (F 0 ). The search stops as soon as the result of the comparison (in step 310) is negative. In this case, the value sought corresponds to the last current value at the end of the sub-steps 311, 309 and 310, for which the comparison was positive, that is to say the value before decrementation. This amounts to incrementing the current value again by 0.1 (in step 312).
Les figures 4 et 5 illustrent une mise en œuvre préférée de l'étape de recherche d'une valeur de fréquence, respectivement au centième près, et au millième près, dans l'unité considérée. La mise en œuvre est similaire à celle présentée en référence à la figure 3, et ne sera ainsi pas décrite dans le détail. Figures 4 and 5 illustrate a preferred implementation of the step of finding a frequency value, respectively to the nearest hundredth, and to the nearest thousandth, in the unit in question. The implementation is similar to that presented with reference to Figure 3, and thus will not be described in detail.
En effet, les étapes 401 à 407 d'une part, et 408 à 412 d'autre part, de la figure 4 sont en tous points identiques aux étapes 301 à 307, et 308 à 312 décrites ci-dessus, à ceci près que les valeurs ^10-i et les Pas d'incrémentation/décrémentation de 0,1 ont été remplacés par des valeurs F 2 et des pas d'incrémentation/décrémentation de 0,01. De même, les étapes 501 à 507 d'une part, et 508 à 512 d'autre part, de la figure 5 sont en tous points identiques aux étapes 301 à 307, et 308 à 312 décrites ci-dessus, à ceci près que les valeurs ^10-i et les pas d'incrémentation/décrémentation de 0,1 ont été remplacés par des valeurs Q_3 et des pas d'incrémentation/décrémentation de 0,001. Un exemple d'application du procédé selon l'invention, et plus spécifiquement des étapes selon les figures 1, 3, 4 et 5 va à présent être décrit en référence aux figures 6 et 7. Sur la figure 6, on a représenté partiellement des échantillons obtenus à l'issue de l'étape 100 (figure 1) à
partir d'un signal provenant d'une source de tension . Indeed, the steps 401 to 407 on the one hand, and 408 to 412 on the other hand, of FIG. 4 are in all respects identical to the steps 301 to 307, and 308 to 312 described above, except that values ^ -i 10 and P as increment / decrement of 0.1 were replaced by values F 2 and no increment / decrement of 0.01. Likewise, the steps 501 to 507 on the one hand, and 508 to 512 on the other hand, of FIG. 5 are in all respects identical to the steps 301 to 307 and 308 to 312 described above, except that values ^ -i 10 and no increment / decrement of 0.1 were replaced by Q_ values 3 and no increment / decrement of 0.001. An example of application of the method according to the invention, and more specifically the steps according to FIGS. 1, 3, 4 and 5, will now be described with reference to FIGS. 6 and 7. FIG. samples obtained at the end of step 100 (Figure 1) at from a signal from a voltage source.
Le signal s(t) considéré est par exemple un signal de tension d'une phase d'un réseau de distribution électrique reçu entre une phase et un neutre d'un compteur 1 représenté schématiquement sur la figure 7. Il pourrait s'agir également de n'importe quel courant de phase et de neutre du réseau de distribution d'énergie électrique. Le compteur comporte des moyens de conversions analogique-numérique 10 et des moyens logiciels 11 permettant, sous le contrôle d'un microcontrôleur 12, de mettre en œuvre le procédé de détermination de la fréquence fondamentale du signal s(t) au millième d'Hertz près. Dans cet exemple, 5000 échantillons à une période d'échantillonnage de 256 psecondes sont disponibles en sortie des moyens de conversion 10 à l'issue de l'étape 100. L'étape d'estimation grossière 200, faite par la méthode de la détection des passages à zéro, permet d'obtenir une première valeur estimée de 51,3 Hz. L'unité considérée est donc ici le Hertz. The signal s (t) considered is for example a voltage signal of a phase of an electrical distribution network received between a phase and a neutral of a counter 1 shown schematically in FIG. 7. It could also be any phase and neutral current of the electrical power distribution network. The counter comprises analogue-digital conversion means 10 and software means 11 making it possible, under the control of a microcontroller 12, to implement the method for determining the fundamental frequency of the signal s (t) to the thousandth of a Hertz near. In this example, 5000 samples at a sampling period of 256 ps are available at the output of the converting means at the end of step 100. The coarse estimation step 200, made by the detection method zero crossings, makes it possible to obtain a first estimated value of 51.3 Hz. The unit considered here is therefore the Hertz.
Le tableau ci-dessous indique les résultats intermédiaires trouvés, conformément aux étapes des figures 3, 4 et 5 : The table below shows the intermediate results found, according to the steps in Figures 3, 4 and 5:
La comparaison des amplitudes (étape 303) donnant ici un résultat négatif, on continue selon l'étape 308. The comparison of the amplitudes (step 303) giving a negative result here, is continued according to step 308.
La valeur de 51,3 Hz constitue donc la meilleure estimation au dixième d'Hertz près que l'on peut obtenir. La recherche d'une estimation au centième d'Hertz près est alors effectuée :
Etape 401 Etape 402 The value of 51.3 Hz is therefore the best estimate to the nearest tenth of Hertz that can be obtained. The search for an estimate to the hundredth of Hertz near is then carried out: Step 401 Step 402
F0 =51,30 Hz A( i¾) = 1477,54 LSB F 0 = 51.30 Hz A (i) = 1477.54 LSB
F =51,31 Hz A( ) = 1477,97 LSB F = 51.31 Hz A () = 1477.97 LSB
10 1 10 1
La comparaison des amplitudes (étape 403) donnant ici un résultat positif, on continue selon l'étape 404 : The comparison of the amplitudes (step 403) gives here a positive result, one continues according to step 404:
La comparaison des amplitudes (étape 406) donnant ici un résultat positif, on réitère les étapes 404 à 405 : Comparing the amplitudes (step 406) giving a positive result here, steps 404 to 405 are repeated:
La valeur de 51,32 Hz constitue donc la meilleure estimation au centième d'Hertz près que l'on peut obtenir. La recherche d'une estimation au millième d'Hertz près est alors effectuée : The value of 51.32 Hz is therefore the best estimate to the hundredth of Hertz that can be obtained. The search for an estimate to the nearest thousandth of a Hertz is then carried out:
La comparaison des amplitudes (étape 403) donnant ici un résultat négatif, on continue selon l'étape 408.
Etape 408 Etape 409 The comparison of the amplitudes (step 403) giving here a negative result, is continued according to step 408. Step 408 Step 409
F0 =51,320 Hz A( 0) = 1478,0926 LSB F 0 = 51.320 Hz A ( 0 ) = 1478.0926 LSB
F , =51,319 Hz A( F, = 51.319 Hz A (
10 5 = 1478,0948 LSB 10 5 = 1478.0948 LSB
La comparaison des amplitudes (étape 410) donnant ici un résultat positif, on continue selon l'étape 411 : The comparison of the amplitudes (step 410) giving here a positive result, one continues according to the step 411:
La valeur de 51,319 Hz constitue donc la meilleure estimation au millième d'Hertz près que l'on peut obtenir. The value of 51.319 Hz is therefore the best estimate to the nearest thousandth of a Hertz that can be obtained.
Ainsi, grâce à l'invention, il est possible de déterminer une fréquence fondamentale voisine de 50 Hz ou de 60 Hz avec une précision relative de 20 ppm, ou plus généralement, toute fréquence fondamentale avec une précision du millième d'unité.
Thus, thanks to the invention, it is possible to determine a fundamental frequency close to 50 Hz or 60 Hz with a relative accuracy of 20 ppm, or more generally, any fundamental frequency with a precision of one thousandth of a unit.