版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

dongming198

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1421.7
帖子: 122
在線: 79.5小時
蟲號: 3411768
注冊: 2014-09-12
性別: GG
專業(yè): 內(nèi)流流體力學(xué)

[求助] LU分解法求解線性方程組FORTRAN程序代碼已有1人參與

如題，請詳細(xì)注釋一下，最好給個算例，謝謝！
另外，LU分解法對線性方程組的系數(shù)矩陣有要求嗎？
我主要是想要一個通用的求解線性方程組程序，其它迭代方法也可以的。
可以用網(wǎng)盤

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有179人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有72人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

相場法fortran源代碼已經(jīng)有166人回復(fù)
有限差分法化之后化簡為線性方程組已經(jīng)有50人回復(fù)
FORTRAN90標(biāo)準(zhǔn)函數(shù)及IMSL標(biāo)準(zhǔn)函數(shù)庫已經(jīng)有166人回復(fù)
跪求求解線性方程組的c語言程序已經(jīng)有3人回復(fù)
用fortran編的用共軛梯度法解決線性方程組的問題。已經(jīng)有5人回復(fù)
fortran real(kind=16) 已經(jīng)有10人回復(fù)
關(guān)于應(yīng)用fortran編寫高斯消去法程序來求解方程組的問題已經(jīng)有5人回復(fù)
如何求解帶有固定參數(shù)的超定方程組？已經(jīng)有21人回復(fù)
用fortran程序遺傳算法解非線性方程組已經(jīng)有7人回復(fù)
這個方程組怎么解啊已經(jīng)有14人回復(fù)
fortran中計算線性積分已經(jīng)有7人回復(fù)
關(guān)于高等代數(shù)線性方程組同解的問題已經(jīng)有5人回復(fù)
求fortran解非線性方程組已經(jīng)有4人回復(fù)
大家都用什么fortran編輯器，推薦一下，最好是能加亮或修改語句顏色的已經(jīng)有15人回復(fù)
【求助】fortran求解矩陣已經(jīng)有7人回復(fù)
【求助】如何寫FORTRAN程序?qū)崿F(xiàn)求平均最近鄰距離已經(jīng)有34人回復(fù)
【求助】用fortran怎么描述下面這樣的情況，我是新手，請大家?guī)兔?/a> 已經(jīng)有8人回復(fù)

【求助】用最小二乘法求解線性方程組的Fortran代碼已經(jīng)有4人回復(fù)
【求助】用fortran求解大型線性方程組時出現(xiàn)的錯誤【已解決】已經(jīng)有11人回復(fù)
【求助】求解非線性方程已經(jīng)有8人回復(fù)

1樓 2015-05-30 17:59:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
計算強帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
dongming198: 金幣+3, ★有幫助, 這個我有的，不過驗證的時候不正確。現(xiàn)在不需要這個啦，還是很感謝啦 2015-06-04 20:48:30
ben_ladeng: 金幣+2, 謝謝指教 2015-06-05 01:00:03

http://blog.chinaunix.net/uid-479984-id-2114914.html

CODE:

Author:RainMan

!2009-12-20

!LU分解法解線性方程組AX=B

!該程序從文件輸入，格式為

!矩陣的維數(shù)是 3

!A矩陣

!4 2 1

!8 7 4

!16 17 15

!B矩陣

!11 34 95 

PROGRAM solve_equation

    REAL,DIMENSION(:,:),ALLOCATABLE::A

    REAL,DIMENSION(:),ALLOCATABLE::B

    CHARACTER*(80):: inputfile = 'dengjin.txt'

    INTEGER DIM

    CALL GET_DIM(inputfile,DIM)

    ALLOCATE(A(DIM,DIM),B(DIM))

    CALL INPUT_MATRIX(A,B,DIM,inputfile)

    CALL LU(A,DIM)

    CALL solve(A,B,B,DIM)

    PRINT*,B

END

SUBROUTINE GET_DIM(inputfile,n)

    CHARACTER*(*) inputfile

    CHARACTER*(80) str

    INTEGER n

    OPEN(10,FILE = inputfile)

    READ(10,*)str,n

    CLOSE(10)

END

SUBROUTINE INPUT_MATRIX(A,B,DIM,inputfile)

    CHARACTER*(*) inputfile

    INTEGER DIM

    INTEGER A(DIM,DIM),B(DIM)

    OPEN(20,FILE=inputfile)

    READ(20,*)

    READ(20,*)

    READ(20,*) ((A(I,J),J=1,DIM),I=1,DIM)

    READ(20,*)

    READ(20,*) (B(I),I=1,DIM)

    CLOSE(20)

END

SUBROUTINE LU(A,DIM)

    INTEGER DIM

    REAL A(DIM,DIM)

    DO I=2,DIM

        DO J=I,DIM

            A(J,I-1) = A(J,I-1)*1.0/A(I-1,I-1)

            m = A(J,I-1)

            DO K=I,DIM

                A(J,K) = A(J,K)-m*A(I-1,K)

            ENDDO

        ENDDO

    ENDDO

END

SUBROUTINE solve(A,B,X,DIM)

    INTEGER:: DIM

    REAL A(DIM,DIM),B(DIM),X(DIM)

!----------LY = B 

    DO I=2,DIM

        DO J=1,I-1

            B(I) = B(I) - A(I,J)*B(J)

        ENDDO

    ENDDO

    B(DIM) = B(DIM)/A(DIM,DIM)

!---------UX=Y Y stores in B,same as X

    DO I=DIM-1,1,-1

        DO J=I+1,DIM

            B(I) = B(I) - A(I,J)*B(J)

        ENDDO

        B(I) = B(I)/A(I,I)

    ENDDO

END

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2樓2015-06-03 14:27:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 dongming198 的主題更新