Tìm Định thức và ma trận nghịch đảo

[conrongchautien]

Bài hoàn chỉnh đây:

/* DINH THUC VA MA TRAN NGHICH DAO */ #include<stdio.h>#include<conio.h> void main(){int i,j,k,m,n,na,dau=1;double det,x,tam;double a[50][50],b[50][50];clrscr();printf("\n CHAO MUNG DEN VOI CHUONG TRINH TINH DINH THUC VA MA TRAN NGHICH DAO");printf("\n");printf("\n NHAP CAP MA TRAN:"); scanf("%d",&n);na=2*n;//so cot ma tran chuyen tiep.//Nhap ma tran A.-------------------------------for(i=0;i<n;i++)    for(j=0;j<n;j++)        {        printf("a[%d][%d]=",i+1,j+1);        scanf("%lf",&a[i][j]);        } //Tao ma tran lien ket.---------------------------for(i=0;i<n;i++)    for(j=n;j<na;j++)        {        if(i==(j-n)) a[i][j]=1;        else a[i][j]=0;        } //Ham giai phuong trinh.--------------------//Dung phep BDSC dua Aij ve ma tran cheo.for(j=0;j<n;j++)    {    for(i=n-1;i>=j;i--)        {        if(a[i][j]==0) continue;//Neu phan tu a bang 0 thi tang i.        if((i>j)&&(a[i][j]!=0))             {            k=i-1;            while((k>=j)&&(a[k][j]==0)) k--;//Xet cac phan tu a khac 0 tren cung cot j            if(k<j)                {                //Doi cho dong thu i va j cho ca ma tran mo rong.                for(m=0;m<na;m++)                    {                    tam=a[i][m];                    a[i][m]=a[j][m];                    a[j][m]=tam;                    }                dau=-dau;                }            if((k>=j)&&(k>=0))                {                x=-a[i][j]/a[k][j];                //Nhan dong thu k cho x roi cong vao dong thu i.                for(m=0;m<na;m++)                    a[i][m]+=a[k][m]*x;                }            }        }    }for(j=n-1;j>=0;j--)    {    for(i=0;i<=j;i++)        {        if(a[i][j]==0) continue;        if((i<j)&&(a[i][j]!=0))            {            k=i+1;            while((k<=j)&&(a[k][j]==0)) k++;            if((k<=j)&&(k<na))                {                x=-a[i][j]/a[k][j];                //Nhan dong thu k voi x roi cong vao dong thu i.                for(m=0;m<na;m++)                    a[i][m]+=a[k][m]*x;                }            }        }    }//Tinh dinh thuc.------------------------------------det=1;for(i=0;i<n;i++)    for(j=0;j<n;j++)        if(i==j) det*=a[i][j];printf("\n Det=%5.2lf",det*dau);//Tinh ma tran nghich dao.-----------------------if(det!=0)    {    for(i=0;i<n;i++)        for(j=n;j<na;j++)            {            a[i][j]=a[i][j]/a[i][i];            }     for(i=0;i<n;i++)        for(j=n;j<na;j++)            {            b[i][j-n]=a[i][j];            }    //In ra ma tran nghich dao.-----------------    printf("\n MA TRAN NGHICH DAO CUA A:\n");    for(i=0;i<n;i++)        {        for(j=0;j<n;j++)            printf("%5.1lf",b[i][j]);        printf("\n");        }    }else printf("\n KHONG CO MA TRAN NGHICH DAO.");printf("\n GOOBYE.CHUC MOI NGUOI VUI VE!");getch();}//6:01 PM 4/22/2010

Cám ơn thầy onminh rất nhiều.

[hoan_it_cgh]

các dajka bài này vẫn có chỗ sai
ma trận là
0 2
1 3
định thức =-2 chứ ko phải là 2
còn ma trận nghịch đảo thì chuẩn rồi .mọi ng Xem lại giúp e với

[conrongchautien]

các dajka bài này vẫn có chỗ sai
ma trận là
0 2
1 3
định thức =-2 chứ ko phải là 2
còn ma trận nghịch đảo thì chuẩn rồi .mọi ng Xem lại giúp e với

