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

時間 2021-12-25 01:29:00

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樓:匿名使用者

希望能幫到你

c語言程式設計求任意對稱正定矩陣的逆

ifndef matrix dot h define matrix dot h template void swap t a,t b class cmatrix 返回矩陣行數 int getcol const 返回矩陣列數 int rowswap int x,int y 行交換,成功返回1,否則0 ...

分塊矩陣求逆矩陣有哪些公式,分塊矩陣的逆矩陣怎麼求

喵小採 1 a00bxa 1 00b 1 aa 1 00 a0 0b 1 0a 1 0b 00 bb 1 對於加法,相容要求兩個矩陣按同樣的方式分塊 而對於乘法,在矩陣a與矩陣b相乘時,對b的一個分塊方式,a可以有幾種分塊方式與之相容,這時便要考慮哪種分塊方式使運算更加簡便。擴充套件資料 一 總結應...

一階矩陣的逆矩陣怎麼求,二階矩陣逆矩陣的公式是哪個

一階矩陣的行列式就是其元素值 不需要證明,就是定義 其逆矩陣的元素值就是他元素值的倒數 也不需要證明,a a 1 就可以看出。到底應該怎麼樣去求逆矩陣才好呢?親,您好,我是你的答題小老師,正在為你整理題目的答案,請您稍等片刻。親,最簡單的辦法是用增廣矩陣。如果要求逆的矩,陣是a,則對增廣矩陣 ae ...