Discussion :

[Howyd]

Bonjour,
je cherche à créer la fonction en bas de la page 12 (Thiele).
http://physinfo.org/Acc_Conv/Acc_Conv_Part3.pdf

J'ai tous mes coefficients rho (doubles) stockés dans la première ligne d'une matrice carrée rho[] de taille n. Mes z_i sont dans un vecteur x. Je ne vois pas comment générer la fraction, le but final étant de pouvoir la réduire, la tracer, la comparer avec la fonction qu'on lui demande d'approximer, etc...
J'ai essayé avec une boucle for qui part de l'étage le plus en bas et de remonter mais je n'ai pas abouti.
On ne se moque pas:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
g=zeros(n);%vecteur qui contiendrait les étages successifs
syms t;
 
for ii=n:-1:1
    if ii==n
        g(ii)=(t-x(n)); %calcul de l'étage le plus en bas
     elseif ii>1
         g(ii)=(t-x(ii)) / (rho(1,ii+1)-rho(1,ii-1)+g(ii+1));%on remonte
     elseif ii==1
         g(ii)=y(1)+ (t-x(1))/rho(1,2)+g(ii+1);%cas particulier du premier étage
    end
 
end

Si quelqu'un a la solution/piste, merci.
Howyd

[Invité]

Bonjour,

Lorsque tu suis un algorithme, prend les mêmes notations dans ton code, ça aide à suivre
Pour la fonction f, je te conseille d'utiliser une fonction récursive qui prendrait en paramètre ρ, z, les zi, ainsi que le numéro n de l'étage courant.
Le dénominateur devient alors un appel à f à l'étage n+1