bpfault.asv

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

tic, % 开始计时
[P,T,R,S1,S2,S,meanp,stdp,meant,stdt]=nninit; % BP网络初始化
disp_fqre=1000; max_epoch=20000;err_goal=0.00005;lr=0.001;
 TP=[disp_fqre  max_epoch err_goal lr];
[W1,B1]=rands(S1,R);
 [W2,B2]=rands(S2,S1);
[W1,B1,W2,B2,epochS,errors]=trainbp(W1,B1,'tansig',W2,B2,'purelin',P,T,TP);
% 仿真结果
TT=simuff(P,W1,B1,'tansig',W2,B2,'purelin')
E=sum((T-TT).^2)./10;
E=sqrt(E)
e1=(T-TT)./T
c=poststd(TT,meant,stdt)
toc  % 结束计时 
%前面的网络训练过程可以是脚本文件运行
%后面的待预测的数据必须每次从命令窗口输入才是每次得到同一的预测结果 只有这样才是每次利用前面已经训练好的网络
%检验数据： 
%##########################################
 % xk1=input('     please input .....xk==')
  %pnew=xk1';
  %pnew=trastd(pnew,meanp,stdp);%将数据通过自标准化再用来输入网络
  %an=simuff(pnew,W1,B1,'tansig',W2,B2,'purelin');%网络模拟出来的值
  %anew=poststd(an,meant,stdt)%模拟出的值再线性返回样本预测值
  %tk1=input('     please input .....tk==')
  %e2=(tk1-anew)./tk1
  %xk2=input('     please input .....xk==')
  %data
  data2
  pnew=xk';
  pnew=trastd(pnew,meanp,stdp);%将数据通过自标准化再用来输入网络
  an=simuff(pnew,W1,B1,'tansig',W2,B2,'purelin');%网络模拟出来的值
  anew=poststd(an,meant,stdt)%模拟出的值再线性返回样本预测值
  %tk2=input('     please input .....tk==')
  tk=tk'
  e2=(tk-anew)./tk