sinx的泰勒展开估计值计算

废话不多说，接下来我要介绍我最近研究的sinx泰勒展开计算

首先呢，我们要知道sinx的通项公式

大概就是这个样子

然后呢我们要构思sinx的计算，总体思路呢就是用for循环计算加和每一项的值。而且这里有几个点要注意：

1.首先就是比较能够想到的，就是什么时候要用浮点数什么时候不用，尤其是定义函数的时候很容易就int定义了，所以这里要想清楚来；

2.其次就是在通项公式里面的，每隔一项要变换正负，接下来的代码中我会用一个计数的ret表示，初始ret会为1，每加和一项都要令ret=-ret变换正负号这样可以省去这个比较麻烦的表示方法；

3.其次就是要注意分子和分母的循环次数不一样，换句话说，如果不定义函数直接把分子和分母放在同一个for循环里面的话，你会发现第n项分子是从1到2n-1，有2n-1个，分子呢只有n个，这样的话计算不出准确的值了

构思的话就先定义一个阶乘函数

int jiecheng(int n)//定义阶乘函数
{
	int i = 1;
	int result = 1;
	for (i = 1; i <= n; i++)
	{
		result *= i;
	}
	return result;
}

这样我们就可以用jiecheng（i）表示从1乘到i，就不会出现上面说的第三点，同时也很好地简化了代码 

其次就是定义sinx函数了，这样我们就可以在主函数里面简简单单地引用，而且代码整体可以更加层次分明

double Talorsinx(double x, int n)//定义sinx求值函数
{
	double sinx = 0.0;
	int i = 1;
	int ret = 1;//计数器变换每项的正负号

	for (i = 1; i <= 2 * n - 1; i += 2)
	{
		sinx += ret * pow(x, i) / jiecheng(i);
		ret = -ret;//交替变换正负
	}
	return sinx;
}

我们规定x为弧度制，n为第n项（也就是整数），那么就可以在主函数里面直接使用scanf写入x和n就好了

int main()
{
	int n = 0;
	double x = 0.0;
	printf("请输入sinx的x(用弧度制表示)：\n");
	scanf("%lf", &x);
	printf("请输入项数n(n为整数)\n");
	scanf("%d", &n);
	double sinx = Talorsinx(x, n);
	printf("sin(%.4lf)计算值约为%.6lf\n", x, sinx);
	return 0;
}

在最后输出结果的时候我使用了%.4lf和%.6lf就是表示最后一位的最高位，否则会出现inf的警告（大概就是小数位太高了），其实这个精确数字也可以自行修改。

最后全部的代码就在这里了

#define _CRT_SECURE_NO_WARNINGS 1
#include <stdio.h>
#include <math.h>
int jiecheng(int n)//定义阶乘函数
{
	int i = 1;
	int result = 1;
	for (i = 1; i <= n; i++)
	{
		result *= i;
	}
	return result;
}


double Talorsinx(double x, int n)//定义sinx求值函数
{
	double sinx = 0.0;
	int i = 1;
	int ret = 1;//计数器变换每项的正负号

	for (i = 1; i <= 2 * n - 1; i += 2)
	{
		sinx += ret * pow(x, i) / jiecheng(i);
		ret = -ret;//交替变换正负
	}
	return sinx;
}


int main()
{
	int n = 0;
	double x = 0.0;
	printf("请输入sinx的x(用弧度制表示)：\n");
	scanf("%lf", &x);
	printf("请输入项数n(n为整数)\n");
	scanf("%d", &n);
	double sinx = Talorsinx(x, n);
	printf("sin(%.4lf)计算值约为%.6lf\n", x, sinx);
	return 0;
}

ok，今天的分享就到这里了。下次再来做一个更不错的程序。