#include "stdio.h"
#include "math.h"
#define MAX 201
#define NMLT 5
#define PI 3.1415926
#define alph 0.1
int method=2;//Œë�·•]‰¿‚ÌŽè–@
//double inpRes[NMLT]={2,-0.2, 0.5,-0.1,0.05};
//double inpRes[NMLT]={2,-0.2, 0.5,0,0};
double inpRes[NMLT]={1.5,0.3, 0.1,0.05,0.00};//–¢’mƒVƒXƒeƒ€�M�†�¶�¬—pƒCƒ“ƒpƒ‹ƒX‰ž“š
double Xsig[MAX], Ysig[MAX], Out[MAX];//Xsig:“ü—Í, Ysig:�o—Í, Out:ƒtƒBƒ‹ƒ^�o—Í
double w[NMLT]={0,0,0,0,0};//‹tƒtƒBƒ‹ƒ^—p�æŽZ’l
double Convolution(double w[], int i){//�ô‚Ý�ž‚݉‰ŽZ
	double V=0;
	for(int ii=i, k=0;k<NMLT;k++, ii--)	if(ii>=0) V+=w[k]*Xsig[ii];
	return V;
}
double MultX(int i){// <X|X>(X‚ÆX‚Ì“]’u�s—ñ‚Æ‚Ì“à�Ï)
	double V=0;
    for(int ii=i, k=0;k<NMLT;k++, ii--)if(ii>=0) V+=Xsig[ii]*Xsig[ii];
	return V;
}
int main(void){
	double err;	printf("No,X,Y,Out,w0,w1,w2,w3,w4,err");
	for(int i=0;i<MAX;i++){
		//Xsig[i]=sin(i*PI/20);
		Xsig[i]=sin(i*PI/20)+cos(i*PI/40)-cos(i*3*PI/30)/3;//“ü—͂̃Vƒ~ƒ…ƒŒ�[ƒVƒ‡ƒ“
		Ysig[i]=Convolution(inpRes,i);//–¢’mƒVƒXƒeƒ€‚̃Vƒ~ƒ…ƒŒ�[ƒVƒ‡ƒ“
		Out[i]=Convolution(w,i);      //�o—Í‚ÌŒvŽZ
		double XX=MultX(i),ep2=0;     //Œë�·ŒvŽZ‚ÆWi‚Ì�ÄŒvŽZ
		if(method==2){err=Ysig[i]-Out[i]; ep2+=err*err;}//Œë�·•]‰¿�@Žè–@‚Q
		if(XX>1E-12){
			for(int ii=i,k=0;k<NMLT;k++,ii--){
				if(method==1) {err=Ysig[ii]-Out[ii];ep2+=err*err;} //Œë�·•]‰¿Žè–@‚P
				w[k]+=alph*err*Xsig[ii]/XX;
			}
		}
		Out[i]=Convolution(w,i);//�o—Í�ÄŒvŽZ
		if(method==1) ep2/=NMLT;//Œë�·’l‚Ì•]‰¿
		ep2=sqrt(ep2);
		printf("\n %d, %7.4lf, %7.4lf, %7.4lf",i,Xsig[i],Ysig[i],Out[i]);
		for(int k=0;k<NMLT;k++) printf(",%7.4lf",w[k]);
		printf(",%7.4lf",ep2);
	}
	getchar();
	return 1;
}