fortran語言矩陣求逆,跪求一個Fortran語言下的大型稀疏矩陣求逆和矩陣乘積的快速演算法,請各位大俠幫忙啊

2021-12-25 01:38:35 字數 2268 閱讀 2528

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 1 4。求A的逆矩陣求A的特徵值和特徵向量

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的逆矩陣,A的轉置求逆為什麼等於A的求逆的轉置矩陣??

使用初等行變換的方法求逆矩陣 即用行變換把矩陣 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...