%定義輸入信號及目標信號
Time=1:80;
p1=sin(1:20);
p2=sin(1:20)*2;
t1=ones(1,20);
t2=ones(1,20)*2;
p=[p1 p2 p1 p2];
t=[t1 t2 t1 t2];
Pseq=con2seq(p);
Tseq=con2seq(t);
%繪出輸入信號及目標信號曲線
figure;
plot(Time,cat(2,Pseq{:}),'--',Time,cat(2,Tseq{:}));
xlabel('時間');
ylabel('輸入及目標');
title('輸入信號及目標信號');
%開始生成網絡
R=1; % 一個輸入元件
S2=1; % 第二層一個輸出神經元
S1=10;% 在第一層中包含10個遞歸神經元單元
net=newelm([-2 2],[S1 S2],{'tansig','purelin'},'traingdx');
%訓練網絡
net.trainParam.epochs = 500; %設置訓練時間為500個時間單位
[net,tr]=train(net,Pseq,Tseq);
%開始測試網絡性能
a = sim(net,Pseq);
%繪出輸出信號及目標信號曲線
figure;
plot(Time,cat(2,a{:}),Time,cat(2,Tseq{:}),'--');
xlabel('時間');
ylabel('輸出及目標');
title('輸出信號及目標信號');