- 在线时间
- 84 小时
- 最后登录
- 2014-10-4
- 注册时间
- 2013-7-16
- 听众数
- 8
- 收听数
- 1
- 能力
- 0 分
- 体力
- 949 点
- 威望
- 0 点
- 阅读权限
- 40
- 积分
- 372
- 相册
- 1
- 日志
- 0
- 记录
- 0
- 帖子
- 172
- 主题
- 23
- 精华
- 0
- 分享
- 0
- 好友
- 35
升级   24% TA的每日心情 | 开心 2014-10-4 18:19 |
---|
签到天数: 102 天 [LV.6]常住居民II
- 自我介绍
- 本人是在校学生,即将参加有史以来的第一场国赛,希望能在这里获得帮助。
 |
本帖最后由 数学不烦 于 2013-9-9 21:56 编辑
# }5 ^( l0 Y9 R& p5 _) b. T) g* w: G; t; ~9 R C
各位朋友,有一个问题请教你们。 我用BP神经网络做预测的时候(下面是我的代码)遇到一个问题。 首先我尝试在代码里面不把数据归一化,这样的出来的结果很差(R=0.46).但我把数据作归一化 与 反归一化 之后 确实R提高了很多(R=0.99),这时问题就来了,请看最下面的测试结果代码,那里只有输入(没有输出,因为是预测),premnmx 归一化要有输出的变量T,但我只有输入数据,没有输出数据,因此反归一的结果(最后的预测结果)也无法做到。那么我的预测就不是做不了吗? 很纠结,查了很多资料了,还是不知道怎么解决。, l# T" w* j2 B/ x
国赛准备到了,希望热心的朋友懂的不吝赐教,我把问题写得很认真,希望各位能帮我解决这个问题,谢谢了。9 N; a+ B+ l5 n5 ], t
: O! E/ Y+ Z* w/ }% xclear/ \6 S2 X# G# ]: i* X! J
clc) \; [- C9 u- ?2 e `( j
%输入数据--以一列作为一个输入数据
! }% G/ F+ z- \9 L7 I0 L% sx1=[227.5 286.9 236.3 300.6 243.0];
3 H6 D2 d& o# A3 Y! x* j4 P7 Zx2=[286.9 236.3 300.6 243.0 354.1];
; j+ m* G8 r: B% g. y) |5 B: B: }8 k1 Qx3=[236.3 300.6 243.0 354.1 294.1];* p& ^, n$ d# p! j/ {6 l
x4=[300.6 243.0 354.1 294.1 361.4];
/ \" ^# e8 [. f& \x5=[243.0 354.1 294.1 361.4 333.7];
0 l/ ]0 ~! j7 c, j3 _/ `1 ^x6=[354.1 294.1 361.4 333.7 426.4];% z9 u% }: z# I% X" I
x7=[294.1 361.4 333.7 426.4 434.4];
3 _3 d, P6 t5 q" i' Xx8=[361.4 333.7 426.4 434.4 449.9];
" L7 C1 F+ y7 O3 I. }x9=[333.7 426.4 434.4 449.9 431.9];
6 K# q( t/ T; }0 T9 M0 ^$ R
. O$ I Y: H0 s( C% IP=[x1;x2;x3;x4;x5;x6;x7;x8;x9]';
* a7 B! m0 w" ]. ]7 q' t%输出数据--以一列作为一个输出数据
- ~7 Q6 v7 X# Qy1=[354.1 294.1 361.4 333.7 426.4];7 v( S% s, ` p/ u5 u X) O
y2=[294.1 361.4 333.7 426.4 434.4];$ b% z! j7 e5 t' I6 w
y3=[361.4 333.7 426.4 434.4 449.9];
! F5 d/ W; D: L/ Yy4=[333.7 426.4 434.4 449.9 431.9];
" x3 i8 A- e- i" Ay5=[426.4 434.4 449.9 431.9 454.0];2 T9 T" l$ e$ x7 X
y6=[434.4 449.9 431.9 454.0 655.6];
# P) C* Y* |+ v+ W- M: X6 Ey7=[449.9 431.9 454.0 655.6 907.0];+ X# f) H' ?2 U! m1 p! b- \
y8=[431.9 454.0 655.6 907.0 1112.5];
* R. Z3 q# v9 @& C! ry9=[454.0 655.6 907.0 1112.5 1401.1];% \$ L" |- [: ~
6 c; k T% `$ z/ Z5 j. m1 `) W6 nT=[y1;y2;y3;y4;y5;y6;y7;y8;y9]';
9 i2 k7 z( U! Y9 L7 s" m
, o }0 i% S+ v0 S' Q9 E%归一化
" {8 x5 D6 M4 [% [P,minp,maxp,T,mint,maxt] = premnmx(P,T);
|( s/ Q& W# G. U" ?7 b- `5 U& C/ y2 h+ y%建立网络
' Y6 Z' A, A& ~6 \- j%其中[10,3]代表有一个隐含层,里面有10个神经元;“3”代表有三个输出数据
7 H u: p7 z: l, t4 e. f/ c p%{'tansig','purelin'}与[10,3]对应,意思为隐含层 与 输出层
% j, A2 `! o: n% u" @%'traingdm' 表示调整权重,缩小误差的数学函数
' M/ d2 x9 ?/ k% W%↓创建一个新的前向神经网络 # ^* a/ T3 |2 y4 F
net = newff(minmax(P),[6,5],{'tansig','purelin'},'traingdm')7 K) ]6 V5 ]; f5 l
. p4 V9 S7 N& O" l) A
% 当前输入层权值和阈值
% d3 }( N3 z p4 h% ^; `' W DinputWeights = net.IW(1,1)+ n: N" o; l2 d& V$ U
inputbias = net.b(1)% R+ z8 A# D% l- o
* j" W; a5 {2 y" o9 ~
% 当前网络层权值和阈值
$ Q0 j$ {# [# u7 a6 G: E8 tlayerWeights = net.LW(2,1)
8 f& X7 Y, L0 _" A1 _5 v5 Zlayerbias = net.b(2)
$ Q7 f% m: J U: B* ^' ~7 i2 o6 ]9 G1 K. I* E5 X# H
% 设置训练参数
6 u* K1 i ~) xnet.trainParam.show = 50;; f7 U5 J/ U3 v; B; y0 a$ h
net.trainParam.lr = 0.05;5 U. @' {4 S$ N" k. ~& J5 p
net.trainParam.mc = 0.9;
$ d, z/ Z. W" z' nnet.trainParam.epochs = 10000; %每计算10000次出一次结果9 q7 L' C" m* R: w) W- i, j0 P
net.trainParam.goal = 1e-6; %期望达到的误差(认为是合适的误差)3 Z- p- O- d# q& ]" _9 y# W
8 s9 S7 d( ?% D4 Q$ d
% 调用 TRAINGDM 算法训练 BP 网络
7 o2 r2 z, k% f" u# n% l3 t[net,tr] = train(net,P,T);
( c' [$ a B8 h. w7 u1 r' L: M5 F5 W3 [& E' C8 n8 ` v: C6 R
% 对 BP 网络进行仿真
! w' X8 w; C/ T. [A = sim(net,P); %A为输出结果* t; F, \# m0 @/ v$ N, U. h
! Q$ n3 `6 |# Q, V
%反归一化. [: h0 H6 h( q$ D
% A = postmnmx(A,mint,maxt);7 Q+ V9 |" @; v7 y/ F
3 W- k D3 ?9 X- r8 g/ j W4 D+ Y
% 计算仿真误差 8 _" R$ ~* d+ h3 b. j( n
E = T-A; %原来的输出 - 训练输出 = 误差
. d& ^- t3 z; n2 }5 gMSE = mse(E) %输出误差4 ]2 ~+ G2 @8 n! p; b+ ?* g
0 j& W7 W- w( R" W' D7 u+ C! H9 C%下面是输入数据,测试结果
' `9 s+ o- R4 c7 P9 m" Rx=[454.0 655.6 907.0 1112.5 1401.1]';4 _6 V/ y6 H& u' v
Y = sim(net,x)1 n9 {% b z3 K7 T. x
|
zan
|