Transformation en Z

Transformation en Z

Transformée en Z

La transformée en Z est un outil mathématique de traitement du signal, qui est l'équivalent discret de la transformée de Laplace.

Elle est utilisée entre autres pour le calcul de filtres numériques à réponse impulsionnelle infinie et en automatique pour modéliser des systèmes dynamiques de manière discrète.

La transformée de Fourier discrète est un cas particulier de la transformée en Z.

Sommaire

Définition

Sa définition mathématique est la suivante : la transformation en Z est une application qui transforme une suite s (définie sur les entiers) en une fonction S d'une variable complexe nommée z, telle que

S(z) = \mathcal{Z}\{s(n)\} =\sum_{n=-\infty}^{+\infty}s(n)z^{-n},\quad z \in \{z\in\mathbb{C}|\sum_{n=-\infty}^{+\infty}s(n)z^{-n} \quad converge\}

La variable n représente en général le temps discrétisé, la variable complexe z ne représente rien de particulier, il s'agit d'une création purement abstraite. Lorsqu'on travaille sur s(n) on dit que l'on est dans le domaine temporel, lorsqu'on travaille sur S(z) le domaine est appelé fréquentiel par analogie avec la transformée de Fourier, mais strictement parlant on est dans un domaine parfaitement abstrait.

Si \forall n<0,\ s(n)=0, on parle de signal causal. Inversement, si \forall n>0,\ s(n)=0, on parle de signal anti-causal.

Existence de la transformée en Z

Le domaine de convergence est le sous-ensemble de \mathbb{C} dans lequel la série converge.
Autrement dit, le domaine de convergence de la transformée en z de la suite (x_{n})_{n\in\mathbb{Z}} est l'ensemble :

\left\{z\in\mathbb{C} | \sum_{n=-\infty}^{\infty}x_{n}z^{-n} \quad\mathrm{existe}\right \}

On l'appelle également couronne de convergence. En effet, en posant z=\rho e^{i\theta}~, il vient :

|S(z)|=\left| \sum_{n=-\infty}^{\infty}x_{n}z^{-n}\right|\leqslant \sum_{n=-\infty}^{\infty}\left|x_{n}\right|\rho^{-n}

Donc S(z) existe si xn a une croissance au plus exponentielle, auquel cas le domaine de convergence est compris dans une couronne :

- de petit rayon le majorant de la base du côté des n négatifs
- de grand rayon le majorant de la base du côté des n positifs

Dans toute la suite de l'article, les transformées en Z ne seront valables que dans ce domaine de convergence sans que cela soit reprécisé.

Propriétés de la transformée en Z

Linéarité

La transformée en Z d'une combinaison linéaire de deux signaux est la combinaison linéaire des transformées en Z de chaque signal.

\mathcal{Z}\{a_1 x_1(n) + a_2 x_2(n)\} = a_1 \mathcal{Z}\{x_1(n)\} + a_2 \mathcal{Z}\{x_2(n)\} \

Décalage temporel

Le décalage temporel d'un signal de k échantillons se traduit par la multiplication de la transformée en Z du signal par z−k.

\mathcal{Z}\{x(n-k)\} = z^{-k}\mathcal{Z}\{x(n)\} \

Convolution

La transformée en Z d'un produit de convolution est le produit des transformées en Z

\mathcal{Z}\{x(n) * y(n)\} = \mathcal{Z}\{x(n)\} \mathcal{Z}\{y(n)\} \

Multiplication par une exponentielle

\mathcal{Z}\{a^{n}x(n)\} = X\left(\frac{z}{a}\right)

Multiplication par la variable d'évolution

De façon générale :

\mathcal{Z}\{n^{k}x(n)\} = \left(-z \frac{\mathrm{d} }{\mathrm{d}z}\right)^{k}\mathcal{Z}\{x(n)\}\

où \textstyle\left(-z \frac{\mathrm{d} }{\mathrm{d}z}\right)^{k}\mathcal{Z}\{x(n)\} signifie que l'on applique k fois à \mathcal{Z}\{x(n)\} l'opérateur \textstyle -z\frac{\mathrm{d} }{\mathrm{d}z}

Si l'on écrit cette formule au rang k=1, on obtient la formule de dérivation :

