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

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

, (2.1)

где х и у - входная и выходная величины. Если в данное уравнение вместо x(t) и y(t) подставить функции X(s) и Y(s) комплексного переменного s такие, что

и , (2.2)

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

a₂ s² Y(s) + a₁ s Y(s) + a₀ Y(s) = b₁ X(s) + b₀ X(s).

Такой переход от ДУ к алгебраическому уравнению называется преобразованием Лапласа, формулы (2.2) соответственно формулами преобразования Лапласа, а полученное уравнение - операторным уравнением.

Новые функции X(s) и Y(s) называются изображениями x(t) и y(t) по Лапласу, тогда как x(t) и y(t) являются оригиналами по отношению к X(s) и Y(s).

Переход от одной модели к другой достаточно прост и заключается в замене знаков дифференциалов на операторы sⁿ, знаков интегралов на множители , а самих x(t) и y(t) - изображениями X(s) и Y(s).

Таблица 1.1 - Преобразования Лапласа

Оригинал x(t)	Изображение X(s)
d-функция

t
t²
tⁿ
e^-^a^t
a^.x(t)	a^.X(s)

x(t - a)	X(s)^.e^-^a^s
	s^n.X(s)

Таблица 1.2 - Формулы обратного преобразования Лапласа (дополнение)

Изображение X(s)	Оригинал x(t)
	a Î R, M Î R (a и М - действительные числа)	M^.e^-^a^t
a = a + j^. w M = C + j^.D (a и М – комплексные числа)	2^.e^a^*t.[C^.cos(w^.t) - D^.sin(w^.t)] для пары комплексных корней

Для обратного перехода от операторного уравнения к функциям от времени используется метод обратного преобразования Лапласа. Общая формула обратного преобразования Лапласа:

, (2.3)

где f(t) - оригинал, F(jw) - изображение при s = jw, j - мнимая единица, w - частота.

Эта формула достаточно сложна, поэтому были разработаны специальные таблицы (см. таблицы 1.1 и 1.2), в которые сведены наиболее часто встречающиеся функции F(s) и их оригиналы f(t). Они позволяют отказаться от прямого использования формулы (2.3). Более полные таблицы преобразований Лапласа можно найти, например, в [22, 23].

Существует несколько теорем преобразования Лапласа.

Теорема 1. Теорема линейности. Изображение суммы функций равно сумме изображений, то есть, если f₁ имеет изображение F₁(s) (или более кратко f₁ « F₁(s) ), f₂ « F₂(s) и т.д., то

a₁^.f₁ + a₂^.f₂ + … + a_n^.f_n « a₁^.F₁(s) + a₂^.F₂(s) + … + a_n^.F_n(s).

Теорема 2. Теорема дифференцирования. Если f(t) имеет изображение F(s), то при нулевых начальных условиях (т.е. при f(0) = 0, f’(0) = 0 и т.д.) производные f(t) будут иметь изображения:

f’(t) « s^.F(s) – для первой производной,

f ”(t) « s^2.F(s) – для второй производной,

f⁽ⁿ⁾(t) « sⁿ^.F(s) – для n-й производной.

При ненулевых начальных условиях:

f’(t) « s^.F(s) – f(0) – для первой производной,

f ”(t) « s^2.F(s) – s^.f(0) – f’(0) – для второй производной,

f⁽ⁿ⁾(t) « s^n.F(s) – s^n-1.f(0) - s^n-2.f’(0) - … - f^(n-1)(0) – для n-й.

Теорема 3. Теорема смещения.

f(t)^.e^a^×^t « F(s - a).

Например, если 1(t) « (см. таблицу 1.1), то 1^.e^a^×^t « .

Теорема 4. Теорема запаздывания.

f(t - t) « F(s)^.e^-^t^×^s,

где t - запаздывание по времени.

Например, если 1(t) « , то 1(t - t) « .

Теорема 5. Теорема интегрирования.

Теорема 6. О начальных и конечных значениях.

где f(0) – начальное значение функции (при t = 0),

f_уст – конечное (значение в установившемся режиме).

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

единичное ступенчатое воздействие имеет изображение X(s) = ,

дельта-функция X(s) = 1,

