windspeed=xlsread('windspeeddata');
past_data=windspeed(:,3);
y1=past_data(1:50);
for i=1:25
b(:,i)=y1((i-1)+1:i*1+25);
end
for i=1:25;
sys=armax(b(:,i),[1 2]);
end
sys1=armax(b(:,1),[1 2]);
sys2=armax(b(:,2),[1 2]);
sys3=armax(b(:,3),[1 2]);
sys4=armax(b(:,4),[1 2]);
sys5=armax(b(:,5),[1 2]);
sys6=armax(b(:,6),[1 2]);
.....
for i=1:25;
eval(['sys' num2str(i) '=armax(b(:,i),[1 2]);'])
end