Работа посвящена моделированию характеристик излучения докритических рентгеновских пульсаров с сильным магнитным полем. При соударении аккрецирующего вещества с поверхностью нейтронной звезды порождается излучение в рентгеновском диапазоне, причём давление этого излучения настолько велико, что оно влияет на динамику аккрецирующей плазмы. Поэтому задачи о нахождении характеристик излучения рентгеновского пульсара и гидродинамики падающего вещества должны решаться самосогласованным образом. Расчёты осложняются тем, что сильное магнитное поле искажает амплитуды элементарных процессов взаимодействия излучения и вещества. В данной работе произведено самосогласованное радиационно-гидродинамическое моделирование потока аккрецирующей плазмы на полюса нейтронной звезды. Рассмотрен случай докритических светимостей, учтено комптоновское рассеяние в сильном магнитном поле с сечением, содержащим основной циклотронный резонанс. В таких предположениях получены спектры и поляризация излучения аккреционной колонки.