百科问答小站 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++中,“?:”这个符号叫什么名字? 
  C语言中逗号表达式存在的意义是什么? 
  视频RAM与IO内存映射的问题,大家可以指导一下吗? 
  C 语言用 换行后就无法再回到上一行了吗? 
  自己开发的Windows7 64位驱动程序需要签名,怎么解决? 
  计算机存储浮点数的时候,为什么不使用十进制来存储呢? 
  为什么很多新型编程语言都抛弃了 C 语言风格的 for 语句? 
  DOS 界面是如何制作的? 
  既然有 memcpy_s 这种安全实现,为啥不禁用危险的 memcpy 或者更新 memcpy 源码? 

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





© 2024-12-26 - tinynew.org. All Rights Reserved.
© 2024-12-26 - tinynew.org. 保留所有权利