Mình vừa check lại rồi. Đúng là do mình sai sót. mình sửa lại cho chính xác rồi. Bạn click vào đây xem lại nha

[onminh]

Bài hoàn chỉnh đây:
Cám ơn thầy onminh rất nhiều.

Em sửa làm giảm đi cái hay. Sao không viết dau = -dau sau mỗi lần đổi hàng mà lại viết if(dau%2!=0) det*=-1;. Củ chuối wá. Hix...

[conrongchautien]

Em sửa làm giảm đi cái hay. Sao không viết dau = -dau sau mỗi lần đổi hàng mà lại viết if(dau%2!=0) det*=-1;. Củ chuối wá. Hix...

Em sửa lại rồi. Tại nghĩ sao làm vậy thôi, chưa bik tối ưu
http://forums.congdongcviet.com/show...3&postcount=21
Cám ơn thầy!

[hoan_it_cgh]

#include<conio.h>
#include<stdio.h>
#include<iostream.h>
#include<math.h>
#include<iomanip.h>
#include <dos.h>
#include <math.h>
#include <stdlib.h>
#include <ctype.h>

class MT
{ private:
 int n;
 double a[10][10];
 public:

 //---------------TOAN TU CONG-------------
MT MT::operator+(const MT &right)
{//if ( n!=right.n) exit (0) ;
 MT c;
  for (int i=1;i<=n;i++)
     for (int j=1;j<=n;j++)
       {     c.a[i][j] = a[i][j] + right.a[i][j] ;  }
     c.n=n;
 return c ;
}

 //-----------------TOAN TU TRU---------------
MT MT::operator-(const MT &b)
 {
  MT d;
 if (n!=b.n) cout<<"Kh«ng tÝnh ®­îc ";
else
for(int i=1;i<=n;i++)
  for(int j=1;j<=n;j++)
  { d.a[i][j]=a[i][j]-b.a[i][j];}
  d.n=n;
  return d;
  }
//-----------------TOAN TU NHAN--------------
MT MT::operator*(const MT &t)
{ MT c;
  for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
      for(int k=1;k<=n;k++)
    {  c.a[i][j]+=a[i][k]*t.a[k][j];}
     c.n=n;
     return c;
     }


