Анализ электрической цепи

Министерство Образования Российской Федерации Хабаровский Государственный Технический Университет

Курсовая работа

по Электротехнике Выполнил: ст-т гр. ПО-22

Бернадцкий С. Л Проверил: д. ф.-м. н., профессор Кузьменко А.П.

Хабаровск, 2004 г.

1. Анализ цепи во временной области методом переменных состояния при постоянных воздействиях

Требуется:

1.1. Составить уравнения состояния цепи для t ??0.

1.2. Найти точные решения уравнений состояния.

1.3. Найти решения уравнений состояния, используя по выбору студента один из численных методов. Вид решаемых уравнений:

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

Схема соответствующая варианту 3, имеет вид:

e (t) = const = E = 5 В.

I (t) = I*1(t) = 2*1(t) А.

1.1 Определим независимые начальные условия для исследуемой электрической цепи

Моменту времени до коммутации (отсутствует источник тока, катушка вырождается в провод, на месте конденсатора — обрыв соединения, так как протекает постоянный ток; источник тока играет роль ключа: I (t)=2*1(t)) соответствует схема:

Используя 1-ый и 2-ой законы Кирхгофа, составим систему уравнений для вышеуказанной схемы (направление обхода контура выберем в соответствии движению часовой стрелки):

Откуда находим:

iL (0) = e (t)/R1=E/R1=5/1000=5*10−3 А.

UC (0) = 0 В.

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

Таким образом, получаем матрицу независимых начальных условий:

Используя 1-ый и 2-ой законы Кирхгофа, составим систему уравнений для исходной схемы:

В соответствии с методом переменных состояния, перепишем систему в дифференциальном виде. Для C элемента запишем ЗТК (), для L элемента ЗНК ().

Полученная система имеет вид:

Решая данную систему относительно неизвестных переменных состояния (iL, UC):

1. Выразим из 1-го уравнения

подставим в остальные уравнения Очевидно те неизвестные, которые не интересуют нас и которые были уже подставлены в другие уравнения могут быть исключены из системы.

2. Выразим из 2-го уравнения

подставим в остальные уравнения

3. Выразим из 3-го уравнения

подставим в остальные уравнения, упростим

4. Выразим из 4-го уравнения

подставим в остальные уравнения, упростим

5. Выразим из 5-го уравнения

подставим в остальные уравнения, упростим

6. Выразим из 6-го уравнения

подставим в остальные уравнения, упростим

7. Выразим из 7-го уравнения, подставим в остальные уравнения, упростим

Откуда:

1.2 Найти точные решения уравнений состояния

Запишем полученные уравнения в матричной форме вида:

(Заметим, что число элементов данной матрицы равно числу реактивных элементов в исследуемой электрической цепи),

матрица соединений, которая содержит элементы, связывающие iL, UC,

(из условия)

— матрица, учитывающая внешние источники.

Таким образом, получаем:

Записываем решение для переменных состояния через экспоненциальную матричную функцию:

Так как, в нашем случае, матрица F = const, то вид уравнения упрощается:

где — матрица обратная матрице A

Найдем собственные значения матрицы A:

Откуда

Разложим экспоненциальную матричную функцию в ряд Тейлора:

Число элементов разложения равно числу переменных состояния (=2).

Находим через найденные выше собственные значения:

Откуда С учетом, получим:

;

Получаем функции состояния:

Построим графики полученных зависимостей:

Рис. 1.1. График зависимости тока на катушке от времени

Рис. 1.2. График зависимости напряжения на конденсаторе от времени.

1.3 Решение численным методом

В качестве метода решения дифференциальных уравнений используем метод Эйлера:

.

Решим дифференциальные уравнения:

Результаты численного решения:

Графики соответствующие полученным данным имеют вид:

График зависимости тока на катушке от времени (синим пунктиром показан график соответствующий численному решению) Рис. 1.3 Графики зависимости тока на катушке от времени.

