Discussion :

[JeaMachine]

Bonsoir,

Voilà, j'aimerais exposer mon problème à un public un peu plus compétent. J'ai préalablement cherché une solution à ce dernier soucis sur la toile sans trouver de solution.
J'ai édité une fonction permettant de calculer le nombre de racines évidentes d'un polynôme d'un ordre quelconque. Dans cette fonction, j'utilise "polyder" pour dériver le polynôme. Or, lorsque j'exécute ma fonction dans Matlab, j'obtiens rapidement l'erreur suivante :

??? Maximum recursion limit of 500 reached. Use set(0,'RecursionLimit',N)
to change the limit. Be aware that exceeding your available stack space can
crash MATLAB and/or your computer.

Error in ==> polyder

Apparemment, l'espace en mémoire vive que j'autorise à Matlab d'utiliser est rapidement saturé et c'est ce qui provoque l'erreur. Néanmoins, j'ai essayé d'augmenter la limite de récursion à la limite du plantage de Matlab sans succès.
J'aimerais donc savoir s'il existe un moyen pour utiliser ma fonction dans des conditions optimales.
Merci de vos réponses.
Cordialement.

JeaMachine.

[Jerome Briot]

Il faudrait que tu nous montres ton code

Je ne comprends pas bien pourquoi tu appelles POLYDER récursivement...

[sigma10]

Bonsoir

voila j'ai le même problème :

Maximum recursion limit of 500 reached. Use set(0,'RecursionLimit',N)
to change the limit.  Be aware that exceeding your available stack space can
crash MATLAB and/or your computer.

Error in ==> ondelettes 

mais où je vais poser mon code car il contient des fonctions ; comment faire ??

j'attends votre réponse

[Invité]

Bonjour,

Et si tu montrais les différentes fonctions? Je ne vois pas le problème
As-tu dans un premier temps essayé d'appliquer ce que te dis le message d'erreur?

[shaiHulud]

Bonjour,

J'imagine que polyeval dérive le polynôme P, et s'appele récursivement sur le polynôme dérivé pour calculer son tableau de variation. polyeval est donc censé faire exactement degrès(P) appels récursifs.

Donc déja, si tu obtiens ton Maximum recursion limit of 500 reached sur un polynome de degrès< 500, c'est qu'il doit y avoir une erreur dans le code. Tu peux la détecter en vérifiant:
- que le degrès de P diminue strictement à chaque appel récursif.
- que la fonction ne fait pas d'appel récursifs sur (au moins) un polynome de degrès -infini, 0 ou 1.

Dans un cadre plus général (pour sigma10 par exemple), il est possible de dérécursifier les appels récursifs terminaux. Au lieu d'écrire une fonction du type

Code :

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

1
2
3
4
5
6
7
8
9
 
function x= recurs1(n)
% Calcule la suite récursive x(1)= 0.5 et x(n+1) = x(n)^2
% Utilise récursivité terminale
if n> 1
    x= recurs1(n-1)^2;
else
    x= 1;
end

Il suffit de faire (ce code est illustratif et pourrait etre vectorisé/optimisé)

Code :

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

1
2
3
4
5
6
7
8
9
 
function x= recurs1(n)
% Calcule la suite récursive x(1)= 0.5 et x(n+1) = x(n)^2
% Sans utiliser récursivité terminale
allx= 1;
for k=2:n
     allx(k)= allx(k-1)^2;
end
x= allx(end);

Le résultat est le même mais la récursion a disparue (et il serait même possible de renvoyer toutes les valeurs d'un coup, et non juste la dernière)

A priori, le code d'estimation par ondelette est récursif terminal (ou équivalent) car il suffit d'estimer à fréquence/2 le résidu de l'estimation des premiers niveaux. Cette solution doit donc marcher pour sigma10.