symbolic dynamic models for highly varying power system loads

79
% A program to forecast a signal for load modeling application
***********************************************************
clear
% Reading data from the file.
Points=5000;
fid=fopen('h:\Ib.dat','rt');
if(fid==(-1))
disp('Error in opening the given input data file');
end
n=1;
while n~=(Points+1)
a(n)=fscanf(fid,'%f',1);
n=n+1;
end
% Calculation of the scaling function
N=10;
minx=min(a);
maxx=max(a);
mat1=[minx 1;maxx 1];
mat2=[1.01;N+0.99];
mat3=inv(mat1)*mat2;
A=mat3(1);
B=mat3(2);
aa1=fix(A*a+B);
x=aa1(:,1:3000); % Historical data set
I_hist=x'-B; % Historical data set, without any shift, for % the
purpose of RMS value comparison.
I=aa1';