"проинтегрируем по одной степени свободы"
По степени свободы хорошо работает в случаях чистых возбуждений по интересующей степени свободы:

В случае же возбуждения по второй оно даёт неправильный результат, очевидный невооружённым глазом:

Состояние, очевидно, симметрично вдоль

, но на одномерном отображении оно выглядит как дважды возбуждённое состояние.
Чтобы сгладить осцилляции по второй оси, я попробовал интегрировать, используя автокорреляционную функцию вдоль второй координаты в качестве веса при усреднении. И в этом случае ситуация исправляется для второго случая:

Но портится для некоторых вариантов первого:

Собственно, вот функция:

а вот её автокорреляция:

Может есть какой-то способ с этими случаями побороться одним средством?