研究生期间的模式识别课程。
具体源代码：

#include <iostream>
#include<stdlib.h>
#include <math.h>
using namespace std;
int s_count;
int par_count;

class sample // 样本类s
{
private:
	int par[3];
	char type; // 样本类别 M F I 三类

public:
	void setpar(int *a)
	{
		//	par = new int[2 + 1]; 
		for (int i = 0; i < 2 + 1; i++)
			par[i] = a[i];
	}
	sample()
	{}
	void setchar(char c)
	{
		type = c;
	}
	void getpar(int *s)
	{
		for (int i = 0; i < 2 + 1; i++)
		{
			s[i] = par[i];
		}
	}
	char gettype()
	{
		return type;
	}
	~sample()
	{}
};
class vector // 向量类
{
private:
	int par[3];
public:

	vector() {}
	int getpar(int i)
	{
		return par[i];
	}
	void getpar(int *a)
	{
		for (int i = 0; i < 2 + 1; i++)
		{
			a[i] = par[i];
		}
	}
	void setpar(int *a)
	{
		//par = new int[2+1];
		for (int i = 0; i < 2 + 1; i++)
		{
			par[i] = a[i];
		}
	}
	~vector()
	{}
};
int martix_multiply(vector w, sample x)
{
	int end = 0;
	/*int *v;
	v = new int[2 + 1];
	int *s;
	s = new int[2 + 1];*/
	int v[3];
	int s[3];
	w.getpar(v);
	x.getpar(s);
	for (int i = 0; i < 2 + 1; i++)
	{
		end += v[i] * s[i];
	}
	return end;
	//delete[]v;
	//delete[]s; 
}
void martix_sub(vector &w, sample x)
{
	/*int *v;
	v = new int[2 + 1];
	int *s;
	s = new int[2 + 1];*/
	int v[3];
	int s[3];
	w.getpar(v);
	x.getpar(s);
	for (int i = 0; i < 2 + 1; i++)
	{
		v[i] = v[i] - s[i];
	}
	w.setpar(v);
	//delete []v;
	//delete []s;
}
void martix_add(vector &w, sample x)
{
	/*int *v;
	v = new int[2 + 1];
	int *s;
	s = new int[2 + 1];*/
	int v[3];
	int s[3];
	w.getpar(v);
	x.getpar(s);
	for (int i = 0; i < 2 + 1; i++)
	{
		v[i] = v[i] + s[i];
	}
	w.setpar(v);
	//delete[]v;
	//delete[]s; 
}
struct s_v
{
	sample s[100];
	vector v[3];
};
int main()
{
	s_v as;
	cin >> s_count;  // s_count : count of sample 
	cin >> par_count; // 2:count of parmeter
	int *temp;
	temp = new int[2 + 1];//三维的一个数组
	int i;
	int j;
	for (i = 0; i < s_count; i++)
	{
		for (j = 0; j < 2; j++)
		{
			cin >> temp[j];
		}//循环为每个样本中的二维点赋值
		temp[2] = 1;
		as.s[i].setpar(temp);
	}
	for (i = 0; i < 3; i++)
	{
		temp[i] = 0;
	}
	for (i = 0; i < 3; i++)
	{
		as.v[i].setpar(temp);
	}
	int d[3] = { 0,0,0 };
	int times = 0; //迭代次数
	bool charge[3];
	for (i = 0; i < 3; i++)
	{
		charge[i] = false;
	}
	while (charge[0] == false || charge[1] == false || charge[2] == false)
	{
		i = times % s_count;
		d[i % 3] = martix_multiply(as.v[i % 3], as.s[i]);
		d[(i + 1) % 3] = martix_multiply(as.v[(i + 1) % 3], as.s[i]);
		d[(i + 2) % 3] = martix_multiply(as.v[(i + 2) % 3], as.s[i]);
		if (d[i] > d[(i + 1) % 3] && d[i] > d[(i + 2) % 3])
		{
			charge[i] = true;
		}
		else if (d[i] <= d[(i + 1) % 3] && d[i] <= d[(i + 2) % 3])
		{
			martix_add(as.v[i % 3], as.s[i]);
			martix_sub(as.v[(i + 1) % 3], as.s[i]);
			martix_sub(as.v[(i + 2) % 3], as.s[i]);
			charge[0] = false;
			charge[1] = false;
			charge[2] = false;
		}
		else if (d[i] > d[(i + 1) % 3 && d[i] <= d[(i + 2) % 3]])
		{
			martix_add(as.v[i % 3], as.s[i]);
			martix_sub(as.v[(i + 2) % 3], as.s[i]);
			charge[0] = false;
			charge[1] = false;
			charge[2] = false;
		}
		else if (d[i] > d[i + 2] % 3 && d[i] <= d[(i + 1) % 3])
		{
			martix_add(as.v[i % 3], as.s[i]);
			martix_sub(as.v[(i + 1) % 3], as.s[i]);
			charge[0] = false;
			charge[1] = false;
			charge[2] = false;
		}
		times++;
	}
	for (i = 0; i < 3; i++)
	{
		cout << "[";
		for (j = 0; j < 2 + 1; j++)
		{
			if (j == 2)
			{
				cout << as.v[i].getpar(j);
			}
			else
				cout << as.v[i].getpar(j) << ",";
		}
		cout << "]" << endl;
	}
	cout << times << endl;
	cin >> i;
}

例子：
w₁: (-1,-1)^T , w₂: (0,0)^T , w₃: (1,1)^T
运行结果：