Discussion :

[ricololo]

Bonjour,
dans le programme suivant (nommé SearchSun.m), je fais 2 boucles (ih et iphi) dans lesquels j'utilisent la fonction fzero.
Pour certaines valeurs de ih et iphi, une erreur est détectée (au niveau de la fonction fzero). J'aimerai que mon programme ignore cette erreur et passe à la valeur de iphi suivante.
Est-ce possible?
Voici le bout de code:

Code :

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

1
2
3
4
5
6
7
 
for ih=32:nH 
    for iphi=1:nfi
     sqq=@(h) hSo(h,om(ih),sinB(iphi),R,R1,R2,RV);
     hphi = fzero(sqq,H(ih),options);
    end
end

et le message d'erreur:

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
Error in ==> refrac at 16
a = [-2.7622e-12   8.7257e-10  -1.0089e-07   4.9801e-06  ... 
 
??? Output argument "dydz" (and maybe others) not assigned during call to refrac.m
 
Error in ==> diffTheta at 16
[nu dnudz]=refrac(z/1000);
 
Error in ==> AngleOde>@(theta,r)diffTheta(theta,r,K,R) at 43
   [thetas,rs] = ode45(@(theta,r) diffTheta(theta,r,K,R),theta,yi,options);   
 
Error in ==> ode45 at 324
    f(:,2) = feval(odeFcn,t+hA(1),y+f*hB(:,1),odeArgs{:});
 
Error in ==> AngleOde at 43
   [thetas,rs] = ode45(@(theta,r) diffTheta(theta,r,K,R),theta,yi,options);    
 
Error in ==> hSo at 13
refr = AngleOde(h,0);
 
Error in ==> @(h)hSo(h,om(ih),sinB(iphi),R,R1,R2,RV)
 
 
Error in ==> fzero at 369
        a = x - dx;  fa = FunFcn(a,varargin{:});
 
Error in ==> SearchSun at 77
     hphi = fzero(sqq,H(ih),options);

[ricololo]

Je peux utiliser "try ... catch ... end".

[Jerome Briot]

En fait, toute la question est de savoir si ces erreurs sont critiques dans la résolution de ton problème.

Si ce n'est pas le cas, tu peux bien entendu utilise TRY-CATCH-END

Tu peux aussi mettre un avertissement :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
for 
 
    try 
 
        fzero(...)
 
    catch
 
        warning(...)
 
    end
 
end