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

Материал из sawiki
Перейти к навигации Перейти к поиску

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

Прямое и обратное преобразование Лапласа некоторых функций:

  • $$x(t) \supset \frac{1}{p}$$
  • $$x(t)t^n \supset \frac{n!}{p^{n+1}}$$
  • $$x(t)t^{\alpha}e^{\beta t} \supset \frac{n!}{(p-\beta)^{n+1}}, \beta \in \mathbb{C}$$
  • $$x(t)t^{\alpha}\cos{\omega t} \supset \frac{\varGamma (\alpha + 1)}{2}\left(\frac{1}{(p-\omega i)^{\alpha+1}} + \frac{1}{(p+\omega i)^{\alpha+1}} \right)$$
  • $$x(t)t^{\alpha}\sin{\omega t} \supset \frac{\varGamma (\alpha + 1)}{2}\left(\frac{1}{(p-\omega i)^{\alpha+1}} - \frac{1}{(p+\omega i)^{\alpha+1}} \right)$$
  • $$x(t)\sin{\omega t} \supset \frac{\omega}{p^2 + \omega^2}$$
  • $$x(t)\cos{\omega t} \supset \frac{p}{p^2 + \omega^2}$$
  • $$x(t)e^{\beta t}\sin{\omega t} \supset \frac{\omega}{(p-\beta)^2 + \omega^2}$$
  • $$x(t)e^{\beta t}\cos{\omega t} \supset \frac{p-\beta}{(p-\beta)^2 + \omega^2}$$
  • $$x(t)e^{\beta t} \supset \frac{1}{p-\beta}$$
  • $$x(t)\operatorname{sh}{\omega t} \supset \frac{\omega}{p^2 - \omega^2}$$
  • $$x(t)\operatorname{ch}{\omega t} \supset \frac{p}{p^2 - \omega^2}$$
  • $$x(t)e^{\beta t}\operatorname{sh}{\omega t} \supset \frac{\omega}{(p-\beta)^2 - \omega^2}$$
  • $$x(t)e^{\beta t}\operatorname{ch}{\omega t} \supset \frac{p-\beta}{(p-\beta)^2 - \omega^2}$$,
    $$\text{где } x(t)$$ - функция Хевисайда.

Свойства:

  • $$f(at) \supset \frac{1}{a}F(\frac{p}{a}) \, , \, a > 0$$
  • $$\frac{1}{b}\,f(\frac{t}{b}) \supset F(pb) \, , \, b > 0 \,;\, b=\frac{1}{a}$$
  • $$x(t-a)f(t-a) \supset e^{-ap}F(p)$$
  • $$x(t)f(t+a) \supset e^{ap}\left(F(p)-\int\limits_0^a f(t)e^{-pt} dt \right) $$
  • $$x(t)e^{\beta t} \supset F(p+\beta)$$
Дифференцируемость:
  • $$f^{(k)}(t) \supset p^{k}F(p) - p^0f^{(k-1)}(0) - pf^{(k-2)}(0) - \ldots - p^{k-1}f(0)$$
  • $$x(t)(-t)^k \supset F^{(k)}(p)$$
Интегрируемость:
  • $$\int\limits_0^t f(\tau) d\tau \supset \frac{F(p)}{p}$$
  • $$\frac{f(t)}{t} \supset \int\limits_0^{+\infty}F(z) dz$$

Примеры:

Нахождение фундаментальной матрицы:

