数学建模社区-数学中国
标题:
求助
[打印本页]
作者:
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]
$ r, q) \" t/ I, y
t=[724.570000000000,746.620000000000,778.270000000000,800.800000000000,827.750000000000,871.100000000000,912.370000000000,954.280000000000,995.010000000000,1037.20000000000]
, Z5 i* u6 u2 q
[pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化
1 U! a" P0 h+ M% }: w4 f! A& {9 y- ?
/ n& h$ p/ {. ?( k! j1 c
net=newff(minmax(pn),[9,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络
" L9 b* _/ H H# W) V
net.trainParam.show=2000; % 训练网络
7 k5 _9 n i: p, A' ~& V
net.trainParam.lr=0.01;
( E. q; q1 Q% c3 D0 `2 P
net.trainParam.epochs=10000;
4 b8 A3 F9 N6 x+ k5 z! E# [
net.trainParam.goal=1e-5;
+ _8 q2 M G6 |( B3 }
[net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络
B+ l/ @# A- G
$ ?- m+ d8 r; j/ D% {1 Q2 k- M
an=sim(net,pn); %对BP网络进行仿真
( X7 x2 p/ M* {1 L6 V
a=postmnmx(an,mint,maxt); %还原数据
0 }& }5 L' Q0 U7 f1 y
x=2001:2010;
2 p2 x% R* ]. k: g# n' ?
newk=a(1,
;
$ k' r' N$ }7 _
plot(x,newk,'r-o',x,t,'b--+')
6 x1 P. n% l4 u* m
legend('网络输出值/万人','实际人口数/万人');
1 j D1 B9 h/ h9 M
xlabel('年份');
6 T+ S: @- _( J7 E2 b3 X3 e
ylabel('人口数');
4 x9 h9 t$ y) b/ Y5 U5 m/ h/ p
title('预测人口数对比表');
+ E* Q+ t% y8 P& Y' I/ K) U# X, e
- x5 R/ T8 Z+ Z
pnew1=[1079.38044723977,1125.65457635540,1173.91252408743,1224.23933874207,1276.72371473226,1331.45814888992,1388.53910347951,1448.06717620014,1510.14727747579,1574.88881534617]
. ~/ G) G, I! T
pnew=pnew1
! X& q( z/ C. o. d" E9 ?+ U! v
pnewn=tramnmx(pnew,minp,maxp);
5 h' T* f ^; g* k
anewn=sim(net,pnewn); %对BP网络进行仿真
5 o: b- [" U$ B a& t
anew=postmnmx(anewn,mint,maxt); %还原数据
/ h( q3 v. Q8 U. m
y=anew'
, C5 t8 c; u7 c1 l& ^; M; X7 n
, L# p0 z: [& W. _4 k
- d- F5 \# V7 A7 e" l( _
预测的结果怎么都一样呀？》？？求助
作者:
wssl103050
时间:
2012-6-19 22:57
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
;