数学建模社区-数学中国
标题:
求助
[打印本页]
作者:
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]
4 {+ |$ T S8 O7 C5 o3 T, d& M
t=[724.570000000000,746.620000000000,778.270000000000,800.800000000000,827.750000000000,871.100000000000,912.370000000000,954.280000000000,995.010000000000,1037.20000000000]
c6 t0 L4 C6 E( V7 t
[pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化
2 Q9 h! F1 v" a9 z0 E
1 ?7 I; G* T( s3 g
net=newff(minmax(pn),[9,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络
! J2 s" g- h# t3 y: u h* z. o
net.trainParam.show=2000; % 训练网络
' |) z. z' u1 ]4 ~2 s; n" {
net.trainParam.lr=0.01;
0 C4 G/ L, L% _, t/ R
net.trainParam.epochs=10000;
, V" O2 G7 I( i+ T: W+ i
net.trainParam.goal=1e-5;
3 K1 [! i; a3 D/ ^" _, `
[net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络
2 }+ r/ t/ D$ p2 K4 I6 a
" n7 j4 \, O5 I4 |% S$ V& |
an=sim(net,pn); %对BP网络进行仿真
/ n" D! e; X3 W% c5 c
a=postmnmx(an,mint,maxt); %还原数据
! h% B" l+ ~2 F, `$ S+ [
x=2001:2010;
0 Y$ X: m/ h( P- S1 R
newk=a(1,
;
9 R* D) ]' s. I
plot(x,newk,'r-o',x,t,'b--+')
3 ?7 Z# a q: q9 q7 K$ n
legend('网络输出值/万人','实际人口数/万人');
" [7 u4 C& l) t6 ]
xlabel('年份');
, _( u$ ]% M" y% F) U) ^- @
ylabel('人口数');
) h/ ~) s! n5 G
title('预测人口数对比表');
' |. e4 w2 P" R0 |2 ^& z- p
: L1 H( |( v2 C" H! v
pnew1=[1079.38044723977,1125.65457635540,1173.91252408743,1224.23933874207,1276.72371473226,1331.45814888992,1388.53910347951,1448.06717620014,1510.14727747579,1574.88881534617]
- b8 R+ Q0 u' _/ k" { K8 v3 G
pnew=pnew1
! f- z% E, e* y" T0 n
pnewn=tramnmx(pnew,minp,maxp);
$ ~& V" L {+ |% D8 ~ l3 b
anewn=sim(net,pnewn); %对BP网络进行仿真
! C9 \: D: `1 b; K& [$ Z& m
anew=postmnmx(anewn,mint,maxt); %还原数据
6 O1 B9 c+ t! d5 g
y=anew'
7 c7 t, [1 K9 T* c
% W$ L* S3 D, J7 {. G
/ J8 `! d- l( L6 I* D' h
预测的结果怎么都一样呀？》？？求助
作者:
wssl103050
时间:
2012-6-19 22:57
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
;9 R* D) ]' s. I