\( \dot X = AX\, , \, X \in \mathbb{R}^{n}, \,A \in \mathbb{R}^{n \times n} \)
X(t) - фундаментальная матрица, если: \(\left\{ \begin{array}{rcl} \frac{dX}{dt}=AX \\ X(0) = I \end{array} \right. \)
Пусть $$X(t) \leftrightarrow Y(p)$$, тогда преобразуем систему по Лапласу $$\Rightarrow$$
$$\Rightarrow PY-I=AY \Rightarrow Y = (PI-A)^{-1}$$

Пример:

\(\left\{ \begin{array}{} \dot X = \begin{pmatrix} 1 & 2 \\ -3 & 1 \end{pmatrix}X \\ X(0) = I \end{array} \right.,\quad Y = \begin{pmatrix} p-1 & -2 \\ 3 & p-1 \end{pmatrix}^{-1} = \frac{1}{(p-1)^2+6}\begin{pmatrix} p-1 & 2 \\ -3 & p-1 \end{pmatrix} \)
От каждого элемента матрицы Y нужно вычислить преобразование Лапласа, и получим:
\( X(t) = \begin{pmatrix} e^t\cos{\sqrt{6}t} & \frac{2}{\sqrt{6}}e^t\sin{\sqrt{6}t} \\ -\frac{3}{\sqrt{6}}e^t\sin{\sqrt{6}t} & e^t\cos{\sqrt{6}t} \end{pmatrix} \)

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

Неоднородное дифференциальное уравнение:

\( \left\{ \begin{array}{} x''(t) + x(t) = 2\cos{t} \\ x(0) = 0, \, x'(0) = -1 \end{array} \right. \\ x(t) \supset X(p) \Rightarrow p^2X + 1 + X = 2\frac{p}{p^2+1} \\ X(p^2+1) = \frac{2p}{p^2+1}-1 \\ X(p) = \frac{2p}{(p^2+1)^2} - \frac{1}{(p^2+1)} \subset x(t) = t\sin{t} - \sin{t} \)

Система линейных неоднородных уравнений:

\( \left\{ \begin{array}{} x'(t) -y'(t) - 2x +2y = 1-2t \\ x''(t) + 2y'(t) + x = 0\\ x(0) = y(0) = x'(0) = 0 \end{array} \right. \\ \text{Преобразуем оба уравнение по Лапласу:} \\ \left\{ \begin{array}{} pX - pY -2(X-Y) = \frac{1}{p}-\frac{2}{p^2} \\ p^2X + 2pY + X = 0 \end{array} \right. \Leftrightarrow \left\{ \begin{array}{} X-Y = \frac{1}{p^2} \\ p^2X + 2pY + X = 0 \end{array} \right. \\ p^2(Y+\frac{1}{p^2}) + 2pY + Y + \frac{1}{p^2} = 0 \\ Y(p) = -\frac{1}{(p+1)^2} - \frac{1}{p^2(p+1)} = -\frac{2}{(p+1)^2} - \frac{1}{p^2} + \frac{2}{p} - \frac{2}{p+1} \\ X(p) = -\frac{2}{(p+1)^2} + \frac{2}{p} - \frac{2}{p+1} \\ \text{И получим: } \\ X(p) \subset x(t) = -2te^{-t} + 2 - 2e^{-t} \\ Y(p) \subset y(t) = -2te^{-t} - t + 2 - 2e^{-t} \)

Расчеты в электрических цепях:

Пусть силе тока соответствует изображение \(i(t) \supset I(p)\).

  • Мгновенные значения напряжения и тока связаны законом Ома: \(u(t)=Ri(t) \Rightarrow \)

Изображение напряжения выглядит следующим образом: \(u(t) \supset U(t) = RI(p).\)

  • Мгновенные значения напряжения и тока индуктивности связаны как: \(u_L(t) = L\frac{di(t)}{dt} \Rightarrow \)

Воспользовавшись свойством дифференцируемости, найдем изображение напряжения: \(u_L(t) \supset U_L(t) = pLI(p) - Li(0).\)

  • Мгновенные значения напряжения и тока емкости связаны как: \(u_С(t) = \frac1C\int_0^t i(\tau)d\tau \Rightarrow \)

Воспользовавшись свойством интегрируемости, найдем изображение напряжения: \(u_C(t) \supset U_C(t) = \frac{I}{Cp}.\)
ЭДС:

  • постоянный ток: \(e(t) = e_0 \supset E(p) = \frac{e_0}{p} \)
  • переменный ток: \(e(t) = e_0\sin\omega t \supset E(p) = \frac{e_0 \omega}{p^2+\omega^2} \)
Пример:
Последовательный колебательный контур.

Пусть имеется цепь, показанная на схеме справа, и известно, что \(i(0) = 0.\)
Тогда по второму закону Киргофа имеем:
\[L\frac{di(t)}{dt} + Ri(t) + \frac1C\int_0^t i(\tau)d\tau = e(t).\]
Применяя к каждому слагаемому соответствующее преобразование, получим:
\[pLI(p) + RI(p) + \frac{I(p)}{Cp} = E(p).\]
\(Z(p) = pL + R + \frac{1}{Cp} -\) импеданс
\(Y(p) = \frac{1}{Z(p)} -\) адмитанс