Попов В.П. Основы теории цепей (1985) (1092095), страница 73
Текст из файла (страница 73)
Данные о координатах и значениях ненулевых элементов компонентных и топологических матриц представляют в виде совокупности одномерных массивов, называемых с п и с к а м и. Дополнительно объем памяти ЭВМ можно сэкономить за счет тога, что ненулевые элементы топологических матриц могут принимать значения только + 1 или — 1.
В связи с тем чта элементы компонентных и топологических матриц численно равны коэффициентам компонентных (7.3) или (7.4) и топологнческих (1.46), (1.62) уравнений, запись в память ЭВМ элементов этих матриц можно рассматривать как занесение в память ЭВМ коэффициентов основной системы уравнений электрического равновесия цепи, а формирование компонентных и топологических матриц равносильно, следовательно, формированию основной системы уравнений электрического равновесия цепи в матричной форме. При использовании принятого в этой главе топологического описания основная система уравнений электрического равновесия цепи, содержащей р ветвей, включает в себя 2р уравнений.
Количество уравнений электрического равновесия может быть уменьшено за счет исключения нз основной системы уравнений зависимых токов и напряжений с помощью методов узловых напряжений и контурных токов. Методы узловых напряжений и контурных токов Известно, что метод узловых напряжений основан на том, что напряжения всех ветвей произвольной электрической цепи могут быть- выражены через ее узловые напряжения.
В матричной форме зависимость напряжений ветвей от узловых напряжений определяется соотношением, получившим название у з л о в о г о п р е о б р а з о в ани я: и — А'н!О (7.5) Здесь А' — транспоиированная матрица узлов; н, и„ вЂ” матрицы- столбцы мгновенных значений напряжений ветвей и узловых напряжений соответственно ° ФЭФФ пример тпь Убедимся в справедливости аютношения (7.5) на примере цепи, сяеми потирай приведена на рис. 7.8, а. Митрица умов этой цена, соотвеасствуюи!ая ее графу, ипображгннол~у на рис 7.3, б, имеет вид — ! О О О 1 О О О 1 ΠΠ— 1 А -.: 1 Π— 1 — 1 О О 1 О О О О 0 — 1 О Подспитляя транспонированную матрицу О 1 узлов О --1 О О О 1 О О О О 1 ΠΠΠΠ— 1 1 — 1 О О О 1 — ! ΠΠΠ— 1 в выражение (7.Б), получаем — 1 О О О иь иьв О 1 О О иэв О О 1 О О О 1 — 1 О О 1 — 1 ΠΠ— 1 "вь иь игь — ию и„— и,„ и„— и„ вЂ” ит — иьь иь ив э и Аналогичное соотношение может быть получено непосредственно из рас.
смотрения схемы цеии (рис 7.3, а). Для формирования системы узловых уравнений воспользуемся компонентным уравнением пепи в форме у. умножая матрицу узлов на каждое из слагаемых, входящих в выражение (7.4), получаем А! — А!у + Ач(н — Ат(нв. 335 Левая часть этого уравнения представляет собой матричную запись уравнения баланса токов (!.46), поэтому А!у + Аа'и — А а'пу = О. Выражая напряжения ветвей через узловые напряжения (7.6) и выполняя преобразования, получаем систему узловых уравнений цепи в матричной форме А'а А<и(о — — А (Уи — ! ) (7.6) или т(((1)н(о = )(о.
(7.7) Из выражений (7.6), (7.7) видно, что матрица узловых проводимостей исследуемой цепи и матрица-столбец узловых токов У«л) = АУА' и )(о = А (Уиу — !у) (7.8) могут быть получены с помощью простых алгебраических операций над матрицей узлов А и компонентными матрицами 1', и, ), Поскольку компонентные и топологические матрицы цепи содержат значительное количество нулевых элементов, прн формировании матрицы узловых проводимостей У(и1 н матрицы-столбца узловых токов )<о используют специальные алгоритмы, учитывающие разреженность матриц А, У, пв, 1 и исключающие тривиальные операции над нулевыми элементамй.
й(ожно убедиться, что система узловых уравнений цепи для мгновенных значений (7.7) и система узловых уравнений той же цепи для комплексных действующих значений (4.22) имеет одинаковую структуру и могут быть получены одна из другой с помощью таких же преобразований, которые необходимы для взаимного преобразования компонентных уравнений для мгновенных значений и компонентных уравнений для комплексных действующих значений. й т О 0 0 0 0 О о )7 о о о 0 0 — 100 О 1 0 0 010 0 — ! 1 0 )та О )7;т 0 0 0 0 0 0 О 0 0 0 0 О вС 0 0 <ваа) х О 0 Ф-з) ' 0 0 0 0 0 О 0 0 0 О О 0 0 0 1 0 0 — 1 — 1 000 — 1 0 0 1 336 ° ФФФФ пример 7.3.
Псоолазул выраасвнал (7.3), сформируем матрацу узловых проводимостей п матрицу-столбец узловых токов цепи, схема которой приведена на рис. 7.3, а. Матрица узлов А н компонентные матрицы цепи У. нв, 1в были получены в примерах 7.1 и 7.2. Подставляя эти матрицы в <7лз), находим о Π— ! — о — 1 1 е(1 ( — +++к)! ( — + — +еС); 1 1 5/-е — 100 О ! 0 0 -Я;0 О 010 0 — 1 1 0 001 0 0 — 1 — 1 000 — ! 0 0 1 — е! Й„ — !00 0 ! 0 0 010 0 — ! 1 0 ! — !е (0) 0 001 0 0 — 1 — 1 000 — 1 0 0 1 — !7 (0) е!й,— 1, (О) !а(0) и (0) 1 — 1~ (О) 0 0 0 0 0 ! 1, (о) 0 (7 (о) О 0 о о 1 — 1 о о о 0 0 0 0 0 я;~ о О атее 0 0 0 О 0 О 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Я;~ О 0 О 0 (еЦ) еО О 0 0 еСО О О О (я(.~)-1 Нетрудно убедиться, что аналогичные выраасения могут быть получены и с применением рассмотренного в ел.
4 илгоритми, используемого в ручных методах форл!ирования уривненид электрического равновесия. Сформировать уравнения электрического равновесия цепи с применением ЭВМ можно и методом контурных токов. Используя компо. нентное уравнение цепи в форме У и учитывая, что токи всех ветвей цепи могут быть выражены через токи главных ветвей с помощью соотношения 1= В'!Нм называемогоконтурным преобразов а н и е м, нетрудно прийти к системе контурных уравнений цепи в матричной форме Х<!г! гп = еср Здесь с!г — матрица-столбец контурных токов; У<!!> = ВУВ'— матрица контурных сопротивлений; еп — — В (Х!, — и„'! — матрица- столбец контурных э.
д. с., В и В' — матрица главных контуров и трапспонированная матрица главных контуров рассматриваемой цепи. В связи с тем что компонентные матрицы У и У содержат в общем случае операторы дифференцирования з и интегрирования з-', узловые и контурные уравнения цепи для мгновенных значений наряду с производными содержат интегралы от неизвестных функций времени. Метод переменных состояния Наличие интегралов в уравнениях электрического равновесия цепи, составленных методами узловых напряжений и контурных токов, значительно затрудняет решение этих уравнений и в течение длительного времени ограничивало возможности применения данных методов при машинном анализе цепей. Интегралы, входящие в уравнения электрического равновесия, могут быть устранены путем дифференцирования, однако при этом повышаегся порядок соответсгвующих уравнений, что также является нежелательным. Поэтому представляет интерес попытка составить уравнения электрического равновесия таким образом, чтобы они вообще не содержали интегралов.
Интегралы в уравнениях электрического равновесия возникакп только тогда, когда напряжение емкости выражают через ток ис = ис(0)+ — д! гс с1г ! Г. с,! о или ток нндуктивности через напряжение ! р гг =!с !!у!+ — д! ис с!г. о Если в качестве независимых переменных выбрать не контурные токи или узловые напряжения, а напряжения емкостей и токи индуктивностей, то уравнения электрического равновесия цепи не будут ~держать интегралов от неизвестных функций времени. Такие уравнения называются уравнениями состояния цепи, а независимые переменные (токи индуктивностей и напряжения емкотей) — переменными состояния.
Такое название отажает тот факт, что именно токи индуктивностей и напряжения ем„~стей определяют запасы энергии в реактивных элементах и, следовательно, характеризуют энергетическое состояние цепи. Рассмотрим методику формирования уравнений состояния на примере простейшей последовательной )сЕС-цепи (см. рис, 2.20, а), основная система уравнений электрического равновесия которой имеет вид аис ес =-С вЂ”; ю е'и иь=-Е— ш иа+ис+ис=-е; (7.9) 1я = 1ь:=1с=1,' Выбирая в качестве независимых переменных (переменных состояния) напряжение емкости и ток индуктивности, выразим остальные переменные, входящие в эти уравнения, через ие и (х: ""с )71ь + 1- — + ис = е; 1с =- С вЂ” ' Н а1 (7,10) или дХЫ1 =- аХ + дп,.
(7.! 1) Ггь 1 Ге Здесь Х = — вектор переменных состояния; и, =- И ~0 à — Я(Е; — 17Е1 Г 1/Е; 0 1 тор внешних воздействий; а =-. ~ ~ 11С; 0 ~' ~0; 0 ~ — мат- рицы, элементы которых определяются параметрами пассивных элементов цепи. ЗЗЗ Получена система уравнений электрического равновесия, в которых в качестве неизвестных фигурируют напряжение емкости и ток индуктивностн.
Разрешим уравнения (7.10) относительно производных: й,/Ж = — — Жь!Š— иеЯ. + еЯЕ; Лис/й = 1,!С и представим полученную систему уравнений в матричной форме Выражение (7.11) является стандартной формой записи уравнений состояния цепи, не содержащей зависимых источников энергии, Очевидно, что число независимых уравнений, составляемых по методу переменных состояния, будет равно числу независимо включенных реактивных элементов, т. е. порядку сложности цепи. Если исследуемая цепь содержит топологические вырождения, к которым относятся емкостные контуры и индуктивные сечения, то система уравнений электрического равновесия цепи наряду с дифференциальными уравнениями (7 11) будет содержать алгебраические уравнения, составленные на основании второго или первого законов Кирхгафа и отражающие связь между напряжениями емкостей или токами индуктивнастей, входящих в соответствующие контуры или сечения.