why matlab Make a mistake in calculating sin and cos while teta in wiling to zero ?
5 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
shahin hashemi
el 26 de En. de 2018
Comentada: Walter Roberson
el 26 de En. de 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 comentarios
Respuesta aceptada
Walter Roberson
el 26 de En. de 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 comentarios
Walter Roberson
el 26 de En. de 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.)
Ver también
Categorías
Más información sobre Symbolic Math Toolbox en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!