График зависимости напряжения на конденсаторе от времени (синим пунктиром показан график соответствующий численному решению).

Рис. 2 Графики зависимости напряжения на конденсаторе от времени.

2. Анализ цепи операторным методом при апериодическом воздействии

В данной части работы требуется:

1. Найти нули и полюсы функции передачи и нанести их на плоскость комплексной частоты.

2. Найти переходную и импульсную характеристики для выходного напряжения или тока.

3. Определить изображение по Лапласу входного импульса.

4. Найти напряжение или ток на выходе цепи, используя HU (p) или HI (p) соответственно.

5. Построить на одном графике переходную и импульсную характеристики цепи, на другом — входной и выходной сигналы.

Источник тока отсутствует; входное напряжение имеет зависимость следующего вида:

Параметры для данного варианта:

Исследуемая схема имеет вид:

1. В соответствии с законами Кирхгофа, составим уравнения для данной схемы:

2. Запишем систему в операторном виде, произведя следующие замены:

=>

3. Определим функцию передачи.

Передаточные (системные) функции цепи могут быть определены как отношение выходной величины ко входной. В зависимости от того какие величины входят в определение передаточной функции различают: передаточные функции по напряжению, по току, передаточные сопротивления и проводимости. Функция передачи по току может быть представлена в виде: HU (p) = UН (p)/Uвх (p), где UН (p) и Uвх (p) операторные изображения выходного и входного сигналов, соответственно.

где

Из составленной выше системы найдем ток, используя аналогичные расчеты из 1-ой части:

Тогда:

4. Найдем нули и полюса функции передачи.

Нули:

Полюса:

Нетрудно заметить, что полюсы передаточной функции p1,2 совпадают с собственными значениями ??, 2?? матрицы A. Это может быть дополнительным способом проверки правильности нахождения передаточной функции цепи. Наиболее наглядным способом охарактеризовать передаточную функцию является графическое расположение ее полюсов и нулей на комплексной плоскости, называемое диаграммой полюсов-нулей. Тип используемых элементов, а также структура цепи ограничивают области комплексной плоскости в которых могут располагаться нули и полюсы. В линейной пассивной цепи с потерями (с резистивными элементами) полюсы передаточной функции лежат в левой полуплоскости. Только при этом условии свободные составляющие токов и напряжений затухают. При отсутствии потерь (резистивных элементов) все корни знаменателя будут чисто мнимыми. Нули передаточной функции, корни числителя, при учете потерь могут располагаться в любой части комплексной плоскости. Их положение не связано с характером изменения во времени свободных составляющих токов и напряжений. Отсутствие нулей передаточной функции на мнимой оси физически означает, что при любой частоте гармонического напряжения на входе цепи на выходе будет какое-то напряжение. При отсутствии резистивных элементов все корни числителя передаточной функции (так же как и знаменателя) находятся на мнимой оси. Передаточные функции, полюса которых не лежат в правой полуплоскости комплексной плоскости, называются устойчивыми.

Знание передаточной функции цепи HU (p) позволяет определить переходную h1(t) и импульсную h?(t)?характеристики цепи.

5. Найдем переходную и импульсную характеристики.

Переходная характеристика цепи представляет собой реакцию цепи на воздействие единичной ступенчатой функции (функции Хэвисайда 1(t), функции включения) и может быть найдена как обратное преобразование Лапласа от HU (p)/p:

Переходная характеристика введена в основном по двум причинам:

1. Если определена данная характеристика, то возможно определить реакцию системы при любой форме внешнего воздействия (посредством интеграла Дюамеля).

2. Единичное ступенчатое воздействие скачкообразное, и поэтому является «тяжелым» для любой системы. Следовательно знать реакцию системы именно при таком воздействии. Иные, более плавные, воздействия будут для системы «легче».

Обратное преобразование Лапласа проводим, используя теорему о разложении:

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

где — корни уравнения M (p)=0

