%[w,b]=solvelin(p,t)
%Linear system identification

clc
pause
clc
time=0:0.025:5;
X=sin(sin(time).*time*10);

P=delaysig(X,0,2);
T=filter([1 0.5 -1.5],1,X);
pause
clc
figure
plot(time,X)
alabel('Time','Input Signal','Input Signal to the System')
pause
clc
figure
plot(time,T)
alabel('Time','Output Signal','Output Signal to the System')
pause
clc

[w,b]=solvelin(P,T)
pause
clc
a=simulin(P,w,b);
figure
plot(time,a,time,T,'+')
alabel('Time','Network output_System Output +','Network and System Output Signals')
pause
clc

e=T-a;
figure
plot(time,e,[min(time) max(time)],[0 0],':r')
alabel('Time','Error','Error Signal')
echo off
disp('End of NN_ID')