OpenEdv-开源电子网

 找回密码
 立即注册
正点原子全套STM32/Linux/FPGA开发资料,上千讲STM32视频教程免费下载...
查看: 9589|回复: 4

发个自己写的sin cos函数,精度 10^-11

[复制链接]

49

主题

340

帖子

0

精华

论坛元老

Rank: 8Rank: 8

积分
5246
金钱
5246
注册时间
2012-8-25
在线时间
1024 小时
发表于 2018-4-20 08:06:34 | 显示全部楼层 |阅读模式
[mw_shl_code=c,true]#define  Pi  3.1415926535897932384626433832795
static inline double  power_d( double x )
{
    return ( x * x );
}


double   sin_Taylor( double x)
{
        double  x_2= power_d(x);
//        return (x*(1.0-x_2*(1.0/6-x_2*(1.0/120- x_2*(1.0/5040  )))));
        //return (x*(1.0-x_2*(1.0/6-x_2*(1.0/120- x_2*(1.0/5040-x_2*(1.0/362880  ))))));
        return (x*(1.0-x_2*(1.0/6-x_2*(1.0/120- x_2*(1.0/5040-x_2*(1.0/362880 -x_2/39916800  ))))));       
}

double   cos_Taylor( double x)
{
        double  x_2= power_d(x);
//        return (1.0-x_2*(1.0/2-x_2*(1.0/24- x_2*(1.0/720-x_2*1.0/40320))));       
        return (1.0-x_2*(1.0/2-x_2*(1.0/24- x_2*(1.0/720-x_2*(1.0/40320 -x_2/3628800  )))));       
}


double  sin_q(double x) // 周期是 2*Pi
{
        while(x> Pi) {x-=2*Pi;}
           while(x<-Pi) {x+=2*Pi;}
    int fuhao=1;
        if(x<0)
        {
        x=-x;
        fuhao=-1;
        }
       
        if(x>Pi/2)
        {
        x= Pi-x;
        }
        if(x< Pi*51/180)//Pi/4)
        {
         x= sin_Taylor(x);       
        }
        else
        {       
         x=        cos_Taylor(x-Pi/2);
        }
        if(fuhao==-1)
        {
                x=-x;
        }
        return x;
}

double  cos_q(double x)
{
        return (sin_q(x + Pi/2) );
}
[/mw_shl_code]
正点原子逻辑分析仪DL16劲爆上市
回复

使用道具 举报

32

主题

883

帖子

0

精华

论坛元老

Rank: 8Rank: 8

积分
3990
金钱
3990
注册时间
2015-11-14
在线时间
541 小时
发表于 2018-11-10 14:29:32 | 显示全部楼层
这个运算速度回怎样呢,会不会比math库的快
回复 支持 反对

使用道具 举报

0

主题

1

帖子

0

精华

新手上路

积分
39
金钱
39
注册时间
2016-9-8
在线时间
7 小时
发表于 2018-11-25 21:11:20 | 显示全部楼层
参考下,谢谢
回复 支持 反对

使用道具 举报

2

主题

10

帖子

0

精华

中级会员

Rank: 3Rank: 3

积分
298
金钱
298
注册时间
2021-1-8
在线时间
37 小时
发表于 2021-5-27 11:45:30 | 显示全部楼层
Electronic 发表于 2018-11-10 14:29
这个运算速度回怎样呢,会不会比math库的快

      #include "math.h"
          float ss;
          float angle=0.72;
          ss=sin(angle*125);
          talab[0]=ss;
大佬向你请教一下  为什么我仿真的ss始终等于0啊  
回复 支持 反对

使用道具 举报

1

主题

89

帖子

0

精华

金牌会员

Rank: 6Rank: 6

积分
2915
金钱
2915
注册时间
2019-4-3
在线时间
235 小时
发表于 2021-7-20 12:57:13 | 显示全部楼层
,首先,感谢大佬分享。
小的用c-free 5.0 pro试了下,未能通过编译 。
若要如何,全凭自己!
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则



关闭

原子哥极力推荐上一条 /2 下一条

正点原子公众号

QQ|手机版|OpenEdv-开源电子网 ( 粤ICP备12000418号-1 )

GMT+8, 2024-11-22 17:26

Powered by OpenEdv-开源电子网

© 2001-2030 OpenEdv-开源电子网

快速回复 返回顶部 返回列表