三角函数sin cos模拟测试

背景： 裸奔程序时，有些函数库报错，干脆自己实现库函数。

优点：速度快，适合嵌入式、图像处理方面。

缺点：不适用于误差要求比较高的场合。

源码：

#include <stdio.h>
#include <math.h>

#define abs(x)	(x < 0 ? (0 - x) : x)
#define Pi 3.1415926535898
#define DEC (Pi/180)

float sine(float x) 

{

	 const float B = 1.2732395447; 
	 const float C = -0.4052847346;
	 const float P = 0.225;//0.2310792853;//0.225;
	 float y = 0;

	 while(abs(x) > Pi)
	 {
		if(x > 0)
			x = -(x - Pi);
		else
			x = -(x + Pi);
	 }
	
	
	 y = B * x + C * x * abs(x);
	 y = P * (y * abs(y) - y) + y;
	 return (float)y;
 }


float cose(float x) 
{
  const float Q = 1.5707963268;
  const float PI =3.1415926536;
  x += Q; 

 return( sine(x));
}

int main(void)
{

	int i;
	for(i = -360; i <= 360; i+=10)
	{
	//	printf("sin(%d):%8f\tsine(%d):%8f\n", i, sin(i * DEC), i, sine(i * DEC));
		
		printf("cos(%d):%8f\tcose(%d):%8f\n", i, cos(i * DEC), i, cose(i * DEC));
	
	}

}

结果：

$ ./test
cos(-360):1.000000      cose(-360):1.000000
cos(-350):0.984808      cose(-350):0.984911
cos(-340):0.939693      cose(-340):0.940055
cos(-330):0.866025      cose(-330):0.866666
cos(-320):0.766044      cose(-320):0.766804
cos(-310):0.642788      cose(-310):0.643347
cos(-300):0.500000      cose(-300):0.500000
cos(-290):0.342020      cose(-290):0.341289
cos(-280):0.173648      cose(-280):0.172565
cos(-270):0.000000      cose(-270):-0.000000
cos(-260):-0.173648     cose(-260):-0.172565
cos(-250):-0.342020     cose(-250):-0.341289
cos(-240):-0.500000     cose(-240):-0.500000
cos(-230):-0.642788     cose(-230):-0.643347
cos(-220):-0.766044     cose(-220):-0.766804
cos(-210):-0.866025     cose(-210):-0.866667
cos(-200):-0.939693     cose(-200):-0.940055
cos(-190):-0.984808     cose(-190):-0.984911
cos(-180):-1.000000     cose(-180):-1.000000
cos(-170):-0.984808     cose(-170):-0.984911
cos(-160):-0.939693     cose(-160):-0.940055
cos(-150):-0.866025     cose(-150):-0.866667
cos(-140):-0.766044     cose(-140):-0.766804
cos(-130):-0.642788     cose(-130):-0.643347
cos(-120):-0.500000     cose(-120):-0.500000
cos(-110):-0.342020     cose(-110):-0.341289
cos(-100):-0.173648     cose(-100):-0.172565
cos(-90):-0.000000      cose(-90):0.000000
cos(-80):0.173648       cose(-80):0.172565
cos(-70):0.342020       cose(-70):0.341289
cos(-60):0.500000       cose(-60):0.500000
cos(-50):0.642788       cose(-50):0.643347
cos(-40):0.766044       cose(-40):0.766804
cos(-30):0.866025       cose(-30):0.866667
cos(-20):0.939693       cose(-20):0.940055
cos(-10):0.984808       cose(-10):0.984911
cos(0):1.000000 cose(0):1.000000
cos(10):0.984808        cose(10):0.984911
cos(20):0.939693        cose(20):0.940055
cos(30):0.866025        cose(30):0.866667
cos(40):0.766044        cose(40):0.766804
cos(50):0.642788        cose(50):0.643347
cos(60):0.500000        cose(60):0.500000
cos(70):0.342020        cose(70):0.341289
cos(80):0.173648        cose(80):0.172565
cos(90):-0.000000       cose(90):-0.000000
cos(100):-0.173648      cose(100):-0.172565
cos(110):-0.342020      cose(110):-0.341289
cos(120):-0.500000      cose(120):-0.500000
cos(130):-0.642788      cose(130):-0.643347
cos(140):-0.766044      cose(140):-0.766804
cos(150):-0.866025      cose(150):-0.866667
cos(160):-0.939693      cose(160):-0.940055
cos(170):-0.984808      cose(170):-0.984911
cos(180):-1.000000      cose(180):-1.000000
cos(190):-0.984808      cose(190):-0.984911
cos(200):-0.939693      cose(200):-0.940055
cos(210):-0.866025      cose(210):-0.866667
cos(220):-0.766044      cose(220):-0.766804
cos(230):-0.642788      cose(230):-0.643347
cos(240):-0.500000      cose(240):-0.500000
cos(250):-0.342020      cose(250):-0.341289
cos(260):-0.173648      cose(260):-0.172565
cos(270):0.000000       cose(270):0.000000
cos(280):0.173648       cose(280):0.172565
cos(290):0.342020       cose(290):0.341289
cos(300):0.500000       cose(300):0.500000
cos(310):0.642788       cose(310):0.643347
cos(320):0.766044       cose(320):0.766804
cos(330):0.866025       cose(330):0.866667
cos(340):0.939693       cose(340):0.940055
cos(350):0.984808       cose(350):0.984911
cos(360):1.000000       cose(360):1.000000