Обеспечение точности и устойчивости компьютерного моделирования радиоэлектронных устройств является важной задачей при их проектировании. Наибольшие трудности при компьютерном моделировании радиоэлектронных устройств возникают в случае анализа автоколебательных цепей, поскольку математические модели таких цепей могут быть жесткими и осциллирующими одновременно. Целью данной работы является обоснование эффективного численного метода решения обыкновенных дифференциальных уравнений, имеющего более высокую точность анализа во временной области различных видов автоколебательных цепей по сравнению с известными методами. Предлагаемый метод является гибридным методом и базируется на известных методах Гира и трапеций, использующихся в пакетах схемотехнического моделирования радиоэлектронных устройств. Для оценки точности предлагаемого метода и известных методов использовалась обобщенная модель автоколебательной цепи, для которой было определено аналитическое решение в установившемся режиме работы. Точность численного решения определялась на основе анализа погрешностей оценки основных параметров колебательного процесса – амплитуды и частоты колебаний. Сравнительный анализ погрешностей оценки амплитуды и частоты колебаний автогенераторов в различных режимах работы показывает высокую эффективность предлагаемого гибридного метода как для анализа генераторов гармонических колебаний, так и для анализа генераторов релаксационных колебаний. Дальнейшее повышение точности гибридного метода возможно за счет применения неявных методов Рунге-Кутты высокого порядка точности, а именно подклассов Радо IIA и Лобатто IIIA, которые обладают L- и P-устойчивостью соответственно. Следует отметить, что с ростом порядка точности неявных методов Рунге-Кутты возрастает вычислительная трудоемкость данных методов, но для подклассов Радо IIA и Лобатто IIIA рост вычислительной трудоемкости будет минимальным.
