Discussion :

[elglantosimpatico]

Bonjour,

J'ai une question toute simple dont je n'ai pas trouvé la réponse sur la FAQ.
Je fais un calcul sur matlab or, selon les paramètres d'entrée, je peux avoir de gros problème numériques.

Voici le calcul :

Code :

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

1
2
b = 2*n^2;
resultat = b*(-pi*csc(l*pi)/(b*gamma(-l))+(b^(-1+l/2)*exp(-b/2)*mfun('WhittakerM',l/2,(l+1)/2,b))/(-1 - l));

je plotte ça pour l allant de 1 à 10 sans problème pour n=2. Cependant, lorsque j'augmente n, des erreurs calcul apparaissent.

Je rentre la même fonction sous Maple et les mêmes erreurs apparaissent. Cependant, si dans Maple je rentre :

Code :

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

Digits:=30:

alors les erreurs disparaissent.
Il y a-t-il un moyen de faire quelque chose de similaire sous matlab?

[soft001]

Tu peux regarder la fonction :

Code :

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

format

[Jean Dumoncel]

Bonjour,

format modifie le format d'affichage des données, il ne modifie pas la précision utilisée pour les calculs :

format -
Set display format for output

Maple est un logiciel de calcul formel, à la différence de Matlab qui fait du calcul numérique. Je ne sais pas si il existe une solution en utilisant les fonctions de base de Matlab, mais tu peux essayer la toolbox de calcul symbolique de matlab.

[soft001]

Bonjour,

format modifie le format d'affichage des données, il ne modifie pas la précision utilisée pour les calculs :


Maple est un logiciel de calcul formel, à la différence de Matlab qui fait du calcul numérique. Je ne sais pas si il existe une solution en utilisant les fonctions de base de Matlab, mais tu peux essayer la toolbox de calcul symbolique de matlab.

Oui je sais bien, mais personnellement j'ai supposé que l'erreur est d'ordre de 10^-4, alors en regardant le résultat à la 4ème chiffre après la virgule (et vu que le format short est mis par défaut) ce terme là sera arrondit et par suite on va penser qu'on a une erreur

[Jean Dumoncel]

Effectivement le posteur n'a pas précisé l'ordre de grandeur de l'erreur...

[elglantosimpatico]

Effectivement le posteur n'a pas précisé l'ordre de grandeur de l'erreur...

Comment puis-je savoir cela?