Originally posted by brqhl_ing at 2010-08-09 18:33:26:

!响应系统的时间序列写出
         
        delta=((z(1,i)-z(4,i))**2.0+(z(2,i)-z(5,i))**2.0+
     &(z(3,i)-z(6,i))**2.0)**0.5

Originally posted by sxf2012 at 2010-08-11 14:59:02:
这儿的z（1,i）中的z是个函数么？没有定义！！！

Originally posted by sxf2012 at 2010-08-11 19:07:20:
那需要另写一个子函数呀，：误差函数子函数function Z(*,*)

FUNCTION MERRF(X)
        DOUBLE PRECISION MERRF,X,MGAM2
        IF (X.GE.0.0) THEN
          MERRF=MGAM2(0.5D0,X*X)
        ELSE
          MERRF=-MGAM2(0.5D0,X*X)
        END IF
        RETURN
        END


        FUNCTION MGAM2(A,X)
        DOUBLE PRECISION MGAM2,A,X
        DOUBLE PRECISION MGAM1,P,Q,D,S,S1,P0,Q0,P1,Q1,QQ
        IF ((A.LE.0.0).OR.(X.LT.0.0)) THEN
          IF (A.LE.0.0) THEN
            WRITE(*,*) '  ERR** A<=0!'
          END IF
          IF (X.LT.0.0) THEN
            WRITE(*,*) '  ERR** X<0!'
          END IF
          MGAM2=-1.0
        END IF
        IF (X+1.0.EQ.1.0) THEN
          MGAM2=0.0
          RETURN
        END IF
        IF (X.GT.1.0D+35) THEN
          MGAM2=1.0
          RETURN
        END IF
        Q=LOG(X)
        Q=A*Q
        QQ=EXP(Q)
        IF (X.LT.1.0+A) THEN
          P=A
          D=1.0/A
          S=D
          DO 10 N=1,100
            P=1.0+P
            D=D*X/P
            S=S+D
            IF (ABS(D).LT.ABS(S)*1.0D-07) THEN
              S=S*EXP(-X)*QQ/MGAM1(A)
              MGAM2=S
              RETURN
            END IF
10          CONTINUE
        ELSE
          S=1.0/X
          P0=0.0
          P1=1.0
          Q0=1.0
          Q1=X
          DO 20 N=1,100
            P0=P1+(N-A)*P0
            Q0=Q1+(N-A)*Q0
            P=X*P0+N*P1
            Q=X*Q0+N*Q1
            IF (ABS(Q)+1.0.NE.1.0) THEN
              S1=P/Q
              P1=P
              Q1=Q
              IF (ABS((S1-S)/S1).LT.1.0D-07) THEN
                S=S1*EXP(-X)*QQ/MGAM1(A)
                MGAM2=1.0-S
                RETURN
              END IF
              S=S1
            END IF
            P1=P
            Q1=Q
20          CONTINUE
        END IF
        WRITE(*,*) ' A too large !'
        S=1.0-S*EXP(-X)*QQ/MGAM1(A)
        MGAM2=S
        RETURN
        END



Originally posted by sxf2012 at 2010-08-12 19:09:39:
从徐士良那本书上看到的：把你程序中的Z(*,*)改成MGAM2(*,*)就该可以了吧，不知道第一个子函数MERRF(X),你那上边用不用得到？
[code]
FUNCTION MERRF(X)
        DOUBLE PRECISION MERRF,X,MGAM2
        IF (X.GE.0.0) THE ...

Originally posted by brqhl_ing at 2010-08-22 18:49:12:

这个是解决那个误差函数的吗？作为一个子函数来调用的吗》？整体的调试了没有？最近家里有点事！

Originally posted by sxf2012 at 2010-08-28 13:41:19:


我不知道你的误差函数是不是通常所指 的那个:
erf(x)=\frac{2}{\pi}\int ^x_0{exp{-t^2}dt}