\mathcal{Z}\{nx(n)\} = -z \frac{\mathrm{d} }{\mathrm{d}z}X(z)\

Théorème de la valeur initiale

Soit x(n)\, un signal causal et X(z)\, sa transformée en Z. Alors :

x(0) = \lim_{n \to 0}x(n)=\lim_{z \to +\infty}X(z)

Théorème de la valeur finale

Soit x(n)\, un signal causal et X(z)\, sa transformée en Z. Alors lorsque la limite existe, on peut écrire :

\lim_{n \to +\infty}x(n)=\lim_{z \to 1}(z-1)X(z)

Transformée en Z inverse

La transformée en Z inverse est donné par :

x(n) = Z^{-1} \{X(z) \}= \frac{1}{2 \pi j} \oint_{C} X(z) z^{n-1}\mathrm dz \

où C est un chemin fermé parcouru dans le sens inverse des aiguilles d'une montre et appartenant entièrement au domaine de convergence.

En pratique, ce calcul s'effectue souvent à l'aide du théorème des résidus et la formule devient dans le cas d'un signal causal :

x(n) = \sum_{ z_i={\rm p\hat{o}les\; de\; } z^{n-1}X(z) } \operatorname{Res}\{z^{n-1}X(z)\}_{z=z_i}\,


Autres méthodes d'inversion

D'autres méthodes d'inversion pour passer de X(z) à x(n) sont: la lecture à l'envers de la table des transformées usuelles; l'application des règles de décalage, de combinaisons linéaires, de produit de convolution. En désespoir de cause, on peut toujours essayer de procéder par identification en donnant à z k+1 valeurs numériques et en recherchant les coefficients x(0) à x(k) qui sont solutions d'un système de k+1 équations linéaires à k+1 inconnues. Ou bien essayer de trouver un développement de Taylor ou Maclaurin de la fonction à inverser. Un cas particulier favorable se présente lorsque la fonction X(z) est une fraction rationnelle. En effet lorsque :X(z) = \frac{P(z)}{Q(z)} , P et Q étant deux polynômes en 1/z, on peut effectuer la division jusqu'au degré de précision souhaité, et l'on obtient directement les valeurs numériques des coefficients x(n), n variant de 0 à m. En l'occurrence on adopte plutôt dans ce cas la notation H(z) = {NUM(z)}/{DENOM(z)}\ . La raison en est que, pour les systèmes discrets ou échantillonnés, la fonction de transfert s'écrit h(n) et sa transformée en Z se présente souvent sous cette forme de quotient entre une sortie (en z) et une entrée (en z): H(z) = {NUM(z)}/{DENOM(z)}\ . Un exemple concret pour illustrer cette démarche:


Quotient de polynômes en z, approximation numérique.

Attention, cette méthode est purement numérique, elle ne fournit pas l'expression analytique de la série inverse. Dans cet exemple, H(z) est le rapport de deux polynômes en 1/z. Le numérateur ressemble à la multiplication par 2 du dénominateur décalé de 1 période, mais on choisit des valeurs numériques un peu inexactes pour éviter un parfait quotient égal à 2/z. Le numérateur, de puissance 11, est une expression de la forme :

\textstyle \scriptstyle NUM(z)=num_0 + num_1 (1/z)^1 + num_2 (1/z)^2 + \cdots + num_{11} (1/z)^{11}
NUM(z) = 0 + 0(1 / z)1 + 2,3 * (1 / z)2 + 4,22 * (1 / z)3 + 6,2 * (1 / z)4 + 8,21 * (1 / z)5 + 10,2 * (1 / z)6 + 12,2 * (1 / z)7 + 12,22 * (1 / z)8 + 12,4 * (1 / z)9 + 12,4 * (1 / z)10 + 12,4 * (1 / z)11


Le dénominateur, de puissance 10, est:

DENOM(z) = 0 + 1,1 * (1 / z)1 + 2,1 * (1 / z)2 + 3,1 * (1 / z)3 + 4,1 * (1 / z)4 + 5,1 * (1 / z)5 + 6,1 * (1 / z)6 + 6,1 * (1 / z)7 + 6,2 * (1 / z)8 + 6,2 * (1 / z)9 + 6,2 * (1 / z)10

Ici la division des polynômes ne "tombe pas juste", nous nous contentons d'une approximation du quotient Q(z), de la forme