Обратное преобразование Лапласа в результате дает следующую зависимость:

(Данное выражение легко может быть проверено на крайних точках временного интервала (при и):)

Импульсная характеристика:

Импульсная характеристика цепи h?(t)?представляет собой реакцию цепи на воздействие единичной импульсной функции?? t) и может быть найдена как обратное преобразование Лапласа от передаточной функции:

Дельта функция (или функция Дирака) определяется как и представляет собой предельный случай импульса очень большого значения и очень малой продолжительности, когда его длительность стремится к нулю, но площадь остается равной единице:

Обратное преобразование Лапласа в результате дает следующую зависимость:

6. Определим изображение входного импульса.

Запишем зависимость входного напряжения от времени для указанной в условии формы сигнала:

Или в операторном виде:

Преобразование Лапласа для 1-ой части выражения табличное, для второй части распишем:

Найденная зависимость потребуется при нахождении напряжения на выходе цепи.

7. Найдем напряжение на выходе цепи.

Так как нас интересует временная зависимость напряжения на выходе цепи, используем для ее нахождения обратное преобразование Лапласа:

Обратное преобразование Лапласа в результате дает следующую зависимость:

Полученная зависимость имеет вид:

8. Построим графики, соответствующие полученным результатам.

Рис. 2 .1. Переходная характеристика электрический цепь анализ воздействие Рис. 2.2. Импульсная характеристика Рис. 2.3. Входное и выходное напряжения

3. Анализ цепи частотным методом при апериодическом воздействии

Требуется:

2.3.1. Найти и построить амплитудно-фазовую (АФХ), амплитудно-частотную (АЧХ) и фазо-частотную (ФЧХ) характеристики функции передачи HU (jw).

2.3.2. Определить полосу пропускания цепи по уровню 0,707 |H (j?)|макс.

2.3.3. Найти и построить амплитудный и фазовый спектры входного сигнала. Определить ширину спектра входного сигнала по уровню 0,1|F (j?)|макс.

2.3.4. Сопоставляя спектры входного сигнала с частотными характеристиками цепи, дать предварительные заключения об ожидаемых искажениях сигнала на выходе цепи. Сверить эти качественные оценки с сигналом на выходе, полученным в п. 2.2.5.

2.3.5. Найти и построить амплитудный и фазовый спектры выходного сигнала.

2.3.6. Определить выходной сигнал по вещественной или мнимой частотной характеристике, используя приближённый метод Гиллемина.

2.3.1. Найти и построить амплитудно-фазовую (АФХ), амплитудно-частотную (АЧХ) и фазо-частотную (ФЧХ) характеристики функции передачи HU (jw).

Для анализа цепи используем ранее полученную функцию передачи

Для исследования частотных характеристик воспользуемся тем фактом, что необходимое нам преобразование Фурье есть частный случай преобразования Лапласа, когда действительная часть оператора Лапласа, т. е.. Таким образом, необходимо произвести подстановку в функцию передачи p=jw:

Амплитудно-частотная характеристика (АЧХ) является одной из самых важных характеристик любой цепи и позволяет исследовать искажения вносимые цепью в спектр входного сигнала. Наличие частотно — зависимых элементов (L и C) в исследуемой цепи приводит к неравномерному изменению составляющих спектра входного сигнала. Наиболее простой способ получения АЧХ цепи — это замена в выражении для HU (p) операторной переменной p на мнимую частоту jw и нахождение модуля полученной комплексной функции частоты: |HU (jw)|.

Рис. 3.1 АЧХ схемы.

ФЧХ схемы определяется как аргумент передаточной функции при w=0.

Рис. 3.2. ФЧХ схемы.

