鱼C论坛

 找回密码
 立即注册
查看: 1063|回复: 4

[已解决]编写一个程序

[复制链接]
发表于 2021-1-12 11:52:54 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
题目描述:编写程序,基于两个数组T={1,2,3,4,5}和S={12.1,32.4,45.3,36.7,78.9},拟合出最高次幂为5的多项式,并在4.6处插值。
烦请大佬指点
最佳答案
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[j][i] > E || _matrix[j][i] < -E)
                                {
                                        _matrix[j][k] /= _matrix[j][i];
                                }
                                
                                if(j > i)
                                {
                                        _matrix[j][k] -= _matrix[i][k];
                                }
                        }
                }
        }
        
        //回代
        for(i=n-1; i>=0 ; i--)
        {
                temp = 0.0;
                
                for(j=n-1 ; j>i; j--)
                {
                        temp += result[j] * _matrix[i][j];
                }
                
                result[i] = _matrix[i][n] - temp;
        }
}

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

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

        Guass(matrix,result,5);

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

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

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

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

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

        free(matrix);
        
        return 0;
}
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

发表于 2021-1-12 13:24:35 | 显示全部楼层
不是很明白什么意思
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 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[j][i] > E || _matrix[j][i] < -E)
                                {
                                        _matrix[j][k] /= _matrix[j][i];
                                }
                                
                                if(j > i)
                                {
                                        _matrix[j][k] -= _matrix[i][k];
                                }
                        }
                }
        }
        
        //回代
        for(i=n-1; i>=0 ; i--)
        {
                temp = 0.0;
                
                for(j=n-1 ; j>i; j--)
                {
                        temp += result[j] * _matrix[i][j];
                }
                
                result[i] = _matrix[i][n] - temp;
        }
}

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

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

        Guass(matrix,result,5);

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

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

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

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

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

        free(matrix);
        
        return 0;
}
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-1-12 14:06:49 | 显示全部楼层
运行结果如图:
未命名.JPG
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-1-12 23:19:00 | 显示全部楼层
这种表述也行?
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

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

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2025-1-12 03:55

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

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