1樓:匿名使用者
沒有必要自己編吧。
1 呼叫imsl/mkl等現成的數學庫,比如.i.就可以直接求逆。
2 直接copy徐士良那本< 2樓:心隨神飛 ! aa為原矩陣,b為存放aa的逆矩陣,n為矩陣aa的維數subroutine nizhen(aa,b,n)integer n,i,j,k real:: aa(n,n),b(n,n),a(n,n)a=aa do i=1,n b(i,i)=1 enddo do i=1,n b(i,:)=b(i,:)/a(i,i) a(i,i:n)=a(i,i:n)/a(i,i)do j=i+1,n do k=1,n b(j,k)=b(j,k)-b(i,k)*a(j,i)enddo a(j,i:n)=a(j,i:n)-a(i,i:n)*a(j,i)enddo enddo do i=n,1,-1 do j=i-1,1,-1 do k=1,n b(j,k)=b(j,k)-b(i,k)*a(j,i)enddo enddo enddoend 跪求一個fortran語言下的大型稀疏矩陣求逆和矩陣乘積的快速演算法,請各位大俠幫忙啊
20 3樓:冽十三 可以參照 徐士良 fortran常用演算法程式集,我記得裡面有你需要的。 跪求fortran上矩陣相乘.矩陣求逆的程式 4樓: cc 矩陣乘法範例 c by perng 1997/9/17program matmul_demo implicit none integer n parameter(n=3) integer a(n,n) ! matrix ainteger b(n,n) ! matrix binteger c(n,n) ! matrix cdata b /1,2,3,4,5,6,7,8,9/data c /9,8,7,6,5,4,3,2,1/call matmul(a,b,n,n,c,n,n)write(*,*) 'matrix a:' call output(a,n) stop endc c 輸出矩陣的子程式 csubroutine output(a,n)implicit none integer n,a(n,n) integer i,j character for*20 data for /'(??(1x,i3))'/c 用字串來設定輸出格式 write( for(2:3), '(i2)' ) ndo i=1,n write( *, fmt=for ) (a(i,j),j=1,n)end do return endc c 矩陣乘法的子程式 csubroutine matmul(a,b,br,bc,c,cr,cc) implicit none integer br ! row of matrix binteger bc ! column of matrix b integer b(br,bc) ! matrix binteger cr ! row of matrix cinteger cc ! column of matrix c integer c(cr,cc) ! matrix cinteger a(br,cc) ! matrix ainteger i,j,k ! 迴圈的計數器! bc若不等於cr, 這兩個矩陣無法相乘if ( bc .ne. cr ) thenwrite(*,*) 'matrix size error!' stop end if do i=1,br do j=1,cc a(i,j)=0 do k=1,bc a(i,j)=a(i,j)+b(i,k)*c(k,j)end do end do end do returnend 5樓:匿名使用者 希望能幫到你 a 1 2 3 1 3 1 6 1 6 a e 1 2 1 4 1 4 2 2 5 6 2 3 所以a的特徵值為 2,3.a 2e x 0 的基礎解係為 a1 2,1 t.所以a的屬於特徵值2的特徵向量為 k1a1 k1 2,1 t,其中k1為任意非零常數.a 3e x 0 的基礎解係為 a2 1,... 矩陣的加號逆是因為矩陣不是n n的,所以不能求出標準逆,加號逆就是 a a a 也就是廣義逆 其中a 表示a的轉置 如果a是m n的矩陣 加號逆 叫做矩陣的廣義逆。傳統的矩陣求逆,必須滿足兩個專條件 屬 1.矩陣是方陣 2.矩陣可逆 這樣的矩陣才能求逆。不滿足這兩個條件的時候,矩陣就沒有逆,但是有所... 使用初等行變換的方法求逆矩陣 即用行變換把矩陣 a,e 化成 e,b 的形式,那麼b就等於a的逆在這裡 a,e 1 2 3 1 0 0 2 2 1 0 1 0 3 4 3 0 0 1 r2 2r1,r3 3r1 1 2 3 1 0 0 0 2 5 2 1 0 0 2 6 3 0 1 r1 r2,r2...已知矩陣A 1 2 1 4。求A的逆矩陣求A的特徵值和特徵向量
矩陣的加號逆是什麼意思?如何求矩陣的加號逆
求A的逆矩陣,A的轉置求逆為什麼等於A的求逆的轉置矩陣??