数学建模社区-数学中国
标题:
求助
[打印本页]
作者:
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 ^, F. [- u/ J% O2 |
t=[724.570000000000,746.620000000000,778.270000000000,800.800000000000,827.750000000000,871.100000000000,912.370000000000,954.280000000000,995.010000000000,1037.20000000000]
6 V. q( Y( J9 |* [
[pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化
! z: d2 S2 s5 T7 O# R( S# W# ?
8 p6 n4 S# ?: Q1 p+ I4 f
net=newff(minmax(pn),[9,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络
# A; A, p4 T) m/ u( v' L
net.trainParam.show=2000; % 训练网络
z* e9 N3 ^8 c( o/ _
net.trainParam.lr=0.01;
- o( Q) m* Q: s! F# {6 ^
net.trainParam.epochs=10000;
9 d) Z0 |4 `" ]! a3 Z
net.trainParam.goal=1e-5;
/ n* T, p" h. y! \ H; B4 w
[net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络
/ v/ [1 i' B+ ^4 p* B& B3 ^% ~! w9 Y3 R; w( L
" o3 p1 ~9 ]" L3 k! q, Q
an=sim(net,pn); %对BP网络进行仿真
" I! J7 |9 L+ w6 s6 |& y, p# |
a=postmnmx(an,mint,maxt); %还原数据
7 B9 y T9 _0 Y5 ]* s
x=2001:2010;
6 e) u6 m; w+ L/ u6 Y7 y& w
newk=a(1,
;
0 R# }$ p1 v! X( r% d
plot(x,newk,'r-o',x,t,'b--+')
X" F* c' {/ K* @' |
legend('网络输出值/万人','实际人口数/万人');
8 m: j2 M4 t( d4 @2 g- |9 Q
xlabel('年份');
% N V+ W. N0 }4 T' S. I9 r+ {3 }
ylabel('人口数');
, N% D m7 [2 B7 D; O. j k
title('预测人口数对比表');
) J" N# z: }5 R' ^
+ c5 B/ K! ]4 c. _
pnew1=[1079.38044723977,1125.65457635540,1173.91252408743,1224.23933874207,1276.72371473226,1331.45814888992,1388.53910347951,1448.06717620014,1510.14727747579,1574.88881534617]
0 E: u, o3 V" [7 U9 P1 g" H
pnew=pnew1
2 N, d. F8 \9 |: T, N' p0 p
pnewn=tramnmx(pnew,minp,maxp);
: N5 Q5 L q. ]2 f$ C
anewn=sim(net,pnewn); %对BP网络进行仿真
% E6 S+ o6 a) B5 \
anew=postmnmx(anewn,mint,maxt); %还原数据
o' R' f- X/ u f
y=anew'
; {; D h* I6 |) R
- I5 V" `6 x+ a! g0 o# {9 g
6 i! _2 d* y. b) l. d- s" T1 Z& L
预测的结果怎么都一样呀？》？？求助
作者:
wssl103050
时间:
2012-6-19 22:57
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
;0 R# }$ p1 v! X( r% d