www.gusucode.com > nnet 案例源码 matlab代码程序 > nnet/GsTimeSeriesExample.m
%% GsTimeSeries % %% load data % Copyright 2015 The MathWorks, Inc. load ph_dataset inputSeries = phInputs; targetSeries = phTargets; %% create network inputDelays = 1:4; feedbackDelays = 1:4; hiddenLayerSize = 10; net = narxnet(inputDelays,feedbackDelays,hiddenLayerSize); %% prepare data [inputs,inputStates,layerStates,targets] = ... preparets(net,inputSeries,{},targetSeries); %% data division net.divideParam.trainRatio = 70/100; net.divideParam.valRatio = 15/100; net.divideParam.testRatio = 15/100; %% train network [net,tr] = train(net,inputs,targets,inputStates,layerStates); %% test network outputs = net(inputs,inputStates,layerStates); errors = gsubtract(targets,outputs); performance = perform(net,targets,outputs) %% view network view(net) %% plot performance figure, plotperform(tr) %% close loop netc = closeloop(net); netc.name = [net.name ' - Closed Loop']; view(netc) [xc,xic,aic,tc] = preparets(netc,inputSeries,{},targetSeries); yc = netc(xc,xic,aic); perfc = perform(netc,tc,yc) %% remove delay nets = removedelay(net); nets.name = [net.name ' - Predict One Step Ahead']; view(nets) [xs,xis,ais,ts] = preparets(nets,inputSeries,{},targetSeries); ys = nets(xs,xis,ais); earlyPredictPerformance = perform(nets,ts,ys)