uybar(mu,:) = filter(1,[1, -Az(mu,mu)],T*uexcite(mu,:));