气象常用程序包-相当给力！

☆涅垚YTZ☆

果断分享了啊

Sun

好东西  分享学习

cuijxin

第三章 多元统计分析里有一小节3.4.2.2多级判别的Fisher准则概念 我想把这个算法转化成C++实现，但是不是气象专业，高数 统计等基础知识也放下很久了 有没有人可以帮我分析一下这段Fortran代码啊 主要是算完D阵的特征值和特征向量之后

ALLOCATE(AA(NP))
        ALLOCATE(X2(NP))
        ALLOCATE(X22(NP))
        NP1=0
        DO I=1,NP
          AA1=1
          DO J=I,NP
            AA(J)=AA1/(1+DT(J))
            AA1=AA(J)
          END DO
          X2(I)=-(N-1-0.5*(M+G))*LOG(AA(I))
          N0=(M+1-I)*(G-I)        !检验的自由度
          ALF=0.05   !检验的置信度
          H=-0.005         !检验时的积分步长
          NN=(1.-ALF)/ABS(H)
          CALL EULER1(H,NN,N0,X22(I))  !哑元依次为积分步长、积分步数、自由度及χ2值
          IF(X2(I)>X22(I))THEN
            NP1=I
          ELSE
            GOTO 30
          END IF
        END DO
30        CONTINUE
        WRITE(12,'("NP1=",I2)')NP1
        WRITE(12,'(5X," χ2检验 ")')
        WRITE(12,'("X2=",<NP1>F8.4)')(X2(I),I=1,NP1)
        WRITE(12,'("X22="<NP1>F8.4)')(X22(I),I=1,NP1)
!   经判别,有NP1个判别函数显著
!        计算每一样品与各组判别函数重心的距离
        ALLOCATE(DS(N,G,NP1))
        DO L=1,NP1
          DO I=1,N
            DO J=1,G
              DO K=1,M
                XX1(K)=X(I,K)-XV(K,J)
              END DO
                  XX2=0
                  XX3=0
                  DO K=1,M
                    XX2=XX2+XX1(K)*V(K,L)
                    XX3=XX3+V(K,L)*XX1(K)
              END DO
                  DS(I,J,L)=XX2*XX3
            END DO
          END DO
        END DO
        ALLOCATE(KGJ(N,NP1))
        DO L=1,NP1
          DO I=1,N
            DMIN=1.0E+30
            DO J=1,G
              IF(DS(I,J,L)<DMIN)THEN
                    DMIN=DS(I,J,L)
                    KGJ(I,L)=J
                  ENDIF
            END DO
          END DO
    END DO
        WRITE(12,'(4X,"样品与判别函数重心的距离")')
        DO L=1,NP1
            DO J=1,G
            WRITE(12,'(" 组别=",I2,<N>F6.3)')J,(DS(I,J,L),I=1,N)
          END DO
          WRITE(12,'("最后分组",<N>I6)')(KGJ(I,L),I=1,N)
        END DO
        DEALLOCATE(AA)
        DEALLOCATE(X2)
        DEALLOCATE(X22)
        DEALLOCATE(DS)
        DEALLOCATE(KGJ)
        END

cuijxin

还有程序里的 F函数 欧拉函数 还有GAMMA函数最好也能解释一下啊

远山

好东西，顶一下

majunlu

果断分享了啊

永远的爱

我是刚刚注册的心用户，在此发现好多很好的资料，非常感谢。

永远的爱

怎么不能购买，总是提示要登录啊！可是我是已经登录了啊！

jodia

正好需要使用谱分析，谢谢分享~~

cnphoenix

给力，谢谢楼主