static void Main(string[] args)
{
System.Console.WriteLine("Inicio do Programa");
// parametros para o inicio da execução
int execucoes = 3;
int indDataSet = 3;
int populationSize = 100;
int iterations = 30000;
int selectionMethod = 2;
int functionsSet = 1;
int geneticMethod = 0;
//Thread workerThread = null;
bool needToStop = false;
String dataset = "";
WriteData saida = new WriteData();
//String nomearq = "Experimento" + ((geneticMethod == 0) ? "_GP_" : "_GEP_");
//saida.escreveArquivo("DataSet" + "\t" + "Metodo" + "\t" + "Execucao" + "\t" + "MMRE" + "\t" + "PRED(25)" + "\r\n", nomearq, 0);
//para executar os experimentos com a GP e GEP conforme parametros acima
for (int k = 0; k <= geneticMethod; k++)
{
String nomearq = "Experimento" + ((k == 0) ? "_GP_" : "_GEP_");
saida.escreveArquivo("DataSet" + "\t" + "Metodo" + "\t" + "Execucao" + "\t" + "MMRE" + "\t" + "PRED(25)" + "\r\n", nomearq, 0);
//para definir o dataset a ser utilizado conforme parametros acima
for (int j = 1; j <= indDataSet; j++)
{
switch (j)
{
case 1: dataset = "abran"; break;
case 2: dataset = "albrecht"; break;
case 3: dataset = "bailey"; break;
case 4: dataset = "belady"; break;
case 5: dataset = "boehm"; break;
case 6: dataset = "heiat"; break;
case 7: dataset = "jss"; break;
case 8: dataset = "kemerer"; break;
case 9: dataset = "miyaza"; break;
case 10: dataset = "shepp"; break;
case 11: dataset = "desharn"; break;
case 12: dataset = "kitchen"; break;
default: dataset = "teste"; break;
}
String caminho = @"D:\Projetos C\Instancias\";
String nomeArq = caminho + dataset + ".dat";
double[,] data = null;
// leitura e impressoa dos dados do dataset
data = (new LoadData()).leituraArquivo(nomeArq);
for (int i = 0; i < data.GetLength(0); i++)
{
System.Console.WriteLine(String.Format("{0:N}", data[i, 0]) + ", " + String.Format("{0:N}", data[i, 1]));
}
// quantidade de execuções a serem realizadas definidas nos parametros acima
for (int i = 1; i <= execucoes; i++)
{
// Criação do objeto solução
Solution solution = new Solution(data, populationSize, iterations, selectionMethod, functionsSet, k, needToStop, dataset, i);
// Criação do objeto do thread. Não inicia o thread.
//workerThread = new Thread(new ThreadStart(solution.SearchSolution));
// Inicia o thread.
//workerThread.Start();
System.Console.WriteLine("Iniciando a execução: " + i.ToString() + "\r\n");
solution.SearchSolution();
}
}
}
System.Console.WriteLine("Fim do Programa");
}
public void SearchSolution()
{
WriteData saida = new WriteData();
// create fitness function
MMREFitness fitness = new MMREFitness(data);
//SymbolicRegressionFitness fitness = new SymbolicRegressionFitness(data, new double[] { 1, 2, 3, 5, 7, 9});
// create gene function
IGPGene gene = (functionsSet == 0) ? (IGPGene)new SimpleGeneFunction(data.GetLength(0)) : (IGPGene)new ExpressionGeneFunction(data.GetLength(0));
// create population
Population population = new Population(populationSize, (geneticMethod == 0) ?
(IChromosome)new GPTreeChromosome(gene) :
(IChromosome)new GEPChromosome(gene, 15),
fitness,
(selectionMethod == 0) ? (ISelectionMethod)new EliteSelection() :
(selectionMethod == 1) ? (ISelectionMethod)new RankSelection() :
(ISelectionMethod)new RouletteWheelSelection()
);
// iterations ou numero de gerações a serem utilizadas
int i = 1;
// solution array
double[,] solution = new double[data.GetLength(0), 2];
double[] input = new double[data.GetLength(0)];
double[] output = new double[data.GetLength(0)];
//double[] input = new double[1];
// Alexander - Define o valor de entrada X com base no dataset
for (int j = 0; j < data.GetLength(0); j++)
{
input[j] = data[j, 1];
}
// Alexander - Define se saida sera baseada na função ou nos dados existentes no dataset
for (int j = 0; j < data.GetLength(0); j++)
{
if (useFunction)
output[j] = func(input[j]);
else
output[j] = data[j, 0];
System.Console.WriteLine(String.Format("{0:N}", input[j]) + ", " + String.Format("{0:N}", output[j]));
}
// loop
while (!needToStop)
{
System.Console.WriteLine("Geração: " + i.ToString());
hits = 0;
if (i <= 1) {
//Grava a população gerada e seu fitness em cada rodada
nomearq = "Populacao" + ((geneticMethod == 0) ? "_GP_" : "_GEP_") + dataset + "_" + executions.ToString() + "_" + "Geracao_" + i.ToString() + "_";
saida.escreveArquivo("Geracao: " + i.ToString() + "\r\n" + population.toString(), nomearq, 1);
nomearq = "AnalisedeTempo" + ((geneticMethod == 0) ? "_GP_" : "_GEP_") + dataset + "_" + executions.ToString() + "_" + iterations.ToString();
saida.escreveArquivo("DataSet" + "\t" + "Metodo" + "\t" + "Geracao" + "\t" + "MMRE" + "\r\n", nomearq, 0);
}
// run one epoch of genetic algorithm
population.ActualGeneration = i;
//population.RunEpoch();
//population.FindBestChromosome();
try
{
// get best solution
population.FindBestChromosome();
string bestFunction = population.BestChromosome.ToString().Trim();
System.Console.WriteLine("Função: " + RPN2Infix.PostfixToInfix(bestFunction));
npred = 0;
sum = 0.0;
error = 0.0;
pred = 0.00;
double result = 0.0;
double resultgerado = 0.0;
// calculate best function
for (int j = 0; j < data.GetLength(0); j++)
{
double HIT_LEVEL = 0.01;
double PROBABLY_ZERO = 1.11E-15;
double BIG_NUMBER = 1.0e15;
System.Console.WriteLine("Valor de entrada: " + input[j]);
resultgerado = PolishGerExpression.Evaluate(bestFunction, input, j);
// fitness (atual - estimado) / atual
result = Math.Abs((output[j] - resultgerado) / output[j]);
if (!(result < BIG_NUMBER)) // *NOT* (input.x >= BIG_NUMBER)
result = BIG_NUMBER;
else if (result < PROBABLY_ZERO) // slightly off
result = 0.0;
if (result <= HIT_LEVEL) hits++; // whatever!
//Somatorio do erro
sum += result;
//Somatoria do pred com 25%
if ((resultgerado <= output[j] * 1.25) && (resultgerado >= output[j] * 0.75)) npred++;
//impressao dos dados gerados e esperados
System.Console.WriteLine(" Valor gerado: " + resultgerado + " resultado esperado: " + data[j, 0] + " erro relativo: " + result);
}
// calculate error MMRE
error = sum / data.GetLength(0);
//calculate Pred(25%)
pred = (((double) npred) / data.GetLength(0));
if (error < bestMMRE)
{
bestMMRE = error;
nomearq = "AnalisedeTempo" + ((geneticMethod == 0) ? "_GP_" : "_GEP_") + dataset + "_" + executions.ToString() + "_" + iterations.ToString();
saida.escreveArquivo(dataset + "\t" + ((geneticMethod == 0) ? "GP" : "GEP") + "\t" + i.ToString() + "\t" + bestMMRE.ToString() + "\r\n", nomearq, 0);
}
System.Console.WriteLine("Erro acumulado: " + sum);
System.Console.WriteLine("Erro Médio: " + error);
System.Console.WriteLine("Hits: " + hits);
System.Console.WriteLine("");
}
catch (Exception e)
{
System.Console.WriteLine(e.Message);
}
// increase current iteration/geracao
i++;
//
if ((iterations != 0) && (i > iterations))
//break;
needToStop = true;
else
population.RunEpoch();
}
// show solution
//System.Console.WriteLine(population.BestChromosome.ToString());
string expressao = population.BestChromosome.ToString().Trim();
string expressaosubst = PolishGerExpression.SubstituteVariables(expressao, input);
string expressaosimpl = PolishGerExpression.SimplifyExpression(expressaosubst);
string resultado = "";
resultado = "Expressão NPR gerada: " + expressao + "\r\n";
resultado = resultado + "Expressão formatada 1: " + (RPN2Infix.PostfixToInfix(expressao)) + "\r\n";
resultado = resultado + "Expressão formatada 2: " + (RPN2Infix.Parse(expressao.Replace(",", "."))) + "\r\n";
//resultado = resultado + "Expressão NPR com substitução: " + expressaosubst + "\r\n";
//resultado = resultado + "Expressão com substitução formatada: " + RPN2Infix.PostfixToInfix(expressaosubst) + "\r\n";
//resultado = resultado + "Expressão NPR com substitução simplificada: " + expressaosimpl + "\r\n";
resultado = resultado + "Expressão com substitução simplificada formatada: " + RPN2Infix.PostfixToInfix(expressaosimpl) + "\r\n";
resultado = resultado + "Resultado para a expressão: " + "\r\n"; // +testeexpressao + "\r\n";
resultado = resultado + "Erro acumulado: " + sum.ToString() + "\r\n";
resultado = resultado + "Erro Médio: " + error.ToString() + "\r\n";
resultado = resultado + "Pred(25): " + pred.ToString() + "\r\n";
resultado = resultado + "Hits: " + hits.ToString() + "\r\n";
resultado = resultado + "Geração: " + population.BestGeneration + "\r\n";
nomearq = "Resultado" + ((geneticMethod == 0) ? "_GP_" : "_GEP_") + dataset + "_" + executions.ToString() + "_";
saida.escreveArquivo(resultado, nomearq, 1);
nomearq = "Experimento" + ((geneticMethod == 0) ? "_GP_" : "_GEP_");
saida.escreveArquivo(dataset + "\t" + ((geneticMethod == 0) ? "GP" : "GEP") + "\t" + executions.ToString() + "\t" + error.ToString() + "\t" + pred.ToString() + "\r\n", nomearq, 0);
//System.Console.WriteLine("Fim do Programa");
}
public void SearchSolution()
{
WriteData saida = new WriteData();
// create fitness function
MMREFitness fitness = new MMREFitness(data);
//SymbolicRegressionFitness fitness = new SymbolicRegressionFitness(data, new double[] { 1, 2, 3, 5, 7, 9});
// create gene function
IGPGene gene = (functionsSet == 0) ? (IGPGene) new SimpleGeneFunction(data.GetLength(0)) : (IGPGene) new ExpressionGeneFunction(data.GetLength(0));
// create population
Population population = new Population(populationSize, (geneticMethod == 0) ?
(IChromosome) new GPTreeChromosome(gene) :
(IChromosome) new GEPChromosome(gene, 15),
fitness,
(selectionMethod == 0) ? (ISelectionMethod) new EliteSelection() :
(selectionMethod == 1) ? (ISelectionMethod) new RankSelection() :
(ISelectionMethod) new RouletteWheelSelection()
);
// iterations ou numero de gerações a serem utilizadas
int i = 1;
// solution array
double[,] solution = new double[data.GetLength(0), 2];
double[] input = new double[data.GetLength(0)];
double[] output = new double[data.GetLength(0)];
//double[] input = new double[1];
// Alexander - Define o valor de entrada X com base no dataset
for (int j = 0; j < data.GetLength(0); j++)
{
input[j] = data[j, 1];
}
// Alexander - Define se saida sera baseada na função ou nos dados existentes no dataset
for (int j = 0; j < data.GetLength(0); j++)
{
if (useFunction)
{
output[j] = func(input[j]);
}
else
{
output[j] = data[j, 0];
}
System.Console.WriteLine(String.Format("{0:N}", input[j]) + ", " + String.Format("{0:N}", output[j]));
}
// loop
while (!needToStop)
{
System.Console.WriteLine("Geração: " + i.ToString());
hits = 0;
if (i <= 1)
{
//Grava a população gerada e seu fitness em cada rodada
nomearq = "Populacao" + ((geneticMethod == 0) ? "_GP_" : "_GEP_") + dataset + "_" + executions.ToString() + "_" + "Geracao_" + i.ToString() + "_";
saida.escreveArquivo("Geracao: " + i.ToString() + "\r\n" + population.toString(), nomearq, 1);
nomearq = "AnalisedeTempo" + ((geneticMethod == 0) ? "_GP_" : "_GEP_") + dataset + "_" + executions.ToString() + "_" + iterations.ToString();
saida.escreveArquivo("DataSet" + "\t" + "Metodo" + "\t" + "Geracao" + "\t" + "MMRE" + "\r\n", nomearq, 0);
}
// run one epoch of genetic algorithm
population.ActualGeneration = i;
//population.RunEpoch();
//population.FindBestChromosome();
try
{
// get best solution
population.FindBestChromosome();
string bestFunction = population.BestChromosome.ToString().Trim();
System.Console.WriteLine("Função: " + RPN2Infix.PostfixToInfix(bestFunction));
npred = 0;
sum = 0.0;
error = 0.0;
pred = 0.00;
double result = 0.0;
double resultgerado = 0.0;
// calculate best function
for (int j = 0; j < data.GetLength(0); j++)
{
double HIT_LEVEL = 0.01;
double PROBABLY_ZERO = 1.11E-15;
double BIG_NUMBER = 1.0e15;
System.Console.WriteLine("Valor de entrada: " + input[j]);
resultgerado = PolishGerExpression.Evaluate(bestFunction, input, j);
// fitness (atual - estimado) / atual
result = Math.Abs((output[j] - resultgerado) / output[j]);
if (!(result < BIG_NUMBER)) // *NOT* (input.x >= BIG_NUMBER)
{
result = BIG_NUMBER;
}
else if (result < PROBABLY_ZERO) // slightly off
{
result = 0.0;
}
if (result <= HIT_LEVEL)
{
hits++; // whatever!
}
//Somatorio do erro
sum += result;
//Somatoria do pred com 25%
if ((resultgerado <= output[j] * 1.25) && (resultgerado >= output[j] * 0.75))
{
npred++;
}
//impressao dos dados gerados e esperados
System.Console.WriteLine(" Valor gerado: " + resultgerado + " resultado esperado: " + data[j, 0] + " erro relativo: " + result);
}
// calculate error MMRE
error = sum / data.GetLength(0);
//calculate Pred(25%)
pred = (((double)npred) / data.GetLength(0));
if (error < bestMMRE)
{
bestMMRE = error;
nomearq = "AnalisedeTempo" + ((geneticMethod == 0) ? "_GP_" : "_GEP_") + dataset + "_" + executions.ToString() + "_" + iterations.ToString();
saida.escreveArquivo(dataset + "\t" + ((geneticMethod == 0) ? "GP" : "GEP") + "\t" + i.ToString() + "\t" + bestMMRE.ToString() + "\r\n", nomearq, 0);
}
System.Console.WriteLine("Erro acumulado: " + sum);
System.Console.WriteLine("Erro Médio: " + error);
System.Console.WriteLine("Hits: " + hits);
System.Console.WriteLine("");
}
catch (Exception e)
{
System.Console.WriteLine(e.Message);
}
// increase current iteration/geracao
i++;
//
if ((iterations != 0) && (i > iterations))
{
//break;
needToStop = true;
}
else
{
population.RunEpoch();
}
}
// show solution
//System.Console.WriteLine(population.BestChromosome.ToString());
string expressao = population.BestChromosome.ToString().Trim();
string expressaosubst = PolishGerExpression.SubstituteVariables(expressao, input);
string expressaosimpl = PolishGerExpression.SimplifyExpression(expressaosubst);
string resultado = "";
resultado = "Expressão NPR gerada: " + expressao + "\r\n";
resultado = resultado + "Expressão formatada 1: " + (RPN2Infix.PostfixToInfix(expressao)) + "\r\n";
resultado = resultado + "Expressão formatada 2: " + (RPN2Infix.Parse(expressao.Replace(",", "."))) + "\r\n";
//resultado = resultado + "Expressão NPR com substitução: " + expressaosubst + "\r\n";
//resultado = resultado + "Expressão com substitução formatada: " + RPN2Infix.PostfixToInfix(expressaosubst) + "\r\n";
//resultado = resultado + "Expressão NPR com substitução simplificada: " + expressaosimpl + "\r\n";
resultado = resultado + "Expressão com substitução simplificada formatada: " + RPN2Infix.PostfixToInfix(expressaosimpl) + "\r\n";
resultado = resultado + "Resultado para a expressão: " + "\r\n"; // +testeexpressao + "\r\n";
resultado = resultado + "Erro acumulado: " + sum.ToString() + "\r\n";
resultado = resultado + "Erro Médio: " + error.ToString() + "\r\n";
resultado = resultado + "Pred(25): " + pred.ToString() + "\r\n";
resultado = resultado + "Hits: " + hits.ToString() + "\r\n";
resultado = resultado + "Geração: " + population.BestGeneration + "\r\n";
nomearq = "Resultado" + ((geneticMethod == 0) ? "_GP_" : "_GEP_") + dataset + "_" + executions.ToString() + "_";
saida.escreveArquivo(resultado, nomearq, 1);
nomearq = "Experimento" + ((geneticMethod == 0) ? "_GP_" : "_GEP_");
saida.escreveArquivo(dataset + "\t" + ((geneticMethod == 0) ? "GP" : "GEP") + "\t" + executions.ToString() + "\t" + error.ToString() + "\t" + pred.ToString() + "\r\n", nomearq, 0);
//System.Console.WriteLine("Fim do Programa");
}
static void Main(string[] args)
{
System.Console.WriteLine("Inicio do Programa");
// parametros para o inicio da execução
int execucoes = 3;
int indDataSet = 3;
int populationSize = 100;
int iterations = 30000;
int selectionMethod = 2;
int functionsSet = 1;
int geneticMethod = 0;
//Thread workerThread = null;
bool needToStop = false;
String dataset = "";
WriteData saida = new WriteData();
//String nomearq = "Experimento" + ((geneticMethod == 0) ? "_GP_" : "_GEP_");
//saida.escreveArquivo("DataSet" + "\t" + "Metodo" + "\t" + "Execucao" + "\t" + "MMRE" + "\t" + "PRED(25)" + "\r\n", nomearq, 0);
//para executar os experimentos com a GP e GEP conforme parametros acima
for (int k = 0; k <= geneticMethod; k++)
{
String nomearq = "Experimento" + ((k == 0) ? "_GP_" : "_GEP_");
saida.escreveArquivo("DataSet" + "\t" + "Metodo" + "\t" + "Execucao" + "\t" + "MMRE" + "\t" + "PRED(25)" + "\r\n", nomearq, 0);
//para definir o dataset a ser utilizado conforme parametros acima
for (int j = 1; j <= indDataSet; j++)
{
switch (j)
{
case 1: dataset = "abran"; break;
case 2: dataset = "albrecht"; break;
case 3: dataset = "bailey"; break;
case 4: dataset = "belady"; break;
case 5: dataset = "boehm"; break;
case 6: dataset = "heiat"; break;
case 7: dataset = "jss"; break;
case 8: dataset = "kemerer"; break;
case 9: dataset = "miyaza"; break;
case 10: dataset = "shepp"; break;
case 11: dataset = "desharn"; break;
case 12: dataset = "kitchen"; break;
default: dataset = "teste"; break;
}
String caminho = @"D:\Projetos C\Instancias\";
String nomeArq = caminho + dataset + ".dat";
double[,] data = null;
// leitura e impressoa dos dados do dataset
data = (new LoadData()).leituraArquivo(nomeArq);
for (int i = 0; i < data.GetLength(0); i++)
{
System.Console.WriteLine(String.Format("{0:N}", data[i, 0]) + ", " + String.Format("{0:N}", data[i, 1]));
}
// quantidade de execuções a serem realizadas definidas nos parametros acima
for (int i = 1; i <= execucoes; i++)
{
// Criação do objeto solução
Solution solution = new Solution(data, populationSize, iterations, selectionMethod, functionsSet, k, needToStop, dataset, i);
// Criação do objeto do thread. Não inicia o thread.
//workerThread = new Thread(new ThreadStart(solution.SearchSolution));
// Inicia o thread.
//workerThread.Start();
System.Console.WriteLine("Iniciando a execução: " + i.ToString() + "\r\n");
solution.SearchSolution();
}
}
}
System.Console.WriteLine("Fim do Programa");
}