- 在线时间
- 28 小时
- 最后登录
- 2013-2-22
- 注册时间
- 2012-4-29
- 听众数
- 4
- 收听数
- 0
- 能力
- 0 分
- 体力
- 385 点
- 威望
- 0 点
- 阅读权限
- 30
- 积分
- 152
- 相册
- 0
- 日志
- 0
- 记录
- 2
- 帖子
- 72
- 主题
- 5
- 精华
- 0
- 分享
- 0
- 好友
- 28
升级   26% TA的每日心情 | 开心 2013-2-22 18:11 |
|---|
签到天数: 54 天 [LV.5]常住居民I
- 自我介绍
- 对建模很感兴趣
 群组: 数学建模培训课堂1 群组: 数学建摸协会 |
p=[724.570000000000,739.774223126620,771.489090673942,804.563606600594,839.056060404793,875.027240514377,912.540541418428,951.662075391767,992.460789009143,1035.00858465455]9 C h% B$ a1 |! ?, ]/ B T5 q
t=[724.570000000000,746.620000000000,778.270000000000,800.800000000000,827.750000000000,871.100000000000,912.370000000000,954.280000000000,995.010000000000,1037.20000000000]
; O2 z7 F8 M! J5 R [pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化; x; G$ a- n4 [- P0 c$ t! J
' V: u. P& D5 q
net=newff(minmax(pn),[9,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络
4 K3 H0 U6 W6 L, o1 G6 xnet.trainParam.show=2000; % 训练网络' q: f- Y( [- J, b) s5 f
net.trainParam.lr=0.01;9 o, P6 r6 r4 _6 P1 |
net.trainParam.epochs=10000;1 ^5 W# s* R+ M' Z9 E7 v1 |8 k) f
net.trainParam.goal=1e-5;
( Z* L' P$ ^/ k8 X) g[net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络
7 I1 T, ~- e& p) ]" l0 s; n, [' T0 X# j
an=sim(net,pn); %对BP网络进行仿真0 i. [& N& r5 H9 ^. h
a=postmnmx(an,mint,maxt); %还原数据
: c% ]' }, q V) P- [; `+ ux=2001:2010;
2 B* q) _' `- b, m$ Bnewk=a(1, ;6 h( R; x% `& |
plot(x,newk,'r-o',x,t,'b--+')
1 y7 j7 `, K% j2 V& G- m7 \legend('网络输出值/万人','实际人口数/万人');
9 u L) t; s% w) S- yxlabel('年份');
9 _/ Q" Q, F5 y9 ]4 }" Z% s. a$ Tylabel('人口数');! V. }2 p5 ?/ g, V: t3 W% G
title('预测人口数对比表'); d0 \3 E3 g3 N5 o9 a9 Y/ |2 r- m
* @' L' k0 P" E; p3 e$ kpnew1=[1079.38044723977,1125.65457635540,1173.91252408743,1224.23933874207,1276.72371473226,1331.45814888992,1388.53910347951,1448.06717620014,1510.14727747579,1574.88881534617]
. w6 @2 G1 E ^4 m4 R4 ypnew=pnew1! N- W! W; W" B* z, M h
pnewn=tramnmx(pnew,minp,maxp);1 v" B* }/ H, o9 |5 G' U5 G5 Q. j
anewn=sim(net,pnewn); %对BP网络进行仿真
% \3 u$ H+ q% _3 J5 J+ f8 kanew=postmnmx(anewn,mint,maxt); %还原数据
' o4 x. y+ `( xy=anew'
: ?% z& h" M$ g1 }% L0 G1 U% F3 y% @# B* k
6 V$ F& I6 K# J7 Z9 O预测的结果怎么都一样呀?》??求助 |
zan
|