数学表达式中sinx泰勒近似的实现

今天咱们来看这样的一道题目，他是以某个数学表达式为背景，并且给了sinx的泰勒展开式让我们对sinx进行近似并且用这个值来计算上面的数学表达式，最后输出这个表达式的值

题目描述

这道题很明显是考察对于函数的声明和调用，主函数要用if函数分为两种情况分别对应这两种表达式，而我们还要在主函数之外定义一个新函数sin(double x)以此来对表达式里面出现过的sinx进行近似值的求解。

但是题目还要求我们将计算精度定在10^-6，那我们又该如何实现呢？大家可以发现泰勒展开式是在进行求和累加，那么毋庸置疑要实现这样的操作循环结构是少不了的。而所谓的泰勒展开本质就是在对函数进行拟合，累加的项越多，它的拟合效果就越好，如图👇

所以累加的项数越多，它的精度就会越高，而加的数的绝对值就会越来越小，以此来逐渐靠近sinx的真实值。那么我们应该如何保证10^-6的精度呢，应该累加几项应该循环几次呢？事实上，我们只要保证在加第i项的时候，这个第i项的绝对值小于10^-6就可以了，但是为了保证精度，我们要将判断条件定为第i项的值小于10^-7。而知道了结束循环的条件，我们要用的是do-while循环保证他至少执行一步，使sinx是有值的。

大家观察表达式会发现k的表达式里面有开方，那么就要用到sqrt()这个函数，这个函数我在“洛谷 P5708 【深基2.习2】三角形面积“ 有过介绍，这里我就不多赘述了。 除此，大家在泰勒展开求和形式的式子里面也发现了阶乘和n次幂的结构，阶乘不用我多说，大家只要用for循环或者while循环就能实现，而n次幂就要引用pow(int x,int y)这个函数了，他表示x^y，而且x、y都可以是整数或者浮点数。 最后还有判断时候要将最后一项的绝对值和10^-7进行比较，fabs(int x)就是对x取绝对值。而值得注意的是，这三个函数都在头文件#include <math.h>里面，所以，记得加上对应的头文件！！！

接下来我们看看代码是怎样的👇

C++代码

#include <bits/stdc++.h>        //万能头文件
#include <math.h>            //数学函数对应的头文件
using namespace std;
double sin(double x);        //对函数的声明
int main()
{
    double r,s,k;
    cin >> r >> s ;            //输入对应的r和s
    if(r*r<=s*s){
        k=sqrt(sin(r)*sin(r)+sin(s)*sin(s));        //第一个表达式
    }
    else{
        k=sin(r*s)*0.5;            //第二个表达式
    }
    cout << k ;               //输出答案即可
    return 0;
}
double sin(double x)            //输入值返回值均为double型
{
    int i=1,j,m=1;            //i用来表示n； j用来求阶乘；  m作为泰勒展开的某一项
    double a,sum=0;
    do{
        for(j=2*i-2;j<=2*i-1;j++){
            if(j<=0)
            j=1;
            m*=j;
        }
        a=pow(-1,i-1)/m*pow(x,2*i-1);
        sum+=a;
        i++;            //求一个项数n+1，即i+1
    }while(fabs(a)>=1e-7);        //用1e-7表示10^-7
    return sum;
}

这就是这道水题的全部内容，相对还是比较清晰明了的，希望对大家有帮助~