sim函数

cichishia

我用matlab建了一个神经网络，newff实用全部缺省项查了sim传递函数是tansig和purelinsim的结果就是和手工计算结果不一
( ?6 d& K5 c$ f, Ciw=net.iw{1,1}
lw=net.lw{2,1}
. Z# [% b5 t8 U, o) U5 E# qb1=net.b{1}
b2=net.b{2}
/ ^5 w& d9 H1 X$ T5 \4 Hh=tansig(iw*input+b1)
output=purelin(lw*h+b2)
% K$ A1 U* r5 N) [+ E用了以上的指令
求教大神
& x3 @9 N4 U, N: K/ S

shuddkk

帮你顶一下

zzz.dan

%% 清空环境变量
: r5 L% s* z/ O2 dclear all
+ N/ o4 t+ n0 g9 O0 S! h' w/ \' `$ S- Bclc
) C' E0 `  F, u%% 训练数据
' r" a) u5 ^* ~* b" W2 BP = [0 1 2 3 4 5 6 7 8 9 10];
3 W0 D* x5 K: d6 MT = [0 1 2 3 4 3 2 1 2 3 4];
0 ]& D4 \& ]3 C: K- g%% 创建网络
  u4 X" e" Q0 y; s2 u$ L' z' R4 wnet = newff(minmax(P),[5 1],{'tansig','purelin'},'trainlm');
& ~* V# s7 Y( \: u7 [% 设置训练参数
0 Q5 ~9 t4 O0 N7 c( `: W' E2 A$ a  h: pnet.trainParam.epochs = 500;
net.trainParam.show=10;
# D6 Y' n$ g7 E0 w8 ~net.trainParam.goal=1e-2;
net.trainParam.lr=0.1;
%% 训练网络
net = train(net,P,T);
: E% O8 l6 _, A%% 记录权值和阈值
w1=net.iw{1};
5 b( T9 ^+ D) F/ [2 r/ i6 U3 Aw2=net.lw{2};
% D/ x2 r" ?5 \2 a0 a) A; jb1=net.b{1};
6 t: y. _$ q! x, J1 @. N0 }b2=net.b{2};
%% 仿真测试
Y1= sim(net,P)
# u0 S* F: ?7 V% \%% 手动计算一
5 F& n5 c& a1 z8 ma=tansig(w1*P,b1);
Y2=purelin(w2*a,b2)
, Q. O- a# ~6 f: Q; \+ X%% 手动计算二
b=tansig(w1*P+repmat(b1,1,size(P,2)));
Y3=purelin(w2*b+repmat(b2,1,size(P,2)))
/ O5 V( M" y( V8 G, f* j; o%% 判断是否相等
3 U( P0 J$ W7 r( ~+ k8 G, Iflag=isequal(Y1,Y2) && isequal(Y1,Y3)

zzz.dan

运行结果：
Y1 =
    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590

) ~& w8 ?- U/ d# L2 uY2 =
7 [& Y' c, c: L) B    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590

8 y; L; N* O2 H& X+ ~Y3 =
1 H7 h, c" e9 U" E" J    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590

  v  D9 R7 n* ?  \; zflag =
* {$ W0 Z' u3 f5 D, K     1

llbnmo