matlab : calcul norme dans L2

Avec matlab, ça ne prend que 10 caractères : si X est le vecteur, il suffit de faire

sqrt(X.*X)

(NB : .* désigne le produit élément par élément de deux matrices de même taille)

Ah oui, j'ai oublié de sommer les composantes : il faut faire

sqrt(sum(X.*X))

Euh... Oui. Tu as vraiment besoin de moi pour remplacer X par (u1 - u2) dans une formule?

Ah bon? Il faudra que tu m'expliques pourquoi.

Bonjour,
Je suis trés serieux, parceque j'ai passé pas mal du temps en cherchant la bonne réponse, mais j'était faché contre la réponse de Luca quand j'ai demandé a lui le détail.. en tout cas le norm dasn $L^2 $ s'est un integral alors il faut tout d'abord connaitre la propriété du vecteur qu'on veux le normé, et comme mon vecteur viens de la solution d'une EDP en utilisant method EF alors je peut calculer l'integral untilisant des method approché sur chaque triangle du maillage, et en fin je somme sur le domain.

Jaaour, non tu n'es pas sérieux. Pendant 4 ou 5 messages tu as laissé les gens croire que tu cherchais comment calculer la norme d'un vecteur, et Lucas t'a répondu la bonne réponse, alors que ce n'est pas ça que tu cherches. Et oui les mots ont un sens et ta question se pose correctement sous la forme suivante :

"si $X$ est le vecteur obtenu par la méthode des éléments finis, qui représente la fonction $u$ solution du problème dans une base de fonctions, comment calculer la norme $L^2$ de $u$ en fonction de $X$ ?"

Avoue que ça n'a {\it rien à voir} avec la norme d'un vecteur, qui est effectivement $||X||=\sqrt{\sum_i x_i^2}$. Relis ton premier message ; comment voulais-tu que quelqu'un devine la question que tu te posais vraiment avec ce que tu as écrit ? Pourquoi n'avoir pas directement été précis sur ta demande ?

Pour ta question, il y a en tout premier lieu une question de maths, avant l'implémentation en Matlab : si les $w_i$ sont tes fonctions de base, et $X=(x_i)$ le vecteur des coordonnées de $u$, alors $u=\sum_i x_i w_i$ et donc
$$||u||_{L^2}^2=\langle \sum_i x_i w_i, \sum_j x_j w_j \rangle_{L^2}=\sum_{i,j} x_i x_j \langle w_i,w_j \rangle_{L^2}=X'QX$$ où $X'$ est la transposée de $X$ et $Q$ la matrice carrée dont le $(i,j)$-ième coefficient est le produit scalaire $L^2$ : $q_{ij}=\langle w_i,w_j \rangle_{L^2}=\int_{\Omega} w_iw_j$, c'est-à-dire la matrice du produit scalaire $L^2$ dans la base $(w_i)$. La première chose à faire est de calculer la matrice $Q$, à la main (normalement tu l'as essentiellement déjà fait si ton EDP comporte un terme en $u$). Une fois ceci fait tu auras pour tout $u$, de coordonnées $X$ : $||u||_{L^2}=\sqrt{X'QX}$ ce qui en Matlab doit s'écrire {\bf sqrt(X'*Q*X)}.

Avoir la volonté, c'est bien ! C'est déjà ça. Après, tu peux commencer par lire un des premiers messages de ce fil. Ca date de quatre ans, mais c'est sans doute toujours d'actualité.