数学建模社区-数学中国
标题:
求助
[打印本页]
作者:
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]
i' P/ u; n! H: x+ ?: `& p# K9 P6 K$ h
t=[724.570000000000,746.620000000000,778.270000000000,800.800000000000,827.750000000000,871.100000000000,912.370000000000,954.280000000000,995.010000000000,1037.20000000000]
' c8 h+ ^# I) E; w4 a5 j* d8 ~8 \
[pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化
# F+ U6 ?9 E. D; n8 ]
. y- v0 j' s$ z, E
net=newff(minmax(pn),[9,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络
V( v2 Z' t2 _% L
net.trainParam.show=2000; % 训练网络
' Y8 X6 y) f9 d1 p) _, r
net.trainParam.lr=0.01;
a5 U7 |& j* r. I4 L5 y
net.trainParam.epochs=10000;
" H! Z& C8 |1 Y
net.trainParam.goal=1e-5;
% ~+ L% `! A! ~9 p6 K" T. {
[net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络
, ^- Y2 \2 w K* M) u+ ~
9 q) h$ b' C1 p& @. k' r) W5 ^% {
an=sim(net,pn); %对BP网络进行仿真
- a3 H9 o5 G5 E& ]! O3 F3 c% }5 O/ I
a=postmnmx(an,mint,maxt); %还原数据
" h3 {/ z3 R3 y% p: `
x=2001:2010;
9 X2 d& R2 B \
newk=a(1,
;
; S; s- i* i7 U; o( m1 _+ W9 H
plot(x,newk,'r-o',x,t,'b--+')
) D: g( n! g8 H4 G* y% o, Z% |
legend('网络输出值/万人','实际人口数/万人');
( i: q5 z: b$ s0 b
xlabel('年份');
" h' [* i0 ]3 I2 _/ t- r0 D
ylabel('人口数');
" T; j: @6 f$ m2 v6 B
title('预测人口数对比表');
, Z. M+ s$ a+ @5 z8 l9 _
0 u5 Q! o4 s+ h* s2 g0 v1 _
pnew1=[1079.38044723977,1125.65457635540,1173.91252408743,1224.23933874207,1276.72371473226,1331.45814888992,1388.53910347951,1448.06717620014,1510.14727747579,1574.88881534617]
; x8 k; R: Q8 {
pnew=pnew1
! X" w ~0 f" F9 k0 h k
pnewn=tramnmx(pnew,minp,maxp);
9 T* U+ \7 F! N' k& T; c
anewn=sim(net,pnewn); %对BP网络进行仿真
- B9 \, x% {( o
anew=postmnmx(anewn,mint,maxt); %还原数据
9 |% g$ q( V4 s- w' m* a, ? B
y=anew'
* `0 U3 Q' n$ G
5 x; y& H% {1 S2 V& P- N
, }/ r& @2 ~8 I) b( H7 g! n: K8 {
预测的结果怎么都一样呀？》？？求助
作者:
wssl103050
时间:
2012-6-19 22:57
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
;