Файл: Голембо, З. Б. Алгоритмизация и программирование электротехнических задач на электронных цифровых вычислительных машинах учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 94

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Матрица проводимостей ветвей Y также квадратная матрица порядка р. Элементы, расположенные по ее диагонали, равны проводимостям ветвей. Из принципа взаимности следует, что Следовательно, матрица проводимостей будет симмет-

рична

YU

Y i 2

•••

У\р

У21

У22

•••

 

Y2p

Y \ p

Y2p

 

Y

"V

 

 

 

 

б. Основные уравнения цепи

Для произвольно взятой k-й ветви

где 11к — напряжение, приложенное к концам ветви; Eh — э. д. с. Аналогичные уравнения можно записать для всех р^ветвей цепи

U г + е2 = Z 2 / 2 ;

Up + ep = Zp/p.

Данная система уравнений в матричной форме^будет иметь вид:

U + e = ZI.

(6.1)

Ut+e,

и2г

Ui+ei

Up + ep

Уравнение для токов ветвей, сходящихся в каком-либо узле,, следующее:

Ф ( / ) Е Л = О,

где k—ветви рассматриваемого узла.

Количество таких уравнений равно количеству узлов цепи q. (количество линейно независимых уравнений q — 1). Так как произведение матрицы соединений П на любой р-мерный вектор,, определяющий те или иные величины ветвей, соответствует сум­ мированию этих величин по узлам, то систему уравнений, полу-

79'


ченных на основании первого закона Кирхгофа, можно записать в матричной форме:

П/ = 0.

(6.2)

где

k — ветви

рассматриваемого контура.

 

Учитывая систему независимых контуров, определенную матри­

цей

Г, и применяя к каждому контуру этой системы

второй за­

кон

Кирхгофа,

имеем р — q + 1 линейно независимых

уравнений.

Эти уравнения

в матричном представлении имеют следующий вид:

 

 

rtu = utr = o.

(6.3)

Индекс t указывает на то, что соответствующая матрица транс­ понирована.

Уравнение (6.3) определяет взаимосвязь между напряжениями ветвей сложной электрической цепи.

Контурные токи 1к связаны с токами ветвей / соотношением (6.4)

в. Уравнения цепи, составленные по методу контурных токов

Если рассматривать цепь как совокупность замкнутых конту­ ров, отдельные части которых принадлежат двум или большему чис­ лу смежных контуров, а остальные (образующие внешний контур) — только одному контуру, то можно выбрать для каждого контура некоторый ток. Выбор контурных токов может происходить раз­ лично, но при этом необходимо помнить, что токи должны быть не­ зависимыми и допускать существование токов в каждой отдельной ветви цепи. Величина и направление контурного тока определяется величиной и направлением действительных токов в участках схе­ мы: действительные токи в участках, принадлежащих только одному контуру, равны соответствующим контурным токам, а действитель­ ные токи в участках, принадлежащих нескольким контурам, явля­ ются алгебраическими суммами соответствующих контурных токов. Переход к уравнениям, записанным относительно контурных токов, может быть реализован с помощью уравнений, составленных отно­ сительно напряжений и токов ветвей (направление токов в ветвях принимается совпадающим с направлением ветвей). Рассмотрим подробно эту операцию перехода. Величины, относящиеся к ветвям схемы, обозначим большими буквами, к контурам — малыми бук­ вами .

Соотношение между параметрами в ветви электрической цепи определяется в матричной форме зависимостью вида (6.1)

U + E = Z I .

80


Для того чтобы исключить U, применим к этому уравнению транспонированную матрицу Г,, при этом э. д. с. в ветвях обра­ щаются в соответствующие э. д. с. в контурах. Учитывая напря­ жение по контуру, равное нулю Г, U = 0, и уравнения связи между токами в ветвях и контурными токами I = Г1к , получаем

 

