请教，哪儿有sin，cos函数的实现

scutan

请教一下，sin，cos函数的实现。在网上搜了一下没有搜到相关的源码。谢谢

reiase

泰勒展开或者贝赛耳曲线？

aaaaa5aa

/******************************************************************
*
* File Name:       pw_sin.h
* Description:     The sin() function map. Using it to calculate a
*                                         angle's sine value and cosine value.
*                                         e.g.
*                                                 angle d=30 degrees,
*                                                 sin(d) = (sin_map[30] >> SIN_MAP_SHIFT)
*                                                 result: sin(30) = 0.5
*                  
* Author :         Jiang Peifu
* Date:            2009.04.27
*/
#ifndef PW_SIN_H
#define PW_SIN_H
/* -------------------------------------------------------------------------- */
#define SIN_MAP_SHIFT 15
#define DEGREE_MIN 0
#define DEGREE_MAX 90
#define SIN_MAP_SIZE (DEGREE_MAX - DEGREE_MIN + 1)
#define rsin(_r_, _d_) ((_r_ * sin_map[_d_] + 0x4000) >> SIN_MAP_SHIFT)
#define rcos(_r_, _d_) ((_r_ * sin_map[DEGREE_MAX-_d_] + 0x4000) >> SIN_MAP_SHIFT)
int sin_map[SIN_MAP_SIZE] = {
    0,
    572,
    1144,
    1715,
    2286,
    2856,
    3425,
    3993,
    4560,
    5126,
    5690,
    6252,
    6813,
    7371,
    7927,
    8481,
    9032,
    9580,
    10126,
    10668,
    11207,
    11743,
    12275,
    12803,
    13328,
    13848,
    14365,
    14876,
    15384,
    15886,
    16384,
    16877,
    17364,
    17847,
    18324,
    18795,
    19261,
    19720,
    20174,
    20622,
    21063,
    21498,
    21926,
    22348,
    22763,
    23170,
    23571,
    23965,
    24351,
    24730,
    25102,
    25466,
    25822,
    26170,
    26510,
    26842,
    27166,
    27482,
    27789,
    28088,
    28378,
    28660,
    28932,
    29197,
    29452,
    29698,
    29935,
    30163,
    30382,
    30592,
    30792,
    30983,
    31164,
    31336,
    31499,
    31651,
    31795,
    31928,
    32052,
    32166,
    32270,
    32365,
    32449,
    32524,
    32588,
    32643,
    32688,
    32723,
    32748,
    32763,
    32768,
};
#endif /* PW_SIN_H */

看下这个

aaaaa5aa

文件i_sin.h中，保存了0-90度角的sin值向左位移15后的结果。对每个角度的三角函数做运算时，以角度为小标在sin_map表中查找相应的sin值（位移后的结果）。之所以要对sin值向左位移，是保证可以用整数来表达sin值。位移的大小（本文中为15），标称了sin值取值的精度，越大，sin值取值越精确。

VIP_fuck

libc里边没有吗？

aaaaa5aa

那个下面没有，其实可以去math目录下看看

aaaaa5aa

我记得南开100题中有一道，也许对你有用

题目21：下列程序的功能是：利用以下所示的简单迭代方法求方程：cos(x)-x=0的一个实根。xn+1=cos(xn)
迭代步骤如下：
（1）取X1初值为0.0；
（2）X0=X1，把X1的值赋给X0；
（3）X1=COS（X0），求出一个新的X1；
（4）若X0-X1绝对值小0.000001,执行步骤(5),否则执行步骤(2);
(5)所求X1就是方程cos(X)-X=0的一个实根,作为函数值返回。

scutan

原帖由 aaaaa5aa 于 2009-7-9 13:23 发表
文件i_sin.h中，保存了0-90度角的sin值向左位移15后的结果。对每个角度的三角函数做运算时，以角度为小标在sin_map表中查找相应的sin值（位移后的结果）。之所以要对sin值向左位移，是保证可以用整数来表达sin值 ...

先谢谢，有没有那种不用查表的方式呢？而且精度要比较高点的。

因为我现在用的环境自带的sin等函数都只是单精度的，使得最终结果差别较，所以我想能否找一个源码来自己去调用。

cheng_lai_shun

LZ可以自己写个！
先把sin、cos用级数表示出来！

ideawu

用google codesearch: http://www.google.com/codesearch ... =rr&cs_r=lang:c
各种语言的版本: http://www.google.com/codesearch ... sin&sbtn=Search