Приложения преобразования Лапласа: различия между версиями
Miron1 (обсуждение | вклад) |
Miron1 (обсуждение | вклад) |
||
(не показаны 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} & | + | 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} | \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)} -\) адмитанс