Численное решение дискретного BHH-уравнения в нормальном случае

Известно, что решение полулинейного матричного уравнения $X-A\overline XB=C$ можно свести к решению классического уравнения Стейна. Нормальный случай означает, что коэффициенты левой части полученного уравнения суть нормальные матрицы. Предлагается способ решения исходного полулинейного уравнения в нормальном случае, позволяющий для уравнений порядка $n=3000$ почти вдвое сократить время вычислений по сравнению с библиотечной функцией dlyap, решающей уравнения Стейна в системе Matlab.