嵌入式

编程与应用
posts - 14, comments - 1, trackbacks - 0, articles - 0
  C++博客 :: 首页 :: 新随笔 :: 联系 :: 聚合  :: 管理

高斯约旦求矩阵的逆

Posted on 2010-04-17 10:55 陈显锋 阅读(774) 评论(1)  编辑 收藏 引用
#include<iostream>
using namespace std;

#define n 3
int main()
{
 int i,j,k;double m;
 double a[n][n],E[n][n];

 //单位矩阵E[n][n]
 for(i=0;i<n;i++){
     for(j=0;j<n;j++){
         if(i==j)
             E[i][j]=1;
         else
             E[i][j]=0;
     }
}
 
//输入矩阵a[n][n]
 cout<<"输入矩阵："<<endl;
 for(i=0;i<n;i++)
  for(j=0;j<n;j++)
   cin>>a[i][j];

//上三角变换
 for(k=0;k<n-1;k++){
  for(i=k+1;i<n;i++){
   m=a[i][k]/a[k][k];
   for(j=0;j<n;j++){
    a[i][j]=a[i][j]-m*a[k][j];
    E[i][j]=E[i][j]-m*E[k][j];
   }
  }
 }
 //下三角变换
 for(k=n-1;k>0;k--){
     for(i=k-1;i>=0;i--){
         m=a[i][k]/a[k][k];
         for(j=0;j<n;j++){
             a[i][j]=a[i][j]-m*a[k][j];
             E[i][j]=E[i][j]-m*E[k][j];
         }
     }
 }
//单位矩阵……
 for(i=0;i<n;i++)
     for(j=0;j<n;j++){
         E[i][j]=E[i][j]/a[i][i];
     }

 //输出矩阵的逆
 cout<<"矩阵的逆是："<<endl;
for(i=0;i<n;i++){
    for(j=0,k=0;j<n;j++){
        cout<<E[i][j]<<"\t";k++;
        if(k==n)
            cout<<endl;
    }
}
}




Feedback

# re: 高斯约旦求矩阵的逆  回复  更多评论   

2010-08-17 09:07 by zhan
这不是高斯约旦算法

只有注册用户登录后才能发表评论。
网站导航: 博客园   IT新闻   BlogJava   知识库   博问   管理