guiyihua.m

这是对一维时间序列进行归一化计算的matlab程序

r=[13.43 13.7  20 26.16	32.21 35.46	39.41 42.99	46.7;
   17.03 11.74 16.28 21.16 23.91 25.7 31.43 37.29 42.94;
   10.39 12.88  21.7 25.27	31.66 37.09	41.86 47.57	54.99;
   234.19 238.64 243.18	248	252 255.69 	258.3 269.26 296.01;
   35392 44422 50599 53744 60935 69051 75870 80851 93356]; 
r=r';
p=r(1:5,:)
t=r(6,:);
pmax=5*p(:,5);
pmin=0.5*p(:,1);
tmax=5*t(5);
tmin=0.5*t(1);
for i=1:5   
    p(:,i)=(p(:,i)-pmin)./(pmax-pmin)  ;
    T(1,i)=(t(1,i)-tmin)./(tmax-tmin) ;
end
    eg=0.02 ;
    sc=1;
    net=newrb(p,T,eg,sc) ;
    c=sim(net,p);
    P=[144.63;46.71;42.94;54.99;296.01] ;
    P=P./pmax;
    a=sim(net,P) ;
    b=tmax.*a ;
    P1=[167.77;50.17;52.86;63.75;296.98];
    P1=(P1-pmin)./(pmax-pmin) ;
    a1=sim(net,P1) ;
    b1=(tmax-tmin).*a1+tmin;
    P2=[194.61;53.9;65.09;73.95;309.07];
    P2=(P2-pmin)./(pmax-pmin) ;
    a2=sim(net,P2);
    b2=(tmax-tmin).*a2+tmin ;
    P3=[225.75;57.91;80.15;85.78;322.14] ;
    P3=(P3-pmin)./(pmax-pmin);
    a3=sim(net,P3) ;
    b3=(tmax-tmin).*a3+tmin ;
    P4=[261.87;62.21;98.68;99.51;336.2];
    P4=(P4-pmin)./(pmax-pmin);
    a4=sim(net,P4);
    b4=(tmax-tmin).*a4+tmin;
    P5=[303.77;66.83;121.51;115.43;351.25];
    P5=(P5-pmin)./(pmax-pmin) ;
    a5=sim(net,P5);
    b5=(tmax-tmin).*a5+tmin; 
    P6=[352.38;71.8;149.61;133.88;367.28];
    P6=(P6-pmin)./(pmax-pmin) ;
    a6=sim(net,P6) ;
    b6=(tmax-tmin).*a6+tmin ;
    d=c-T;
    Wt=[b,b1,b2,b3,b4,b5,b6];
    P=[P1,P2,P3,P4,P5,P6]