应用了归一化的预测

「爱情、让人受尽委屈。」 2022-05-16 10:29 125阅读 0赞

应用了归一化的预测，在归一化的过程中使用了premnmx和postmnmx，并在最后给出了这两个函数的应用情况。

%应用了归一化的预测 clc; clear; originalData=rands(20,14); inputSampledata=originalData'; outputData=rands(20,1); outputSampledata=outputData'; % creating neural network and setting trainging parameters gwwnet=newff(minmax(inputSampledata),\[4,1\],\{'tansig','purelin'\},'traingdm'); gwwnet.trainParam.show = 50; gwwnet.trainParam.lr = 0.05; gwwnet.trainParam.epochs = 5000; gwwnet.trainParam.goal = 1e-3; \[input,mininput,maxinput,output,minoutput,maxoutput\] = premnmx(inputSampledata,outputSampledata); % Preprocess data so that minimum is -1 and maximum is 1 % % Examples % % Here is the code to normalize a given data set so that the inputs % and targets will fall in the range \[-1,1\]. % p = \[-10 -7.5 -5 -2.5 0 2.5 5 7.5 10\]; % t = \[0 7.07 -10 -7.07 0 7.07 10 7.07 0\]; % \[pn,minp,maxp,tn,mint,maxt\] = premnmx(p,t); % % If you just want to normalize the input, \[pn,minp,maxp\] = premnmx(p); % % Algorithm % pn = 2\*(p-minp)/(maxp-minp) - 1; \[gwwnet,tr\]=train(gwwnet,input,output); y=sim(gwwnet,input); %data offset (converting the network output data to it original unit) nnoutput = postmnmx(y,minoutput,maxoutput); % % Postprocess data that has been preprocessed by premnmx % Syntax % % \[P,T\] = postmnmx(PN,minp,maxp,TN,mint,maxt) % \[p\] = postmnmx(PN,minp,maxp) % Description % % postmnmx postprocesses the network training set % that was preprocessed by premnmx. It converts the data % back into unnormalized units. % postmnmx takes these inputs, % PN -- R x Q matrix of normalized input vectors % minp -- R x 1 vector containing minimums for each P % maxp -- R x 1 vector containing maximums for each P % TN -- S x Q matrix of normalized target vectors % mint -- S x 1 vector containing minimums for each T % maxt -- S x 1 vector containing maximums for each T % and returns, % P -- R x Q matrix of input (column) vectors % T -- R x Q matrix of target vectors % Examples % % In this example we normalize a set of training data with premnmx, % create and train a network using the normalized data, simulate the network, % unnormalize the output of the network using postmnmx, % and perform a linear regression between % the network outputs (unnormalized) and the targets % to check the quality of the network training. % p = \[-0.92 0.73 -0.47 0.74 0.29; -0.08 0.86 -0.67 -0.52 0.93\]; % t = \[-0.08 3.4 -0.82 0.69 3.1\]; % \[pn,minp,maxp,tn,mint,maxt\] = premnmx(p,t); % net = newff(minmax(pn),\[5 1\],\{'tansig' 'purelin'\},'trainlm'); % net = train(net,pn,tn); % an = sim(net,pn); % \[a\] = postmnmx(an,mint,maxt); % \[m,b,r\] = postreg(a,t); % % Algorithm % p = 0.5(pn+1)\*(maxp-minp) + minp; %plot time=1:1:20; plot(time,outputSampledata,'-v',time,nnoutput,'o'); legend('actual output','NN output'); xlabel('time');ylabel('Learning fitting curve'); %scenario1 forecasting process column=10; for i=1:column; SceInput(1,i)=inputSampledata(1,20)\*(1.0464^i); SceInput(2,i)=inputSampledata(2,20)\*(1.0631^i); SceInput(3,i)=inputSampledata(3,20)\*(1.0872^i); SceInput(4,i)=inputSampledata(4,20)\*(1.2044^i); SceInput(5,i)=inputSampledata(5,20)\*(1.2326^i); SceInput(6,i)=inputSampledata(6,20)\*(1.0605^i); SceInput(7,i)=2\*(1.01^i); SceInput(8,i)=42\*(1.02^i); SceInput(9,i)=inputSampledata(9,20)\*(1.1426^i); SceInput(10,i)=inputSampledata(10,20)\*(1.017^i); SceInput(11,i)=inputSampledata(11,20)\*(1.0205^i); SceInput(12,i)=inputSampledata(12,20)\*(1.1336^i); SceInput(13,i)=inputSampledata(13,20)\*(1.1599^i); SceInput(14,i)=inputSampledata(14,20)\*(1.1783^i); end for j=1:20; for i=1:14; recalldata(i,j)=inputSampledata(i,j); end end for j=21:30; for i=1:14; recalldata(i,j)=SceInput(i,j-20) end end \[alterinput,mininput,maxinput\] = premnmx(recalldata); %training fvalue=sim(gwwnet,alterinput); %data offset (converting the network output data to it original unit) forecastvalue = postmnmx(fvalue,minoutput,maxoutput); %plot % waitforbuttonpress; % % Wait for key or mouse button press % Syntax % k = waitforbuttonpress % % Description % % k = waitforbuttonpress blocks the caller's execution stream % until the function detects that the user has pressed a mouse button % or a key while the figure window is active. The function returns 0 % if it detects a mouse button press 1 if it detects a key press % % Example % % These statements display text in the Command Window % when the user either clicks a mouse button or % types a key in the figure window: % w = waitforbuttonpress; % if w == 0 % disp('Button press') % else % disp('Key press') % end % clf; figure; time1=1:1:30; time2=1:1:20; plot(time1,forecastvalue,'r-',time2,outputSampledata,'-'); legend('预测曲线','实际曲线'); title('曲线'); xlabel('时间');ylabel('量'); %%%%%%%premnmx及postmnmx使用说明 % clc; % clear; % p = \[1 2 3 4; % 5 6 8 9 \]; % \[p1 pmin pmax\]=premnmx(p); % p2=postmnmx(p1,pmin,pmax); % p1 % pmin % pmax % p2 % % p1 = % % -1.0000 -0.3333 0.3333 1.0000 % -1.0000 -0.5000 0.5000 1.0000 % % % pmin = % % 1 % 5 % % % pmax = % % 4 % 9 % % % p2 = % % 1 2 3 4 % 5 6 8 9