Discussion :

[momo70]

je suis stagiarie,en master , je bosse pour faire un progiciel de la gestion de la serre agricole.
le seul language de programmation qui m'est disponible et que j'ai des notions c'est Matlab:
mon probleme; se resume comme suit:
pour caracteriser le climat exterieur, il me fallait programmer:
mes variables sont:
la latitue LA:
un tableau des mois de janvier a decembre donc il varie de 1a 12 ,T1=[1:12]
chaque mois lui est affecter son nbre de jour cad
T2=[31,28,31,30,31,30,31,31,30,31,30,31],
je veux faire une boucle qui calcul pour chaque jour de mois .
la boucle que j'ai fait ne me donne pas pour tous les jours de l'année.
voila mon la boucle

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
16
17
18
19
20
21
22
23
24
25
26
27
clear;
la=23,45 
% 'declinaison'
mois = [1:12];
jour = [ 31,28,31,30,31,30,31,31,30,31,30,31]
for =1:12,
for k=1:12,
for m=1: jour(k);
% 'declinaison'
% Dl=23.45*cos((30*mois+jour-202)*pi/180)*pi/180,(la formule declinaison)
 
Dl=23.45*cos((30*l+k-202)*pi/180)*pi/180,
 
%’la tangente’
th = (( 1 - ( - tan( Dl) * tan ( la * 3.1416 / 180 )) ^2) ^.5) / (- tan ( la * 3.1416 / 180 )* tan ( Dl)),
%’angle horaire’
ah=atan(th);
%’durée de jour’
if ah==0,d1=12,end
if ah <0, d1=2*(pi+ah)*12/pi
else d1=2*(ah*12/pi),
%’durée de nuit
d0=24-d1,disp([Dl,th,ah,d1,d0]);
 
end 
end
end

j'attend impatiement votre aide pour me debloquer et avancer
recevez tous mes remerciement

[Jerome Briot]

Déjà, il faut programmer plus proprement. Il faut utiliser l'indentation pour rendre le code lisible. J'ai déjà remarqué deux problèmes (en rouge) :

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
16
17
18
19
20
21
22
23
24
25
26
27
28
29
clear;
la=23,45 
% 'declinaison'
mois = [1:12];
jour = [ 31,28,31,30,31,30,31,31,30,31,30,31];
for l=1:12,
   for k=1:12,
      for m=1: jour(k);
         % 'declinaison'
         % Dl=23.45*cos((30*mois+jour-202)*pi/180)*pi/180,(la formule declinaison)
         Dl=23.45*cos((30*l+k-202)*pi/180)*pi/180;
         %’la tangente’
         th = (( 1 - ( - tan( Dl) * tan ( la * 3.1416 / 180 )) ^2) ^.5) / (- tan ( la * 3.1416 / 180 )* tan ( Dl));
         %’angle horaire’
         ah=atan(th);
         %’durée de jour’
         if ah==0
            d1=12;
         end
         if ah <0
            d1=2*(pi+ah)*12/pi;
         else 
            d1=2*(ah*12/pi);
            %’durée de nuit
            d0=24-d1,disp([Dl,th,ah,d1,d0]);			
         end         
      end        
   end    
end

[momo70]

vous avez raison sur la probleme de l'indentation, mais je l'ai ecris sur la messagerie, en effet je n'ai pas respecter cette regle.
mon probleme sur matlab je respecte cette regle d'indentation, mais je me demande sur la formulation de la boucle qui ne me calcule pas pour toutes les valeurs introduites en l'occurence que les mois varie de 1 a 12 et les jours de 1 a 365.
car je dois avoir ces fonctions pour chaque jour.
Dl=23.45*cos((30*mois+jour-202)*pi/180)*pi/180
il me calcule juste au 12eme jour du mois
or que je veux que le mois varie selon son nombre de jour.
exemple pour janvier mois=1, je dois avoir 31 valeurs de ces differentes fonction
pour fevrier le mois c le numero2 je dois avoir 28 valeurs etc.
j'espere que vous avez saisi le probleme.

toutes mes remerciements

[Jerome Briot]

Je crois que j'y suis :

Code :

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

1
2
3
4
5
6
7
8
mois = [1:12];
jour = [31,28,31,30,31,30,31,31,30,31,30,31]
 
for m=1:numel(mois)
   for n=1:jour(m)
      Dl=23.45*cos((30*m+n-202)*pi/180)*pi/180;
   end
end