smo

2020-01-07 23:58:00 阅读：376 来源： 互联网

还在更新中


import numpy as np
import random

# build auxiliary function

class optstruct:
    def __init__(self, datax, datay, C, toler):
        self.x = datax
        self.label = datay
        self.c = C
        self.tol = toler
        self.m = datax.shape[0]  # how many samples are there
        self.alphas = np.mat(np.zeros((self.m, 1)))  # initialize alphas
        self.b = 0  # initialize b
        self.ercache = np.mat(np.zeros((self.m, 2)))

def calcek(os,k):
    fxk=float(np.multiply(os.alphas,os.label).T*(os.x*os.x[k,:].T))+os.b # fx=wx+b,w=alphasi*yi*xi sum from n
    Ek=fxk-float(os.label[k])  # Ek=yi-labeli
    return Ek



def pickJ(i,os,Ei):
    maxk=-1; maxdeltaE=0;Ej=0
    os.ercache[i]=[1,Ei]
    valerchache=np.nonzero(os.ercache[:,0].A)[0] # work for .A matrix convert array  # index value
    if (len(valerchache))>1: # if having value , do it
        for k in valerchache: # traversal index
            if k==i:
                continue
            Ek=calcek(os,k)
            deltas=abs(Ei-Ek)
            if deltas>maxdeltaE:  # heuristic algorithm
                maxk=k
                maxdeltaE=deltas
                Ej=Ek
        return maxk,Ej
    else:
        j=i
        while(j==i):
            j=int(random.uniform(0,os.m))
        Ej=calcek(os,j)
        return j,Ej
def updateEk(os,k):
    Ek=calcek(os,k)
    os.ercache[k]=[1,Ek]

标签：calcek,Ek,self,smo,np,os,ercache
来源： https://www.cnblogs.com/tangjunjun/p/12164363.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

smo