Заметим, что в рассматриваемом графе нет неправильных раз мещении; поэтому в уравнениях:
— токи в емкостных ветвях выражаются как алгебраические суммы токов резистивных ветвей;
— напряжение резистивных ветвей (и, следовательно, токи ветвей) выражаются через напряжения емкостных ветвей и ис точников питания.
Таким образом, в этом случае могут быть сразу записаны уравнения переменных состояния в нормальной форме (П.11).
При наличии неправильных размещений задача получения уравнений в нормальной форме оказывается более сложной. Дей ствительно, если в дереве оказалась ветвь, соответствующая ре зистивной ветви схемы (как, например, в графе рис. П.Зе), то в правую часть уравнений напряжений, составленных по ЛІ-мат- риде, войдут не только напряжения на емкостях и напряжения источников, но и напряжение на активном сопротивлении; полу чение уравнений переменных состояния в нормальной форме свя зано с дополнительным преобразованием: необходимо выразить напряжение на активном сопротивлении через емкостные.
Аналогично, если ветвь графа, соответствующая емкостной ветви схемы, не оказалась в выбранном дереве, в правых частях системы уравнений для токов окажется емкостный ток и для по лучения нормальной формы уравнений необходимо выразить этот ток через токи резисторов и выполнить соответствующие преобра зования.
На базе описанной выше формализованной методики записи уравнений нетрудно разработать программу автоматического по лучения математической модели схемы па ЦВМ [37]. Однако при наличии неправильных размещений получение моделей по этой ме тодике становится сложной и трудоемкой задачей: в этих слу чаях более предпочтительно получение математической модели схемы в матричной форме [39].
П.5. АЛГОРИТМЫ РАСЧЕТА
Пусть тем или иным способом получена математическая мо дель схемы — система уравнений в нормальной форме, например система (П.11). Для решения этих уравнений на ЦВМ следует вы брать некоторый численный метод (например, Эйлера, Рунге — Кут та, Адамса или их модификации).
Выберем метод Эйлера: приближенное решение дифферен
циального |
уравнения |
— f (z) производится |
по алгоритму |
z n + i = Z n |
- \ - h f ( z n ) 1 где h |
— шаг интегрирования, п — номер шага. |
Запишем алгоритм решения системы (П.11) в виде |
|
«и (я + |
1) = |
и„ («) + |
АФ, [«п (п), иа (я)] |
I |
|
«і2(я + |
1) = |
Иі2.(я) + |
АФ[цп (п), щ2 (п)] |
J |