%定義輸入信號并繪出其曲線
time=0:0.025:5;
X=sin(sin(time).*time*10);
plot(time,X);
title('輸入信號T');
xlabel('時間');
ylabel('輸入信號');
figure;
%定義系統(tǒng)線性變換函數(shù),繪出系統(tǒng)輸出曲線
T=X*2+0.8;
plot(time,T);
title('系統(tǒng)輸出T');
xlabel('時間');
ylabel('系統(tǒng)輸出');
%定義網(wǎng)絡(luò)輸入
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));
%建立網(wǎng)絡(luò)
net=newlind(P,T);
%測試網(wǎng)絡(luò)
a=sim(net,P);
%繪出網(wǎng)絡(luò)輸出a與系統(tǒng)輸出T
figure;
plot(time,a,'+',time,T,'--');
title('網(wǎng)絡(luò)輸出a與系統(tǒng)輸出T');
xlabel('時間');
ylabel('系統(tǒng)輸出--  網(wǎng)絡(luò)輸出+');
%計算誤差,并繪出其曲線
e=T-a;
figure;
plot(time,e);
title('輸出誤差');
xlabel('時間');
ylabel('誤差');