Produit de matrices indexées

[einstein123]

Bonjour,

Mon problème est le suivant:

j'ai définit une suite de matrice qui représente réellement une fonction à double variable comme suit:

%Definition des fonctions de bases pour la plaque Wi(x,y)%
    
for i=1:n;
    for j=1:n;
        for kk=1:length(x);
            for jj=1:length(y);
            w0nxy(kk,jj,i,j)=w0n(kk,i)*w0n(jj,j);
            end
        end
    end
end

%Passage de l'indice double à l'indice simple pour calcul de tenseurs%
for i=1:n;
    for j=1:n;
        for kk=1:length(x);
            for jj=1:length(y);
                k=(3*(i-1))+j;
                w0nxyI(kk,jj,k)=w0nxy(kk,jj,i,j);
            end
        end
    end
end

Donc kk c'est x, jj c'est y...i c'est l'ordre de la fonction utilisée j également...
Je suis après passé à un indice unique k pour faciliter la manipulation donc j'ai maintenant la suite de matrice (ou fonctions à double variables): w0nxyI(kk,jj,k) c'est une matrice ou fonction x y d'indice k.


Maintenant ce que je veux réaliser c'est le produit de ces matrices (produit matriciel de la ième et j ème matrice) pour différents indices c'est à dire pour i et j donnés je fais:

w0nxyI(:,i)*w0nwyI(:,j) pour avoir une matrice M(i,j).

Je veux juste avoir une idée comment formuler cette manip, car vraiment je suis bloqué...

Merci pour vos aides précieuses....
-----

[lou_ibmix_xi]

Si j'ai bien compris, pour la matrice (k, n), il te suffit de faire:

w0nxyI(:, :, k) * w0nxyI(:, :, n)

1 remarque en passant, évite d'utiliser "i" comme indice, c'est le symbole pour l'imaginaire pur...

[einstein123]

Merci énormément pour la réponse, il me reste un petit autre blocage si cela ne vous gêne pas cher ami...
C'est que après avoir fait le produit comme suit:

for i=1:n;
    for j=1:n;
        Mn(:,:,i,j)=w0nxyI(:,:,i)*w0nxyI(:,:,j);
    end
end

Donc maintenant je veux faire l'intégration de la fonction Mn(:,:,i,j) sur dx et dy....
J'ai déjà fait celà pour une seule variable et ça a marché maintenant je fait ça

for i=1:n;
    for j=1:n;
        Mn(:,:,i,j)=w0nxyI(:,:,i)*w0nxyI(:,:,j);
    end
end

Mn(:,:,i,j) est une fonction de x et y (discrétisé par kk et jj), j'ai tapé la commande suivante:

for i=1:n;
    for j=1:n;
        Mn(:,:,i,j)=w0nxyI(:,:,i)*w0nxyI(:,:,j);
        M(i,j)=trapz(0:dx:1, 0:dy:1, Mn(:,:,i,j));
    end
end

et ça me donne l'erreur suivante:

??? Error using ==> permute
ORDER contains an invalid permutation index

Error in ==> trapz at 39
y = permute(y,perm);

Error in ==> Plaque at 146
M(i,j)=trapz(0:dx:1, 0:dy:1, Mn(:,:,i,j));

Merci à l'avance l'ami...

[Dlzlogic]

Si j'ai bien compris, pour la matrice (k, n), il te suffit de faire:

w0nxyI(:, :, k) * w0nxyI(:, :, n)

1 remarque en passant, évite d'utiliser "i" comme indice, c'est le symbole pour l'imaginaire pur...

Juste une question, le symbole 'i' serait donc en quelque sorte un mot-clé pour matlab ?

[A voir en vidéo sur Futura]

[lou_ibmix_xi]

Juste une question, le symbole 'i' serait donc en quelque sorte un mot-clé pour matlab ?

i*i

et ça me donne l'erreur suivante:

??? Error using ==> permute
ORDER contains an invalid permutation index

Error in ==> trapz at 39
y = permute(y,perm);

Error in ==> Plaque at 146
M(i,j)=trapz(0:dx:1, 0:dy:1, Mn(:,:,i,j));
Merci à l'avance l'ami...

Je ne peux pas plus t'aider que l'aide en ligne. Chez moi (octave), "help trapz", ne parle que de deux arguments, et tu en donnes trois.