sim函数

cichishia

我用matlab建了一个神经网络，newff实用全部缺省项查了sim传递函数是tansig和purelinsim的结果就是和手工计算结果不一
- c6 z* Y; M/ t) z  _iw=net.iw{1,1}
lw=net.lw{2,1}
! X9 I+ n) A, A; s1 Nb1=net.b{1}
7 w9 {+ K- i' C3 R: sb2=net.b{2}
/ t) b) L9 D  Z7 Z1 w! oh=tansig(iw*input+b1)
output=purelin(lw*h+b2)
- n0 l* }% ~* E. `& }用了以上的指令
6 e/ G9 i5 v3 J求教大神
- \9 s4 u1 ~/ N% c6 j2 s

shuddkk

帮你顶一下

zzz.dan

%% 清空环境变量
) y% n' L0 p% l. {7 ?clear all
clc
%% 训练数据
1 R: u  b- V  |- ~! rP = [0 1 2 3 4 5 6 7 8 9 10];
T = [0 1 2 3 4 3 2 1 2 3 4];
%% 创建网络
/ j) k, n9 E. y* f1 l6 Inet = newff(minmax(P),[5 1],{'tansig','purelin'},'trainlm');
2 S, |" `9 U; O6 h' x% 设置训练参数
3 o8 @: ~: X5 A8 ynet.trainParam.epochs = 500;
net.trainParam.show=10;
: }9 D& z! C( u$ Z7 Z' [( y0 Gnet.trainParam.goal=1e-2;
net.trainParam.lr=0.1;
3 k7 r. H- m5 ~. Y%% 训练网络
  n8 O8 b6 e7 [( bnet = train(net,P,T);
%% 记录权值和阈值
, F. Q% P! C! V. d3 j0 Ew1=net.iw{1};
! ]+ j" r! ^. Nw2=net.lw{2};
  i! s* b, I! s$ w  Vb1=net.b{1};
1 _( b- o/ ?, [5 Yb2=net.b{2};
& _! W9 c8 |( c' m+ l2 X6 k6 h' ~%% 仿真测试
" `$ p3 o5 B7 X4 iY1= sim(net,P)
, V7 h+ m* f4 Y* v5 }/ z%% 手动计算一
a=tansig(w1*P,b1);
0 F2 V! h7 ]! x5 |1 s5 A- |Y2=purelin(w2*a,b2)
- w: _+ B' _5 Z2 _7 X%% 手动计算二
* p* H4 ]* O9 ?) U& H% }) U3 {) z. eb=tansig(w1*P+repmat(b1,1,size(P,2)));
Y3=purelin(w2*b+repmat(b2,1,size(P,2)))
# D/ g# J0 F4 C%% 判断是否相等
flag=isequal(Y1,Y2) && isequal(Y1,Y3)

zzz.dan

运行结果：
! D  n& P  i" RY1 =
+ a% y) N& y0 q$ X    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590
/ n/ w0 r% F! S& B/ s  _2 j$ f3 z
Y2 =
/ ~/ P" e" t3 D6 ~8 N    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590
* ]& K, v; ?9 |2 z3 m% c  A
Y3 =
  K/ K0 _! F& Q    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590
# e- }8 h8 @7 `6 j( X+ g
flag =
     1

llbnmo