plantv.asv

氧乐果控制过程

clf reset
pausetime=0.1
load datapv.mat
k1=0.000024015;in=in1*k1;
out=out1/2^15
x11=in(1:50:4000)
x21=out(1:50:4000)
p=[x11';x21']
T=x21'
plot(p,T,'+')
title('training vectors');xlabel('input vector p');ylabel('target vector t');pause
[R,Q]=size(p)
[S2,Q]=size(T)
s1=5
[w10,b10]=nwtan(s1,r)
w20=rands(s2,s1)*0.5
b20=rands(s2,1)*0.5
plot(p,purelin(w30*tansig(w10*p,b10),b20))
pause2(pausetime)
hold on
h=get(gca,'children')
h=h(1)

disp_freq=10
max_epoch=2000
err_goal=0.001
lr=0.01
lr_inc=1.05
lr_dec=0.7
err_ratio=1.04

TP=[disp_freq max_epoch err_goal lr lr_inc lr_dec err_ratio]
flops(0)
[w1,b1,w2,b2,epoch,TP]=trainbpa(w10,b10,'tansig',w20,b20,'purelin',p,T,TP)
totalflops=flops;
delete(h)
h=plot(p,purelin(w2*tansig(w1*p2,b1