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