Амплитуднофазочастотная характеристика цепи (годограф) связывает воедино изменение коэффициента передачи (в нашем случае, по напряжению — |HU (jw)| и фазового сдвига между выходным и входным напряжением во всем диапазоне частот. Годограф включает сведения которые содержатся как в АЧХ, так и в ФЧХ.

Годограф является параметрической кривой, параметром которой является частота w. Длина вектора, проведенного из начала координат к какойлибо точке годографа соответствует абсолютному значению передаточной функции на этой частоте |HU (jw)|, а угол между ним и положительным направлением вещественной оси — аргументу передаточной функции??? arg (HU (jw)). На рис. 7 представлен годограф для рассматриваемой цепи. Нулевой частоте соответствует точка с координатой 0 на вещественной оси, очень большой (в пределе бесконечной) частоте соответствует точка с координатой 0.2641 на вещественной оси. На этих граничных частотах влияние реактивных элементов на фазовый сдвиг отсутствует.

Рис. 3.3. АФХ цепи.

Определить полосу пропускания цепи по уровню 0,707 |H (j?)|макс.

По построенной характеристике может быть определена полоса пропускания. Полосой пропускания цепи называют диапазон частот для которых коэффициент передачи не более чем в ?2 отличается от его максимального значения. Это же соответствует снижению уровня сигнала на 3 дБ. Для рассматриваемой цепи максимальное значение передаточной функции составляет |HU (jw)|max = 0.2641. Границе полосы пропускания соответствует значение передаточной функции |HU (jw)|max/??2 = 0.707|HU (jw)|max =0.187. Это значение достигается на частоте w=97 551.6407 c-1. Таким образом полоса пропускания равна? w=[97 551.6407,?]. Если основные гармоники сигнала лежат в этой полосе частот, то не происходит искажения формы сигнала.

Найти и построить амплитудный и фазовый спектры входного сигнала. Определить ширину спектра входного сигнала по уровню 0,1|U (j?)|макс.

АЧХ входного сигнала находится как модуль от спектральной характеристики U (j?), которую получаем на основе ранее найденного изображения сигнала.

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

Построим амплитудно-частотную характеристику входного сигнала .

Рис. 3.4. АЧХ входного сигнала.

Определяя по уровню 0,1|U (j?)|макс. ширину спектра, получим? w=2.703*105 с-1

Ограничивая спектр сигнала определенной по уровню 0,1|U (j?)|макс. шириной спектра? w, мы учитываем?? W?w/Wt) · 100% ?? 96?? от полной энергии Wt сигнала. Это следует из использования теоремы Рейли для расчета данного отношения (

)).

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

ФЧХ входного сигнала определяется как аргумент от спектральной характеристики:

Рис. 3.5. ФЧХ входного сигнала.

Найти и построить амплитудный и фазовый спектры выходного сигнала.

Используя известную функцию передачи можем найти АЧХ и ФЧХ выходного сигнала.

АЧХ входного сигнала находится как модуль от спектральной характеристики U (j?), умноженный на модуль функции передачи Рис. 3.6. АЧХ выходного сигнала.

ФЧХ входного сигнала находится как сумма аргументов спектральной характеристики U (j?), и функции передачи .

Рис. 3.7. ФЧХ выходного сигнала.

Определить выходной сигнал по вещественной или мнимой частотной характеристике, используя приближённый метод Гиллемина.

Метод Гиллемина является одним из методов позволяющих восстановить функцию времени (какой — либо сигнал) по известной вещественной (или мнимой) частотной характеристике. Метод основан на такой аппроксимации, когда аппроксимирующая частотную характеристику функция либо ее производные состоят из последовательности бесконечно коротких импульсов. Последовательность бесконечно коротких импульсов представляет собой заданную функцию в так называемой квантованной форме. Погрешность метода преимущественно связана со ступенчатым характером аппроксимирующей функции. Уменьшение этой погрешности требует увеличения общего числа членов в аппроксимации. Исходная частотная характеристика аппроксимируется кусочнолинейным образом, после чего два последовательных дифференцирования позволяют свести аппроксимирующую функцию к последовательности бесконечно коротких импульсов. Окончательное выражение для искомой функции времени f (t) полученной по вещественной частотной характеристике имеет вид:

