距离判别源代码清单详解.pdfVIP

距离判别源代码清单距离判别源代码清单 距离判别源代码清单距离判别源代码清单 Option Base 1 ****************************************************************************** ************** 函数部分 ****************************************************************************** ************** Sub putout(A() As Double, m%, n%) m 为矩阵的行，n 为矩阵的列 ****************************************************************************** ********** 功能 ： 输出矩阵 参数 ： a 为待输出的矩阵；m 为a 矩阵的行；n 为a 矩阵的列 返回值： 无返回值 ****************************************************************************** ********** For i = 1 To m For j = 1 To n Print Tab(j * 20); A(i, j); Format(a(i, j), 0.0000); Next j Print Next i End Sub Sub pingjun(A() As Double, m%, qd%, zd%, pj() As Double) ****************************************************************************** ********** 功能 ： 求m 列矩阵从第qd 列到第zd 列的平均 参数 ： a 为待求矩阵；m 为矩阵的行数；qd 为起始列下标；zd 终止列下标 返回值： pj 返回最终结果，是一m 为列向量 ****************************************************************************** ********** ReDim pj(m) As Double For i = 1 To m For j = qd To zd pj(i) = pj(i) + A(j, i) Next j pj(i) = pj(i) / (zd - qd + 1) Next i End Sub Sub licha(A() As Double, m%, qd%, zd%, lc() As Double) ****************************************************************************** ********** 功能 ： 求m 列矩阵从第qd 列到第zd 列的离差 参数 ： a 为待求矩阵；m 为矩阵的行数；qd 为起始列下标；zd 终止列下标 返回值： lc 返回最终结果，是一m 为列向量 ****************************************************************************** ********** ReDim lc(m, m) As Double 原来的lc 矩阵可能已存入非零数据，所以要重新声明使其数据均 为零 Dim pj1() As Double Call pingjun(A(), m, qd, zd, pj1()) For r = 1 To m For t = 1 To m For j = qd To zd sum = sum + (A(j, t) - pj1(t)) * (A(j, r) - pj1(r)) Next j lc(t, r) = sum sum = 0 Next t, r End Sub Sub chengfa(n() As Double,