Manipuler des nombres complexes (et quaternions, et autres objets) avec Maxima

Publié le 18 décembre 2012 par raf

Publié dans Algèbre, Informatique et mathématiques Laisser un commentaire

Le but de cet article est de présenter comment manipuler des quaternions avec le logiciel Maxima.

Pour commencer, il faut charger le package atensor qui va nous permettre de réaliser tous les calculs:

load(atensor);

Ensuite, on va effectuer les deux commandes suivantes:

init_atensor(clifford, 0, 0, 2);
q(a,b,c,d):=a+b.v[1]+c.v[2]+d.v[1].v[2];

La première indique que l’on veut travailler avec les quaternions. Il s’agit en fait de la définition de la forme quadratique qui définit l’algèbre de Clifford. Le premier 0 représente la dimension positive de l’algèbre, le second la dimension dégénérée et la dernière la dimension négative. Ainsi, ici, on considère la forme quadratique $Q(x_1,x_2)=-x_1^2-x_2^2$ définie sur $\mathbb R^2$.
La deuxième ligne nous donne un moyen de définir facilement des quaternions, via leurs composantes.

Maintenant, on peut commencer à travailler avec les quaternions. Les deux vecteurs de base sont stockés dans un vecteur v. Ainsi, i correspond à v[1] et j correspond à v[2]. Le vecteur k, quand a lui, correspond à v[1].v[2]. Par exemple, si je veux définir le nombre complexe $z=1+2i$, je fais

z:1+2.v[1];

Si je veux le multiplier avec le quaternion $3+4i+5j$, je fais:

(1+2*v[1]).(3+4*v[1]+5*v[2]);

et le résultat est
$2\,\left( 5\,\left( {v}_{1} . {v}_{2}\right) +4\,\left( {v}_{1} . {v}_{1}\right) +3\,{v}_{1}\right) +5\,{v}_{2}+4\,{v}_{1}+3$
Si on veut un résultat plus joli, on peut utiliser la commande atensimp:

atensimp((1+2*v[1]).(3+4*v[1]+5*v[2]));

qui nous renvoie
$
2\,\left( 5\,\left( {v}_{1} . {v}_{2}\right) +3\,{v}_{1}-4\right) +5\,{v}_{2}+4\,{v}_{1}+3
$
En combinant cela avec la fonction expand

expand(atensimp((1+2*v[1]).(3+4*v[1]+5*v[2]));)

on trouve finalement
$
10\,\left( {v}_{1} . {v}_{2}\right) +5\,{v}_{2}+10\,{v}_{1}-5
$
qui correspond au quaternion $10k+5j+10i-5$.

On peut définir quelques commandes pour nous simplifier la vie:

conj(x):=subst([v[1]=-v[1],v[2]=-v[2],v[1] . v[2]=-v[1] . v[2]] , x);
norm(x):=sqrt(expand(atensimp(x.conj(x))));
inverse(x):=conj(x).expand(atensimp(x.conj(x)))^(-1);

Ainsi, je peux calculer l’inverse du quaternion $10k+5j+10i-5$:

inverse(q(-5,10,5,10))

On peut aussi faire des calculs plus généraux, comme par exemple multiplier les quaternions $q(a,b,c,d)$ ainsi que $q(e,f,g,h)$. Si vous essayez de le faire via la commande

expand(atensimp(q(a,b,c,d).q(e,f,g,h)));

, vous verrez que le résultat est plutôt… horrible (tellement horrible que je vous l’éparnge ;-))! Le problème est que nos variables a, b, …, h ne sont pas considérées comme réelles. On corrige le problème avec le code suivant:

declare([a,b,c,d,e,f,g,h],scalar);
expand(atensimp(q(a,b,c,d).q(e,f,g,h)));

Remarque

Comme indiqué dans l’introduction, on peut travailler avec une forme plus générale de l’algèbre des quaternions en considérant l’une des algèbres de Clifford. Dans ce cas, il faudrait modifier quelques petites choses:

L’application de conjugaison n’est pas la bonne.
Tous les éléments non-nuls ne sont pas forcément inversibles (il faut se restreindre au groupe de Clifford).

Les algèbres de Clifford – I

Publié le 7 décembre 2012 par raf

Publié dans Algèbre, Géométrie hyperbolique Laisser un commentaire

Motivation

A propos de la classification des isométries de l’espace hyperbolique, on a deux résultats très connus:

Premièrement, les isométries de l’espace hyperbolique $H^{n+1}$ proviennent exactement, via l’extension de Poincaré, de transformations de Möbius de l’espace $\mathbb R^n$.
Deuxièmement, les isométries de l’espace $H^2$ (respectivement $H^3$) sont données par $\mathrm{PGL}(2,\mathbb R)$ (respectivement $\mathrm{PGL}(2,\mathbb C)$). De même, les isométries qui préservent l’orientation de l’espace $H^2$ (respectivement $H^3$) sont données par $\mathrm{PSL}(2,\mathbb R)$ (respectivement $\mathrm{PSL}(2,\mathbb C)$).

Les algèbres de Clifford fournissent une manière de généraliser cela aux dimensions supérieures.

Conventions et remarques

Le terme algèbre désigne pour moi une algèbre associative, unitaire (si, si) et de dimension finie.
Mon but est de présenter cette construction dans les grandes lignes. Je ne vais pas donner tous les détails car bien des textes le font mieux que moi (voir les références).

Construction – Idée de base

