Приложения преобразования Лапласа: различия между версиями

Материал из sawiki
Перейти к навигации Перейти к поиску
 
(не показаны 32 промежуточные версии 2 участников)
Строка 1: Строка 1:
Приведем основные формулы и свойтсва, используемые для решения дифференциальных уравнений и вычисления матричных экспоненциалов с помощью преобразования Лапласа:
+
Приведем основные формулы и свойтсва, используемые для решения дифференциальных уравнений и вычисления матричных экспоненциалов с помощью [https://sawiki.cs.msu.ru/index.php/Преобразование_Лапласа преобразования Лапласа]:
===Прямое и обратное преобразование Лапласа некоторых функций:===  
+
==Прямое и обратное преобразование Лапласа некоторых функций:==  
 
* $$x(t) \supset \frac{1}{p}$$
 
* $$x(t) \supset \frac{1}{p}$$
 
* $$x(t)t^n \supset \frac{n!}{p^{n+1}}$$
 
* $$x(t)t^n \supset \frac{n!}{p^{n+1}}$$
Строка 14: Строка 14:
 
* $$x(t)\operatorname{ch}{\omega t} \supset \frac{p}{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{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}$$
+
* $$x(t)e^{\beta t}\operatorname{ch}{\omega t} \supset \frac{p-\beta}{(p-\beta)^2 - \omega^2}$$, <br> $$\text{где } x(t)$$ - [https://ru.wikipedia.org/wiki/Функция_Хевисайда функция Хевисайда].
  
====Свойства:====
+
==Свойства:==
 
*$$f(at) \supset \frac{1}{a}F(\frac{p}{a}) \, , \, a > 0$$
 
*$$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}$$
 
*$$\frac{1}{b}\,f(\frac{t}{b}) \supset F(pb) \, , \, b > 0 \,;\, b=\frac{1}{a}$$
Строка 30: Строка 30:
 
*$$\int\limits_0^t f(\tau) d\tau \supset \frac{F(p)}{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$$
 
*$$\frac{f(t)}{t} \supset \int\limits_0^{+\infty}F(z) dz$$
===Примеры:===
+
 
 +
==Примеры:==
 
====Нахождение фундаментальной матрицы:====
 
====Нахождение фундаментальной матрицы:====
 
\(
 
\(
Строка 62: Строка 63:
 
\)
 
\)
 
<br />
 
<br />
От каждого элемента матрицы Y нужно вычислить преобразование Лапласа, и получим:
+
От каждого элемента матрицы Y нужно вычислить [https://sawiki.cs.msu.ru/index.php/Преобразование_Лапласа преобразование Лапласа], и получим:
 
<br />
 
<br />
 
\(
 
\(
 
X(t) = \begin{pmatrix}
 
X(t) = \begin{pmatrix}
e^t\cos{\sqrt{6}t} & -2e^t\sin{\sqrt{6}t} \\
+
e^t\cos{\sqrt{6}t} & \frac{2}{\sqrt{6}}e^t\sin{\sqrt{6}t} \\
-3e^t\sin{\sqrt{6}t} & e^t\cos{\sqrt{6}t}
+
-\frac{3}{\sqrt{6}}e^t\sin{\sqrt{6}t} & e^t\cos{\sqrt{6}t}
 
\end{pmatrix}
 
\end{pmatrix}
 
\)
 
\)
Строка 88: Строка 89:
 
\left\{ \begin{array}{}
 
\left\{ \begin{array}{}
 
x'(t) -y'(t) - 2x +2y = 1-2t \\  
 
x'(t) -y'(t) - 2x +2y = 1-2t \\  
x''(t) + 2y'(t) + x = \\
+
x''(t) + 2y'(t) + x = 0\\
 
x(0) = y(0) = x'(0) = 0
 
x(0) = y(0) = x'(0) = 0
 
\end{array}
 
\end{array}
Строка 110: Строка 111:
 
Y(p) \subset y(t) = -2te^{-t} - t + 2 - 2e^{-t}  
 
Y(p) \subset y(t) = -2te^{-t} - t + 2 - 2e^{-t}  
 
\)
 
\)
 +
 +
====Расчеты в электрических цепях:====
 +
Пусть силе тока соответствует изображение \(i(t) \supset I(p)\).<br>
 +
* Мгновенные значения напряжения и тока связаны законом Ома: \(u(t)=Ri(t) \Rightarrow \) <br>
 +
Изображение напряжения выглядит следующим образом: \(u(t) \supset U(t) = RI(p).\)
 +
* Мгновенные значения напряжения и тока индуктивности связаны как: \(u_L(t) = L\frac{di(t)}{dt} \Rightarrow \) <br>
 +
Воспользовавшись свойством дифференцируемости, найдем изображение напряжения: \(u_L(t) \supset U_L(t) = pLI(p) - Li(0).\)
 +
* Мгновенные значения напряжения и тока емкости связаны как: \(u_С(t) = \frac1C\int_0^t i(\tau)d\tau \Rightarrow \) <br>
 +
Воспользовавшись свойством интегрируемости, найдем изображение напряжения: \(u_C(t) \supset U_C(t) = \frac{I}{Cp}.\) <br>
 +
'''ЭДС''':
 +
* постоянный ток: \(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} \)
 +
=====Пример:=====
 +
[[File:pcc.png|frame|right|Последовательный колебательный контур.]]
 +
Пусть имеется цепь, показанная на схеме справа, и известно, что \(i(0) = 0.\) <br>
 +
Тогда по второму закону Киргофа имеем: <br>
 +
\[L\frac{di(t)}{dt} + Ri(t) + \frac1C\int_0^t i(\tau)d\tau = e(t).\]<br>
 +
Применяя к каждому слагаемому соответствующее преобразование, получим: <br>
 +
\[pLI(p) + RI(p) + \frac{I(p)}{Cp} = E(p).\] <br>
 +
\(Z(p) = pL + R + \frac{1}{Cp} -\)  '''импеданс''' <br>
 +
\(Y(p) = \frac{1}{Z(p)} -\)  '''адмитанс'''

Текущая версия на 21:29, 29 декабря 2020

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

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

  • $$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)} -\) адмитанс