荔园在线

荔园之美,在春之萌芽,在夏之绽放,在秋之收获,在冬之沉淀

[回到开始] [上一篇][下一篇]


发信人: 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软件 网络书店