• 当前位置:首页>>C语言>>C语言编程实例>>学数学也要学点偷懒的技术
  • 学数学也要学点偷懒的技术
  •     我做这个微小软件其实是想在做线性代数作业偷懒,她很快帮你解出N阶行列式的值,N阶线性方程的解,还能进行矩阵的相乘。我拿它出来其实是让大家在学习其他课程是也能学C语言。以后我还会对它升级。大家多多留意喔!

     

    #include<stdio.h>
    #define N  10
    #define M 10
    #define S 10
     
    void njie_valu();
    void njief_value();
    float Dvalue(float a[N][N],int n);
    void juzhench();

    void juzhench()
    {
    int i,j,k,m,n,s,s1;
    float a[M][S],b[S][N],c[M][N];
    printf("矩阵A是m*n矩阵。\n") ;
    printf("请分别输入m,s的值:");
    scanf("%d,%d",&m,&s);
    printf("矩阵B是s*n矩阵\n");
    printf("请分别输入s,n的值: ");
    scanf("%d,%d",&s1,&n);
    if(s!=s1)  {
       printf("两次输入的s的值不一样,不能进行运算!\n");
       exit(0);
      }
    for(i=0;i<m;i++) {
        for(j=0;j<n;j++)
          c[i][j]=0;
       }
    printf("请输入矩阵A\n");
    for(i=0;i<m;i++)  {
      for(j=0;j<s;j++)
         scanf("%f",&a[i][j]);
       }
    printf("请输入矩阵B:\n");
    for(i=0;i<s;i++) {
       for(j=0;j<n;j++)
         scanf("%f",&b[i][j]);
       }
    for(i=0;i<m;i++) {
      for(j=0;j<n;j++) {
         for(k=0;k<s;k++)
           c[i][j]=c[i][j]+a[i][k]*b[k][j];
          }
      }
    printf("相乘后的结果是:\n");
    for(i=0;i<m;i++) {
       for(j=0;j<n;j++)
          printf("%f   ",c[i][j]);
       printf("\n");
      }
    getch();
    }

    void njie_valu()
    {
    int i,j,n,k,h=0;
    float x[N],a[N][N],D=1.0;
    printf("说明:输入的n阶行列式,必须保证主对角线上的前n-1个元素都不为零!\n");
    printf("请输入是几阶行列式:\n");
    scanf("%d",&n);
    printf("请以行列式的正确顺序输入个元素:\n");
    for(i=0;i<n;i++){
         for(j=0;j<n;j++)
            scanf("%f",&a[i][j]);
      }
    for(i=0,j=i;i<n-1,j<n-1;i++,j++)
        if(a[i][j]==0)
             h++;
    if(h==n-1) {
         printf("请确保主对角线上前n-1个元素都不为零,请您调整!☆\n");
         printf("");
         exit(0);
       }
    for(k=0;k<n-1;k++)  {
        for(i=k+1;i<n;i++)
           x[i]=a[i][k]/a[k][k];
        for(i=k+1;i<n;i++)  {
             for(j=0;j<n;j++)
               a[i][j]=a[i][j]-a[k][j]*x[i];
             }
       }
    for(i=0;i<n;i++)  {
        for(j=0;j<n;j++)
            printf("%f   ",a[i][j]);
        printf("\n");
       }
    printf("\n");
    for(i=0,j=0;i<n,j<n;i++,j++)
        D=D*a[i][j];
    printf("D= %f",D);
    getch();
    }

     

    void njief_value()
    {
    float D;
    float Dvalue(float [N][N],int n);
    int i,j,k,n,h,m;
    float a[N][N],x[N],b[N],D1[N],a1[N][N],vable[N];
    printf("说明::输入的n阶行列式,必须保证n阶线性方程的系数及系数和常数项构成的行列式主对角线上的前n-1个元素都不为零!\n");
    printf("请输入是几阶线性方程: ");
      scanf("%d",&n);
    printf("请以行列式的正确顺序输入个元素:\n");
    for(i=0;i<n;i++){
         for(j=0;j<n;j++)
            scanf("%f",&a[i][j]);
      }
    for(i=0;i<n;i++){
      for(j=0;j<n;j++)
        printf("%f  ",a[i][j]);
      printf("\n");
      }
    printf("按正确顺序输入常数项:\n");
    for(i=0;i<n;i++)
      scanf("%f",&b[i]);
    for(i=0;i<n;i++)  {
       for(j=0;j<n;j++)
         a1[i][j]=a[i][j];
       }
    D=Dvalue(a1,n);
    printf("D= %f",D);
    getch();
    for(i=0;i<n;i++)  {
       for(j=0;j<n;j++)
         a1[i][j]=a[i][j];
       }
    for(i=0;i<n;i++) {
         for(k=0;k<n;k++)
      a1[k][i]=b[k];
         for(m=0;m<n;m++) {
     for(h=0;h<n;h++)
       printf("%f  ",a1[m][h]);
           printf("\n");
           }
          printf("\n");
         D1[i]=Dvalue(a1,n);
         printf("D%d= %f\n",i+1,D1[i]);
         getch();
         for(h=0;h<n;h++) {
            for(m=0;m<n;m++)
             a1[h][m]=a[h][m];
           }
       }
     for(i=0;i<n;i++)
        vable[i]=D1[i]/D;
     for(i=0;i<n;i++){
        printf("x%d=%f",i+1,vable[i]);
        printf("\n");
       }
    getch();
     }

    [1] [2] 下一页