%定義輸入信號并繪出其曲線
time=0:0.025:5;
X=sin(sin(time).*time*10);
plot(time,X);
title('輸入信號T');
xlabel('時間');
ylabel('輸入信號');
figure;
%定義系統線性變換函數,繪出系統輸出曲線
T=X*2+0.8;
plot(time,T);
title('系統輸出T');
xlabel('時間');
ylabel('系統輸出');
%定義網絡輸入
Q=size(X,2);
P=zeros(3,Q);
P(1,1:Q)=X(1,1:Q);
P(2,2:Q)=X(1,1:(Q-1));
P(3,3:Q)=X(1,1:(Q-2));
%建立網絡
net=newlind(P,T);
%測試網絡
a=sim(net,P);
%繪出網絡輸出a與系統輸出T
figure;
plot(time,a,'+',time,T,'--');
title('網絡輸出a與系統輸出T');
xlabel('時間');
ylabel('系統輸出--  網絡輸出+');
%計算誤差,并繪出其曲線
e=T-a;
figure;
plot(time,e);
title('輸出誤差');
xlabel('時間');
ylabel('誤差');