Calcul de fonction de transfert en z
Calculez rapidement une fonction de transfert discrète de type IIR, vérifiez la stabilité via les pôles, obtenez le gain statique, visualisez la réponse fréquentielle et interprétez le comportement de votre système en domaine z.
Calculateur interactif H(z)
Le calculateur utilise la forme standard : H(z) = (b0 + b1 z^-1 + b2 z^-2) / (1 + a1 z^-1 + a2 z^-2).
Cette valeur détermine l’échelle de la réponse fréquentielle jusqu’à fs/2.
Coefficients du numérateur
Coefficients du dénominateur
Comprendre le calcul de fonction de transfert en z
Le calcul de fonction de transfert en z est une étape fondamentale en traitement numérique du signal, en automatique discrète, en filtrage digital et dans la modélisation des systèmes échantillonnés. Dès qu’un système est observé ou commandé à intervalles réguliers, l’analyse dans le domaine z devient plus pertinente que l’analyse purement continue dans le domaine de Laplace. Une fonction de transfert en z, notée H(z), permet de relier l’entrée discrète x[n] à la sortie discrète y[n] à travers une relation algébrique compacte. Elle synthétise le comportement dynamique du système, sa stabilité, son gain, sa sélectivité fréquentielle et sa sensibilité à l’échantillonnage.
Dans la pratique, le calcul de H(z) intervient dans de nombreux cas : conversion d’un correcteur analogique vers un correcteur numérique, conception de filtres passe-bas ou passe-bande, estimation de la réponse d’un algorithme récursif, simulation de modèles physiques discrétisés, ou encore vérification de la stabilité d’un système commandé en boucle fermée. Le grand intérêt du domaine z est qu’il transforme des équations aux différences en expressions rationnelles. Cela facilite le calcul des pôles, des zéros et de la réponse fréquentielle.
Définition générale
Une fonction de transfert discrète est souvent écrite sous la forme :
H(z) = Y(z) / X(z) = (b0 + b1 z^-1 + b2 z^-2 + … + bm z^-m) / (1 + a1 z^-1 + a2 z^-2 + … + an z^-n)
Cette forme correspond directement à l’équation aux différences :
y[n] = -a1 y[n-1] – a2 y[n-2] – … + b0 x[n] + b1 x[n-1] + b2 x[n-2] + …
Autrement dit, le numérateur décrit l’influence des échantillons d’entrée présents et passés, tandis que le dénominateur encode la mémoire interne du système via la rétroaction. Les systèmes sans rétroaction sont appelés FIR, alors que les systèmes avec rétroaction sont souvent IIR.
Pourquoi le domaine z est central en systèmes discrets
Le domaine z joue pour les systèmes numériques un rôle comparable à celui de Laplace pour les systèmes continus. Il permet d’étudier la réponse à des séquences, la convergence, les modes internes et les caractéristiques fréquentielles. Il devient particulièrement utile dès qu’un système est échantillonné, car les outils continus ne traduisent pas toujours correctement les effets de l’échantillonnage, du maintien d’ordre zéro ou des délais numériques.
- Analyse de stabilité : un système discret est stable si tous ses pôles sont strictement à l’intérieur du cercle unité.
- Réponse fréquentielle : on l’obtient en évaluant H(z) sur le cercle unité, c’est-à-dire pour z = e^(jω).
- Implémentation : la structure algébrique de H(z) se traduit directement en code embarqué, DSP ou microcontrôleur.
- Conception numérique : filtres audio, capteurs intelligents, commande de moteurs, observateurs, asservissements et télécommunications utilisent tous cette représentation.
Méthode de calcul étape par étape
1. Écrire l’équation aux différences
La première étape consiste à poser correctement la relation entre entrée et sortie. Par exemple, si l’on a :
y[n] – 1.143 y[n-1] + 0.4128 y[n-2] = 0.0675 x[n] + 0.1349 x[n-1] + 0.0675 x[n-2]
alors la fonction de transfert s’écrit :
H(z) = (0.0675 + 0.1349 z^-1 + 0.0675 z^-2) / (1 – 1.143 z^-1 + 0.4128 z^-2)
Le calculateur ci-dessus reprend exactement cette convention d’écriture. Cela vous permet de passer très vite d’un jeu de coefficients à une représentation exploitable.
2. Identifier les zéros
Les zéros sont les racines du numérateur. Ils correspondent aux fréquences ou modes pour lesquels la sortie est atténuée ou annulée. Dans le cadre d’un filtre numérique, leur position dans le plan z détermine souvent les bandes rejetées. Des zéros proches de z = -1 influencent fortement la zone haute fréquence, tandis que des zéros proches de z = 1 affectent le comportement à basse fréquence.
3. Identifier les pôles
Les pôles sont les racines du dénominateur. Ce sont eux qui gouvernent la dynamique interne, la résonance, le temps de réponse et surtout la stabilité. Dans le domaine z, le critère est simple mais capital : si tous les pôles ont un module inférieur à 1, le système est BIBO stable. Si un pôle atteint le cercle unité, le système devient marginal. Si un pôle se trouve à l’extérieur, la réponse peut diverger.
4. Calculer le gain statique
Le gain en continu, souvent appelé gain DC, se calcule en évaluant H(z) pour z = 1. Cette opération revient à sommer les coefficients du numérateur puis à les diviser par la somme adaptée des coefficients du dénominateur. Le gain DC est très utile en automatique numérique, car il renseigne sur la fidélité de la transmission des composantes lentes ou constantes.
5. Tracer la réponse fréquentielle
Pour étudier le comportement spectral, on remplace z par e^(jω), avec ω variant entre 0 et π rad/échantillon. On obtient alors la magnitude et la phase. Dans un contexte d’ingénierie, cette courbe permet de distinguer un filtre passe-bas, passe-haut, coupe-bande ou résonant. Le graphique fourni par le calculateur affiche le module en dB de 0 à fs/2, ce qui est généralement la forme la plus exploitable sur le plan pratique.
Interprétation physique des coefficients
Les coefficients b0, b1, b2, etc., représentent la contribution directe de l’entrée courante et des entrées passées. Plus leur répartition est symétrique, plus on rencontre souvent un comportement de filtrage régulier, notamment dans certaines structures de filtres passe-bas. Les coefficients a1, a2, etc., quant à eux, contrôlent la mémoire et la réinjection interne. Une légère variation de ces coefficients peut déplacer significativement les pôles et donc changer la stabilité ou la sélectivité.
En électronique numérique, en commande industrielle et en audio embarqué, il est courant de quantifier les coefficients avec une précision limitée. C’est pourquoi il est important de recalculer les pôles après arrondi. Un filtre stable en virgule flottante peut devenir moins robuste, voire instable, après implantation en précision fixe si les pôles se rapprochent trop du cercle unité.
| Fréquence d’échantillonnage | Fréquence de Nyquist | Usage réel fréquent | Impact sur le calcul en z |
|---|---|---|---|
| 8 000 Hz | 4 000 Hz | Téléphonie numérique classique | Résolution limitée dans l’aigu, suffisante pour la voix. |
| 44 100 Hz | 22 050 Hz | Audio CD | Très courant pour filtres audio et égalisation. |
| 48 000 Hz | 24 000 Hz | Audio pro, vidéo, interfaces multimédia | Standard robuste pour traitements temps réel. |
| 1 000 Hz | 500 Hz | Commande numérique et capteurs industriels | Fréquent pour asservissements de vitesse ou température rapide. |
| 100 Hz | 50 Hz | Supervision lente, instrumentation | Adéquat pour systèmes à dynamique lente. |
Stabilité et localisation des pôles
Le point clé du calcul de fonction de transfert en z reste souvent la stabilité. En continu, on observe les pôles dans le demi-plan gauche. En discret, le critère change : le cercle unité devient la frontière. Si un pôle est proche de 1, la réponse est lente et peut présenter une longue traîne. S’il est proche de -1, des oscillations alternées apparaissent. Deux pôles complexes conjugués proches du cercle unité créent des résonances marquées et un amortissement faible.
- Calculez les racines du dénominateur.
- Évaluez leur module.
- Vérifiez que chaque module est strictement inférieur à 1.
- Contrôlez ensuite la marge pratique de stabilité, surtout en implantation fixe.
Une bonne pratique industrielle consiste à éviter des pôles trop proches du cercle unité lorsque la quantification, les retards logiciels ou les saturations risquent de dégrader le comportement. Dans les applications de contrôle embarqué, on préfère souvent une stabilité non seulement théorique mais aussi robuste.
Différence entre fonction de transfert continue et discrète
La fonction de transfert continue H(s) et la fonction de transfert discrète H(z) décrivent des systèmes dynamiques, mais elles ne vivent pas dans le même espace mathématique. Lors d’une discrétisation, un système continu est converti en système échantillonné. La relation entre s et z dépend alors de la méthode utilisée : maintien d’ordre zéro, transformation bilinéaire, invariance impulsionnelle ou autres techniques plus spécialisées.
| Critère | Domaine continu | Domaine z | Conséquence pratique |
|---|---|---|---|
| Variable complexe | s = σ + jω | z = re^(jω) | Le plan d’analyse et les critères de stabilité changent. |
| Stabilité | Pôles dans le demi-plan gauche | Pôles dans le cercle unité | Une discrétisation mal choisie peut altérer les marges. |
| Réponse fréquentielle | Évaluation sur jω | Évaluation sur le cercle unité | La fréquence est bornée entre 0 et fs/2. |
| Mise en oeuvre | Équations différentielles | Équations aux différences | Le code temps réel découle naturellement de H(z). |
Influence de la précision numérique
Le calcul de fonction de transfert en z ne s’arrête pas à la formule théorique. En pratique, la précision machine influence directement les performances. Plus le nombre de bits est réduit, plus les coefficients et les états internes sont quantifiés. Cela peut déplacer les pôles, ajouter du bruit ou provoquer de la saturation. En audio numérique, en contrôle moteur ou en instrumentation embarquée, ces effets sont bien réels.
La dynamique théorique d’un codage uniforme est souvent estimée par la relation 6,02N + 1,76 dB, où N est le nombre de bits. Cela donne des repères concrets pour l’implantation des filtres :
| Résolution | Dynamique théorique approximative | Usage courant | Effet possible sur H(z) |
|---|---|---|---|
| 8 bits | 49,92 dB | Microcontrôleurs simples, capteurs bas coût | Quantification visible, bruit et imprécision plus élevés. |
| 12 bits | 74,00 dB | Instrumentation industrielle | Compromis correct pour de nombreuses boucles de commande. |
| 16 bits | 98,08 dB | Audio, DSP embarqué | Très bon niveau pour la majorité des filtres numériques. |
| 24 bits | 146,24 dB | Audio studio, calcul haute précision | Grande marge pour les traitements sensibles. |
Exemple d’interprétation rapide
Supposons que vous calculiez une fonction avec un gain DC proche de 1, des pôles complexes conjugués de module 0,64 et des zéros près de z = -1. Vous pouvez déjà conclure plusieurs choses : le système transmet correctement les basses fréquences, la dynamique est amortie et stable, et les hautes fréquences seront plutôt atténuées. Si la courbe de magnitude montre une coupure progressive sans pic résonant important, vous êtes probablement face à un filtre passe-bas bien amorti. À l’inverse, si les pôles s’approchent du cercle unité avec une phase qui change vite, un pic de résonance apparaîtra souvent autour de la fréquence propre discrète.
Erreurs courantes lors du calcul
- Confondre les signes des coefficients du dénominateur : l’écriture normalisée doit être cohérente entre équation aux différences et H(z).
- Oublier que la stabilité se juge dans le cercle unité : appliquer un critère continu en discret conduit à des erreurs.
- Utiliser une fréquence d’échantillonnage trop faible : le repliement spectral et une mauvaise représentation dynamique deviennent alors probables.
- Négliger la quantification : surtout pour les filtres d’ordre élevé ou les applications embarquées.
- Interpréter uniquement le gain : il faut aussi examiner la phase, les pôles et la bande utile.
Quand utiliser ce calculateur
Ce calculateur convient particulièrement lorsque vous avez déjà les coefficients d’un filtre ou d’un système discret et que vous souhaitez : vérifier la stabilité, obtenir la forme analytique de H(z), estimer le gain statique, tracer une réponse fréquentielle rapide et identifier les pôles et les zéros principaux. Il est idéal pour l’enseignement, la vérification de prototypes et l’analyse technique de premiers et deuxièmes ordres. Pour des systèmes d’ordre supérieur, la logique reste identique, mais il est préférable d’utiliser ensuite des outils de synthèse plus avancés.
Ressources académiques et institutionnelles recommandées
Pour approfondir le calcul de fonction de transfert en z, consultez ces sources reconnues : MIT OpenCourseWare, Stanford University, University of Illinois.
Conclusion
Le calcul de fonction de transfert en z constitue l’un des langages les plus puissants de l’ingénierie numérique moderne. Il sert autant à analyser qu’à concevoir. En partant de quelques coefficients, vous pouvez déterminer le comportement global d’un système discret : stabilité, filtrage, résonance, gain et sensibilité fréquentielle. Bien maîtrisé, ce calcul permet de passer d’une simple équation aux différences à une lecture complète du système. Le calculateur interactif ci-dessus vous donne précisément cette passerelle, avec une interprétation immédiate et une visualisation concrète de la réponse en fréquence.