\sum_{n \geq 0} q_n(1/z)^n
jusqu'à la puissance 10 :
Q(z) = 0 + 2,090909 * (1 / z)1 − 0,155372 * (1 / z)2 + 0,040421 * (1 / z)3 + 0,0309047 * (1 / z)4 − 0,015368 * (1 / z)5 + 0,007694 * (1 / z)6 + 0,101526 * (1 / z)7 − 0,176646 * (1 / z)8 + 0,061258 * (1 / z)9 + 0,015904 * (1 / z)10

Le reste R(z) de cette division incomplète est:

R(z) = 0 + 0 * (1 / z)1 + 0 * (1 / z)2 + 0 * (1 / z)3 + 0 * (1 / z)4 + 0 * (1 / z)5 + 0 * (1 / z)6 + 0 * (1 / z)7 + 0 * (1 / z)8 + 0 * (1 / z)9 + 0 * (1 / z)10 + 0 * (1 / z)11 + 0,550806 * (1 / z)12 − 0,413006 * (1 / z)13 − 0,063683 * (1 / z)14 + 0,040876 * (1 / z)15 − 0,052647 * (1 / z)16 − 0,011071 * (1 / z)17 + 0,616793 * (1 / z)18 − 0,478404 * (1 / z)19 − 0,098602(1 / z)20

On peut vérifier sur un tableur (du commerce ou open-source) ou à la main que ces polynômes répondent bien à la définition de la division euclidienne: H(z) = NUM(z)/DENOM(z)= Q(z)+ R(z)/DENOM(z). On suppose que le reste est négligeable par rapport aux coefficients du quotient. Les schémas de ces divers polynômes peuvent être visualisés sur un tableur comme suit.

ZtransfoNumDenomQuotient.PNG
ZtransfoPolynReste.PNG


Par curiosité on peut afficher la réponse impulsionnelle de l'approximation Q(z) de H(z). De même on peut afficher la réponse indicielle de Q(z) à un échelon de Heavyside.

Ztransfo2reponses.PNG

Si nous nous contentions d'une approximation moins précise de H(z) par le quotient Q(z), de la forme

\sum_{n \geq 0} q_n(1/z)^n
jusqu'à la puissance 5 par exemple :
\textstyle \scriptstyle Q(z)=0 + 2,090909* (1/z)^1 -0,155372* (1/z)^2 + 0,040421* (1/z)^3 + 0,0309047* (1/z)^4 -0,015368 *(1/z)^5 + 0
, nous obtiendrions des courbes de réponse lègèrement différentes, beaucoup moins précises (imprécision 6 fois plus forte environ). Le choix du degré d'approximation, autrement dit du meilleur compromis entre la précision et la lourdeur des calculs, est dicté par l'examen concret du problème spécifique que l'on traite.


Procédé par identification approximative des coeffs de X(z).

Pour passer de X(z) à x(n) , si aucune méthode ne semble déboucher, en désespoir de cause on peut toujours essayer de procéder par identification en donnant à z k+1 valeurs numériques et en recherchant les coefficients x(0) à x(k) qui sont solutions d'un système de k+1 équations linéaires à k+1 inconnues. Exemple:


Utilisation des fractions rationnelles, exemple de la fonction de transfert de la Suite de Fibonacci.

La Suite de Fibonacci peut être vue comme une suite f(n)= 0, 1 , 1, 2, 3, 5, 8, 13 etc. qui dépend donc de la suite "rampe" r(n)= 0, 1 ,2 ,3 , 4 etc. Ces deux suites peuvent faire l'objet de calculs utilisant l'outil de la transformée en Z. Nous établirons d'abord deux identités remarquables issues de considérations banales. Premièrement, la suite "rampe unitaire", r, qui n'est rien d'autre que la suite des entiers naturels, r0 = 0,r1 = 1,r2 = 2, 3, 4 , 5 , 6 etc est telle qu'un nombre est égal à la moyenne entre son précédent et son suivant, soit

0 1 2 3 4 5 6 7 8 9
0-0-2-4-6 -8 -10 -12 etc
0+0+0 1 2 3 4 5 etc

sauf au début avec 0 et 1, tous les termes de cette combinaison ont un résultat nul.

rn − 1 + 2rn + rn + 1 = 0

