Главная Случайная страница Контакты | Мы поможем в написании вашей работы! | ||
|
В четвертой главе была рассмотрена процедура составления дифференциальных уравнений путем приведения исходной -цепи к резистивной форме. Возвратимся теперь к электрической цепи, изображенной на рис. 5.11 а) и остановимся на решении дифференциальных уравнений, используя матричный метод анализа переходных процессов.
Уравнения динамики цепи имеют вид:
, (5.54)
где – ЭДС источника электроэнергии,
(5.55)
Поведение -цепи зависит от собственных чисел матрицы , которые мы можем изменять путем вариации значений активных сопротивлений, индуктивности и емкости.
Рассмотрим несколько частных случаев.
5.5.1. Сопротивления ,
Элементы матриц уравнения (5.23) принимают следующие значения:
, (5.56)
Электрическая цепь становится неразветвленной, состоящей из последовательно соединенных и (рис. 5.11 б)).
Поскольку активные сопротивления в цепи отсутствуют, динамическая система становится консервативной, и теоретически переходный процесс должен протекать без рассеяния (потерь) энергии.
До замыкания ключа ток через индуктивность и напряжение на емкости равны нулю. Следовательно, вектор начальных условий .
Для определения характеристического уравнения системы (5.22) составим определитель (детерминант) матрицы:
(5.57)
Поскольку собственные значения матрицы численно совпадают с корнями характеристического уравнения (5.25), то они равны: , где – собственная частота цепи (рис. 5.11 б)).
Для удобства проведения математических преобразований воспользуемся операторным методом. Применение операторного метода позволяет в области изображений (частотной области) обращаться с матричным дифференциальным уравнением как с алгебраическим.
Процедура решения состоит в следующем. Сначала с помощью таблиц преобразований Лапласа осуществляется операция перевода дифференциального уравнения из области оригиналов (временная область) в область изображений (частотная область). Затем выполняются решения алгебраического уравнения в частотной области относительно вектора состояния, записанного в операторной форме. На последнем этапе также с помощью таблиц преобразований Лапласа производится переход от алгебраического уравнения в частотной области к уравнению во временной области (переход из области изображений в область оригиналов). Все данные таблиц преобразований Лапласа в среде MatLAB можно найти в приложении Symbolic Toolbox - пакете символьной математики [20].
При нулевых начальных условиях переход в область изображений для слагаемых уравнения (5.22) будет заключаться в использовании следующих данных:
, , , (5.58)
где – оператор Лапласа, есть знак соответствия между выражениями во временной области (слева от знака) и частотной области (справа от знака).
С учетом (5.26) получим операторную форму записи матричного дифференциального уравнения (5.22):
(5.59)
Обращаясь с формулой (5.27) как с алгебраическим матричным уравнением, мы вправе записать: .
Его решение относительно вектора переменных состояния в операторной форме:
(5.60)
С учетом (5.24) уравнение (5.28) примет вид:
(5.61)
Обратная матрица в уравнении (5.29) определится с помощью стандартной процедуры инверсии: , где – присоединенная для матрица, – определитель матрицы . В рассматриваемом случае инверсная матрица для :
.
Таким образом, в операторной форме ток через индуктивность:
(5.62)
Напомним, что согласно таблицам преобразований Лапласа, выполняется условие соответствия:
(5.63)
Воспользуемся этим соответствием для перехода из области изображений в область оригиналов. Току через индуктивность, определенному формулой (5.30), во временной области соответствует выражение:
, (5.64)
где – амплитуда тока. Таким образом, ток в электрической цепи изменяется по закону синуса.
Напряжение на емкости в операторной форме, согласно (5.29), равно:
. (5.65)
Из правил выполнения преобразований Лапласа известно, что интегрирование функции в области оригиналов (при нулевых начальных условиях) соответствует делению на оператор этой функции в области изображений. Обращаясь вновь к уравнению (5.33) и соответствию (5.31), мы видим, что в знаменателе (5.33) содержится оператор . Следовательно, выражение в частотной области соответствует интегралу правой части (5.31) во временной области.
Поэтому для напряжения на емкости как функции времени имеем: , где – постоянная интегрирования, определяемая из начальных условий. Действительно, напряжение . Тогда постоянная интегрирования .
Подставляя это выражение в записанное выше выражение для , мы получим:
. (5.66)
Напряжение на индуктивности определим по формуле:
. (5.67)
Для наглядного представления кривых, определяемых формулами (5.32), (5.34) и (5.35), построим графики с помощью MatLAB и выведем их на печать. Составим программу, представленную файлом .
Программа содержит текстовые комментарии, которые начинаются со знака «%». Затем задается вектор времени , изменяющийся от 0 до 6 (при шаге дискретности ). Принятые расчетные значения: ; . Построение графиков выполняется оператором . По окончании построения графиков оператором задается пауза. Обратите внимание на то, что все вычисления производятся в процессе выполнения оператора .
Файл
%File “sah36.m”.
%Plots for electrical system.
%LC-circuit with DC source
Im=0.5; E=2;
x=0:0.02*pi:6.0*pi;
subplot(221);
plot(x, Im*sin(x), x, E*(1-cos(x)), x, E*cos(x)), grid
pause
Графические построения приведены на рис.5.12.
Из приведенных графиков следует, что кривая напряжения на емкости располагается выше оси абсцисс (то есть не принимает отрицательных значений). Напряжение имеет постоянную составляющую и составляющую , изменяющуюся с собственной частотой . В любой момент времени сумма напряжений на емкости и индуктивности равна постоянному значению, равному , что отвечает закону Кирхгофа:
Дата публикования: 2014-11-03; Прочитано: 413 | Нарушение авторского права страницы | Мы поможем в написании вашей работы!