百科问答小站 logo
百科问答小站 font logo



求教,atan2函数有什么简易的替代品? 第1页

  

user avatar   zhang-hao-72 网友的相关建议: 
      

自己答一个:后来考虑了一下其实有asin就行了,用了查表+线性插值的方案。查找表只用了0到90度,间隔5.625度的17个点, 再乘以0xffff,如下:

       const unsigned short sintab[] = {     0x0000, 0x1918, 0x31f1, 0x4a50, 0x61f7, 0x78ad, 0x8e39, 0xa267, 0xb504, 0xc5e3,      0xd4da, 0xe1c5, 0xec82, 0xf4f9, 0xfb14, 0xfec3, 0xffff };     

然后查表+线性插值,输入限定在0~0xffff。

       int int_asin(int x) {     int ret, sign=1;     if(x<0) {         x=-x;         sign=-1;        // 处理负数     }     for(int i = 0; i < sizeof(sintab) / sizeof(sintab[0]); i++) {         if(sintab[i] <= x && x <= sintab[i + 1]) {             ret =  (((x - sintab[i]) * 90 / (sintab[i+1] - sintab[i]) + i * 90) ) >> 4;             break;         }     }     return ret * sign; }     

实测算出来的角度值和math库的asin基本相符,个别情况差1度。

       int main(void) {     for(int i = -90; i <= 90; i++) {         int y = (sin(i * 3.141593 / 180.0) * 65535.0);         int x1 = int_asin(y);          float x2 = asin(y / 65536.0) * 180.0 / 3.141593;         printf("%d %d %.3f %d
", i, x1, round(x2), x1-round(x2));     }      return 0; }     




  

相关话题

  C 语言中的 double 类型所能表示的数的范围为什么这么大,不是说只占32(或64)位吗? 
  请问学 C 有必要死磕指针吗? 
  这个图左边这个程序的输出值不太懂为啥是10? 
  Linux 内核的 C 代码风格是怎样的? 
  求教,atan2函数有什么简易的替代品? 
  请问如何知道I2C线上面有没有信息传输? 
  求十亿内所有质数的和,怎么做最快? 
  一条C语言语句不一定是原子操作,但是一个汇编指令是原子操作吗? 
  DOS 界面是如何制作的? 
  C 语言自带函数返回值为指针类型的数组为什么不需要释放内存? 

前一个讨论
c语言如何定义没有返回值的main函数?
下一个讨论
假设可控核聚变实现,地球是否会有危机,比如全球变暖?





© 2024-11-25 - tinynew.org. All Rights Reserved.
© 2024-11-25 - tinynew.org. 保留所有权利