Discussion :

[melane]

Bonjour,

Je suis novice sur Matlab et j'ai quelques soucis de calculs.

j'ai un calcul du type
x=a*(k/w)*(cos(k(u+i))/cos(ki))*sin(kj-wt)

Toutes mes lettres représentent des matrices d'une colonne et pleins de lignes (aux alentours de 10^6) sauf u et t qui sont deux variables.
J'ai attribué un pas de temps à t et un pas de distance à u. Ces deux variables sont donc représentées par des matrices d'une ligne et plusieurs colonnes.

Tout d'abord je suis obligée de décomposer mon calcul car out of memory (en même temps vu le nombre de données...).

Mais mon plus gros problème c'est que du coup, mes variables appliques des dimensions de matrices différentes ce qui m'empeche d'appliqué mon calcul.

Je voudrais faire ce calcul pour chaque temps et distance données.

J'ai tenté les boucles for mais je ne sais pas si j'utilise la bonne fonction ou si je l'ai mal faite, en tout cas ça n'a pas marché.

Je ne sais pas si j'ai été assez claire, mais auriez-vous des idées?

Merci

[shaiHulud]

Bonjour

- Remplace tous tes produits et division * et / par les produits /division termes à termes .* et ./
- Mets tous tes vecteurs avec la même orientation, en transposant les lignes en colonnes ou le contraire (transposé de x <=> x' ; x.' si x est complexe)

Cela suffit en général à éviter les boucles simples (c'est le principe de calcul vectorisé en matlab). Pour des formules plus ardues, on peut préférer une boucle.

NB: cela suppose que les valeurs u(n) t(n) x(n) etc coincident. Si tu cherches à croiser tes variables ([u(1) x(1), u(1) x(2), u(2) x(1)...]), utilise meshgrid, comme indiqué par FLB, et fais toutes les opérations vectorialisées sur des matrices plutôt que des vecteurs

[FLB]

Salut,
Pour chaque couple (u,t) x est un tableau de 10^6 valeurs?
Ce petit lien explique comment va évoluer la mémoire de matlab :
http://matlab.developpez.com/faq/?pa...mation_memoire
Chaque couple (u,t) prendra 1e6*8/1024/1024 = 7.6Mo de mémoire -> tu as vraiment besoin de tout stocker au préalable?

Si tu veux parcourir u pour chaque valeur de t (et vice versa), regarde ce que fait la méthode meshgrid.

[melane]

Bonjour,

Merci pour vos réponses,
VV33D j'ai déjà fait ce que tu me suggères et ça m'a super bien dépanné pour la gestion de la mémoire. Mais cela ne fonctionne pas pour ce problème ci.


FLB je vais tenter de voir avec meshgrid. Je n'y avais pas pensé merci. Je ne suis pas obligée de tout stocker, je peux très bien regarder chaque ligne indépendament des autres mais in fine il va falloir que je les visualise ensemble. Et comme j'ai des tas de fichier à traiter je me suis dis que c'était peut être la meilleur solution.

Je vous tiens au courant!

[melane]

Alors après quelques vérifications j'ai remarqué que j'avais en effet "oublié" quelques points à mes opérations.
Du coup je les ai rajouté et tout est nickel, ça marche avec une boucle et pas besoin de faire un meshgrid, à une exception. Il me fixe les variables sur une valeur (alors que je les avais défini avec des pas).

Par exemple j'ai défini:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
t=(0:5:40)
u=(0:1:31)
for u=u
        for t=t
             x=a.*(k./w).*(cos(k(u+i))./cos(ki)).*sin(kj-wt)
        end
end

et quand mon calcul se termine j'ai une valeur unique:
t=15
u=30
alors que je voudrais voir la variation de chaque ligne de mes tableaux (10^6) avec chaque pas de t et u.

Est-ce que je me suis trompée quelque part ou bien je fais totalement fausse route?

[Jerome Briot]

Essaie ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
ku = 0;
 
for u = 0:1:31
 
    ku = ku+1;
    kt = 1;
 
    for t = 0:5:40
 
        x(ku,kt) = a.*(k./w).*(cos(k(u+i))./cos(ki)).*sin(kj-wt);
        kt = kt+1;
 
    end
 
end

Par curiosité, tu pourrais nous montrer ce que renvoie ceci à la fin de l'exécution du code ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

whos a k w i ki kj wt