symbolic dynamic models for highly varying power system loads

86
x=aa1(:,1:100);
I_hist=x'-B; % Historical data set, without any shift, for % the
purpose of RMS value comparison.
I=aa1';
for t=1:100
II(t)=1*sin(2*pi*60*0.0001*t);
end
% II is a synthetic current source with unity amplitude and % points
separated by the same sampling frequency as in
% historical data set.
IFFT=fft(II); % FFT of the above current signal.
fac=abs(IFFT(2))% Finding a factor to obtain the true FFT.
ph1=angle(IFFT(2)) % Phase angle of the fundamental.
hisFFT=fft(a(1:100)); % FFT of the historical data set.
M=fac*abs(hisFFT(2)) % True magnitude of the fundamental
% component of the historical data set.
ph=angle(hisFFT(2)) % Phase angle of the fundamental
% component of the historical data set.
x11=x'; % Historical data set.
nn=size(x11);
n=nn(1,1);
nw=5;
%calculate dictionary of all words. The dictionary is d
f=0;
%w=word length
for w=1:nw
for s=1:n-w+1
c=zeros(w,1);
c(1:w)=x11(s:s+w-1);