Transformée en Z

Un article de Wikipédia, l'encyclopédie libre.

Vous avez de nouveaux messages (diff^?).

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

1 Définition
2 Existence de la transformée en z
3 Propriétés de la transformée en z
4 Transformée en z inverse
5 Relation avec les autres transformées
- 5.1 Transformée de Laplace
- 5.2 Transformée de Fourier discrète
6 Transformées en z usuelles

[modifier] 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 \mathbb{C}$

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.

[modifier] Existence de la transformée en z

Le domaine de convergence est le sous-ensemble de $\mathbb{C}$ dans lequel la série converge. On l'appelle également couronne de convergence.
Autrement dit, le domaine de convergence est défini par l'ensemble :

$\left\{z : \sum_{n=-\infty}^{\infty}x[n]z^{-n} < \infty\right\}$

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é.

[modifier] Propriétés de la transformée en z

[modifier] Linearité

La transformée en z d'une combinaison linéaire de deux signaux est la combinaison linéaire des transformée 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)\} \$

[modifier] 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)\} \$

[modifier] 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)\} \$

[modifier] Multiplication par une exponentielle

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

[modifier] 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ù $\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 $-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)\$

[modifier] 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)$

[modifier] 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)$

[modifier] 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

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 :

$NUM(z)=num_0 + num_1 (1/z)^1 + num_2 (1/z)^2 + \cdots + num_{11} (1/z)^{11}$

N U M (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:

D E N O M (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.

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.

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 :

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).

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, $r 0 = 0, r 1 = 1, r 2 = 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.

r n - 1 + 2 r n + r n + 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):

$(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 $f n = f n - 1 + f n - 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):

$(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

$\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 à

$\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.

$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 :

$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 :

$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, $z 0$ et $z 1$ 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 $z 0$ et $z 1$ : $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:

$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 - z 0)$ peut se travailler ainsi:

$\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 :

$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]$ .

Puis nous remplaçons z0 et z1 par leur définition d'origine :

$f_n = \frac{1}{\sqrt{5}}\left[\left(\frac{1+\sqrt{5}}{2}\right)^n - \left(\frac{1-\sqrt{5}}{2}\right)^n\right]$ .

C'est la formule de Binet, son auteur est Jacques Philippe Marie Binet.

[modifier] Relation avec les autres transformées

[modifier] 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^{s T} \$

où $T = 1/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

$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)$

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 :

$X_q(s) = \int_{0^-}^{\infty} x_q(t) e^{-s t} \,\mathrm dt$

$\ = \int_{0^-}^{\infty} \sum_{n=-\infty}^{\infty} x[n] \delta(t - n T) e^{-s t} \,\mathrm dt$

$\ = \sum_{n=-\infty}^{\infty} x[n] \int_{0^-}^{\infty} \delta(t - n T) e^{-s t} \,\mathrm dt$

$\ = \sum_{n=-\infty}^{\infty} x[n] e^{-n s T}.$

Ce qui est précisement 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^{s 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(s) = X(z) \Big|_{z=e^{sT}}.$

[modifier] Transformée de Fourier discrète

La transformée en z est la généralisation de la transformée de Fourier en temps discret (TFtD). La TFtD peut être trouvé en évaluant $X(z)\$ en $z=e^{j\omega}\$ ou, en d'autres termes, en l'évaluant sur le cercle unité.

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

[modifier] 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\|\,$