Discussion :

[NELLLY]

Salut
ci-dessous une partie de mon programme

Code :

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

1
2
3
4
5
6
7
8
9
 
PB=(0:0.001:0.999)';
IPB=(1000*PB)+1;
C=cell(length(PB),7);
for i=1:length(PB)
    C{i,1}=zeros(IPB(i),length(p));
    size(C{i,1})%pour voir où ça plante
    pause
end

J'ai 2 messages qui s'affichent:

1. Size vector should be a row vector.

   Je ne sais pas pourquoi il affiche ce message alors que les valeurs de IPB sont entières et malgré ça il affiche ensuite la taille des C{i,1}.
2. Au dessus de i=639 le message
   

   Error using zeros, out of memory.

Que dois-je faire?

[FLB]

Code :

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

1
2
PB=(0:0.001:0.999)';
IPB=(1000*PB)+1;

Le problème est que IPB ne représente pas un entier mais un flottant (Pourquoi 0.3-0.2-0.1 est-il différent de 0 ?)
Utilise pour indexer :

Code :

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

IPB=round((1000*PB)+1);

Ce qui corrigera "Warning: Size vector should be a row vector with integer elements."

Ensuite pour les problèmes de mémoire, il est toujours bon de la préallouer

Code :

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

C{length(PB),1} = zeros(round(IPB(length(PB))),length(p));

Avant ta boucle for. Chez moi ca tourne sans soucis avec length(p)=100

[NELLLY]

Salut
merci pour la solution proposée. La première erreur est évitée cependant la deuxième est toujours présente. Vu que vous avez de l'expérience avec Matlab est-ce -que si j'augmente la capacité de la carte mémoire de mon pc je n'aurais plus de problème de mémoire avec MATLAB?

[Invité]

Plus de problème de mémoire, c'est un peu idyllique... disons moins

Pour N = length(p), tu essaies de stocker

Code :

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

N*sum(IPB)*8 /1024 /1024 /1024 (Go)

Ce qui fait pour N = 100 environ 372 Go !!
Comment déterminer approximativement la quantité de mémoire nécessaire pour le stockage d'une variable ?

As-tu réellement besoin tout cela en même temps ?
Que fais-tu par la suite avec C ?

[NELLLY]

J'ai oublié de noter que

Code :

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

length(p)=1001

car

Code :

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

p=(0:0.001:1)'

Pour chaque valeur de PB j'ai une matrice (1001*1001) et un cout à déterminer.
En premier lieu j'ai essayé d'utiliser une boucle for pour PB

Code :

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

for PB=0:0.001:0.999

et vu que j'ai deux autres variables de décision, la première est entière et la deuxième est réelle. Le temps d'exécution du programme est très grand j'ai attendu 6 heures et le calcul n'était pas encore terminé. Par suite j'ai modifié l'algorithme en traitant PB comme un vecteur et par la suite j'ai utilisé des multidimensional cell arrays pour que le calcul soit plus rapide et j'étais affronté à une mémoire insuffisante.
J'ai trouvé une alternative provisoire, je peux départager le vecteur PB en sous vecteurs et à chaque je cherche une solution et je compare ensuite moi meme las solutions obtenus pour chaque sous vecteur et je choisis la solution minimale. Pour le moment si je mets

Code :

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

1
2
p=(0:0.001:1)';
PB=(0:0.01:0.05)';

je n'ai plus une mémoire insuffisante, il me reste ensuite de boucler sur deux autres variables. et chercher ensuite une deuxième solution pour

Code :

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

PB=(0.06:0.01:0.11)'

etc.
Concernant le langage C je l'ai utilisé auparavant pour d'autres programmes utilisant des boucles for mais ne contenant pas de calcul matriciel et comparé à son analogue en MATLAB leurs temps d'exécution était beaucoup plus rapide. Ce pendant pour mes programmes utilisant du calcul matriciel j'utilise MATLAB car j'ai des multidimensional arrays and cell arrays.
Si vous avez d'autres solutions faites les moi parvenir.