wd2.txt

这是本人编写的有关矩阵数值计算的几个C语言程序,方便大家学习<矩阵计算方法>!

/*算法4.1----p96*/
#include"stdio.h"
#include"math.h"
main()
{
 int i,n,k,j;
 float beta[100],x[100];
Tag_D:
 k=0;
 printf("\n\n请输入Vandermande方阵的阶数n：\n");
 scanf("%d",&n);
 n=n-1;
 printf("\n\n请输入原始数据即b: \n");
 for(i=0;i<=n;i++)
    {
     printf("Input beta[%d]> ",i);
     scanf("%f,",&beta[i]);
    }
 printf("\n\n请输入V方阵的参数即V(x0,x1,...,xn): \n");
 for(i=0;i<=n;i++)
    {
     printf("Input x[%d]> ",i);
     scanf("%f,",&x[i]);

    }
 for(i=0;i<=n;i++)                                             /*检查输入是否有误*/
    for(j=i+1;j<=n;j++)
       if(x[i]==x[j])
         {
          printf("\n\nx[%d]=x[%d],输入有误！\n",i,j);
          goto Tag_D; 
         }
Tag_0:
 for(i=n;i>=k+1;i--)
    beta[i]=(float)(beta[i]-x[k]*beta[i-1]);
 if(k<n-1)
   {
    k=k+1;
    goto Tag_0;
   }
Tag_1:
 for(i=k+1;i<=n;i++)
    beta[i]=(float)(beta[i]/(x[i]-x[i-k-1]));
 for(i=k;i<=n-1;i++)
    beta[i]=beta[i]-beta[i+1];
 if(k>0)
   {
    k=k-1;
    goto Tag_1;
   }
 printf("\n\n输出数据Z: \n");
 for(i=0;i<=n;i++)
    printf("%.3f  ",beta[i]);

}