%定義輸入信號(hào)及目標(biāo)信號(hào)
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);
%繪出輸入信號(hào)及目標(biāo)信號(hào)曲線
figure;
plot(Time,cat(2,Pseq{:}),'--',Time,cat(2,Tseq{:}));
xlabel('時(shí)間');
ylabel('輸入及目標(biāo)');
title('輸入信號(hào)及目標(biāo)信號(hào)');
%開始生成網(wǎng)絡(luò)
R=1; % 一個(gè)輸入元件
S2=1; % 第二層一個(gè)輸出神經(jīng)元
S1=10;% 在第一層中包含10個(gè)遞歸神經(jīng)元單元
net=newelm([-2 2],[S1 S2],{'tansig','purelin'},'traingdx');
%訓(xùn)練網(wǎng)絡(luò)
net.trainParam.epochs = 500; %設(shè)置訓(xùn)練時(shí)間為500個(gè)時(shí)間單位
[net,tr]=train(net,Pseq,Tseq);
%開始測試網(wǎng)絡(luò)性能
a = sim(net,Pseq);
%繪出輸出信號(hào)及目標(biāo)信號(hào)曲線
figure;
plot(Time,cat(2,a{:}),Time,cat(2,Tseq{:}),'--');
xlabel('時(shí)間');
ylabel('輸出及目標(biāo)');
title('輸出信號(hào)及目標(biāo)信號(hào)');