линейное воздействие X(s) = .

Пример. Решение ДУ с использованием преобразований Лапласа.

Допустим, входной сигнал имеет форму единичного ступенчатого воздействия, т.е. x(t) = 1. Тогда изображение входного сигнала, согласно таблице 1.1, имеет вид X(s) = .

Производим преобразование исходного ДУ по Лапласу и подставляем X(s):

s²×Y(s) + 5×s×Y(s) + 6×Y(s) = 2×s×X(s) + 12×X(s),

s²×Y(s) + 5×s×Y(s) + 6×Y(s) = 2×s + 12 ,

Y(s)×(s³ + 5s² + 6s) = 2×s + 12.

Определяется выражение для Y:

Оригинал полученной функции отсутствует в таблице оригиналов и изображений. Для решения задачи его поиска дробь разбивается на сумму простых дробей с учетом того, что знаменатель может быть представлен в виде s(s + 2)(s + 3):

= = - + .

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

y(t) = 2 - 4^.e^-2^t + 2^.e^-3^t. ¨

При решении ДУ с использованием преобразований Лапласа часто встает промежуточная задача разбиения дроби на сумму простых дробей. Существуют два пути решения этой задачи:

- путем решения системы уравнений относительно коэффициентов числителей,

- путем расчета коэффициентов числителей по известным формулам.

Общий алгоритм разбиения дроби на сумму простых дробей:

шаг 1 – определяются корни знаменателя s_i (знаменатель дроби приравниватся к нулю и решается полученное уравнение относительно s);

шаг 2 – каждому корню ставится в соответствие простая дробь вида , где М_i – неизвестный коэффициент; если имеет место кратный корень с кратностью k, то ему ставится в соответствие k дробей вида ;

шаг 3 – определяются коэффициенты M_i по одному из вариантов расчета.

Первый вариант. Определение M_i с помощью системы уравнений.

Все дроби приводятся к одному знаменателю, затем путем сравнения коэффициентов при равных степенях s числителя полученной дроби и числителя исходной определяется система из n уравнений, где n – степень знаменателя (количество корней s_i и коэффициентов M_i). Решение системы относительно M_i дает искомые коэффициенты.

Пример. Декомпозиция дроби из предыдущего примера. В исходной дроби n = 3, поэтому решение уравнения s³ + 5s² + 6s = 0 дает 3 корня: s₀ = 0, s₁ = -2 и s₂ = -3, которым соответствуют знаменатели простых дробей вида s, (s – s₁) = (s + 2) и (s – s₂) = (s + 3). Исходная дробь декомпозируется на три дроби:

= = + + .

Далее дроби приводятся к общему знаменателю:

= .

Сравнивая получившуюся дробь с исходной, можно составить систему из трех уравнений с тремя неизвестными (при 2-й степени s в исходной дроби стоит 0, при 1-й стоит 2, свободный член равен 12):

М₀ + М₁ + М₂ = 0 M₀ = 2

5^.М₀ + 3^.М₁ + 2^.М₂ = 2 à M₁ = -4

6^.М₀ = 12 M₂ = 2

Следовательно, дробь можно представить как сумму трех дробей:

= - + .¨

Второй вариант. Определение коэффициентов M_i по формулам.

Также как и в 1-м варианте необходимо найти корни знаменателя исходной дроби вида . Для определения M_i существуют формулы для каждого вида корней:

- Для нулевого корня s_i = 0 знаменатель исходной дроби можно записать в виде A(s) = s^.A₁(s); тогда коэффициент M_i можно определить как .

- Для ненулевого некратного корня (действительного или комплексного) s_i:

где A’(s) – производная знаменателя по s.

Примечание - Комплексные корни при решении уравнений появляются комплексно-сопряженными парами вида s_i = a_i ± j×w_i , где a_i – действительныя часть корня, w_i – мнимая часть, j – мнимая единица. Поэтому коэффициенты для этих корней также будут комплексно-сопряженными: M_i = c_i ± d_i. То есть достаточно определить коэффициент только для одного корня, для парного корня он будет комплексно-сопряженным.

- Для корня s_i кратности k исходная дробь может быть представлена в виде

;