Новый байесовский метод точнее восстанавливает гамма-спектры
Гамма-астрономия сталкивается с фундаментальной проблемой: наблюдаемые спектры гамма-излучения искажаются из-за конечного энергетического разрешения детекторов. Восстановление истинных спектров — некорректная обратная задача, где прямое обращение усиливает статистические шумы. Существующие методы часто не дают надежных оценок погрешностей.
Ученые из международной группы предложили новое решение — эмпирико-байесовский иерархический метод развертки. Он сохраняет пуассоновскую структуру счета, гарантирует неотрицательность спектра и учитывает фон через совместную likelihood для режимов ON и OFF.
Метод использует априорное распределение для истинного спектра, центрированное на автоматически выбранном опорном спектре Ричардсона-Люси. Адаптивная ширина априора остается широкой в слабо ограниченных областях, не навязывая жестких предположений.
Апостериорный вывод выполняется с помощью семплера No-U-Turn, что позволяет получать одновременные доверительные интервалы для разрешенного спектра. Это дает ученым инструмент для количественной оценки неопределенности — то, чего часто не хватало в классических методах.
Авторы сравнили свой метод с недавним частотным регуляризованным методом максимального правдоподобия. В тестовых сценариях как с высокой, так и с низкой статистикой результаты оказались очень близкими, что подтверждает надежность байесовского подхода.
Новый метод представляет гибкую и расширяемую основу для развертки в гамма-астрономии. В будущем он может быть адаптирован для других областей, где необходимо восстанавливать распределения по данным с искажениями и ограниченной статистикой.
