数学建模社区-数学中国
标题:
求助
[打印本页]
作者:
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]
" V B" U. P4 E3 j% T
t=[724.570000000000,746.620000000000,778.270000000000,800.800000000000,827.750000000000,871.100000000000,912.370000000000,954.280000000000,995.010000000000,1037.20000000000]
2 t4 b; O8 \" O
[pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化
. k' @6 ~' D' T4 L( P K
& k) m4 O( R6 m; k$ M
net=newff(minmax(pn),[9,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络
4 y" U# w# l. K
net.trainParam.show=2000; % 训练网络
4 X" A& |" O0 ]- b$ ~' `0 ~. X
net.trainParam.lr=0.01;
" J$ C6 Z( x: F* x
net.trainParam.epochs=10000;
+ d7 I2 y! n/ k
net.trainParam.goal=1e-5;
" I- A. r6 B' w& e( U5 x6 ]
[net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络
0 o8 v$ a. _ s+ a+ t0 J; P
% b6 H. X2 @0 J2 D* n/ N
an=sim(net,pn); %对BP网络进行仿真
5 ~* p; J$ h, W% n* \# W+ W3 I
a=postmnmx(an,mint,maxt); %还原数据
& J. z7 N; }. E
x=2001:2010;
% Z( q6 r' l/ k, W
newk=a(1,
;
. [: k) _( s# [4 |6 z, R
plot(x,newk,'r-o',x,t,'b--+')
% C! J0 a* x+ O# H5 z }1 y) t4 B
legend('网络输出值/万人','实际人口数/万人');
* m6 _% Q' j+ w9 h8 L$ W& m
xlabel('年份');
5 c6 g7 E$ W% J6 o; G: {
ylabel('人口数');
, f- A6 w. p% E% ]9 d+ r. K
title('预测人口数对比表');
; q+ b- H2 ]; E3 |9 x1 Z
" A/ `8 _: I' Q0 h4 N, y# @- Q) H
pnew1=[1079.38044723977,1125.65457635540,1173.91252408743,1224.23933874207,1276.72371473226,1331.45814888992,1388.53910347951,1448.06717620014,1510.14727747579,1574.88881534617]
0 [/ a. j5 M( r4 ~) ^
pnew=pnew1
: I0 ?9 \3 D* q7 K8 e
pnewn=tramnmx(pnew,minp,maxp);
2 a7 k' s2 J& {9 _) c
anewn=sim(net,pnewn); %对BP网络进行仿真
6 ?/ ~- S1 D/ U _+ d( a9 N3 X
anew=postmnmx(anewn,mint,maxt); %还原数据
- |& Z5 Z7 q d/ Y1 Y
y=anew'
' e2 \8 H5 L$ C# \" [! Q/ h. m
( m4 P5 f3 y' E3 [! c+ c H: N
) F$ P: s4 ?: I8 v
预测的结果怎么都一样呀?》??求助
作者:
wssl103050
时间:
2012-6-19 22:57
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5