bp.m

短期负荷预测通常是指24小时的日负荷预测和168小时的周负荷预测

clear;
 %T = [27,32,20,31,26,23,28,23,22,33,24,30,33,32,25,36,32,35,40,29,22,27,36,31,24,37,23,25,22,30,24,26,29,27,25,27,30,31,25,30,32,30,29,31,28,37,23,25,26,14,21,23,25,26];
 T=[2 5 3 8];
 for i=1:4;
     P(i)=i;
 end
 [r,q]=size(P);
 [s2,q]=size(T);
 s1=5;
 [w1,b1]=rands(s1,r);
 [w2,b2]=rands(s2,s1);
 net.trainParam.epochs =100;
 net.trainParam.goal =0.01;
 Tp=[10 100 0.01 0.01];
 [w1,b1,w2,b2,epochs,errors]=trainbp(w1,b1,'tansig',w2,b2,'purelin',P,T,Tp)

 net = newff([0 1],{w1 w2},{b1 b2},{'tansig' 'purelin'});
 Y2 = sim(net,P);
 hold on;
 plot(P,T,P,Y2,'o')