svm手写数字识别

svm手写数字识别

from numpy import *


def loadDataSet（fileName）:
    dataMat = []
    labelMat = []
    fr = open（fileName）
    for line in fr.readlines（）:
        lineArr = line.strip（）.split（‘\t‘）
        dataMat.append（[float（lineArr[0]） float（lineArr[1]）]）
        labelMat.append（float（lineArr[2]））
    return dataMat labelMat


def selectJrand（i m）:
    j = i
    while （j == i）:
        j = int（random.uniform（0 m））
    return j


# 调整alpha值
def clipAlpha（aj H L）:
    if aj > H:
        aj = H
    if L > aj:
        aj = L
    return aj


def smoSimple（dataMatIn classLabels C toler maxIter）:
    dataMatrix = mat（dataMatIn）
    labelMat = mat（classLabels）.transpose（）
    b = 0
    m n = shape（dataMatrix）
    alphas = mat（zeros（（m 1）））
    iter = 0
    while （iter < maxIter）:
        alphaPairsChanged = 0
        for i in range（m）:
            fXi = float（multiply（alphas labelMat）.T * （dataMatrix * dataMatrix[i :].T）） + b
            Ei = fXi - float（labelMat[i]）
            if （（labelMat[i] * Ei < -toler） and （alphas[i] < C）） or （（labelMat[i] * Ei > toler） and （alphas[i] > 0））:
                j = selectJrand（i m）  # 从0到m选出不等于i的随机数
                fXj = float（multiply（alphas labelMat）.T * （dataMatrix * dataMatrix[j :].T）） + b
                Ej = fXj - float（labelMat[j]）
                alphaIold = alphas[i].copy（）
                alphaJold = alphas[j].copy（）
                if （labelMat[i] != labelMat[j]）:
                    L = max（0 alphas[j] - alphas[i]）
                    H = min（C C + alphas[j] - alphas[i]）
                else:
                    L = max（0 alphas[j] + alphas[i] - C）
                    H = min（C alphas[j] + alphas[i]）
                if L == H:
                    print（“L==H“）
                    continue
                eta = 2.0 * dataMatrix[i :] * dataMatrix[j :].T \
                      - dataMatrix[i :] * dataMatrix[i :].T \
                      - dataMatrix[j :] * dataMatrix[j :].T
                if eta >= 0:
                    print（“eta>=0“）
                    continue
                alphas[j] -= labelMat[j] * （Ei - Ej） / eta
                alphas[j] = clipAlpha（alphas[j] H L）
                if （abs（alphas[j] - alphaJold） < 0.00001）:
                    print（“j not moving enough“）
                    continue
                alphas[i] += labelMat[j] * labelMat[i] * （alphaJold - alphas[j]）
                b1 = b - Ei - labelMat[i] * （alphas[i] - alphaIold） * dataMatrix[i :] * dataMatrix[i :].T - labelMat[
                                                                                                                  j] * （
                                                                                                                  alphas[
                                                                                                                      j] - alphaJold） * dataMatrix[
                          

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-02-10 16:05  test5\
     目录           0  2017-02-10 23:25  test5\.idea\
     目录           0  2017-02-07 15:03  test5\.idea\inspectionProfiles\
     文件         228  2017-02-07 15:03  test5\.idea\inspectionProfiles\profiles_settings.xml
     文件         226  2017-02-07 15:06  test5\.idea\misc.xml
     文件         262  2017-02-07 15:03  test5\.idea\modules.xml
     文件         398  2017-02-07 15:06  test5\.idea\test5.iml
     文件       25948  2017-02-10 23:25  test5\.idea\workspace.xml
     文件       14260  2017-02-10 16:05  test5\svmMLiA.py
     目录           0  2017-02-10 16:04  test5\testDigits\
     文件        1088  2010-10-07 05:35  test5\testDigits\1_0.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_1.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_10.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_11.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_12.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_13.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_14.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_15.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_16.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_17.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_18.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_19.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_2.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_20.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_21.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_22.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_23.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_24.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_25.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_26.txt
     文件        1088  2010-10-07 05:35  test5\testDigits\1_27.txt
............此处省略573个文件信息