仿真程序.m

这是在论文中

clf;
figure(1);
clc
Bininput=[7207 7405 7805 8100 8400 8801 9001 9229 9401 9601 9801 10001 10200 10400 10879 11000 11400 11701 12100 12506]/15000;
Toutput=[-1.3 2.4 8.8 12.7 16.3 20.5 22.3 24.4 25.8 27.4 29.0 30.5 31.9 33.2 36.3 37.0 39.3 41.0 43.1 45.1]/50;
plot(Bininput,Toutput,'b-');
pause
s=2;
[w1,b1,w2,b2]=initff(Bininput,s,'tansig',Toutput,'purelin');
echo off
df=1;
me=1200;
eg=0.0001;
lr=0.05;
tp=[df me eg lr];
figure(2);
[w1,b1,w2,b2,ep,tr]=trainlm(w1,b1,'tansig',w2,b2,'purelin',Bininput,Toutput,tp);
pause
w1
b1
w2
b2
pause
figure(3);
ploterr(tr,eg);
pause
figure(4);
a=simuff(Bininput,w1,b1,'tansig',w2,b2,'purelin');
plot(Bininput*15000,a*50,'r.',Bininput*15000,Toutput*50,'bo')
Bininput=[7238 7405 7501 7800 8092 8401 8701 9101 9501 9901 10101 10300 10501 10800 11100 11501 11801 12100 12404 12506]/15000;
Toutput=[-0.7 2.4 4.1 8.7 12.6 16.3 19.5 23.3 26.7 29.7 31.2 32.6 33.9 35.8 37.6 39.9 41.5 43.1 44.6 45.1]/50;
figure(5);
a=simuff(Bininput,w1,b1,'tansig',w2,b2,'purelin');
plot(Bininput*15000,a*50,'r.',Bininput*15000,Toutput*50,'bo')
echo off