fandoghi_fill=fandoghi_fill(1:2500,1:51);
fandoghi=fandoghi_fill;
ahmad=ahmad_fill;
data_train=[fandoghi(1:158,1:40),ahmad(1:158,1:40)];
data_test=[fandoghi(1:158,41:51),ahmad(1:158,41:51)];
% maximum=max(data_train);
% minimum=min(data_train);
data_train_1=data_train(1,:);
%%%%%%%%%%zero_crossing
% information=data_train;
% [count_zero_crossing,count_zero]=zero_crossing(information);
y=fft(data_train);
mag=abs(y);
magn=mag(1:80,:);
[magni,fstdm]=mapstd(magn);
magni=removeconstantrows(magni);
[magnitude,fpcam]=processpca(magni,0.01);
ph=angle(y);
pha=ph(1:80,:);
[phas,fstdph]=mapstd(pha);
[phase,fpcaph]=processpca(phas,0.01);
inputnetwork=[data_train_1;magnitude;phase];
range=minmax(inputnetwork);
net=newff(range,[30 15 2]);
net.trainParam.epochs =100;
net=train(net,inputnetwork,targettrain);
output=sim(net,inputnetwork);
output_digital=hardlim(output-0.9999);
data_test_1=data_test(1,:);
x=fft(data_test);
amp=abs(x);
ampl=amp(1:80,:);
ampli=mapstd('apply',ampl,fstdm);
amplitude=processpca('apply',ampli,fpcam);
dar=angle(x);
dara=dar(1:80,:);
daraj=mapstd('apply',dara,fstdph);
daraje=processpca('apply',daraj,fpcaph);
testnetwork=[data_test_1;amplitude;daraje];
outputtest=sim(net,testnetwork);
outputtest_digital=hardlim(outputtest-0.9999);