à partir d'un certain rang, ce qui implique la multiplication suivante de polynômes en 1/z, le premier étant le trinôme (1 -2 1):

\textstyle \scriptstyle (1-2 (1/z)^1+1 (1/z)^2) \cdot (0 + 1 (1/z)^1 + 2 \cdot (1/z)^2 + 3 \cdot (1/z)^3 + \cdots + r_{n} (1/z)^n + \cdots )= 1 (1/z)^1

Deuxièmement, la suite de Fibonacci est par définition fn = fn − 1 + fn − 2, à partir d'un certain rang, ce qui s'illustre ainsi

0 1 1 2 3 5 8 13 etc
- 0-1-1-2-3-5-8 etc
----0-1-1-2-3-5-8 etc

Sauf au début avec 0 et 1, tous les termes de cette combinaison ont un résultat nul et ceci implique la multiplication suivante de polynômes en 1/z, le premier étant le trinôme (1,-1,-1):

\textstyle \scriptstyle (1-1 (1/z)^1-1 (1/z)^2) \cdot (0 + 1 (1/z)^1 + 1 \cdot (1/z)^2 + 2 \cdot (1/z)^3 + \cdots + f_{n} (1/z)^n + \cdots )= 1 (1/z)^1

De ces deux multiplications nous pouvons déduire l'égalité de deux fractions

\textstyle \scriptstyle \frac{0 + 1 (1/z)^1 + 1 \cdot (1/z)^2 + 2 \cdot (1/z)^3 + \cdots + f_{n} (1/z)^n )+ \cdots }{0 + 1 (1/z)^1 + 2 \cdot (1/z)^2 + 3 \cdot (1/z)^3 + \cdots + r_{n} (1/z)^n )+ \cdots } = \frac{1-2 (1/z)^1+1 (1/z)^2}{1-1 (1/z)^1-1 (1/z)^2}

Fraction qui est égale finalement à

\textstyle \scriptstyle \frac{0 + 1 (1/z)^1 + 1 \cdot (1/z)^2 + 2 \cdot (1/z)^3 + \cdots + f_{n} (1/z)^n )+ \cdots }{0 + 1 (1/z)^1 + 2 \cdot (1/z)^2 + 3 \cdot (1/z)^3 + \cdots + r_{n} (1/z)^n )+ \cdots } = \frac{(z-1)^2}{+z^2-z-1}

Cette formule que l'on peut exprimer en termes de transformée en Z n'est autre que celle de la fonction de transfert entre la rampe et la suite de Fibonacci.

\textstyle \scriptstyle H(z)=\frac{ F(z) }{R(z) } = \frac{(z-1)^2}{z^2-z-1}

Comme la transformée R de la rampe fait partie des formules usuelles (\frac{z}{(z-1)^2} ), ceci nous donne :

\textstyle \scriptstyle H(z)=\frac{ F(z) }{ \frac{ z}{ (z-1)^2} } = \frac{(z-1)^2}{z^2-z-1}

Après simplifications nous voici dès lors en possession d'une formule simple de la transformée en Z de la suite de Fibonacci :

\textstyle \scriptstyle F(z) = \frac{z}{z^2-z-1}

Nous pourrions nous en contenter, mais si l'on souhaite obtenir une formule analytique de la suite de Fibonacci , f(n), il faut procéder à la transformation inverse. La méthode des fractions rationnelles peut être tentée. Le dénominateur possède deux pôles, z0 et z1 qui sont le nombre d'or :z_0 = \varphi = {1 + \sqrt{5} \over 2} et l'opposé de la section dorée :z_1 = 1-\varphi = {1-\sqrt{5} \over 2}. Pour les calculs rencontrés ci-dessous on se servira des propriétés suivantes de z0 et z1 : z_0-z_1=(2 \cdot z_0-1)=\sqrt{5}, et : (z-z_0) \cdot (z-z_1)=z^2-z-1. La fonction se décompose en fractions rationnelles élémentaires que l'on réécrit un peu:

\textstyle \scriptstyle F(z) = \frac{z}{z^2-z-1}= \frac1\sqrt5 \cdot\left( \frac{z_0}{z-z_0}- \frac{z_1}{z-z_1} \right) = \frac1\sqrt5 \cdot\left( z_0\cdot\frac1{z-z_0} - z_1\cdot\frac1{z-z_1} \right).

