symbolic dynamic models for highly varying power system loads

84
end
end
newwd;
wdlen;
mm=wdlen+1;
emp=any(newwd(:,2:nw));
if(emp==0)
newI(p)=newwd(1);
else
newI(p)=newwd(mm);
end %end of if-else
x=[x,newI(p)];
% New value of current, newI is the last symbol of the new % word,
selected from the mini dictionary.
clear ss k temp nx minx minfr len wdlen min wtfr sm newfr cfr sz newwd
p=p+1;
end % end of while
Ipred=newI'-B % predicted current data set
Iact=I(3001:5000)-B; % expected current
RMSpred=sqrt(mean(Ipred.^2)) %RMS of predicted data
Avgpred=mean(Ipred) %Average value of predicted data
RMSact=sqrt(mean(Iact.^2))%RMS value of actual data
Avgact=mean(Iact) % Average value of actual data
Error_act=((RMSpred-RMSact)/RMSact)*100
RMShist=sqrt(mean(I_hist.^2)) %RMS of historical data
Error_hist=((RMSpred-RMShist)/RMShist)*100
Avghist=mean(I_hist) % Average value of historical data