avr中怎么运用atan(x,y) 函数
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/25 03:37:39
avr中怎么运用atan(x,y) 函数
avr中怎么运用atan(x,y) 函数
avr中怎么运用atan(x,y) 函数
# define PiOv2 (3.14159f / 2.0f)
# define SignBit(f)\x05\x05((*(const unsigned long *)&(f)) >> 31)
float Atan_16bits_Pre( float y,float x )
{
\x05float a,s;
\x05if ( fabs( y ) > fabs( x ) )
\x05{
\x05\x05a = x / y;
\x05\x05s = a * a;
\x05\x05s = - ( ( ( ( ( ( ( ( ( 0.0028662257f * s - 0.0161657367f ) * s + 0.0429096138f ) * s - 0.0752896400f )
\x05\x05\x05* s + 0.1065626393f ) * s - 0.1420889944f ) * s + 0.1999355085f ) * s - 0.3333314528f ) * s ) + 1.0f ) * a;
\x05\x05if ( SignBit( a ) ) {
\x05\x05\x05return s - PiOv2;
\x05\x05} else {
\x05\x05\x05return s + PiOv2;
\x05\x05}
\x05}
\x05else
\x05{
\x05\x05a = y / x;
\x05\x05s = a * a;
\x05\x05return ( ( ( ( ( ( ( ( ( 0.0028662257f * s - 0.0161657367f ) * s + 0.0429096138f ) * s - 0.0752896400f )
\x05\x05\x05* s + 0.1065626393f ) * s - 0.1420889944f ) * s + 0.1999355085f ) * s - 0.3333314528f ) * s ) + 1.0f ) * a;
\x05}
}
atan(x,y)需要CPU浮点指令级支持,AVR中一般没有专用浮点单元,所以atan也是软件实现的.如果你在AVR的库中一时没有找到atan,可以用上面这段程序代替,注意,上面程序只提供了两字节精度支持,大概小数点后3位左右.
上述程序来至id的4号引擎,基本上都是大神写的,你很难找到一段效率比它还高的代码,上述代码在一定情况下和浮点指令的速度差不多.
你看看 math.h