- •Математическое обеспечение сапр.
- •Лекция №2
- •Устойчивость явного и неявного методов временного анализа.
- •Временной анализ с повышенной точностью.
- •Прогноз начальных условий.
- •Временной анализ нелинейных схем.
- •Линеаризация нелинейных схем. Разложим функцию в ряд Тейлора относительно точкиUk:
- •Определение производных в сложных функциях.
- •Линеаризация системы нелинейных уравнений.
- •Лекция №5 Алгоритмические модели схем во временной области.
- •Анализ добротных схем, минуя переходной процесс.
- •Лекция №6
- •Алгоритм определения стационарного режима автогенератора.
- •Лекция №7 Частотный анализ схем
- •Прямой метод формирования матрицы проводимости.
- •Лекция №8 Узловой анализ схем с комплексными коэффициентами.
- •Частотный анализ нелинейных схем. Метод гармонического баланса.
- •Узловое уравнение с нелинейными коэффициентами.
- •Метод решения нелинейных алгебраических уравнений – метод простых итераций.
- •Формирование нелинейных уравнений с2м-полюсником. Формирование 2м-полюсников.
- •Очевидно, что
Временной анализ нелинейных схем.
Любая схема, содержащая нелинейный элемент, называется нелинейной. Функция описывающая математическую модель схем, в общем виде (для линейных схем) имеет вид:
,
где
– переменные составляющие,
– их производные.
Дискретная модель в общем виде имеет вид:
,
то есть в конкретный момент времени,
производная заменяется на функцию
.
Для нелинейных схем:
Ток
в каком-либо узле есть нелинейная функция
,
гдеi
и U
– ток и напряжение в установившемся
режиме.
Линеаризация нелинейных схем. Разложим функцию в ряд Тейлора относительно точкиUk:
,
где
Uk – значение в известной точке, R – члены разложения высших порядков.
Отбросим члены разложения высших порядков получим линейную модель
![]()
Поскольку мы пренебрегли членами высших порядков разложения, то решение не будет точным.
Суть
метода заключается в том, что при
![]()
.
Докажем это предположение:
Пусть
имеется нелинейное уравнение
,
представим его в виде ряда Тейлора:
![]()
выражаем
:
.
Воспользуемся графическим методом решения:

Тогда
,
очевидно,
.
При
приближаемся кx.
Этот метод носит название – метод
Ньютона линеаризации нелинейных
уравнений.
Замечание: метод сходится всегда, если начальные условия заданы близко к исходной точке.
Определение производных в сложных функциях.
Сложная функция, пример:
.
Сложная функция в общем виде:
.
Производная этой функции имеет вид:
, (1)
Найдем величину
:
,
– коэффициентk-
коэффициент по амплитуде,i– коэффициент по времени. (Так как решаем
относительно текущего момента времени,
то
).
В простейшем случае
.
То есть (1) имеет вид:
.
Пример:
,
вычислим производную:
.
Линеаризация системы нелинейных уравнений.
Пусть рассматривается функция вида:
.
Многомерное разложение в ряд Тейлора:
Математическая
модель схемы описывается системой
дифференциальных уравнений:
(2)
Обозначим:
,
а
![]()

То есть математическая модель схемы имеет вид:
(3).
Разложим (3) в ряд Тейлора:
,
–множество
частных производных появившихся при
многомерном разложении в ряд Тейлора.
Более
подробно
:
![]()

Эта матрица производных называется якобиан.
Лекция №5 Алгоритмические модели схем во временной области.
Математическая модель схемы имеет вид:
,
то есть это система дифференциальных
нелинейных уравнений.
Для описания данной модели в компьютере необходимо произвести алгеброизацию (замену производных на линейную функцию), а так же линеаризацию (разложение в ряд, с целью линеаризирование нелинейных функций).
В результате линеаризации (см. прошлую лекцию) получили:
![]()
Это линеаризированная модель схемы. Произведем алгброическое ппреобразование получим решение.
,
Этот метод носит название Ньютона - Рафсона.
Недостаток
метода
алгоритмизации Ньютона - Рафсона –
присутствие в алгоритме операции
обращения матрицы Якоби
.
Представим (1) как:
.
Это метод Ньютона Бройдена - представляет собой систему линейных алгебраических уравнений с правой частью (и в отличие от предыдущего метода никаких обратных матриц!).
Таким образом, этот метод дает решение нелинейной системы дифференциальных уравнений в определенный момент времени, то есть на самом деле формулы имеет вид:

При реализации метода задаются следующие величины:
,
где h
– шаг дискретизации по времени,
-
начальные условия,M
– порядок полинома.Вычисляется:
,
то есть находятся коэффициенты системы
и вычисляются ее правая часть.Решаем систему (2), находим
.Производим контроль результата. (То есть, сверяем точность):
,
где
– число точек на временном интервале,
– модуль, либо среднеквадратичное
значение разности.
Если
получилось, что
допустимое,
то необходимо уменьшить шаг по времени.
(Шаг дискрета по времени, например
).
Если
<
допустимое,
то
принять за
,
и повторить все действия, начиная с
первого.