  void nhap();
  void in();
  void dt_nd(MT &a);
  };
  //================Dinh thuc nghich dao==========
  void MT::dt_nd(MT &u)
{
int i,j,k,m,n,na,dau=0; //Thay vì dau=1;
double det,x,tam;
double b[10][10];
/*clrscr();
printf("\n CHAO MUNG DEN VOI CHUONG TRINH TINH DINH THUC VA MA TRAN NGHICH DAO");
printf("\n");
printf("\n NHAP CAP MA TRAN:"); scanf("%d",&n);*/
na=2*n;//so cot ma tran chuyen tiep.
//Nhap ma tran A.-------------------------------
/*for(i=0;i<n;i++)
    for(j=0;j<n;j++)
    {
    printf("a[%d][%d]=",i+1,j+1);
      //  scanf("%lf",&a[i][j]);
      cin>>a[i][j];
    }
  */
//Tao ma tran lien ket.---------------------------
for(i=0;i<n;i++)
    for(j=n;j<na;j++)
    {
    if(i==(j-n)) a[i][j]=1;
    else a[i][j]=0;
    }

//Ham giai phuong trinh.--------------------
//Dung phep BDSC dua Aij ve ma tran cheo.
for(j=0;j<n;j++)
    {
    for(i=n-1;i>=j;i--)
    {
    if(a[i][j]==0) continue;//Neu phan tu a bang 0 thi tang i.
    if((i>j)&&(a[i][j]!=0))
        {
        k=i-1;
        while((k>=j)&&(a[k][j]==0)) k--;//Xet cac phan tu a khac 0 tren cung cot j
        if(k<j)
        {
        //Doi cho dong thu i va j cho ca ma tran mo rong.
        for(m=0;m<na;m++)
            {
            tam=a[i][m];
            a[i][m]=a[j][m];
            a[j][m]=tam;
            }
        dau++; // thay vì :dau=-dau;
        }
        if((k>=j)&&(k>=0))
        {
        x=-a[i][j]/a[k][j];
        //Nhan dong thu k cho x roi cong vao dong thu i.
        for(m=0;m<na;m++)
            a[i][m]+=a[k][m]*x;
        }
        }
    }
    }
for(j=n-1;j>=0;j--)
    {
    for(i=0;i<=j;i++)
    {
    if(a[i][j]==0) continue;
        if((i<j)&&(a[i][j]!=0)) 
            { 
            k=i+1; 
            while((k<=j)&&(a[k][j]==0)) k++; 
            if((k<=j)&&(k<na)) 
                { 
                x=-a[i][j]/a[k][j]; 
                //Nhan dong thu k voi x roi cong vao dong thu i. 
                for(m=0;m<na;m++) 
                    a[i][m]+=a[k][m]*x; 
                } 
            } 
        } 
    } 
//Tinh dinh thuc.------------------------------------ 
det=1; 
for(i=0;i<n;i++) 
    for(j=0;j<n;j++) 
        if(i==j) det*=a[i][j]; // B? :*dau 
        if(dau%2!=0) det*=-1; //N?u s? l?n d?i ch? 2 dòng là l? thì d?i d?u giá tr? d?nh th?c 
printf("\n Det=%5.2lf",det); 
//Tinh ma tran nghich dao.----------------------- 
if(det!=0) 
    { 
    for(i=0;i<n;i++) 
        for(j=n;j<na;j++) 
            { 
            a[i][j]=a[i][j]/a[i][i]; 
            } 

    for(i=0;i<n;i++) 
        for(j=n;j<na;j++) 
            { 
            b[i][j-n]=a[i][j]; 
            } 
    //In ra ma tran nghich dao.----------------- 
    printf("\n MA TRAN NGHICH DAO CUA A:\n"); 
    for(i=0;i<n;i++) 
        { 
        for(j=0;j<n;j++) 
            printf("%5.1lf",b[i][j]); 
        printf("\n"); 
        } 
    } 
else printf("\n KHONG CO MA TRAN NGHICH DAO."); 
printf("\n GOOBYE.CHUC MOI NGUOI VUI VE!");
return;
getch(); 
} 

   void MT::nhap()
  { cout<<"\n CÊp ma trËn=: ";
  cin>>n;
  for(int i=1;i<=n;i++)
  for(int j=1;j<=n;j++)
  { cout<<"\n PhÇn tö hµng"<<i<<" cét"<<j<<"= ";
  cin>>a[i][j];
  }
  }
void MT::in()
   { for(int i=1;i<=n;i++)
    { cout<<"\n";
      for(int j=1;j<=n;j++)
      cout<<setw(6)<<a[i][j]<<" ";
     }
    }

  void main()
  { MT a,b,h,k,tich;


  a.nhap();
  clrscr();
  cout<<"\nMa trËn A: ";
  a.in();
 b.nhap();
 clrscr();
cout<<"\nMa trËn B: ";
a.in();
b.in();
h=a+b;
k=a-b;
tich=a*b;
cout<<"\nMa trËn tæng: ";
h.in();
cout<<"\nMa trËn hiÖu";
k.in();
cout<<"\nMa trËn tÝch";
tich.in();
a.dt_nd(a);
  getch();
  } 

