Discussion :

[bichou]

bonjour,
je dois programmer une relation de recurrence avec matlab afin de trouver alpha:

Code :

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

1
2
3
for k=2:6 
alpha(k+1) = alpha(k)-(alpha(k)-alpha(k-1))*h(alpha(k))/(h(alpha(k))-h(alpha(k-1))) ; 
end

ou h est la fonction definie par:

Code :

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

h=1+psi(1/alpha)/alpha-sum(V(i))^(alpha)*log(sum(V(i))/sum(V(i))^alpha+log(alpha/N*sum(V(i))^(alpha)/alpha))

mais je ne sais pas comment rentrer cette fonction, vu qu'on me dit qu'il ne connait pas alpha...ce que justement je veux calculer...

pourriez vous m'aider?
merci

[Jerome Briot]

Euhh... tu n'aurais pas besoin d'une première valeur pour alpha ? alpha(1) doit être connu non ?

[bichou]

si j'ai mes deux premieres valeurs.je connais alpha(1) et alpha(2)
mais je ne les ai definies qu'au moment de faire ma recurrence..
j'ai d'abord rentre ma fonction h et c'est la qu'il me dit qu'il ne connait pas alpha et apres j'ai fait ma recurrence.

[Jerome Briot]

Le plus simple ici est d'intégrer h dans la formule... peut importe la longueur finale de celle-ci au final.