Une fraction du type 1 / (z − z0) peut se travailler ainsi:

\textstyle \scriptstyle \frac{ 1}{(z- z_0 )}= \frac{ z}{(z- z_{0})} \cdot \frac{ 1}{z}

La première partie étant la transformée de la formule usuelle exponentielle , z_0^{n} , la seconde partie 1/z étant le retard pur d'un cran. Si bien que la transformée inverse de cette fraction élémentaire est z_0^{n-1}, en appliquant les règles de combinaisons linéaire nous calculons la suite cherchée :

\textstyle \scriptstyle f_n = \frac{1}{\sqrt{5}}\left( z_{0} \cdot z_{0}^{n-1} - z_{1} \cdot z_{1}^{n-1}\right) =\frac{1}{\sqrt{5}}\left( z_{0}^{n} - z_{1}^{n}\right) .

Relation avec les autres transformées

Transformée de Laplace

La transformée en Z est simplement la transformée de Laplace d'un signal discret en effectuant la substitution :

z \equiv e^{p T} \

où T = \tfrac1{f_s} est la période d'échantillonnage (en unité de temps i.e. secondes) et f_s \ est la fréquence d'échantillonnage (en échantillon par seconde ou hertz)

Pour démontrer cela, posons :

q(t) \equiv \sum_{n=-\infty}^{\infty} \delta(t - n T)

qui est un peigne de Dirac et

\begin{align} x_q(t) &\equiv x(t) q(t) = x(t) \sum_{n=-\infty}^{\infty} \delta(t - n T) \\ & = \sum_{n=-\infty}^{\infty} x(n T) \delta(t - n T) = \sum_{n=-\infty}^{\infty} x[n] \delta(t - n T) \end{align}

qui est la représentation en temps continu du signal x(t) \ échantillonné et où x[n] \equiv x(nT) \ sont les échantillons de x(t) \ .

La transformée de Laplace du signal échantillonné x_q(t) \ est :

\begin{align}X_q(s) &= \int_{0^-}^{\infty} x_q(t) e^{-p t} \, \mathrm dt \\ &= \int_{0^-}^{\infty} \sum_{n=-\infty}^{\infty} x[n] \delta(t - n T) e^{-p t} \, \mathrm dt \end{align}

Le passage du symbole d'intégration à l'intérieur de la série est difficile à justifier si le résultat final n'est pas convergent.

\begin{align}X_q(s) &= \sum_{n=-\infty}^{\infty} x[n] \int_{0^-}^{\infty} \delta(t - n T) e^{-p t} \, \mathrm dt \\ &= \sum_{n=-\infty}^{\infty} x[n] e^{-n p T} \end{align}

Ce qui est précisément la définition de la transformée en Z d'un signal discret x[n] \

X(z) = \sum_{n=-\infty}^{\infty} x[n] z^{-n}

avec la substitution z \leftarrow e^{p T} \ .

En comparant les deux dernières équations, nous trouvons la relation entre la transformée en Z et la transformée de Laplace d'un signal discret :

X_q(p) = X(z) \Big|_{z=e^{pT}}.

Transformée de Fourier discrète

La transformée en Z est la généralisation de la Transformée de Fourier discrète (TFD). La TFD peut être trouvé en évaluant X(z) en z = eiω ou, en d'autres termes, en l'évaluant sur le cercle unité.

X(e^{i \omega}) = X(z) \Big|_{z=e^{i \omega}}.

Transformées en Z usuelles

\delta[n] \, représente l'impulsion de Dirac et u[n] \, la fonction de Heaviside appelée aussi échelon unitaire, step en anglais.

  Signal x(n) Transformée en Z X(z) Domaine de convergence
