centerneuron.cpp

基于VC开发的神经网络工具箱

#include "../include/Exception.h"
#include "../include/CenterNeuron.h"

#include <cmath>
#include <cstdio>
using namespace std;
namespace annie
{

CenterNeuron::CenterNeuron(int label, int dimension) : Neuron(label)
{
	_dimension = dimension;
	_activationFunction = gaussian;
	_dActivationFunction = dgaussian;

	_center = new real[dimension];
//	_deltaCenter = new real[dimension];
	for (int i=0;i<dimension;i++)
	{
		_center[i] = random();
//		_deltaCenter[i]=(real)0.0;
	}

	_classHeirarchy.push_back(_CENTER_NEURON_STRING);
}

CenterNeuron::CenterNeuron(int label, VECTOR center) : Neuron(label)
{
	_dimension = center.size();
	setCenter(center);
	_activationFunction = gaussian;
	_dActivationFunction = dgaussian;
}

CenterNeuron::CenterNeuron(int label, int dimension, real center[]) : Neuron(label)
{
	_dimension = dimension;
	_activationFunction = gaussian;
	_dActivationFunction = dgaussian;
	setCenter(center);
}

CenterNeuron::~CenterNeuron()
{	
//	delete[] _deltaCenter;
	delete[] _center;	
}

void
CenterNeuron::setActivationFunction(ActivationFunction f, ActivationFunction df)
{
	_activationFunction = f;
	_dActivationFunction = df;
}

void
CenterNeuron::setCenter(real center[])
{
	int dim = getDimension();
	for (int i=0;i<dim;i++)
		_center[i]=center[i];
}

void
CenterNeuron::setCenter(VECTOR center)
{
	if ((int)center.size() != getDimension())
		throw Exception("CenterNeuron::setCenter() - Center provided has different dimension than that of this neuron.");
	VECTOR::iterator it;
	int i;
	for (i=0, it=center.begin(); it!=center.end(); it++,i++)
		_center[i]=(*it);
}

VECTOR
CenterNeuron::getCenter()
{	
	VECTOR center;
	for (int i=0;i<getDimension();i++)
		center.push_back(_center[i]);
	return center;	
}

void
CenterNeuron::_recacheOutput()
{
	if (_outputCacheValid)
		return;
	if (getInputCount() != getDimension())
		throw Exception("CenterNeuron::_recacheOutput() - Attempt to calculate output even though neuron doesn't have as many inputs as it the center's dimension.");
	_activationCache = 0.0;
	LINKS::iterator i;
	int j;
	for (i = _inputLinks.begin(), j = 0; i!=_inputLinks.end(); i++, j++)
	{
		Link *l = (Link *)(*i);
		real temp = l->getSource()->getOutput() - _center[j];
		_activationCache += (temp*temp);
	}
	_activationCache = (real)sqrt(_activationCache);
	_outputCache = _activationFunction(_activationCache);
	
	_outputCacheValid = true;
}

void
CenterNeuron::_recacheError()
{	_errorCacheValid = true;	}

int
CenterNeuron::getDimension()
{	return _dimension;	}

void
CenterNeuron::connect(Neuron *from)
{
	disconnect(from);
	Link *link = new Link(this,from,1.0);
	_inputLinks.push_back(link);
	from->_outputLinks.push_back(link);
	invalidateOutputCache();
}

//void
//CenterNeuron::calculateNewCenter(real learningRate)
//{
//	LINKS::iterator k1,i;
//	int k;
//	real act,dact,factor;
//	Link *inl,*outl;
//	if (getInputCount() != getDimension())
//	{
//		string error(getClassName());
//		error = error + "::calculateNewCenter() - Setup not valid, number of inputs not equal to dimension of center.";
//		throw Exception(error);
//	}
//	
//	for (k1=_inputLinks.begin(), k=0; k1!=_inputLinks.end(); k1++,k++)
//	{
//		inl=(Link*)(*k1);
//		_deltaCenter[k] = (real)0.0;
//		act = getActivation();
//		dact = (real)(exp(-1*act*act)*2);//_dActivationFunction(act); //
//		factor = (inl->getSource()->getOutput() - _center[k])*dact;///act;
//		for (i=_outputLinks.begin(); i!=_outputLinks.end(); i++)
//		{
//			outl=(Link*)(*i);
//			_deltaCenter[k] += outl->getDestination()->getError() * outl->getWeight();
//		}
//		_deltaCenter[k] = _deltaCenter[k] * factor * learningRate;
//	}
//}
//
//void
//CenterNeuron::update()
//{
//	int dim = getDimension();
//	int i;
//	for (i=0; i<dim; i++)
//	{
//		_center[i]+=_deltaCenter[i];
//		_deltaCenter[i]=(real)0.0;
//	}
//	invalidateOutputCache();
//	invalidateErrorCache();
//}

const char *
CenterNeuron::getClassName()
{
	return _CENTER_NEURON_STRING;
}

string
CenterNeuron::toString()
{
	string ans = Neuron::toString();
	char buffer[1000];
	int i=0;
	sprintf(buffer,"- Center = (");

	for (i=0;i<getDimension()-1;i++)
		sprintf(buffer,"%s%f,",buffer,_center[i]);
	sprintf(buffer,"%s%f)\n",buffer,_center[i]);
	ans = ans + buffer;
	return ans;
}

}; //namespace annie