編程求解逆矩陣

1、精選優(yōu)質(zhì)文檔-傾情為你奉上本程序是運(yùn)用矩陣的初等變換的知識(shí)來求解矩陣的逆矩陣的。1.理論基礎(chǔ)：若存在矩陣P使得矩陣A有PA = E那么 P = A。所以如果我們將（A，E）化成了（E，P），那么，P就是我們要求的A。2.由上面的分析可知，我們思路可分為以下3步：（1）申請(qǐng)一個(gè)新的矩陣，它是（A,E）(2) 將A先對(duì)角化（3）再將A單位化E，這是A后面的部分就是A了。3.編程實(shí)現(xiàn)#include<stdio.h>#include<stdlib.h>float A110,A2102*10,B1010;/原矩陣、加上單位矩陣后的矩陣、所求的逆矩陣void ShangSanJi

2、ao(int n)/化成上三角int k,i,j,l=0,temp;if(A200 = 0)while(A2l0=0&&l<n-1)l+;if(l!=0)/如果00處的數(shù)為0，則交換for(i=0;i<2*n;i+)temp = A2li;A2li = A20i;A20i = temp;for(k=0;k<n-1;k+)for(i=k+1;i<n;i+)A1i = A2ik/A2kk;for(l=k+1;l<n;l+)for(j=0;j<2*n;j+)A2lj = A2lj - A2kj*A1l;void DuiJiaoXian(int n)

3、/將對(duì)角線上的數(shù)化1，這樣方便將原矩陣化成單位矩陣int i,j;for(i=n-1;i>=0;i-)A1i=A2ii;for(i=n-1;i>=0;i-)for(j=0;j<2*n;j+)A2ij = A2ij/A1i;void DanWei(int n)/將原矩陣化成單位矩陣int i,j,l;float temp;for(i=n-1;i>0;i-)for(j=i-1;j>=0;j-)temp = A2ji;for(l=0;l<2*n;l+)A2jl = A2jl - A2il*temp;int main()int n;/矩陣的階數(shù)scanf("

4、;%d",&n);if(n>10)printf("矩陣不再表示的范圍內(nèi)！");exit(-1);int i,j;for(i=0;i<n;i+)for(j=0;j<n;j+)scanf("%f",&A2ij);for(j=n;j<2*n;j+)A2ij=0;for(i=0,j=n;i<n,j<2*n;i+,j+)A2ij = 1;/*for(i=0;i<n;i+)for(j=0;j<2*n;j+)printf("%f ",A2ij);printf("n&

5、quot;);printf("n");*/顯示輸入的矩陣ShangSanJiao(n);/化成上三角/*for(i=0;i<n;i+)for(j=0;j<2*n;j+)printf("%f ",A2ij);printf("n");printf("n");*/顯示上三角的矩陣DuiJiaoXian(n);/對(duì)角化1/*for(i=0;i<n;i+)for(j=0;j<2*n;j+)printf("%f ",A2ij);printf("n");printf("n");*/顯示對(duì)角線化1后的矩陣DanWei(n);/將原矩陣單位化/*for(i=0;i<n;i+)for(j=0;j<2*n;j+)printf("%f ",A2ij);printf("n");printf("n");*/將所求矩陣移到B中for(i=0;i<n;i+)for(j=n;j<2*n;j+)Bij-n =A2ij;/顯示所求的逆矩陣for(i=0;i<