sim函数

cichishia

我用matlab建了一个神经网络，newff实用全部缺省项查了sim传递函数是tansig和purelinsim的结果就是和手工计算结果不一
) k: n. S) T8 O7 h( C: biw=net.iw{1,1}
lw=net.lw{2,1}
b1=net.b{1}
b2=net.b{2}
h=tansig(iw*input+b1)
  S& T( [% S+ T" U. u) coutput=purelin(lw*h+b2)
用了以上的指令
# h+ B, {$ z1 J" n求教大神
6 c4 Z. L) m) }: C: W5 p

shuddkk

帮你顶一下

zzz.dan

%% 清空环境变量
4 K5 b& I- L6 W+ k# c1 H  n8 xclear all
clc
%% 训练数据
P = [0 1 2 3 4 5 6 7 8 9 10];
T = [0 1 2 3 4 3 2 1 2 3 4];
  u* ]; p7 `8 f%% 创建网络
; A+ j$ \! s4 y# V8 Vnet = newff(minmax(P),[5 1],{'tansig','purelin'},'trainlm');
% 设置训练参数
" f8 w8 \$ a0 znet.trainParam.epochs = 500;
net.trainParam.show=10;
' I: q" N: k8 V& _4 B' o  v( ]( J& Hnet.trainParam.goal=1e-2;
net.trainParam.lr=0.1;
%% 训练网络
net = train(net,P,T);
%% 记录权值和阈值
w1=net.iw{1};
w2=net.lw{2};
( k3 ]- ]1 z; A' ?b1=net.b{1};
b2=net.b{2};
0 g! W/ d& C% t5 a$ x%% 仿真测试
4 r; k( E) n$ E- N9 N3 BY1= sim(net,P)
# D" d* F5 @/ {2 w. ~; j%% 手动计算一
a=tansig(w1*P,b1);
Y2=purelin(w2*a,b2)
0 g: a# ]* ?+ \0 n+ O9 |%% 手动计算二
4 b( e1 ]: j2 W% i- L) Mb=tansig(w1*P+repmat(b1,1,size(P,2)));
Y3=purelin(w2*b+repmat(b2,1,size(P,2)))
- E- Z. Q. U; `: d7 ]4 h3 Y2 R%% 判断是否相等
flag=isequal(Y1,Y2) && isequal(Y1,Y3)

zzz.dan

运行结果：
$ i2 c9 W* H3 w5 y( V! }# cY1 =
; o/ r/ `  t2 A, t. H    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590
* N0 A7 E1 Y9 O' P' n6 [$ ~3 }
Y2 =
) i2 _2 \+ o) [# R  K    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590

Y3 =
1 M: _, F7 F8 P: z. i7 v" U5 x    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590
4 P' |# H6 p8 `4 T
& T. z! B/ }" \$ X6 F3 A8 Hflag =
     1

llbnmo