編程求解逆矩陣

1、本程序是運(yùn)用矩陣的初等變換的知識(shí)來(lái)求解矩陣的逆矩陣的。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#includefloat A1 10, A2 10 2*10, B 10 10; /原矩陣、加上單位矩陣后的矩陣、所求的逆矩陣vo i d ShangSanJ i ao (i nt n)/化成上三角i nt k, i, j, I

2、二0, temp;if (A20 0 = 0)while(A2l0=0&ln-1)I+；)if (I !=0)/如果0 0處的數(shù)為0,則交換for (i=0;i 2*n;i+)temp = A2Ii;A2li = A20i;A2 0i = temp;)for (k=0;kn1;k+)for(i=k+1;in;i+)A1i = A2ik/A2k k;for(l=k+1; Kn; I+)for (j=0;j=0;i一一)A1i=A2ii;1for (i 二r)T ; i=0; i -)for (j=0;j0; i-)for (j=i-1;j二0;j-)temp = A2j i;for (l=0;

3、K2*n; I+)A2j I = A2j I - A2i I*temp； )int main()int n;矩陣的階數(shù)scanf(%d, &n);if (n10)printf (H矩陣不再表示的范圍內(nèi)! M); exit (T);)int i, j;for (i=0;in;i+)for (j二0;jn;j+) scanf (%fM, &A2i j);for (j=r); j2*n; j+) A2ij二 0;)for (i=0, j=n; in, j2*n; i+, j+)A2ij = 1;/* for (i =0;i n;i+)for (j=0;j2*n;j+)pr intf (%f H, A

4、2i j);1pr intf Cn”);)printfCXn) ;*/顯示輸入的矩陣ShangSanJiao(n) ;/化成上三角/* for (i =0;i n;i+)for (j=0;j2*n;j+)pr intf (M%f H, A2i j);pr intf Cn);)printfCAn) ；*/顯示上三角的矩陣Dui JiaoXian(n) ;/對(duì)角化 1/* for (i =0;i n;i+)for (j=0;j2*n;j+)pr intf (M%f H, A2i j);pr intf CrT);)printf (HnM) ;*/顯示對(duì)角線化1后的矩陣DanWei (n) ;/將原矩陣單位化/* for (i=0;in;i+)for (j=0;j2*n;j+)pr intf (H%f H, A2i j);pr i rrtf (”n“);pr intf Cn);*/將所求矩陣移到B中for (i =0;i n;i+)for(j=n;j2*n;j+)Bij-n =A2ij;)顯示所求的逆矩