Petit problème pour un script

Tout d'abord bonjour,

Je dois faire un programme sous matlab pour du rayonnement thermique.

Le but du programme est de déterminer des températures.

Mon problème : à chaque boucle j'obtiens un vecteur (1X6) (une colonne de 6 lignes), je souhaiterais stocker ou concaténer chaque vecteur (créés par la boucle) dans une matrice finale.

Ainsi je verrai apparaître dans chaque colonne (de cette matrice) les températures à chaque boucle.

Comment faire ? J'ai des soucis d'indexation !!!

Je sais que ce problème est plutôt fastoche pour vous qui êtes des cracks !!!

En tout cas je vous remercie par avance de votre aide.:D

Montre nous ce que tu as essayé de programmer...

comme avez dit Mr dutmatlab on doit voir un bout de ton travail, mais en tous cas ce code te permet de crée ta matrice.

Code:

---

1 2 3 4 5 6 7 8

res=[]; for p=1:7 for m=1:6 temp(m)=rand(1) ;%calcul de la température end res=[res temp(:)]; end

---

Salut,

Je ne suis pas d'accord avec le code de soft001 car je ne sais pas ce qu'il y aura dans la première colonne de "res" et que la boucle sur "p" ne sert à rien à part une mini perte de temps puisque la ligne "res=..." pourrait très bien être mise dans la boucle sur "m"...

Perso je ferais juste ceci :

Code:

---

1 2 3

for m=1:N temp(:,m)=calcul_température; end

---

C'est à dire que tu vas pouvoir faire N fois le "calcul_température" et tu stockes chaque "calcul_température" dans la colonne "m" de "temp".

A la fin de ton code, tu auras une matrice/tableau "temp"...

Maintenant, je pense que ça dépend surtout de ton "calcul_température", c'est pourquoi le grand Dut te demande de nous montrer ce que tu as fait...

A bientôt

Christophe

Un grand merci pour vos contributions !!!!!!!!!!! :D

donc voici mon petit programme, je vous ai détaillé les différentes variables.

C'est un problème de rayonnement thermique : on considère un cube, un flux (Phi=1320 Watts) arrive avec un angle incident (i) sur deux surfaces du cube appelées (S1) et (S2). Ainsi ce flux se décompose comme :
-sur (S1) : phi*cos(i)
-sur(S2) : phi*sin(i)

Ensuite le programme utilise des formules (méthode de Gebhart : "résolution du pb inverse" (utilisés pour les échanges dans le satellite)).

Bref beaucoup de blabla....

Une fois que ces calculs sont effectués, j'obtiens pour chaque boucle (c'est à dire pour une variation d'angle d'incidence) une température (sous forme de vecteur).

Le truc c'est que je souhaiterais mettre chacun de ces vecteurs (dans le programme c'est T1) dans une matrice=CONCATéNER dans une matrice.

PROGRAMME :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52

%Rayonnement sur un satellite   %Rayonnement sur 2 surfaces. % les données de problème   for i=0:(pi/8):pi/2 eps=[0.8 0.8 0.8 0.8 0.8 0.8]; %epsilon : coeeficient d'émissivité. rho=1-eps; %rho coefficient de reflexion flux=[1320*cos(i) 1320*sin(i) 0 0 0 0]; %flux imposé sur la surface avec un angle i. sigma=5.67e-8; % matrice de forme for i=1:6 for j=1:6 if (i==j) F(i,j)=0; else F(i,j)=0.2; end end end   % calcul Bij for i=1:6 for j=1:6 if (i==j) a(i,i)=1-rho(i)*F(i,i); c(i,i)=eps(i)*F(i,i); else a(i,j)=-rho(j)*F(i,j); c(i,j)=eps(j)*F(j,i); end end end % inversion de la matrice B=inv(a)*c; % calcul des temperatures %element du veteur E for j=1:6 E(j)=flux(j)/sigma; end % elements de la matrice D for i=1:6 for j=1:6 if (i==j) d(i,i)=1+eps(i)-eps(i)*B(i,i); else d(i,j)=-eps(j)*B(i,j); end end end % calcul de T1 T=(inv(d)*E'); Tcube=T.^(1/4) %Vecteur où la ième ligne est la température de (Si) end %T_finale=

---

Voilà, j'espère que vous me comprennez.

Si un point n'est pas assez clair, n'hésitez pas !!!!!

Citation:

---

Envoyé par christophe_halgand

Je ne suis pas d'accord avec le code de soft001

---

tu a raison j’avais la possibilité de pas utiliser les boucles imbriqués et de mieux optimiser mon code, mais concernant la première colonne de res elle doit contenir les températures déjà calculées.

Salut,

Alerte ROUGE, tu utilises plusieurs fois "i" comme variable... Ce n'est pas bien du tout... Déjà que si tu fait "help i", tu verras que matlab ne voit pas les choses comme toi...

Donc évites les "i" pour les boucles et de plus évites d'utiliser une variable dans une boucle sur cette même variable...

Je pense que tu as juste à faire ce genre de chose (je garde "i" pour que tu comprennes par rapport à ce que tu as proposé au dessus) :

Code:

---

1 2 3 4 5 6 7

u=1; for i=0:(pi/8):pi/2 blablabla Tcube(:,u)=T^1/4; u=u+1; end Tfinal=Tcube

---

Pour la forme de la matrice, je te propose cette ligne de commande :

Code:

---

F=0.2*ones(6)-0.2*eye(6);

---

Pour l'élément du vecteur E, je te propose cette ligne de commande :

Code:

---

E=flux(1:6)/sigma

---

Voilà, je pense que ça devrait t'aider...

A bientôt

Christophe

Merci beaucoup toulousain, je suis aussi du coin !!!!

je teste ça tout de suite et je te tiendrai au courant

ps: Allez le ST :D:D

super j'ai suivi vos indications (surtout celles christophe_halgand) et ça marche très bien !!! :D

je vous remercie pour vos contributions :king::king: