%%
%對(duì)程序的說(shuō)明%
%做生成數(shù)（使用0.1做因子）之后13年的數(shù)據(jù)按統(tǒng)計(jì)方法標(biāo)準(zhǔn)化后輸入訓(xùn)練仿真對(duì)比結(jié)果。這個(gè)過(guò)程相當(dāng)于插值計(jì)算%
%輸出結(jié)果請(qǐng)最后五句r1rlplotyy（ktokt）plotyy（kakf）%

%這樣做的缺點(diǎn)是不能檢驗(yàn)預(yù)測(cè)效果的好壞還有一點(diǎn)我對(duì)網(wǎng)絡(luò)初始化的過(guò)程沒(méi)有把握%

clear
%?數(shù)據(jù)輸入
inv=[189.6??244.2??328.5??405.3??487.5??568.6??669.9??801.6??785.9??859.2??979.7??1086.3??1241];
d1=[15.3??16.6??17.2??22.8??37.6??30.8??38.9??50.8??58.2??41.2??55.8??71.7??82.3];
d2=[117.7??140.7??177.8??209.8??247.7??290.8??343.2??337.1??308.1??353.6??350.1??367.2??415.7];
d3=[56.6??86.9??133.5??172.6??202.2??247??287.8??413.7??419.6??464.4??573.8??647.4??742.1];
f=[822.3??959.7??1198.4??1607.5??1997.6??2380.9??2683.8??2798.9??2897.4??3253??3561??3882?4433];
k=1991:1:2003;
%?使用13年數(shù)據(jù)訓(xùn)練得到網(wǎng)絡(luò)仿真計(jì)算相對(duì)誤差
%對(duì)原始數(shù)據(jù)初始化的過(guò)程為t（1）=f（1）;??t（i）=t（i-1）+0.1*f（i）（i=2:13）
%?累加?p1（1）=d1（1）;for?i=2:13?p1（i）=0.1*d1（i）+p1（i-1）;end?%
%?累減?a（1）=v（1）;for?i=2:13?a（i）=（v（i）-v（i-1））*10;end?%
%?先賦予0值
p1=zeros（113）;p2=zeros（113）;p3=zeros（113）;t=zeros（113）;
p1（1）=d1（1）;for?i=2:13?p1（i）=0.*d1（i）+p1（i-1）;end
p2（1）=d2（1）;for?i=2:13?p2（i）=0.1*d2（i）+p2（i-1）;end
p3（1）=d3（1）;for?i=2:13?p3（i）=0.1*d3（i）+p3（i-1）;end
t（1）=f（1）;for?i=2:13?t（i）=0.1*f（i）+t（i-1）;end
b=[p1;p2;p3;t]‘;
p=[（p1-mean（p1））./std（p1）;（p2-mean（p2））./std（p2）;（p3-mean（p3））./std（p3）];
tt=（t-?mean（t））./std（t）;%數(shù)據(jù)標(biāo)準(zhǔn)化處理標(biāo)準(zhǔn)化為網(wǎng)絡(luò)輸入p期望輸出tt

%建立網(wǎng)絡(luò)并訓(xùn)練
net=newff（minmax（p）[41]{‘tansig‘‘purelin‘}‘trainlm‘）;
net.iw{11}=zeros（size（net.iw{11}））+0.5;
net.lw{21}=zeros（size（net.lw{21}））+0.75;
net.b{11}=zeros（size（net.b{11}））+0.5;
net.b{21}=zeros（size（net.b{21}））;
net.trainParam.epochs=3000;
net.trainParam.goal?=0.000005;
net=train（netptt）;
%網(wǎng)絡(luò)仿真得到網(wǎng)絡(luò)輸出，并計(jì)算誤差
tt1=sim（netp）;
to=tt1.*std（t）+mean（t）;%利用標(biāo)準(zhǔn)化的逆變換得到t1的近似值to
a=zeros（113）;a（1）=to（1）;
for?i=2:13?a（i）=（to（i）-to（i-1））*10;?end?%累減得到近似的實(shí)際產(chǎn)值a
r_net=（tt1-tt）./tt;%網(wǎng)絡(luò)相對(duì)誤差
r_add=（to-t）./t;%總產(chǎn)值累加數(shù)相對(duì)誤差
r_real=（a-f）./f;%實(shí)際總產(chǎn)值相對(duì)誤差

%輸出結(jié)果%
r1=[r_net;r_add;r_real]‘%網(wǎng)絡(luò)相對(duì)誤差/累加數(shù)相對(duì)誤差/實(shí)際值相對(duì)誤差
r=r_real./r_add?%誤差放大倍數(shù)
l=[to;t;to-t;a;f;a-f]‘%累加數(shù)計(jì)算值/累加數(shù)/絕對(duì)誤差//計(jì)算的實(shí)際值/實(shí)際值/絕對(duì)誤差
plotyy（ktokt）;%總產(chǎn)值的仿真累加值（to）隨時(shí)間變化曲線與實(shí)際累加值（t）隨時(shí)間變化曲線
%在這里設(shè)置斷點(diǎn)以查看這兩條曲線%
plotyy（kakf）;%總產(chǎn)值的仿真值（a）隨時(shí)間變化曲線與實(shí)際值（f）隨時(shí)間變化曲線

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????2474??2005-05-22?15:06??bp源程序\bp.m

?????文件???????1988??2005-06-13?21:54??bp源程序\bp2.m

?????文件??????28792??2004-08-11?20:11??bp源程序\sourse\ANN.mat

?????文件???????8442??2004-08-12?10:20??bp源程序\sourse\AnnData.txt

?????文件???????1388??2004-08-12?10:32??bp源程序\sourse\Predict.asv

?????文件???????1371??2004-08-12?10:27??bp源程序\sourse\Predict.m

?????文件????????273??2004-08-12?10:34??bp源程序\sourse\說(shuō)明.txt

????..AD...?????????0??2006-03-18?19:39??bp源程序\sourse

?????目錄??????????0??2006-03-18?19:38??bp源程序

-----------?---------??----------?-----??----

????????????????44946????????????????????10