% FF2NARX feedforward to narx network
%    [narx] = ff2narx(net,d)
%     net : feedforward network
%       d : delay lines
%    narx : narx network
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [narx] = ff2narx(net,d)

% Get input range for real inputs
R = net.inputs{1}.range;
n = size(R,1)-length(d);
Pr = R(1:n,:);

% Get layer dimensions
h = zeros(size(net.layers))';
for i = 1:length(h)
   h(i) = net.layers{i}.dimensions;
end

% Create narx network
narx = newnarx(Pr,0,d,h,{'tansig','purelin'});
narx.biasConnect = zeros(size(narx.biasConnect));

% Set narx weights from ffnet weights
IW = net.IW{1};
narx.IW{1} = IW(:,1:n);
narx.LW{2,1} = net.LW{2,1};
narx.LW{1,2} = IW(:,n+1:end);