Discussion :

[JayyyJayyy]

Bonjour,

Je viens d'apprendre récemment qu'il était préférable d'utiliser la vectorisation par rapport aux boucles for...
Néanmoins, j'aurais besoin d'un peu d'aide pour les boucles (imbriquées) suivantes :

Code :

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

1
2
3
4
5
6
7
8
9
10
for k=1:N	
    Bprod =prod(B(1:k-1));
    borneInf = 0;
    for i=1:k-1
        borneInf += x(i) * prod(A(i:k-1))
    end
    borneSup = borneInf + x(k);
    integrale = quad(@(x)wblpdf(x,scale,shape)/(1-wblcdf(x,scale,shape)), borneInf, borneSup);
    result += Bprod * integrale
end

Merci

[duf42]

Bonjour,

Une petite remarque: l'opérateur += n'existe pas en MATLAB donc ton code ne marche pas

Duf

[Invité]

Bonjour,

Une autre remarque, pour k = 1:

>> B(1:k-1)
ans =
   Empty matrix: 0-by-1

Mais prod(...) renvoie 1 tout de même
Sinon je doute que l'on puisse vectoriser ces boucles.

[JayyyJayyy]

Bonsoir,

J'ai effectivement posté dans la catégorie MATLAB mais pour être précis, j'utilise Octave qui offre les mêmes fonctionnalités que MATLAB avec (à mon avis) quelques avantages supplémentaires (dont notamment le "+=").

En tout cas, merci !