跪求fortran上矩陣相乘 矩陣求逆的程式

時間 2021-08-30 09:46:49

1樓:

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

2樓:匿名使用者

希望能幫到你

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

3樓:冽十三

可以參照 徐士良 fortran常用演算法程式集,我記得裡面有你需要的。

如何用 fortran 輸出矩陣啊

4樓:匿名使用者

不是輸出的問題,是初始化的問題。fortran是列優先的陣列儲存策略,你初始化的矩陣就是輸出的樣子。

5樓:匿名使用者

舉例:用baifortran輸出每行10個數的du10*10矩陣:

implicit none

integer,parameter:: n=10integer a(n,n),i,j

open(10,file="out.txt")do i=1,n

do j=1,n

a(i,j)=(i-1)*10+j-1

end do

end do

write(10,"(10i5)")(a(i,:),i=1,n)end

fortran簡介:

fortran源自

zhi於「公式翻譯」dao(英語:formulatranslation)的縮寫,是一版種程式語言權。它是世界上最早出現的計算機高階程式設計語言,廣泛應用於科學和工程計算領域。

fortran語言以其特有的功能在數值、科學和工程計算領域發揮著重要作用。fortran 90之前的版本是人們所知曉的fortran(全部字母大寫),從fortran 90以及以後的版本都寫成fortran(僅有第一個字母大寫)。

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

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

行矩陣的逆矩陣怎麼求,n行1列矩陣怎麼求逆矩陣

雨說情感 1 伴隨矩陣法 如果矩陣a可逆,則 的餘因子矩陣的轉置矩陣。a 0,a 為該矩陣對應的行列式的值 a的伴隨矩陣為 其中aij 1 i jmij稱為aij的代數餘子式。2 初等行變換法 在行階梯矩陣的基礎上,即非零行的第一個非零單元為1,且這些非零單元所在的列其它元素都是0。綜上,行最簡型矩...

用matlab求含有未知矩陣的逆矩陣

你的矩陣也不是個方陣,怎麼求逆?隨便編了幾個數,可以求的,只要存在逆矩陣,但結果很複雜 clear all clc syms s a b c d a s,0,0,a i sqrt b 1 c 1 0,s,a i sqrt b c 0 a,b,c,d a b,c,d,a b inv a 直接計算就行 ...