symbolic dynamic models for highly varying power system loads

90
if(isempty(minx))
for jj=1:f
if(d(jj,2:nw)==0)
minx(kk,:)=d(jj,:);
minfr(kk,:)=frac_occ(jj,:);
len(kk,:)=1;
end
kk=kk+1;
end
end
min=[minx,minfr,len]; % Mini dictionary of all the matching %
symbols/sequences
wtfr=minfr.*(len.^2);% weighted fractional occurrence
sm=sum(wtfr(:,1));
newfr=wtfr/sm; % modifying wtfr so that it add to 1.
cfr=cumsum(newfr); % cumulative fractional occurrence
sz=size(minx);
f2=sz(1);
r=rand;
if (r=cfr(i))&(r