机器学习-RBF高斯核函数处理

 机器学习-RBF高斯核函数处理

SVM高斯核函数-RBF优化 
重要了解数学的部分： 
协方差矩阵，高斯核函数公式。 
个人建议具体的求法还是看下面的核心代码吧，更好理解，反正就我个人而言，烦躁的公式，还不如一段代码来的实际。本来想用Java的一个叫jblas的矩阵包，但是想了想，还是自己动手写一下吧。加深一下自己理解。实现的语言用的是java孪生兄弟Scala。我想应该不难懂。矩阵变换用二位数组将就。 
以下代码建议用scala命令行调试

核心代码

def TransposedMatrix(a:Array[Array[Double]]):Array[Array[Double]]={//返回转置矩阵

    val length=a.length

    val width=a(0).length

    var TransposedM =Array.ofDim[Double](width,length)

    for(i <- 0 to length-1)

      for(j <- 0 to width-1) {

        TransposedM(j)(i) = a(i)(j)

      }

    TransposedM

  }

  def R1(i:Int,j:Int,M:Array[Array[Double]]):Double={//R(i,j)=(第i列-第j列)*[（第i列-第j列）转置]

    var sum:Double=0

    for(x <- 0 to M.length-1)

      sum+=Math.pow((M(x)(i)-M(x)(j)),2)

    sum

  }

  def RowAverage(a:Array[Array[Double]]):Array[Double]={//返回列的均值，返回一个列矩阵

    val length=a.length

    val width=a(0).length

    var b=new Array[Double](width)

    for(i <- 0 to width-1)

        for(j <- 0 to length-1)

          b(i) += a(j)(i)

    for(i <- 0 to width-1)

      b(i)=b(i)/length

    b

  }

  def sumOfRow(a:Array[Array[Double]]):Array[Double]={//返回矩阵列的和，返回一个列矩阵

  val length=a.length

    val width=a(0).length

    var b=new Array[Double](width)

    for(i <- 0 to width-1)

      for(j <- 0 to length-1)

        b(i) += a(j)(i)

    for(i <- 0 to width-1)

      b(i)=b(i)

    b

  }

  def sum(i:Int,j:Int,a:Array[Array[Double]]):Double={//i列乘j列的转置

    var result:Double=0

    for(x<- 0 to a.length-1)

      result +=(a(x)(i)*a(x)(j))

    result

  }

  def cov(a:Array[Array[Double]]):Array[Array[Double]]={//将特征矩阵作为参数，返回协方差矩阵

    val m1=TransposedMatrix(a)

    val m2=RowAverage(m1)

    val m3=datasort(m1,m2)//将矩阵中心化

    val width=m3(0).length

    var b =Array.ofDim[Double](width,width)

    for(i <- 0 to width-1)

      for(j <- 0 to width-1)

        b(i)(j)=sum(i,j,m3)

    b

  }

  def datasort(a:Array[Array[Double]],b:Array[Double]):Array[Array[Double]]={//矩阵中心化，将每列减去列的均值

          for(i <- 0 to a(0).length-1)

            for(j <- 0 to a.length-1)

              a(j)(i) -= b(i)

    a

  }

  def gaussMatrix(a:Array[Array[Double]],delta:Array[Double]):Array[Array[Double]]={//a为特征矩阵，delta为协方差矩阵列之和，返回高斯核函数矩阵

    val b=TransposedMatrix(a)

    val length=b(0).length

    var R =Array.ofDim[Double](length,length)

    for(i <- 0 to length-1)

      for(j <- 0 to length-1)

        R(i)(j)=Math.exp(-R1(i,j,b)/delta(j))

    R

  }

val test=Array(Array(2.0, 8.0), Array(3.0, 6.0), Array(9.0, 2.0))

val test2=cov(test)

val rowOfsum=sumOfRow(res65)

gaussMatrix(test,rowOfsum)

欢迎各位看官大爷批评指教。

感谢下面百度知道回复的朋友，实现的代码段很多得到他的启示。 
http://zhidao.baidu.com/link?url=-u5LznclWQ0LbvEx3DB8sofohyP7nJCWws78TsWBNaDR15rDn-7ENoRealHRIM8W8ycioegl_NGAFzQJ33PbZ90ACQQ7eLf8HgR7DAQUJjS