荔园在线
荔园之美,在春之萌芽,在夏之绽放,在秋之收获,在冬之沉淀
[回到开始]
[上一篇][下一篇]
发信人: Nol (天书), 信区: CACE
标 题: 结构电算
发信站: 荔园晨风BBS站 (Thu Apr 4 13:14:20 2002), 转信
阿f你这次再删我的劳动成果,你就成为人民的公敌啦。嘿嘿~~~
第一题:高斯消去法(考虑A(i,i)=0的情况)
REAL*8 A(N,N),B(N),C(N),D
READ(*,*) N
DO I=1,N
DO J=1,N
READ(*,*) A(I,J)
END DO
READ(*,*) B(I)
END DO
DO I=1,N
IF(A(I,I).EQ.0.0) THEN
DO J=1,N
IF((A(J,I).NE.0.0).AND.(A(I,J).NE.O.0))THEN
DO K=1,N
C(K)=A(I,K)
A(I,K)=A(J,K)
A(J,K)=C(K)
D=B(I)
B(I)=B(J)
B(J)=D
END DO
GOTO 10
END IF
10 END DO
END IF
END DO
C 下面的程序与原来的一样
DO I=1,N
DO J=1+I,N
A(I,J)=A(I,J)/A(I,I)
END DO
B(I)=B(I)/A(I,I)
DO J=1+I,N
DO K=1+I,N
A(J,K)=A(J,K)-A(J,I)*A(I,K)
END DO
B(J)=B(J)-A(J,I)*B(I)
END DO
END DO
DO I=N-1,1
DO J=I+1,N
B(I)=B(I)-A(I,J)*B(J)
END DO
END DO
第二题:高斯消去法矩阵求逆(也考虑A(i,i)=0的情况)
考虑A(i,i)=0,同样可以用第一题的算法,先把矩阵的行调整成没有A(i,i)=0的情
况,然后再用上节课结果的算法运算。懒得写了。
--
※ 修改:·Nol 於 Apr 4 13:14:57 修改本文·[FROM: 192.168.41.27]
※ 来源:·荔园晨风BBS站 bbs.szu.edu.cn·[FROM: 192.168.41.27]
[回到开始]
[上一篇][下一篇]
荔园在线首页 友情链接:深圳大学 深大招生 荔园晨风BBS S-Term软件 网络书店