Soit $V$ un $K$-espace vectoriel ainsi que $Q$ une forme quadratique définie sur $V$. Pour paraphraser Wikipédia, l’algèbre de Clifford $Cl(V,Q)$ est l’algèbre (associative et unitaire) engendrée par $V$ la plus générale qui satisfait $v\cdot v=Q(v)\cdot 1$ pour tout $v$ dans $V$ (je vais probablement omettre bien souvent la multiplication par le $1$ de l’algèbre).

Remarque (pour casser un peu le suspense): le cas qui va nous intéresser le plus est $K=\mathbb R$ ainsi que $Q(x_1,\ldots,x_n)=-{x_1}^2-\ldots-{x_n}^2$. On verra qu’en faisant varier n, on obtiendra successivement les nombres réels, les nombres complexes ainsi que les quaternions de Hammilton.

Construction – Propriété universelle

On peut caractériser l’algèbre de Clifford via une propriété universelle: il s’agit de la donnée d’une algèbre (associative, unitaire) définie sur $K$ ainsi que d’une application $K$-linéaire $i:V\longrightarrow Cl(V,Q)$ telle que $i(v)^2=Q(v)\cdot 1$ ; le tout doit satisfaire la propriété universelle suivante

pour toute $K$-algèbre $A$ et toute application $K$-linéaire $f:V\longrightarrow A$ telle que $f(v)^2=Q(v)\cdot 1_A$, il existe un unique morphisme de $K$-alègbres $\tilde f:Cl(V,Q)\longrightarrow A$ tel que $\tilde f\circ i = f$:

$Algèbre de Clifford - Propriété universelle$

Ainsi, si l’algèbre de Clifford associée à un couple $(V,Q)$ existe, alors la propriété universelle donne:

L’algèbre de Clifford $Cl(V,Q)$ est unique (à isomorphisme près).
Il s’agit d’un foncteur covariant de la catégorie des $K$-espaces vectoriels munis d’une forme quadratique (avec comme morphisme les applications linéaires qui préservent la forme quadratique) dans la catégorie des $K$-algèbres (unitaires et associatives).

Construction explicite

On considère pour commencer l’algèbre tensorielle $T(V)$, c’est-à-dire $T(V)=\bigoplus_{n=0}^\infty V^{\otimes n}$ et la multiplication dans $T(V)$ est donnée par le produit tensoriel. Pour terminer, on a plus qu’à considérer l’idéal $I_Q$ engendré par tous les éléments $v\otimes v-Q(v,v)$ ainsi que le quotient $T(V)/I_Q$. Il est alors facile de vérifier la proposition suivante:

Proposition

Le quotient $T(V)/I_Q$ satisfait la propriété universelle ci-dessus. Ainsi, $Cl(V,Q)=T(V)/I_Q$.

Une base pour l’algèbre de Clifford

Si $e_1,\ldots,e_n$ est une $K$-base de $V$, on vérifie que l’ensemble $\Big\{ e_{i_1}\cdot\ldots\cdot e_{i_k}:1\leq i_1<i_2<\ldots<i_k\leq n\Big\}$ est une base de $T(V)/I_Q$ (le produit vide est autorisé et correspond au 1). Ainsi, l’algèbre de Clifford $Cl(V,Q)$ est une algèbre de dimension $2^{\textrm{dim}_KV}$.

Le cas d’une base orthonormée

Si l’on suppose maintenant $e_1,\ldots,e_n$ est orthonormée, on obtient les identités:

$e_i\cdot e_i=1 $, pour tout i;
$e_i\cdot e_j=-e_j\cdot e_i$ pour tous $i\neq j$.

Dans le prochain article, on verra quelques exemples et où ces algèbres nous mènent…

Pourquoi on fait tout cela, en fait? (et autres considérations)

Publié le 5 décembre 2012 par raf

Publié dans Réflexions générales Laisser un commentaire

Même si je n’ai commencé ma thèse que depuis très peu de temps (mais en fait trop, quand je vois à la vitesse à laquelle j’avance ;-)), j’ai eu l’occasion de suivre quelques exposés. En règle générale (mais j’avoue que cela varie selon mon humeur et d’autres facteurs, comme l’éventuelle présence de slides), je commence l’exposé, optimiste, en prenant des notes. Le temps passe, les définitions élémentaires laissent la place aux définitions spécifiques (et moins élémentaires) et, petit à petit, une question se formule dans ma petite tête:

Pourquoi on fait tout cela, déjà?

A ce moment là, il y a deux possibilités:

L’exposé concerne un domaine que je connais -même si très vaguement- ou bien m’intéresse (pour une raison quelconque).
L’exposé concerne un sujet qui m’est (à peut près) totalement inconnu.

Dans le deuxième cas (et surtout si l’exposé est très technique, et se résume à une suite de lemmes techniques et spécifiques -oui ,oui…-), j’ai très vite tendance à décrocher.

Oui, oui! J’aime quand l’orateur commence son exposé par une motivation, quand il essaie de me faire comprendre pourquoi je pourrais avoir envie de l’écouter! Cela ne me gêne pas -trop- de supporter un exposé très technique si je sais où je vais et pourquoi je le fais mais dans le cas contraire…

Voilà ce que j’ai envie de dire: commençons par une motivation (surtout si l’audience est constituée de non spécialistes)!

Pour le reste, il y a l’article, Comment éviter les mauvais exposés de colloquium, de Michèle Audin.

En ce qui concerne l’écrit:

Tips for authors, de Milne
« How to write mathematics badly », exposé donné par J-P Serre; partie 1/3, partie 2/3, partie 3/3

Maths

Monthly Archives: décembre 2012