Matrix computation of subresultant polynomial remainder sequences