EDA365电子论坛网
标题:
sim函数
[打印本页]
作者:
cichishia
时间:
2020-8-3 10:40
标题:
sim函数
我用matlab建了一个神经网络,newff实用全部缺省项查了sim传递函数是tansig和purelinsim的结果就是和手工计算结果不一
3 B, `: Z* u! I
iw=net.iw{1,1}
* o! R% e7 q. n. x; i& V
lw=net.lw{2,1}
# G+ s _9 \, d, U( U7 q
b1=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" a
clear all
5 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 q
P = [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 _ L
net.trainParam.epochs = 500;
+ o& P2 [6 L; ~; N l
net.trainParam.show=10;
9 }- l: k1 ]( D' U; Z: F% F
net.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. u
w2=net.lw{2};
4 _! r. h# q2 O; f6 B9 U4 a3 s/ e
b1=net.b{1};
& g; q L/ \# t2 t; b2 u0 S! h
b2=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 P
Y1 =
/ 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( N
Y2 =
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* K
Y3 =
- 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; f
flag =
" ~( `. G0 y* J$ \6 P3 G( N/ S
1
6 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