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樓:匿名使用者 希望能幫到你 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 ...c語言程式設計求任意對稱正定矩陣的逆
分塊矩陣求逆矩陣有哪些公式,分塊矩陣的逆矩陣怎麼求
一階矩陣的逆矩陣怎麼求,二階矩陣逆矩陣的公式是哪個