%定義輸入信號并繪出其曲線
time1=0:0.005:4;
time2=4.005:0.005:6;
time=[time1 time2];
X=sin(sin(time*4).*time*8);
plot(time,X);
title('輸入信號X');
xlabel('時間');
ylabel('輸入信號');
%定義系統輸出,繪出曲線
steps1=length(time1);
[T1,state]=filter([1 -0.5],1,X(1:steps1));
steps2=length(time2);
T2=filter([0.9 -0.6],1,X((1:steps2) + steps1),state);
T=[T1 T2];
figure;
plot(time,T);
title('系統輸出T');
xlabel('時間');
ylabel('系統輸出');
%將系統輸入和輸出轉換成序列信號
T=con2seq(T);
P=con2seq(X);
%建立網絡
lr=0.5;
delays=[0 1];
net=newlin(minmax(cat(2,P{:})),1,delays,lr);
%訓練網絡
[net,a,e]=adapt(net,P,T);
%繪出網絡輸出a與系統輸出T
figure;
plot(time,cat(2,a{:}),'+',time,cat(2,T{:}),'--');
title('網絡輸出a與系統輸出T');
xlabel('時間');
ylabel('系統輸出--  網絡輸出+');
%繪出誤差曲線
figure;
plot(time,cat(2,e{:}));
title('輸出誤差');
xlabel('時間');
ylabel('誤差');