cos block error (simulink)

15 次查看(过去 30 天)
기범
기범 2023-1-14
评论: Paul 2023-1-17
why is this happen?
cos(pi/2) = 6.123e-17??

采纳的回答

John D'Errico
John D'Errico 2023-1-14
编辑:John D'Errico 2023-1-14
Welcome to the wonderful, wacky world of floating point arithmetic.
Is the number
format long g
pi/2
ans =
1.5707963267949
EXACTLY pi/2?
cos(pi/2)
ans =
6.12323399573677e-17
Of course not. Pi is an irrational number (even subtly worse, transcendental), so it cannot ever be represented exactly in double precision arithmetic, so not in any finite number of digits. And that means pi/2 as stored is just a tiny bit off from the exact value of pi/2. Consequently, cos(pi/2) is not EXACTLY zero. Close. To within floating point trash. But not exactly so.
If you want MATLAB to work in exact multiples of pi, you can do so, but that requires using the symbolic toolbox.
sym(pi/2)
ans = 
cos(sym(pi/2))
ans = 
0
  5 个评论
기범
기범 2023-1-15
Im using cos for my Quaternion,
so I need an exact value.
6.123e-17 will make error for my Quaternion.
I dont want my system slow down as you said.
Can i set 6.123e-17 to Zero in simulink?
ex) below 1*e-10 to Zero
John D'Errico
John D'Errico 2023-1-16
So you would want a rather poor approximation to the cosine function for small values? That is, you would rather see an error of 1e-10, than an error of 1e-17? Why do I feel that is just a really, really, seriously bad idea? You want an exact value some of the time, but far more often, you are willing to accept a really poor approximation?
You could test to see if the absolute value is less then 1e-16, and if so, set that to zero. At last then you are not making a rather large error.

请先登录,再进行评论。

更多回答(1 个)

Paul
Paul 2023-1-14
Base Simulink doesn't have a block that implements functions cosd or cospi. Maybe some other toolbox does. If not, both of those functions support code generation so you can use the Matlab Function block and call either (or both of them) from there.
  5 个评论
기범
기범 2023-1-17
For your recommend,
Im using MATLAB FUNCTION.
error : 1. The 'sym' class does not support code generation.
2. defined function or variable 'phirad'. The first assignment to a local variable
determines its class.
can u help me?
Paul
Paul 2023-1-17
I never said to use syms. I did say to use cosd, etc. if you want have all your angles defined in degrees. So, if theta, phi, and psi are all in degrees on input, then we have
x = r.*sind(phi).*cosd(theta)
and similar for y and z

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Discontinuities 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by