Использование преобразования Лапласа при расчете переходных процессов

Использование преобразования Лапласа при расчете переходных процессов

Преобразование Лапласа

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

. (4.1)

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

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

Пусть функция задана для положительных значений и равна нулю при , причем она может не удовлетворять условию (4.1). Умножая на , где и , получим, пользуясь выражением (3.1):

(4.2)

На основании (3.2) тогда имеем

.

Отсюда видно, что

(4.3)

Вводя комплексную переменную , запишем вместо (4.2)

. (4.4)

Полученное соотношение, преобразующее вещественную функцию вещественного переменного в функцию комплексного переменного , называется преобразованием Лапласа. Функцию называют изображением функции , а саму функцию называют оригиналом. Символически связь между изображением и оригиналом записывается так:

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

Теперь в ином виде представим выражение (4.3), переходя к переменной , Так как , то , a пределы интегрирования в (4.3) при изменении от до теперь будут соответственно и .Тогда получаем:

. (4.5)

Это соотношение, которое известно в математике как формула обращения Римана-Моллина, иногда называют обратным преобразованием Лапласа.

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

В выражении (4.5) интегрирование ведется по прямой, лежащей на плоскости комплексного переменного и отстоящей от оси мнимой величины на расстоянии . В теории функций комплексного переменного^{**) Хапланов М.Г. Теория функций комплексного переменного. М., 2005}

Романовский П.И. Ряды Фурье, теория поля, аналитические и специальные функции, преобразование Лапласа. М., 2001

Привалов И.И. Введение в теорию функций комплексного переменного. М.-Л., 2008.⁾ доказывается, что интегрирование в выражении (4.5) может быть произведено по замкнутому контуру (рис. 4.1). Контур интегрирования должен охватывать все полюсы подинтегралыюй функции, т.е. точки плоскости комплексного переменного, в которых подинтегральная функция выражения (4.5) обращается в бесконечность. Вычисление интеграла при этом сводится к определению суммы вычетов (обозначаемых буквами ) подинтегральной функции в полюсе, т.е.

 (4.6)

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

Значения комплексного переменного , при которых знаменатель обращается в нуль, а следовательно, подинтегральная функция -- в бесконечность, и будут полюсами под интегральной функции. В данном случае полюсы являются корнями уравнения . Число полюсов подинтегральной функции определяется числом корней этого уравнения.

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



Вычет подинтегралъной функции в полюсе порядка находится по формуле

Если полюс простой, т.е. первого порядка, то вычет в этом полюсе определяется выражением

(4.8)

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

Тогда оригинал по этому изображению находится с помощью (4.5)

или, так как полюс подинтегральной функции , то согласно (4.8) вычет подинтегральной функции в этом полюсе



и на основании (4.6) имеем , т.е. получаем исходную функцию. Символически можно записать:

(4.9)

Пусть теперь задана экспоненциальная функция

Изображение этой функции:

По найденному изображению определим теперь оригинал

Полюсом подинтегральной функции является . Вычет в этом полюсе равен

Таким образом, можем записать:

(4.10)

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

. (4.11)

Для дельта функции изображение имеет вид:

. (4.12)

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

Пусть изображение является рациональной дробью

,

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

. (4.13)

Находя по формуле (4.11) оригинал функции , получим выражение для нахождения по изображению оригинала, т.е. функции :

. (4.14)

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

Выражению (4.14) можно придать иной вид для случая, если один из корней уравнения равен нулю. Обозначим , где - полином, не имеющий нулевых корней. Выделим в выражении (4.13) слагаемое, соответствуют ее нулевому корню. Так как в принятом обозначении для производной функции имеем , то вместо (4.13) получаем:

.

Учитывая выражения (4.9) и (4.11), получаем для нахождения оригинала по заданному изображению выражение

, (4.15)

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

Основные свойства преобразования Лапласа

Рассмотрим теперь некоторые свойства преобразования Лапласа, необходимые для решения многих задач по переходным процессам в радиотехнических цепях.

Если функция , где - постоянная величина, а функция имеет изображение , то

. (4.16)

Если задана сумма функций , то, осуществляя прямое преобразование Лапласа, находим:

, (4.17)

где , .

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

. (4.18)

Изображение производной по времени от функции , т.е. изображение от ; находим, производя преобразование Лапласа:

,

или, так как

является изображением функции ,a , то

(4.19)

