%定義輸入信號(hào)并繪出其曲線
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('輸入信號(hào)X');
xlabel('時(shí)間');
ylabel('輸入信號(hào)');
%定義系統(tǒng)輸出,繪出曲線
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ǒng)輸出T');
xlabel('時(shí)間');
ylabel('系統(tǒng)輸出');
%將系統(tǒng)輸入和輸出轉(zhuǎn)換成序列信號(hào)
T=con2seq(T);
P=con2seq(X);
%建立網(wǎng)絡(luò)
lr=0.5;
delays=[0 1];
net=newlin(minmax(cat(2,P{:})),1,delays,lr);
%訓(xùn)練網(wǎng)絡(luò)
[net,a,e]=adapt(net,P,T);
%繪出網(wǎng)絡(luò)輸出a與系統(tǒng)輸出T
figure;
plot(time,cat(2,a{:}),'+',time,cat(2,T{:}),'--');
title('網(wǎng)絡(luò)輸出a與系統(tǒng)輸出T');
xlabel('時(shí)間');
ylabel('系統(tǒng)輸出--  網(wǎng)絡(luò)輸出+');
%繪出誤差曲線
figure;
plot(time,cat(2,e{:}));
title('輸出誤差');
xlabel('時(shí)間');
ylabel('誤差');