Ну, векторного пространства они не образуют - ну хотя бы потому, что среди них нуля нету, если

ненулевая. А описать алгоритмически можно - я по нему и действовал:
Стартуя с двух подобных матриц

и

составляем ОСЛУ

, находим ФСР и общее решение - получаем подпространство

в пространстве матриц. Однако не все нам годятся, а только обратимые. Каждому невырожденному решению

этой ОСЛУ сопоставляем матрицу

и получаем полное описание всех этих

. Если требуется абы какое

, то берём абы какую невырожденную матрицу, удовлетворяющую ОСЛУ.