EDA365电子论坛网

标题: sim函数 [打印本页]

作者: cichishia    时间: 2020-8-3 10:40
标题: sim函数
我用matlab建了一个神经网络,newff实用全部缺省项查了sim传递函数是tansig和purelinsim的结果就是和手工计算结果不一
3 B, `: Z* u! Iiw=net.iw{1,1}
* o! R% e7 q. n. x; i& Vlw=net.lw{2,1}
# G+ s  _9 \, d, U( U7 qb1=net.b{1}5 Y& a, w- u" p4 c9 Y% k3 Q
b2=net.b{2}  s, {" J* w( l: w
h=tansig(iw*input+b1)' E3 E5 ?0 e, d
output=purelin(lw*h+b2)5 t9 S7 E+ s  ^& B3 i5 U
用了以上的指令
; X( s& `2 D# T4 M/ Z8 d8 @求教大神) }$ t2 Z0 J  G# j3 T* i

作者: shuddkk    时间: 2020-8-3 11:12
帮你顶一下
作者: zzz.dan    时间: 2020-8-6 10:53
%% 清空环境变量
6 o" z  O% z/ t/ H- [4 G7 v" aclear all5 d2 X1 P3 U& r: l% n0 Z! Z% |
clc) G+ L9 ^* Z  \9 E0 U7 Y* T( d- T
%% 训练数据
/ }, G$ S7 h  {' ?* C5 P1 qP = [0 1 2 3 4 5 6 7 8 9 10];2 l" R! X* z( K7 `4 J2 X
T = [0 1 2 3 4 3 2 1 2 3 4];. j1 _& X1 f$ x' j4 l5 C( t; f
%% 创建网络
7 v$ u% I8 {3 R: ^$ ?; u! s( |net = newff(minmax(P),[5 1],{'tansig','purelin'},'trainlm');$ U( r/ B, O7 W( C4 r# t9 a
% 设置训练参数
' X  F  ^5 _  Lnet.trainParam.epochs = 500;+ o& P2 [6 L; ~; N  l
net.trainParam.show=10;
9 }- l: k1 ]( D' U; Z: F% Fnet.trainParam.goal=1e-2;9 c( x* X# y+ D; O& e9 g0 C
net.trainParam.lr=0.1;
( O. N) Z! j& o* K%% 训练网络1 `# Y; U, C& r& [4 M( P- `5 d
net = train(net,P,T);! i2 i) c6 M8 O4 a# V, |$ D
%% 记录权值和阈值) L2 y5 L* i9 h4 k+ u
w1=net.iw{1};
: o" S% |; h" C. uw2=net.lw{2};
4 _! r. h# q2 O; f6 B9 U4 a3 s/ eb1=net.b{1};
& g; q  L/ \# t2 t; b2 u0 S! hb2=net.b{2};' H8 t3 V0 v9 z2 ~8 @' C
%% 仿真测试
( F% T6 `1 {3 k. [0 @Y1= sim(net,P)
& `# Z" T9 a: K5 O- [) x2 n%% 手动计算一0 R6 ^# f, S+ [: ]
a=tansig(w1*P,b1);3 \; ^4 L9 O% B8 g# X! q
Y2=purelin(w2*a,b2)
( e9 b# ~$ b0 k; @) `%% 手动计算二* m! a3 b, L: H6 d& k" |1 }
b=tansig(w1*P+repmat(b1,1,size(P,2)));
  L# N, |+ P2 F" G) D7 {Y3=purelin(w2*b+repmat(b2,1,size(P,2)))
" f- l* i2 O0 w- h9 [+ t( v%% 判断是否相等2 `# @0 [6 N6 {
flag=isequal(Y1,Y2) && isequal(Y1,Y3)
作者: zzz.dan    时间: 2020-8-6 10:53
运行结果:
. t# ~0 Y8 `  i3 [+ L4 PY1 =
/ J; U/ N5 @! H. w1 ^# B    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590
: Y/ x  o$ X: f# ?8 [
3 M, A. p! C5 a# q4 j" l+ K. S7 E( NY2 =1 \4 L2 B  a8 w3 i7 \& b
    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590
0 _- t& D2 {/ ?% h, ^3 n
; n. z7 V+ X- K7 o/ L) Z* KY3 =
- x; N+ b' @) c    0.0012    1.0011    1.9858    3.0039    3.9442    3.1098    1.8481    1.1784    1.8534    3.1175    3.9590
; s0 }, Q, D: }. C6 D3 j
, p0 \# b: M! n0 i; X( l; fflag =
" ~( `. G0 y* J$ \6 P3 G( N/ S     16 K- a: o: A& }. K4 ?3 g& P5 |- x

作者: llbnmo    时间: 2020-8-6 10:57





欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2