wnetsim.m

新的神经网络算法源程序

function  z = wnetsim(th, u, y)
%Wnetsim: Simulation using a WaveNet-ARX model
%
%        z = wnetsim(THETA, u, y0) 
%
% z: Simulated output
%
% THETA: the parameters of a WaveNet-ARX model
% u:     the input signal (column vector),
%        for multi-input systems, u=[u1 u2 .. um], each ui is a column vector.
% y0:    initial values of z, set to zeros by default

%By Qinghua Zhang. March, 1994.


[ur,uc]=size(u);

[t,d,w,c,b, nn] = decompth(th);

na=nn(1); nb=nn(2:uc+1); nk=nn(uc+2:2*uc+1);

[yr, xu] = io2reg([(1:ur)' u], [0 nb nk]);

z = zeros(ur,1);
if nargin>=3, z(1:length(y))=y; end

nmax=max([na+1 nb+nk])-1;

for k=nmax:ur-1
   z(k+1)=wavenet([z((k):(-1):(k+1-na))' xu(k-nmax+1,:)]', t, d, w, c, b);
end