que1hivshengjing.m

signal and system answer

clear
t=[5.095820783341429,4.452483209581492,4.189913281945378,3.876701375187596,3.431933896729329,3.477797999814873,3.381813771892126,3.200390410226402,3.064234861011267,2.952680941706319,2.863186994331568,2.755253288942034,2.662849431311355,2.578974239966391,2.510102974111348,2.454728765291937,2.414402683154488,2.386519203783996,2.379641636977733,2.379072502648362,2.392709208639219,2.417585292565604,2.453317547941533,2.470899075643746,2.551853804278938,2.508943938900803,2.386126781280675,2.401407098847278,2.401921349556324,2.386963801780501,2.416437391239251,2.453469509869026,2.458700421496229,2.4879552391846,2.543260621171543,2.605332646386495,2.674171314829457,2.736242320747744,2.793008756258855,2.894207179694991,3.166003561365927,3.342090379785534,3.40695839646896];
p=1:43;
%figure(1)
%plot(p,t,'r');
%hold on
n=20;
net=newff([0,100],[n 1],{'tansig' 'purelin'},'trains');
net.trainParam.epochs=1600;
net.trainParam.lr=0.05;
net.trainParam.goal=0.1;
[net,tr]=train(net,p,t); 
y=sim(net,p);
figure(1)
plot(p,y);
hold on
figure
pk=41:1:60;
y1=sim(net,pk);
plot(p,t,'b')
hold on
plot(pk,y1,'g')