Может быть, имеются какие-то более простые рецепты построения именно унитарных матриц?
Да, есть.
Унитарные матрицы - это далеко идущее обобщение понятия
матрицы поворота. Например, вот это вот:

- одновременно и матрица поворота, и ортогональная (действительная) матрица, и унитарная (комплексная) матрица. Ещё унитарными матрицами являются такие диагональные, что модули диагональных элементов равны 1. И ещё некоторые:

В случае матриц большего размера, так просто их описать не получится. Формулы становятся очень громоздкими. Примеры можно найти в Википедии в статьях про углы Эйлера, или про матрицу ККМ. Но есть общее правило: унитарная матрица - это экспонента от некоторой анти-эрмитовой. То есть, если у вас есть некая

то можно получить унитарную матрицу

А экспоненту надо вычислять в базисе с.в. исходной матрицы

: там она будет диагональной с экспонентами от диагональных элементов

(Ортогональные - в лучшем случае в некотором базисе блочно-диагональные, с блоками вида

и

)