symbolic dynamic models for highly varying power system loads

91
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
% New value of current, newI is the last symbol of the new % word,
selected from the mini dictionary.
sinus(p)=M*sin((2*pi*60*0.0001*p)+ph);
% Sinus is a sinusoidal source formed by taking the
% fundamental component of the historical data set.
junk(p)=((newI(p)-B)/A)-sinus(p); % amount of
% unwanted components in the new current value.
good(p)=a(100+p)-junk(p); % injecting the
% unwanted components in phase opposition, into the actual % future
current.
x=[x,I(100+p)];
x11=[x11;I(100+p)];
newI(p);
clear ss k n nx minx minfr len min wtfr sm sm1 newfr cfr sz
newwd d1 f2 kount1 frac_occ
p=p+1;
end % end of while
filt1=good'/1000;