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?????????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]?>?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）?????????????????????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個文件信息