1 \delta[n] \, 1\, \mathbb{C}\
2 u[n] \, \frac{1}{1-z^{-1}} |z| > 1\,
3 a^n u[n] \, \frac{1}{1-a z^{-1}} |z| > |a|\,
4 n a^n u[n] \, \frac{az^{-1} }{ (1-a z^{-1})^2 } |z| > |a|\,
5 -a^n u[-n-1] \, \frac{1}{1-a z^{-1}} |z| < |a|\,
6 -n a^n u[-n-1] \, \frac{az^{-1} }{ (1-a z^{-1})^2 } |z| < |a|\,
7 \cos(\omega_0 n) u[n] \, \frac{ 1-z^{-1} \cos(\omega_0) }{ 1-2z^{-1}\cos(\omega_0)+ z^{-2} } |z| >1\,
8 \sin(\omega_0 n) u[n] \, \frac{ z^{-1} \sin(\omega_0) }{ 1-2z^{-1}\cos(\omega_0)+ z^{-2} } |z| >1\,
9 a^n \cos(\omega_0 n) u[n] \, \frac{ 1-a z^{-1} \cos( \omega_0) }{ 1-2az^{-1}\cos(\omega_0)+ a^2 z^{-2} } |z| > |a|\,
10 a^n \sin(\omega_0 n) u[n] \, \frac{ az^{-1} \sin(\omega_0) }{ 1-2az^{-1}\cos(\omega_0)+ a^2 z^{-2} } |z| > |a|\,

Lien interne

Filtre numérique


  • Portail des mathématiques Portail des mathématiques
Ce document provient de « Transform%C3%A9e en Z ».

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Transformation en Z de Wikipédia en français (auteurs)

Игры ⚽ Нужно решить контрольную?

Regardez d'autres dictionnaires:

  • transformation — [ trɑ̃sfɔrmasjɔ̃ ] n. f. • 1375; lat. transformatio 1 ♦ Action de transformer, opération par laquelle on transforme. ⇒ conversion. La transformation des matières premières. Industrie de transformation. « un de ces théoriciens qui ont rêvé [...]… …   Encyclopédie Universelle

  • Transformation — (root transform ) may refer to:Transformation is also referred to as a turn.In science: * Transformation (geometry), in mathematics, as a general term applies to mathematical functions. ** Data transformation (statistics) in statistics. *… …   Wikipedia

  • transformation — is the process by which transactions that are open at, or after, the record date/election transfer date are cancelled and/or replaced by new transactions in accordance with the terms of the reorganisation. Euroclear Clearing and Settlement… …   Financial and business terms

  • Transformation — (lateinisch die Umformung) bezeichnet allgemein die Veränderung der Gestalt, Form oder Struktur. Der Begriff wird in einer Vielzahl von Themengebieten verwendet: Konvertierung (Informatik) in der Informatik die Überführung von Daten in ein… …   Deutsch Wikipedia

  • transformation — Transformation. s. f. v. Changement d une forme en une autre. La transformation des metaux. la transformation fabuleuse de Protée. les Poëtes Grecs sont pleins de transformations fabuleuses …   Dictionnaire de l'Académie française

  • Transformation — Trans for*ma tion, n. [L. transformatio: cf. transformation.] The act of transforming, or the state of being transformed; change of form or condition. Specifically: [1913 Webster] (a) (Biol.) Any change in an organism which alters its general… …   The Collaborative International Dictionary of English

  • transformation — (n.) c.1400, from O.Fr. transformation and directly from L. transformationem (nom. transformatio), noun of action from pp. stem of transformare (see TRANSFORM (Cf. transform)) …   Etymology dictionary

  • transformation — transformation. См. трансформация. (Источник: «Англо русский толковый словарь генетических терминов». Арефьев В.А., Лисовенко Л.А., Москва: Изд во ВНИРО, 1995 г.) …   Молекулярная биология и генетика. Толковый словарь.

  • Transformation — (lat., »Umgestaltung«), in der Mathematik zunächst soviel wie Substitution (s. d.), dann aber auch jedes Verfahren, durch das aus gegebenen Figuren neue Figuren von andrer Lage und andrer Gestalt abgeleitet werden. Eine T. in diesem Sinn erhält… …   Meyers Großes Konversations-Lexikon

  • Transformation — Transformation, gewöhnlich so viel wie Substitution (s.d.). Zwischen Scharen von Transformationen und von Substitutionen besteht jedoch der Unterschied, daß bei ersteren die Parameter stetig veränderliche Größen sind, während sie bei letzteren… …   Lexikon der gesamten Technik

  • Transformation — (lat.), Umbildung; transformieren, umformen, umbilden …   Kleines Konversations-Lexikon

Share the article and excerpts

Direct link
… Do a right-click on the link above
and select “Copy Link”