mọi ng xem hộ em cái này với
lỗi cannot convert MT to int la gì nhi? ko hiểu
Sửa giúp em với ('-")

[onminh]

Vấn đề là ở chỗ này: Hàm thành viên void dt_nd(MT &a) được định nghĩa bên ngoài lớp MT, nhưng thiếu MT::. Thứ hai, tham biến MT &a trùng tên với biến trong lớp. Nếu sử dụng nó thì viết a.a[i][j], còn swr dụng b thì viết b.a[i][j].

[hoan_it_cgh]

em đã sửa rồi đã hết lỗi nhưng nó ko hiện kq ra màn hình
có phải việc truyền đối số sai sót chỗ nào ko ạ
onminh sữa giúp em với

[tnglinh]

hoàn xem lại phần nhập ma trận đi
bài của hoàn vòng "for" đều bắt đầu từ 1
nhưng ở phần tính định thức và ma trận nghich đảo vòng for lại bắt đầu từ 0
sửa lại đi rồi xem kết quả nhé

"lỗi cannot convert MT to int la gì nhi ? ko hiểu" cái nè là ko thể chuyển MT sang kiểu int vì nó đang ở kiểu double

[cdcv282]

Đây, nếu code trực tiếp thì dài dòng và khó hiểu hơn đó, chịu khó đọc nhé. Tính det theo Gauss, còn ma trận nghịch đảo theo Gauss-Jordan đó.

#include<stdio.h>#include<conio.h> void main(){    int i,j,k,n, sign = 1;     double det = 1; double temp;    double a[10][10],b[10][10];// Nhap cap cua ma tran:    printf("\nCap cua ma tran vuon, n = "); scanf("%d",&n); // Nhap ma tran a:    for (i = 0; i < n; i++)        for (j = 0; j < n; j++)            {                printf("a[%d][%d] = ", i, j);                scanf("%lf",&a[i][j]);            } // Tinh dinh thuc cua a    for (i = 0; i < n-1; i++)        {            if (a[i][i] == 0) // Nếu gặp phần tử trên đường chéo chính bằng 0 thì tìm hàng khác để đổi                {                    k = i+1;                    while (k < n && a[k][i] == 0) k++;                    if (k==n) { printf("\ndet(a) = 0"); break;} // Không timg thấy, tức det(a) = 0                    for (j = i; j < n; j++)// Đổi hàng i với hàng k                        {                            temp = a[i][j];                            a[i][j] = a[k][j];                            a[k][j] = temp;                        }                    sign = -sign; // Định thức đổi dấu đó                }            for (j = i+1; j < n; j++)// Biển đổi để các phần tử cùng cột ở hàng dưới bằng 0 đó                {                    temp = -a[j][i]/a[i][i];                    for (k = i+1; k < n; k++)   a[j][k] += temp*a[i][k]; // Nhân hàng i với (-a[j][i]/a[i][i]) rồi cộng vào hàng j                }            det *= a[i][i]; // Tính dần det(a)        }     det *= a[n-1][n-1]; // Nhân với phần tử cuối chéo chính nữa là xong    printf("\ndet(a) = %lf\n",det); //Tinh ma tran nguocif (det != 0){    for (i = 0; i < n; i++) // Tạo ra b là ma trận đơn vị đó        for (j = 0; j < n; j++)        if (i==j) b[i][j] = 1; else b[i][j] = 0;     for (i = 0; i < n; i++) // Xử lý từ hàng đầu đến hàng cuối        {            if (a[i][i] == 0) // Nếu gặp phần tử trên đường chéo chính bằng 0 thì đổi hàng                 {                    k = i+1;                    while (k < n && a[k][i] == 0) k++;                    for (j = 0; j < n; j++) // Đổi hàng đó của a đó, cả với ma trận b nữa                        {                            temp = a[i][j];                            a[i][j] = a[k][j];                            a[k][j] = temp;                            temp = b[i][j];                            b[i][j] = b[k][j];                            b[k][j] = temp;                        }                }             temp = a[i][i];            for (j = i; j < n; j++) a[i][j] /= temp;            for (j = 0; j < n; j++) b[i][j] /= temp;             for (j = i+1; j < n; j++)                {                    temp = -a[j][i];                    for (k = i; k < n; k++) a[j][k] += temp*a[i][k];                     for (k = 0; k < n; k++) b[j][k] += temp*b[i][k];                 }        }// Kết thúc quá trình Gauss     for (i = n-1; i > 0; i--) // Bắt đầu quá trình Jordan        for (j = i-1; j >=0; j--)            {                temp = -a[j][i];                for (k = n-1; k >= i; k--)  a[j][k] += temp*a[i][k];                 for (k = 0; k < n; k++) b[j][k] += temp*b[i][k];             }     for (i = 0; i < n; i++)// In ra ma trận b, bây giờ là ma trận nghich đảo của a        {            for (j = 0; j < n; j++) printf("%8.3lf\t", b[i][j]);            printf("\n");        }}else printf("\nKhong co ma tran nghich dao\n");    getch();}[/i]

Em muốn phát triển bài này theo cách viết Modul (chương trình con). Xin các pro giúp cho...