#include <fstream.h>
#include <iostream.h>
#include <time.h>
#include <stdlib.h>
#include <math.h>
#define   HN    30
#define  N  96
#define   daynumber  20
#define   classifynumber   10
double pricepointer[daynumber+1][N];
double loadpointer[daynumber][N];
double goalprice[N];
double nextdayload[N];
//隱含層神經元個數, 請多次修改其值來進行比較
#define   IN  5
struct sample
{
	double input[IN];
	double teach;
	
};

int unitnumber;							
double maxprice=500;
double maxload=1e5;
double xout[IN-1];                 //輸入層的輸出值
double hn_out[HN];                //隱含層的輸出值
double y;                         //輸出層的輸出值
double w[HN][IN-1];
double w1[HN][IN-1];
double w0[HN][IN-1]; 
double hn_b[HN];                  //隱含層的偏置系數
double v[HN];                     //隱含層到輸出層的權值
double v1[HN];
double v0[HN];
double y_b;                       //輸出層的偏置系數
double hn_delta[HN];              //隱含層的delta
double y_delta;                   //輸出層的delta
double error;                     //全局誤差
double error_1=100.0;
double errlimit=1.0e-10;           //預設的全局誤差極限
double alpha=0.6, beta=0.2;       //預設的學習速率
double a;	
long times=5000;                    //預設的訓練次數上限
int i, j,m;
int loop;

double sumtemp;
double errtemp;
struct sample data[daynumber+1];


/*************************************************************
*					數據的讀取							     *															 *
**************************************************************/

void getteachingdataprice()   
{
	ifstream Fread("電價樣本數據.txt",ios::in );
	//for(int m=0; m<daynumber; m++)
	for(int j=0; j<N; j++)
	{
		//for(int j=0; j<N; j++)
		for(int m=0; m<daynumber; m++)
		{
			Fread>>pricepointer[m][j];
			pricepointer[m][j]=pricepointer[m][j]/maxprice;
		}
	}
	
	
	Fread.close();
}



void getteachingdataload()
{
	ifstream Fread("負荷樣本數據.txt",ios::in);
	//for(int m=0; m<daynumber; m++)
	for(int j=0; j<N; j++)
	{
		//for(int j=0; j<N; j++)
		for(int m=0; m<daynumber; m++)
		{
			Fread>>loadpointer[m][j];
			loadpointer[m][j]=loadpointer[m][j]/maxload;
		}
	}

	Fread.close();
}


void gettestingdataprice() //讀入目標電價  
{
	ifstream Fread("下一日實際電價數據.txt",ios::in );
	for(int j=0; j<N; j++)
	{
		Fread>>goalprice[j];
		goalprice[j]=goalprice[j]/maxprice;
	}	
	Fread.close();
}


void gettestingload()
{
	ifstream Fread("下一日負荷預測數據.txt",ios::in );
	for(int j=0; j<N; j++)
	{
		Fread>>nextdayload[j];
		nextdayload[j]=nextdayload[j]/maxload;
	}
	Fread.close();
}