感知机c/c++代码实现

具体公式百度：参看统计学方法中算法描述实现的

一，原始形式

<pre name="code" class="cpp">#include<iostream>
#include<stdlib.h>
#include<string.h>
#include<stdio.h>
using namespace std;
int main() 
{
	int x1[]={3,4,1};
	int x2[]={3,3,1};
	int y[]={1,1,-1};
	int w[]={0,0};
	int b=0;
	int index=0;
	int n=1;
	while(index<3)
	{
		int temp=y[index]*(x1[index]*w[0]+x2[index]*w[1]+b);
		if(temp<=0)
		{
			w[0]+=n*y[index]*x1[index];
			w[1]+=n*y[index]*x2[index];
			b+=n*y[index];
			index=0;
		}
		else{
			index++;
		}
	}
	cout<<"result:"<<endl;
	cout<<"w:"<<w[0]<<" "<<w[1]<<endl;
	cout<<"b:"<<b;
	getchar();
	return 0;
} 

对偶形式

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
int main()
{
	int x1[]={3,4,1};
	int x2[]={3,3,1};
	int y[]={1,1,-1};
	
	int b=0;
	int a[]={0,0,0}; 
	int G[3][3];
	//计算Gram矩阵
	cout<<"Gram矩阵为:"<<endl; 
	for(int i=0;i<3;i++)
	 for(int j=0;j<3;j++)
	  {
		 G[i][j]=x1[i]*x1[j]+x2[i]*x2[j];
		 cout<<G[i][j]<<"\t";
		 if(j==2)
		  cout<<endl;
	  }
	cout<<"开始计算a,b:"<<endl;
	int index=0;
	int n=1;
	while(index<3)
	{
		int temp=0;
		for(int j=0;j<3;j++)
		{
			temp+=a[j]*y[j]*G[j][index];
		}
		int all=y[index]*(temp+b);
		if(all<=0)
		{
			a[index]+=n;
			b+=n*y[index];
			index=0;
			
		}
		else
		{
			index++;
		}
		
	}
	cout<<"a:"<<a[0]<<" "<<a[1]<<" "<<a[2]<<endl;
	cout<<"b:"<<b;
	
	//计算w和b
	int w[]={0,0};
	for(int i=0;i<3;i++)
	{
		w[0]+=a[i]*y[i]*x1[i];
		w[1]+=a[i]*y[i]*x2[i];
	} 
	b=0;
	for(int i=0;i<3;i++)
	{
		b+=a[i]*y[i];
	}
	cout<<endl<<"w:"<<w[0]<<" "<<w[1]<<endl;
	cout<<"b:"<<b;
	getchar();
	return 0;
}