设为首页收藏本站
切换到窄版

鱼C论坛

 找回密码
 立即注册
鱼C论坛»论坛 › 编程语言专区 › C\C++交流 › 编写一个程序
返回列表 发新帖
查看: 1209|回复: 4

[已解决]编写一个程序

[复制链接]
Intermezzo.L
发表于 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处插值。
烦请大佬指点
最佳答案
月排行榜 / 总排行榜
xieglt
2021-1-12 13:56:32
  1. #include <stdio.h>
  2. #include <stdlib.h>

  4. #define E        (0.000000001)

  6. void Guass(double ** _matrix,double * result,int n)
  7. {
  8.         int i,j,k;
  9.         double temp = 0.0;
  10.        
  11.         for(i=0; i<n; i++)
  12.         {
  13.                 for(j=i ; j<n ; j++)
  14.                 {
  15.                         for(k=i+1; k<=n; k++)
  16.                         {       
  17.                                 if(_matrix[j][i] > E || _matrix[j][i] < -E)
  18.                                 {
  19.                                         _matrix[j][k] /= _matrix[j][i];
  20.                                 }
  21.                                
  22.                                 if(j > i)
  23.                                 {
  24.                                         _matrix[j][k] -= _matrix[i][k];
  25.                                 }
  26.                         }
  27.                 }
  28.         }
  29.        
  30.         //回代
  31.         for(i=n-1; i>=0 ; i--)
  32.         {
  33.                 temp = 0.0;
  34.                
  35.                 for(j=n-1 ; j>i; j--)
  36.                 {
  37.                         temp += result[j] * _matrix[i][j];
  38.                 }
  39.                
  40.                 result[i] = _matrix[i][n] - temp;
  41.         }
  42. }

  44. int main(void)
  45. {
  46.         int i,j;
  47.         double temp,sum;
  48.         double result[5] = {0};
  49.         double ** matrix = (double **)malloc(sizeof(double *) * 5);
  50.        
  51.         for(i=0 ; i<5 ; i++)
  52.         {
  53.                 matrix[i] = (double *)malloc(sizeof(double) * (5 + 1));
  54.                 sum = 1;
  55.                 for(j=0; j<5 ; j++)
  56.                 {
  57.                         sum *= (i+1);
  58.                         matrix[i][j] = sum;
  59.                 }
  60.         }

  62. //        T={1,2,3,4,5}和S={12.1,32.4,45.3,36.7,78.9}
  63.         matrix[0][5] = 12.1;
  64.         matrix[1][5] = 32.4;
  65.         matrix[2][5] = 45.3;
  66.         matrix[3][5] = 36.7;
  67.         matrix[4][5] = 78.9;

  69.         Guass(matrix,result,5);

  71.         printf("f(x) = ");
  72.        
  73.         temp = 1;
  74.         sum = 0;

  76.         for(i=0 ; i<5 ; i++)
  77.         {
  78.                 temp *= 4.6;
  79.                 sum += result[i] * temp;

  81.                 if(i > 0)
  82.                 {
  83.                         if(result[i] > 0)
  84.                         {
  85.                                 printf("+");
  86.                         }
  87.                 }
  88.                 printf("%0.4lf*x^%d",result[i],i+1);
  89.         }

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

  93.         for(i=0 ; i<5 ; i++)
  94.         {
  95.                 free(matrix[i]);
  96.         }

  98.         free(matrix);
  99.        
  100.         return 0;
  101. }

复制代码
小甲鱼最新课程 -> https://ilovefishc.com
回复

使用道具 举报

心驰神往
发表于 2021-1-12 13:24:35 | 显示全部楼层
不是很明白什么意思
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

xieglt
发表于 2021-1-12 13:56:32 | 显示全部楼层    本楼为最佳答案   
  1. #include <stdio.h>
  2. #include <stdlib.h>

  4. #define E        (0.000000001)

  6. void Guass(double ** _matrix,double * result,int n)
  7. {
  8.         int i,j,k;
  9.         double temp = 0.0;
  10.        
  11.         for(i=0; i<n; i++)
  12.         {
  13.                 for(j=i ; j<n ; j++)
  14.                 {
  15.                         for(k=i+1; k<=n; k++)
  16.                         {       
  17.                                 if(_matrix[j][i] > E || _matrix[j][i] < -E)
  18.                                 {
  19.                                         _matrix[j][k] /= _matrix[j][i];
  20.                                 }
  21.                                
  22.                                 if(j > i)
  23.                                 {
  24.                                         _matrix[j][k] -= _matrix[i][k];
  25.                                 }
  26.                         }
  27.                 }
  28.         }
  29.        
  30.         //回代
  31.         for(i=n-1; i>=0 ; i--)
  32.         {
  33.                 temp = 0.0;
  34.                
  35.                 for(j=n-1 ; j>i; j--)
  36.                 {
  37.                         temp += result[j] * _matrix[i][j];
  38.                 }
  39.                
  40.                 result[i] = _matrix[i][n] - temp;
  41.         }
  42. }

  44. int main(void)
  45. {
  46.         int i,j;
  47.         double temp,sum;
  48.         double result[5] = {0};
  49.         double ** matrix = (double **)malloc(sizeof(double *) * 5);
  50.        
  51.         for(i=0 ; i<5 ; i++)
  52.         {
  53.                 matrix[i] = (double *)malloc(sizeof(double) * (5 + 1));
  54.                 sum = 1;
  55.                 for(j=0; j<5 ; j++)
  56.                 {
  57.                         sum *= (i+1);
  58.                         matrix[i][j] = sum;
  59.                 }
  60.         }

  62. //        T={1,2,3,4,5}和S={12.1,32.4,45.3,36.7,78.9}
  63.         matrix[0][5] = 12.1;
  64.         matrix[1][5] = 32.4;
  65.         matrix[2][5] = 45.3;
  66.         matrix[3][5] = 36.7;
  67.         matrix[4][5] = 78.9;

  69.         Guass(matrix,result,5);

  71.         printf("f(x) = ");
  72.        
  73.         temp = 1;
  74.         sum = 0;

  76.         for(i=0 ; i<5 ; i++)
  77.         {
  78.                 temp *= 4.6;
  79.                 sum += result[i] * temp;

  81.                 if(i > 0)
  82.                 {
  83.                         if(result[i] > 0)
  84.                         {
  85.                                 printf("+");
  86.                         }
  87.                 }
  88.                 printf("%0.4lf*x^%d",result[i],i+1);
  89.         }

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

  93.         for(i=0 ; i<5 ; i++)
  94.         {
  95.                 free(matrix[i]);
  96.         }

  98.         free(matrix);
  99.        
  100.         return 0;
  101. }

复制代码
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

xieglt
发表于 2021-1-12 14:06:49 | 显示全部楼层
运行结果如图：
未命名.JPG
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

风过无痕1989
发表于 2021-1-12 23:19:00 | 显示全部楼层
这种表述也行？
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

GMT+8, 2025-7-10 12:59

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

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