why matlab Make a mistake in calculating sin and cos while teta in wiling to zero ?
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
shahin hashemi
le 26 Jan 2018
Commenté : Walter Roberson
le 26 Jan 2018
dear all
as it is obvious base on the Taylor series or l'hopital's rule , answer the below equation is 1/2
Taylor series :
(1-cos(teta))/(teta^2)= 1/2! - (teata^2)/4!
and while teta wiling to zero others terms become zero and the answer is 1/2
l'hopital's rule
(1-cos(teta))/(teta^2)= sin(teta)/(2*teta) =(1/2)*(cos(teta))
and again while teta wiling to zero answer is 1/2
but matlab act different and it gives right answer when you assume that teta is 0.01 or 0.001 but when you assume that teta is 0.00001 matlab give zero answer to this equation
is it my fault or something wrong with matlab and i really appreciated if some one can help me and if let me know if there is code for my problem
0 commentaires
Réponse acceptée
Walter Roberson
le 26 Jan 2018
By teta = .0000001 then cos(teta) becomes 1 to the limit of double precision, considering round-off error. Then 1-cos(teta) is 0, so you get 0 as the result.
This is a limitation in numeric processing with finite precision. If double precision had (for example) 80 bits instead of 53, then the problem would still occur, just with a smaller x.
2 commentaires
Walter Roberson
le 26 Jan 2018
The Symbolic toolbox is the main supported tool for extended precision. You could also look at https://www.mathworks.com/matlabcentral/fileexchange/6446-multiple-precision-toolbox-for-matlab (read the comments on the submission for some corrections.)
Voir également
Catégories
En savoir plus sur Logical dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!