数学建模社区-数学中国
标题:
求助
[打印本页]
作者:
422258511
时间:
2012-6-19 20:47
标题:
求助
p=[724.570000000000,739.774223126620,771.489090673942,804.563606600594,839.056060404793,875.027240514377,912.540541418428,951.662075391767,992.460789009143,1035.00858465455]
9 C/ e$ P3 d: h U; Y
t=[724.570000000000,746.620000000000,778.270000000000,800.800000000000,827.750000000000,871.100000000000,912.370000000000,954.280000000000,995.010000000000,1037.20000000000]
: R, l4 i0 _& o* s6 B8 w
[pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化
n3 m4 f' @5 j: T4 a
. [5 @( a* [) P9 W3 C1 Y
net=newff(minmax(pn),[9,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络
5 _9 q' r* ^: G5 m" T
net.trainParam.show=2000; % 训练网络
' j+ N7 v9 J. B( u1 m
net.trainParam.lr=0.01;
X" u/ J* e8 [0 u: s' F, e
net.trainParam.epochs=10000;
+ x6 B. m& k6 b0 d5 ?
net.trainParam.goal=1e-5;
/ y, F1 Z, b- ^3 B, p
[net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络
$ y: d. ~8 N$ z; q& j3 ]$ _' L
- u# n* K. i2 j0 p/ T7 w/ @) P
an=sim(net,pn); %对BP网络进行仿真
# |+ H/ ^# E" G
a=postmnmx(an,mint,maxt); %还原数据
" @; `& @+ g+ U9 J. ]* }; J0 e# U
x=2001:2010;
/ G @ r0 E5 _8 x
newk=a(1,
;
4 @; {0 p7 `4 T" A
plot(x,newk,'r-o',x,t,'b--+')
' C( ^& S! c$ ]) i: V8 K
legend('网络输出值/万人','实际人口数/万人');
' i: `/ T% C' L. i" @/ S
xlabel('年份');
. X, u% z. ], A* `
ylabel('人口数');
8 Q" ` @) n4 E2 G+ x: r% k
title('预测人口数对比表');
- C% L9 z# P5 _# t- H; w
' A8 `4 ]# p3 U, l
pnew1=[1079.38044723977,1125.65457635540,1173.91252408743,1224.23933874207,1276.72371473226,1331.45814888992,1388.53910347951,1448.06717620014,1510.14727747579,1574.88881534617]
$ O3 S' F/ y C' [& O p5 b
pnew=pnew1
' M7 ?6 x, O9 b( b e* t1 ^
pnewn=tramnmx(pnew,minp,maxp);
9 T% K" ?( p1 x' _- x+ Q
anewn=sim(net,pnewn); %对BP网络进行仿真
! |- c$ z7 _# g7 s. U4 T$ `2 f
anew=postmnmx(anewn,mint,maxt); %还原数据
6 ]+ X4 z! q/ l; |- E. l1 M8 i
y=anew'
5 I* P8 b4 w* ?' m6 e* U
8 t" B& n/ m% N, Z9 Y
( J- P' Q% D6 n- G7 s7 \/ f2 L3 Q
预测的结果怎么都一样呀？》？？求助
作者:
wssl103050
时间:
2012-6-19 22:57
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
;4 @; {0 p7 `4 T" A