Разработана 1.5-мерная программа численного моделирования нелинейных нестационарных процессов в приборах клистронного типа, основанная на нестационарной дискретной теории возбуждения резонаторов Л.А. Вайнштейна и методе «частиц в ячейке» для моделирования динамики электронного пучка. Представлены результаты численного моделирования основных режимов колебаний четырехрезонаторного клистронного генератора с внешней запаздывающей обратной связью. Полученные результаты подтверждают качественную картину поведения генератора, изученную ранее для упрощенных теоретических моделей в виде дифференциальных уравнений с запаздывающим аргументом, а также позволяют определить значения выходных параметров (мощность, КПД, ширина спектра и т.д.), важных с практической точки зрения.