function  Y=simulazARX(AA,BB,U,x0)
%x0 = primi n campioni di y dove n è l'ordine del sistema 
    Y=zeros(size(x0,1),size(U,2));
    Y(:,1:size(x0,2))=x0;
    
    for i=size(x0,2)+1:size(U,2)
        for j=2:size(AA,3)
            Y(:,i)=Y(:,i)+AA(:,:,j)*Y(:,i-j+1)+BB(:,:,j)*U(:,i-j+1);              
        end
    end
     %Y(:,i)=-AA(:,:,2)*Y(:,i-1)-AA(:,:,3)*Y(:,i-2)+BB(:,:,1)*U(:,i)+BB(:,:,2)*U(:,i-1)+BB(:,:,3)*U(:,i-2);  
end

% 
