example4_8.m

matlab6.5辅助小波分析与应用例程 希望对大家能有所帮助

%导入原始测量数据
load choles_all;
%对原始数据进行规范化处理
[pn,meanp,stdp,tn,meant,stdt]=prestd(p,t);
[ptrans,transMat]=prepca(pn,0.001);
[R,Q]=size(ptrans) 
%将原始数据分成几个部分作为不同用途
iitst=2:4:Q;
iival=4:4:Q;
iitr=[1:4:Q 3:4:Q];
vv.P=ptrans(:,iival); 
vv.T=tn(:,iival);
vt.P=ptrans(:,iitst); 
vt.T=tn(:,iitst);
ptr=ptrans(:,iitr); 
ttr=tn(:,iitr);
%建立网络
net=newff(minmax(ptr),[5 3],{'tansig' 'purelin'},'trainlm');
%训练网络
net.trainParam.show=5;
[net,tr]=train(net,ptr,ttr,[],[],vv,vt);
%绘出训练过程中各误差的变化曲线
plot(tr.epoch,tr.perf,'r',tr.epoch,tr.vperf,':g',tr.epoch,tr.tperf,'-.b');
legend('训练','确证','测试',-1);
ylabel('平方误差');
xlabel('时间');
pause;
%得到各个级别的结果的线性回归结果，并绘出曲线
an=sim(net,ptrans);
a=poststd(an,meant,stdt);
for i=1:3
figure(i)
[m(i),b(i),r(i)] = postreg(a(i,:),t(i,:));
end