При нулевых начальных условиях, т.е. при , получим:

.

Аналогичные вычисления позволяют найти изображение для интеграла от функции :

, (4.20)

где

при .

Если имеем две функции и их изображения, т.е. если дано , а , то функция , образованная по формуле

,

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

.

Полученное соотношение носит название теоремы свертывания (называемое также теоремой Бореля) и записывается так:

. (4.21)

Это выражение может иметь и другой вид:

(4.22)

При решении задач и выполнении технических расчетов с использованием преобразования Лапласа приходится находить изображение функции по оригиналу и наоборот. Для облегчения расчетов пользуются обычно справочниками, содержащими таблицы функций и их изображения⁾ Диткин В.А., Кузнецов Н.И. Справочник по операционному исчислению. М.-Л., 2001.

Диткин В.А., Прудников А.П. Интегральные преобразования и операционное исчисление. М., 2001.⁾. В этих справочниках часто приводятся изображения, найденные на основании преобразования Карсона:

, (4.23)

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

Расчет переходных процессов с помощью преобразования Лапласа

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

Рассмотрим уравнение Кирхгофа для последовательного контура, на который действует э.д.с. :

(4.24)

При решении задачи операторным методом от этого уравнения необходимо перейти к равенству, связывающему изображения функций и , преобразуя обе части уравнения (4.24) по Лапласу.

Если известны изображения функций, т.е. и , то на основании рассмотренных выше свойств преобразования Лапласа получаем вместо уравнения (4.24) выражение

. (4.25)

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

Из выражения (4.25) находим изображение искомого тока

,

где - импеданс цепи в операторной форме записи. Как видно, операторная форма записи импеданса, или операторное сопротивление цепи. получается из выражения импеданса цепи путем замены на величину . Зная изображение искомой функции, можно найти ее оригинал, т.е. ток .

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

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

1) Составляются уравнения Кирхгофа для исследуемой цепи.

2) Находятся изображения заданных и искомых функций времени.

3) Записывается операторное уравнение, соответствующее исходному.

4) Решается алгебраически операторное уравнение относительно изображения, искомой функции.

5) По изображению искомой функции находится ее оригинал. Для этого применяется способ нахождения вычетов или теорема разложения.

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

,

то операторное уравнение будет

,

оригинал найдем с помощью обратного преобразования Лапласа, которое вычислим с помощью вычетов

полюсами подинтегральной функции являются значения , где

, , , .

Вычет в первом полюсе

.

Вычет во втором полюсе

.

Тогда искомый ток в контуре

.

Теперь рассмотрим пример, когда начальные условия задачи не равны нулю (рис.4.2). Предварительно ключ разомкнут и к цепи, состоящей из сопротивления и индуктивности , подключен последовательно источник гармонического напряжения . В момент ключ замыкается и цепь оказывается короткозамкнутой. Тогда для данной цепи при уравнение Кирхгофа имеет вид:

.

Вводя изображение для тока , запишем уравнение в операторной форме

,

откуда

.

Оригинал для изображения , т.е. искомый ток, найдем в данном случае не с помощью вычетов, а при помощи теоремы разложения, пользуясь формулой (4.14).Так как в данном примере, как видно из выражения для , , а , то . Корнем уравнения является . Подставляя полученные значения в формулу (4.14),находим ток

.

Значение находится из выражения для тока, протекающего в цепи под действием напряжения . В момент мгновенное значение тока в индуктивности имеет вид:

,

где

.

Тогда полученное выражение для тока , описывающее переходный процесс в цепи, принимает вид:

.

Переходной процесс в цепи, замкнутой накоротко, обусловлен энергией магнитного поля, имеющегося у катушки индуктивности в момент коммутации (при ).

Использованная литература

1 Афанасьев В.П. и др. Теория линейных электрических цепей: Учебное пособие для вузов. - М.: Высш. шк., 1973. - 592 с.

2 Белецкий А.Ф. Теория линейных электрических цепей: Учебное пособие для вузов. - М.: Радио и связь, 1986. - 544 с.

З Гринченко Л.В. Методические указания к курсовой работе по дисциплине “Основы радиоэлектроники”. - Х. 2003. -30 с.

4 3ернов Н.В., Карпов В.Г. Теория радиотехнических цепей - М. Энергия, 1972. - 715 с.

5 Попов В.П.Основы теории цепей: Учеб для вузов. M.: Высш. шк., 1985. - 490с.