%% noisy-damage

% Load data from files
data = load('noisy-damage.txt');

% Set training data to first half of data
trainInputs  = data(1:500,2)';
trainTargets = data(1:500,1)';

% Set performance data to second half of data
inputs = data(501:end, 2)';
targets = data(501:end, 1)';

% Take a few samples
inputsamples  = data(1:10,2)';
targetsamples = data(1:10,1)';

plot(trainInputs, trainTargets, '.blue')
figure(); plot(damage_transform(trainInputs), damage_transform(trainTargets), '.red')


%%% Transformations
trainInputs     = damage_transform(trainInputs);
trainTargets    = damage_transform(trainTargets);
inputs          = damage_transform(inputs);
targets         = damage_transform(targets);
inputsamples    = damage_transform(inputsamples);
targetsamples   = damage_transform(targetsamples);

size(trainInputs)
size(trainTargets)
size(inputs)
size(targets)
size(inputsamples)
size(targetsamples)


% Create new ff network
net = newff(inputsamples, targetsamples, [], {'tansig'}, 'trainlm');

%net.biasConnect(layerNumber) = booleanValue;
net.biasConnect(1) = 1;
net.trainParam.epochs = 500;
[net trainPerformance] = train(net, trainInputs, trainTargets);
trainPerformance

%[Y,Xf,Af,E,perf] = sim(net, trainingInput, [], [], trainingTarget)
%[Y,Xf,Af,E,perf] = sim(net, inputs);


%%%%%%%%%% 
Gives me this output


Warning: Imaginary parts of complex X and/or Y arguments ignored
> In damage at 19

ans =

     1   500


ans =

     1   500


ans =

     1   500


ans =

     1   500


ans =

     1    10


ans =

     1    10

??? Error using ==> plus
Matrix dimensions must agree.

Error in ==> calcperf2 at 163
        N{i,ts} = N{i,ts} + Z{k};

Error in ==> trainlm at 253
[perf,El,trainV.Y,Ac,N,Zb,Zi,Zl] = calcperf2(net,X,trainV.Pd,trainV.Tl,trainV.Ai,Q,TS);

Error in ==> network.train at 219
  [net,tr] = feval(net.trainFcn,net,tr,trainV,valV,testV);

Error in ==> damage at 44
[net trainPerformance] = train(net, trainInputs, trainTargets);