%數(shù)據(jù)文件為choles_all安裝文件C:\MATLAB7\toolbox\nnet\nndemos
load choles_all
%prestd數(shù)據(jù)預(yù)處理
[pn,meanp,stdp,tn,ment,stdt]=prestd(p,t);
%prepca主要成分分析
[ptrans,transMat]=prepca(pn,0.001);
%驗證矩陣大小
[R,Q]=size(ptrans)
%數(shù)據(jù)分組，訓(xùn)練，驗證，測試（等間隔抽取）
iitst=2:4:Q;
iival=4:4:Q;
iitr=[1:4:Q 3:4:Q];
val.P=ptrans(:,iival);val.T=tn(:,iival);
test.P=ptrans(:,iitst);test.T=tn(:,iitst);
ptr=ptrans(:,iitr);ttr=tn(:,iitr);
%網(wǎng)絡(luò)建立
net=newff(minmax(ptr),[5,3],{'tansig' 'purelin'},'trainlm');
%訓(xùn)練
[net,tr]=train(net,ptr,ttr,[],[],val,test);
%顯示訓(xùn)練誤差，驗證誤差，測試誤差
plot(tr.epoch,tr.perf,'-',tr.epoch,tr.vperf,':',tr.epoch,tr.tperf,'- .')
legend('Training','Validation','Test',-1);
ylabel('平方差');
xlabel('時間');
%sim仿真，poststd反規(guī)范化轉(zhuǎn)換
an=sim(net,ptrans);
a=poststd(an,ment,stdt);
%hdl，ldl，vhdl線形回歸
for i=1:3
figure(i)
[m(i),b(i),r(i)]=postreg(a(i,:),t(i,:));
end
%改進(jìn)20個神經(jīng)元
net=newff(minmax(ptr),[20,3],{'tansig' 'purelin'},'trainlm');
[net,tr]=train(net,ptr,ttr,[],[],val,test);
plot(tr.epoch,tr.perf,'-',...
tr.epoch,tr.vperf,':',tr.epoch,tr.tperf,'- .')
legend('Training','Validation','Test',-1);
ylabel('平方差');
xlabel('時間');
%線形回歸
an=sim(net,ptrans);
a=poststd(an,ment,stdt);
for i=1:3
figure(i)
[m(i),b(i),r(i)]=postreg(a(i,:),t(i,:));
end
[m(1),b(1),r(1)]=postreg(a(1,:),t(1,:));