Для создания новых устройств при помощи их математических моделей важное значение имеет возможность представления количественных характеристик, определяющих исследуемые процессы нового устройства, в виде аналитических зависимостей и численных данных. Важнейшими для практики теоретическими проблемами оказались разработка методов учета особенностей протекания электромагнитных процессов в зависимости от электромагнитных свойств сред и их конфигурации, формирование принципов и методов анализа и синтеза математических моделей электротехнических устройств, а также диагностирование процессов в них и управление ими. История создания общих подходов к решению этих проблем — важнейшая часть истории ТЭ. Поскольку ТЭ является мостом между фундаментальной и прикладной науками, она вынуждена была выбирать и разрабатывать такие математические методы, которые позволяют производить аналитические исследования и численные расчеты с необходимой для практики точностью. Возможность получения аналитических решений, как правило, определяется сложностью математической модели. Точность же численных расчетов в рамках данной математической модели определяется возможностями аппаратных средств, используемых для выполнения численных расчетов.
По этим причинам методы, предлагаемые в ТЭ, развивались с развитием новых разделов математики и средств вычислительной техники. Появление ЭВМ и их внедрение в практику, начиная с 1950 г., оказало решающее воздействие на ТЭ.
Практика использования ЭВМ для расчета электрических цепей привела к существенному изменению направления теоретических разработок в ТЭ. Ограниченные возможности существоваших до появления ЭВМ средств численных расчетов стимулировали развитие методов, позволяющих снижать количество уравнений, вынуждали развивать теорию подобия и создавать физические и математические аналоги электрических цепей.
Так, например, в 50-е годы широкое распространение получили расчетные столы для исследования сложных цепей и систем, физические модели электроэнергетических систем (М.П. Костенко, Л.Р. Нейман, В.А. Веников, и др.), аналоговые и цифроаналоговые модели (Г.Е. Пухов, Б.Я. Коган, Н.Е. Кобринский, Г. Ольсон и др.) и цифровые дифференциальные анализаторы (А.В. Каляев) для исследования переходных процессов. Ограниченное математическое обеспечение первого и второго поколения ЭВМ вынудило записывать уравнения электрических цепей в виде системы дифференциальных уравнений, разрешенных относительно первой производной, т.е. в форме уравнений Коши. Поскольку динамические свойства цепей определяются появлением ЭДС при изменении магнитного потока и токов смещения вследствие изменения потока электрического смещения, то уравнения Коши для цепей естественно записать для переменных потокосцеплений Ψ и электрических зарядов Q. Поскольку именно таков был подход к выбору значимых переменных в термодинамике, то и в ТЭ эти переменные были названы переменными состояния. В 1957 г. Т.Р. Башков впервые записал уравнения электрических цепей относительно переменных состояния. Развитие матрично-топологического метода в большой мере определялось необходимостью автоматизировать ввод данных о топологии цепи и формировать на этой основе уравнения состояния (Г. Крон, Р. Рорер, Ф.Х. Брэнин, С.Д. Фенвес, Э. Ку, Д.Р. Рос, И.П. Норенков, В.Н. Ильин и др.).
Численные расчеты электрических цепей на ЭВМ показали, что использование стандартного математического обеспечения недостаточно для эффективного использования вычислительной техники. Специфичными для цепей были проблемы учета разреженности узловых и контурных матриц параметров, большого разброса собственных чисел этих матриц, жесткости системы дифференциальных уравнений; проблемы достижения решения для периодических установившихся режимов, выбора переменных для обеспечения однозначности решения уравнений нелинейных цепей и др. Вклад ученых (Н. Сато, В.Ф. Тинней, Г.Д. Хэчел, Р.К. Брейтон, Т.Д. Эприл, Т.Н. Трик, К.В. Гир, К.Г. Бройден, Ю.В. Ракитский, С.М. Устинов, К.С. Демирчян, П. А. Бутырин и др.), работающих в области ТЭ, при решении этих проблем привел к появлению в прикладной математике новых разделов. Так, например, метод обращения слабозаполненных матриц, предложенный Н. Сато и В.Ф. Тиннеем инициировал создание таких матриц в прикладной математике. Проблемы выбора шага численного интегрирования в случае жесткости системы дифференциальных уравнений были успешно решены Ю.В. Ракитским и его учениками при помощи разработки нового метода численного интегрирования, названного системным. По этой методике можно на каждом новом шаге интегрирования удвоить его без снижения точности интегрирования и нарушения устойчивости численного процесса интегрирования. Метод, названный Ю.В. Ракитским квазистационарностью производной, позволил существенно упростить решение жестких систем уравнений, снизить их порядок. Использование метода нахождения решения для установившегося состояния при помощи операторного метода (К.С. Демирчян, П.А. Бутырин) позволило организовать такой процесс численного решения уравнений состояний, в котором сочетаются все достоинства системных методов с непосредственным расчетом установившегося режима.