鱼C论坛

 找回密码
 立即注册
查看: 1205|回复: 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
  1. #include <stdio.h>
  2. #include <stdlib.h>

  3. #define E        (0.000000001)

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

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

  58. //        T={1,2,3,4,5}和S={12.1,32.4,45.3,36.7,78.9}
  59.         matrix[0][5] = 12.1;
  60.         matrix[1][5] = 32.4;
  61.         matrix[2][5] = 45.3;
  62.         matrix[3][5] = 36.7;
  63.         matrix[4][5] = 78.9;

  64.         Guass(matrix,result,5);

  65.         printf("f(x) = ");
  66.        
  67.         temp = 1;
  68.         sum = 0;

  69.         for(i=0 ; i<5 ; i++)
  70.         {
  71.                 temp *= 4.6;
  72.                 sum += result[i] * temp;

  73.                 if(i > 0)
  74.                 {
  75.                         if(result[i] > 0)
  76.                         {
  77.                                 printf("+");
  78.                         }
  79.                 }
  80.                 printf("%0.4lf*x^%d",result[i],i+1);
  81.         }

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

  83.         for(i=0 ; i<5 ; i++)
  84.         {
  85.                 free(matrix[i]);
  86.         }

  87.         free(matrix);
  88.        
  89.         return 0;
  90. }

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

使用道具 举报

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

使用道具 举报

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

  3. #define E        (0.000000001)

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

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

  58. //        T={1,2,3,4,5}和S={12.1,32.4,45.3,36.7,78.9}
  59.         matrix[0][5] = 12.1;
  60.         matrix[1][5] = 32.4;
  61.         matrix[2][5] = 45.3;
  62.         matrix[3][5] = 36.7;
  63.         matrix[4][5] = 78.9;

  64.         Guass(matrix,result,5);

  65.         printf("f(x) = ");
  66.        
  67.         temp = 1;
  68.         sum = 0;

  69.         for(i=0 ; i<5 ; i++)
  70.         {
  71.                 temp *= 4.6;
  72.                 sum += result[i] * temp;

  73.                 if(i > 0)
  74.                 {
  75.                         if(result[i] > 0)
  76.                         {
  77.                                 printf("+");
  78.                         }
  79.                 }
  80.                 printf("%0.4lf*x^%d",result[i],i+1);
  81.         }

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

  83.         for(i=0 ; i<5 ; i++)
  84.         {
  85.                 free(matrix[i]);
  86.         }

  87.         free(matrix);
  88.        
  89.         return 0;
  90. }

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

使用道具 举报

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

使用道具 举报

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

使用道具 举报

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

本版积分规则

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

GMT+8, 2025-7-6 15:56

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

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