nn_wyq.m

一个很好的船舶仿真源程序

x1=[rands(1,5000)]*pi;
x2=[zeros(1,5000)];
x3=[zeros(1,5000)];
P=[x1;x2;x3];
T=-0.1*(x1)+0.053*(x2).*(x2).*(x2)-0.938*(x2);
[R,Q]=size(P);
S1=12;
[S2,Q]=size(T);
pr=[-pi,pi;-0.01,0.01;-0.01,0.01]
net=newff(pr,[S1,S2]);
[W1,B1]=rands(S1,R);
[W2,B2]=rands(S2,S1);
disp_freq=15;
max_epoch=500;
err_goal=0.002;
lr=0.0001;
lr_inc=1.05;
lr_dec=0.7;
momentum=0.95;
err_ratio=1.04;
TP=[disp_freq max_epoch err_goal lr lr_inc lr_dec momentum err_goal];
[W1,B1,W2,B2,TE,TR]=trainbp(W1,B1,'tansig',W2,B2,'purelin',P,T,TP);
plot(TR,net.trainParam.goal);
grid;
x1=[rands(1,5000)]*pi;
x2=[zeros(1,5000)];
x3=[zeros(1,5000)];
P=[x1;x2;x3];
T=-0.1*(x1)+0.053*(x2).*(x2).*(x2)-0.938*(x2);
%T=1.06*x2.*x2.*x2+88*x3+21.24*x2+20*x1;
%T=k*x3+1.24*dt*k*(x2)+0.053*dt*k*(x2).*(x2).*(x2)+3.4*dt*k*(x3);
y=sim(net,P);
figure(2)
plot(x1,T,'*',P,y);
grid;