OpenEdv-开源电子网

 找回密码
 立即注册

扫一扫,访问微社区

正点原子全套STM32开发资料,上千讲STM32视频教程,RT1052教程免费下载啦...

查看: 560|回复: 0

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

[复制链接]

  离线 

22

主题

208

帖子

0

精华

金牌会员

Rank: 6Rank: 6

积分
1790
金钱
1790
注册时间
2012-8-25
在线时间
455 小时
发表于 2018-4-20 08:06:34 | 显示全部楼层 |阅读模式
正点原子公众号
[C] 纯文本查看 复制代码
#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) );
}

回复

使用道具 举报

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

本版积分规则




关闭

必看:"原子哥”力荐上一条 /1 下一条

QQ|联系我们|手机版|官方淘宝店|微信公众平台|OpenEdv-开源电子网 ( 粤ICP备12000418号-1 )

GMT+8, 2018-8-22 13:41

Powered by OpenEdv-开源电子网

© 2001-2030 OpenEdv-开源电子网

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