Здесь ak — величины бесконечно коротких импульсов, wk — координаты импульсов на частотной оси.

Вещественная частотная характеристика Gн (w) может быть определена как:

Рис. 3.7. Вещественная частотная характеристика выходного сигнала.

Интерполируем полученную зависимость, взяв количество разбиений равным 20.

Диапазон частот ограничим величиной примерно равной ширине спектра — .

Тогда интервал разбиений равен /20 =

Рис. 3.8 Интерполированная зависимость.

Найдем производные первого порядка в точках разбиения и построим полученную зависимость на графике.

Рис. 3.8 График первой производной.

Найдем производные первого порядка в точках разбиения и построим полученную зависимость на графике.

Рис. 3.8 График второй производной.

Для определения выходной зависимости напряжения от времени, на основе полученных данных, воспользуемся формулой

Здесь G''k — величины бесконечно коротких импульсов, wk — координаты импульсов на частотной оси.

Восстановленная временная зависимость имеет вид:

Рис. 3.8 Восстановленная по методу Гиллемина временная зависимость выходного напряжения (пунктиром) и временная зависимость выходного напряжения определенная во 2-ой части работы.

4. Анализ цепи частотным методом при периодическом воздействии

Требуется:

1. Разложить в ряд Фурье заданную периодическую последовательность импульсов и построить её амплитудный и фазовый спектры.

2. Построить на одном графике заданную периодическую последовательность импульсов и её аппроксимацию отрезком ряда Фурье, число гармоник которого определяется шириной амплитудного спектра входного сигнала, найденной в п. 2.3.3.

3. Используя рассчитанные в п. 2.3.1 АЧХ и ФЧХ функции передачи цепи, определить напряжение или ток на выходе цепи в виде отрезка ряда Фурье.

4. Построить напряжение или ток на выходе цепи в виде суммы гармоник найденного отрезка ряда Фурье. Графики по пп. 2.4.2 и 2.4.4 построить на одном масштабе времени и разместить их на одном листе один под другим.

Данные задачи: Um = 10 В, tи = 610−5 c, T = 1810−5 c.

4.1. Разложить в ряд Фурье заданную периодическую последовательность импульсов и построить её амплитудный и фазовый спектры.

Ряд Фурье периодической функции можно представить в виде:

где — амплитуда i-ой гармоники, — начальная фаза.

— основная (главная) частота рассчитывается по формуле:

T — период сигнала; рад/с Далее, в задании необходимо аппроксимировать входную и выходную зависимость, количество гармоник, необходимых для аппроксимации выбираем исходя из того, чтобы частоты этих гармоник помещались в ширину спектра, определенную в 3-ей части.

(с учетом округления в большую сторону) Для нахождения амплитудных составляющих воспользуемся спектральной плотностью сигнала, найденной в 3-ей части курсовой работы.

Для нахождения фазы также воспользуемся ею:

4.2. Построить на одном графике заданную периодическую последовательность импульсов и её аппроксимацию отрезком ряда Фурье, число гармоник которого определяется шириной амплитудного спектра входного сигнала, найденной в п. 2.3.3.

Рис. 4.1. Амплитудный спектр входного сигнала Рис. 4.2. Фазовый спектр входного сигнала

4.3. Используя рассчитанные в п. 2.3.1 АЧХ и ФЧХ функции передачи цепи, определить напряжение или ток на выходе цепи в виде отрезка ряда Фурье.

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

Рис. 4.3. Амплитудный спектр выходного сигнала Рис. 4.4. Фазовый спектр выходного сигнала Рис 4.3. Входной сигнал и его аппроксимация рядом Фурье Таким образом по найденным коэффициентам и значениям начальных фаз, строим функцию зависимости сигнала от времени.

Рис. 4.4. Выходной сигнал и его аппроксимация рядом Фурье