r ( E = r , z r i , .

 

Далее, введя контурные величины ек и z,

получим

 

е,. = Г,Ё;

z = r , Z r .

(6.5)

Окончательно имеем

 

 

 

 

ек =

 

z 1к,

(6-6)

где z по отношению

к 1к есть подобное преобразование, переводя­

щее поле контурных значений токов в контурные э. д. с.

Элементы матрицы Z соответствуют контурным сопротивлениям.

Решение уравнения

(6.6) дает контурные токи:

 

i K z - 'E K = YeK ,

(6-7)

где Y — контурная

проводимость.

 

 

 

Найдя Y, определяем значения токов в ветвях, для чего преоб­

разуем уравнение (6.7) посредством

оператора

Г:

 

Г 1К = Г z"1 ек

=

Г z- 1 Г, Ё.

(6.7а)

Таким образом, формулы (6.5) — (6.7а) показывают связь меж­ ду э. д. с , токами и полными сопротивлениями ветвей.

Операция преобразования матрицы Г„ отражающей в алгебраи­ ческой форме схему соединения ветвей электрической цепи, фак­ тически представляет собой переход от токов в ветвях к контурным

токам.

 

 

 

 

 

 

 

 

Таким образом, чтобы составить

уравнения,

включающие в

себя

контурные

токи, достаточно

построить

матрицу

сопротивле­

ний

ветвей Z,

матрицу

э. д. с. в

ветвях -Е, матрицу

совпадений

Г. Последовательность

умножения

 

матриц

производится по

фор­

мулам (6.5) и (6.6). Контурные токи

определяются

из

(6.6), а

токи

в ветвях—из

(6.7а).

 

 

 

 

 

 

 

При реализации алгоритмов уравнений контурных токов пре­ дусматривается использование стандартных подпрограмм умноже­ ния симметричных матриц, транспонирования этих матриц и реше­ ния системы алгебраических уравнений.

г. Уравнения цепи, составленные по методу узловых напряжений

В уравнениях цепи, составленных по методу узловых напряже­ ний, неизвестными являются узловые напряжения. Приняв за условный нуль напряжение какого-нибудь одного из /г+1 узлов

4—622

81


цепи, определяем по модулю и фазе напряжения всех остальных п узлов. Изменения, которые необходимо внести в уравнения, сос­ тавленные относительно токов и напряжений, связаны с исключе­

нием из указанных

уравнений цепи токов в ветвях.

 

Чтобы исключить I из (6.1), умножим последнее на обратную

матрицу полных сопротивлений Z - 1 .

 

Поскольку Z"1 =

1/Z, то

 

 

Z"1 U + Z " J E = 1.

(6.8)

Умножая (6.8) на П , а также используя (6.2), получаем

урав­

нение

 

 

 

П Z - 1 U + П Z - 1 Е = 0.

(6.9)

Связь между узловыми напряжениями и и напряжениями U, приложенными к ветвям, в матричной форме может быть записана в следующем виде:

откуда

un = U„

 

(6.10)

и = П,и,.

 

(6.11)

 

 

Если напряжение одного из узлов принять равным 0, то вме­

сто матрицы П используется матрица П ' . Выражение

(6.11) при­

мет вид:

 

 

 

 

 

U ^ I b u , .

(6.12)

Уравнение

(6.9) с учетом (6.12) запишется как

 

 

IT Z - 1 П u/ + П' Z"1

Е = 0.

(6.13)

Обозначим

I I ' Z - 1 П' через

Y, тогда

(6.13) можно

представить

в виде уравнения

 

 

 

 

Y u / +

n / Z ~ 1 E = 0.

(6.14)

Матрица Y представляет собой матрицу проводимостей схемы. Столбцовая матрица э. д. с. в ветвях и столбцовая матрица узло­ вых напряжений соответственно

Е, Е2 и2

ип

При составлении уравнений по методу узловых напряжений, пользуясь геометрической теорией цепей, сначала рассматриваем ветви схем и составляем матрицу соединения (условный нуль отсче­ та выбирается заранее). Элементами этой матрицы являются ко­ эффициенты уравнений, составленных по второму закону Кирхго-

82


фа. Затем при помощи матрицы соединения образовываем систему уравнений для узловых напряжений заданной цепи.

Для узла, принятого за условный нуль отсчета, нет необходи­ мости составлять отдельное уравнение, так как в нем автоматически выполняется условие непрерывности тока. Строение системы урав­ нений, сформулированных на основе метода узловых напряжений, предусматривает число уравнений, равное числу узлов цепи без нулевого узла.

При программировании используются стандартные подпрограм­ мы умножения матриц, транспонирования матриц и решения алгеб­ раических уравнений. Вычислительные операции легко програм­ мируются.

д.

Оценка

методов

узлового

напряжения

и

контурных токов

для алгоритмизации

различных

электротехнических

задач

Электротехнические методы контурных токов и узловых напря­ жений содержат в себе много общего, так как уравнения, состав­ ленные по этим методам, аналогичны по своей структуре. Но при­ менительно к определенному классу физических цепей выбор между этими уравнениями имеет принципиальное значение для построе­ ния алгоритма. От него в свою очередь зависит выбор математичес­ кого метода решения.

Совокупность электротехнического и математического методов определяет конструктивные средства построения алгоритма. Вви­ ду этого, естественно, возникает вопрос, какой из алгоритмов будет наиболее полным? Не существует ни одного электротехнического метода, который был бы наиболее эффективным во всех случаях анализа цепей.

Серьезным ограничением для применения метода контурных токов с целью создания алгоритмов при изучении электромагнитных процессов в сложных электрических цепях является трудность вы­ бора наиболее удобной системы независимых контуров. Транспони­ рованная матрица для одной и той же схемы может быть различного вида, определяемого выбором независимых контуров. Если реали­ зацию выбора независимых контуров предусмотреть в конструкции самого алгоритма, то для целого ряда сложных цепей и цепей не­ плоского характера алгоритм становится громоздким и крайне неэкономным (резко возрастают число арифметических операций, количество запоминаемой информации).

Выбор системы независимых контуров связан со всем комплек­ сом сведений об их возможных комбинациях и является задачей с конечным числом вариантов. В каждом отдельном случае принци­ пиально могут быть произведены полный перебор всех имеющихся возможностей и сопоставление их между собой. Такой полный ана­ лиз в сложных цепях с большим числом контуров связан с расчетом значительного числа вариантов, а в ряде случаев ввиду большого числа возможных комбинаций практически нецелесообразен. В дан-

4*

83