感知机算法的对偶形式

4,算法代码

/*************************************************************************
    > File Name: durl.cpp
    > Author:zhangtx
    > Mail: 18510665908@163.com 
    > Created Time: 2018年01月08日 星期一 10时04分10秒
 ************************************************************************/
#include <iomanip>
#include<iostream>
using namespace std;
int a[3]={0,0,0};
int b=0;
int ep=1;
int samples[3][3]={
    {
        3,3,0
    },
    {
        4,3,0
    },
    {
        1,1,0
    }

};
int y[3]={1,1,-1};
int TransSamples[3][3];
int gram[3][3];
void TransposeMatrix(int matrixSrc[][3],int matrixDest[][3],int row,int col)
{
    for(int idxRow=0;idxRow<row;idxRow++)
        for(int idxCol=0;idxCol<col;idxCol++)
            matrixDest[idxCol][idxRow]=matrixSrc[idxRow][idxCol];
}

void MultiMatrix(int matrixLeft[][3],int leftRow,int leftCol,int matrixRight[][3],int rightRow,int rightCol,int result[][3])
{
    for(int idxLeftRow=0;idxLeftRow<leftRow;idxLeftRow++)
    {
        for(int idxRightCol=0;idxRightCol<rightCol;idxRightCol++)
        {

            int sum=0;
            for(int idxRightRow=0;idxRightRow<rightRow;idxRightRow++)
            {
                sum+=matrixLeft[idxLeftRow][idxRightRow]*matrixRight[idxRightRow][idxRightCol];
            }
            result[idxLeftRow][idxRightCol]=sum;
        }
    }
}
void PrintMatrix(int matrix[][3],int row,int col,char *info)
{
    cout<<endl<<info<<endl;
    for(int idxRow=0;idxRow<row;idxRow++)
    {
        cout<<endl;
        for(int idxCol=0;idxCol<col;idxCol++)
        {
            cout<<setw(8)<<matrix[idxRow][idxCol];
        }
    }
    cout<<endl;
}

int w[2]={
    0,0
};

int main(int argc,char *argv[])
{
    PrintMatrix(samples,3,3,"samples");
    TransposeMatrix(samples,TransSamples,3,3);
    PrintMatrix(TransSamples,3,3,"Transposes");
    MultiMatrix(samples,3,3,TransSamples,3,3,gram);
    PrintMatrix(gram,3,3,"result");

    /*train*/
    int iterator=0;
    while(true)
    {
        bool over=true;

        iterator++;
        for(int i=0;i<3;i++)
        {
            int val=0;
            for(int j=0;j<3;j++)
            {
                val+=a[j]*y[j]*gram[j][i];
            }
            
            if (y[i]*(val+b)<=0)
            {
                cout<<" iterator " <<iterator<<"a:";
                for(int k=0;k<3;k++)
                {
                    cout<<setw(8)<<a[k];
                }
                cout<<"   b"<<setw(10)<<b<<endl;
                a[i]=a[i]+ep;
                b=b+y[i];
                over=false;
            }
        }

        if (over==true)
        {

            cout<<" iterator " <<iterator<<"a:";
            for(int i=0;i<3;i++)
            {
                cout<<setw(8)<<a[i];
            }
            cout<<"   b"<<setw(10)<<b<<endl;
            break;
        }
    }

    for(int i=0;i<2;i++)
    {
        for(int j=0;j<3;j++)
            w[i]+=a[j]*y[j]*samples[j][i];
    }
    cout<<endl<<"w="<<w[0]<<"  "<<w[1]<<endl;

    cout<<"the func is "<<w[0]<<"*x1"<<"+"<<w[1]<<"*x2"<<"+("<<b<<")"<<endl;
    return 0;
}