论坛 › C\C++交流 › 编写一个程序

Intermezzo.L 发表于 2021-1-12 11:52:54

编写一个程序

题目描述：编写程序，基于两个数组T={1,2,3,4,5}和S={12.1,32.4,45.3,36.7,78.9}，拟合出最高次幂为5的多项式，并在4.6处插值。
烦请大佬指点

心驰神往 发表于 2021-1-12 13:24:35

不是很明白什么意思

xieglt 发表于 2021-1-12 13:56:32

#include <stdio.h>
#include <stdlib.h>

#define E        (0.000000001)

void Guass(double ** _matrix,double * result,int n)
{
        int i,j,k;
        double temp = 0.0;
       
        for(i=0; i<n; i++)
        {
                for(j=i ; j<n ; j++)
                {
                        for(k=i+1; k<=n; k++)
                        {       
                                if(_matrix > E || _matrix < -E)
                                {
                                        _matrix /= _matrix;
                                }
                               
                                if(j > i)
                                {
                                        _matrix -= _matrix;
                                }
                        }
                }
        }
       
        //回代
        for(i=n-1; i>=0 ; i--)
        {
                temp = 0.0;
               
                for(j=n-1 ; j>i; j--)
                {
                        temp += result * _matrix;
                }
               
                result = _matrix - temp;
        }
}

int main(void)
{
        int i,j;
        double temp,sum;
        double result = {0};
        double ** matrix = (double **)malloc(sizeof(double *) * 5);
       
        for(i=0 ; i<5 ; i++)
        {
                matrix = (double *)malloc(sizeof(double) * (5 + 1));
                sum = 1;
                for(j=0; j<5 ; j++)
                {
                        sum *= (i+1);
                        matrix = sum;
                }
        }

//        T={1,2,3,4,5}和S={12.1,32.4,45.3,36.7,78.9}
        matrix = 12.1;
        matrix = 32.4;
        matrix = 45.3;
        matrix = 36.7;
        matrix = 78.9;

        Guass(matrix,result,5);

        printf("f(x) = ");
       
        temp = 1;
        sum = 0;

        for(i=0 ; i<5 ; i++)
        {
                temp *= 4.6;
                sum += result * temp;

                if(i > 0)
                {
                        if(result > 0)
                        {
                                printf("+");
                        }
                }
                printf("%0.4lf*x^%d",result,i+1);
        }

        printf("\nf(4.6) = %lf\n",sum);

        for(i=0 ; i<5 ; i++)
        {
                free(matrix);
        }

        free(matrix);
       
        return 0;
}

xieglt 发表于 2021-1-12 14:06:49

运行结果如图：

风过无痕1989 发表于 2021-1-12 23:19:00

这种表述也行？

查看完整版本: 编写一个程序