Дискуссионный математический форумМатематический форум

Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 4 часа [ Летнее время ]
MathHelpPlanet.com RSS-лента Математического форума

Часовой пояс: UTC + 4 часа [ Летнее время ]


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

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


Содержание

Основные определения


1. Оригинал — это комплекснозначная функция [math]f(t)[/math] действительного аргумента [math]t[/math], которая удовлетворяет следующим условиям:


а) [math]f(t)=0[/math] при [math]t<0[/math];

б) на любом конечном отрезке [math][a;b]\in[0;+\infty)[/math] функция [math]f(t)[/math] имеет не более чем конечное число точек разрыва первого рода;

в) [math]f(t)f(t)[/math] имеет ограниченный рост, т.е. возрастает не быстрее показательной функции: существуют такие постоянные [math]M>0[/math] и [math]\sigma\geqslant0[/math], что [math]|f(t)|<M\,e^{\sigma t}[/math] при [math]t>0[/math].


Замечания 5.1


1. Величина [math]\sigma_0=\inf\sigma[/math] называется показателем роста функции [math]f(t)[/math]. Для любой ограниченной функции, являющейся оригиналом, можно принять [math]\sigma_0=0[/math].

2. Обозначим [math]f(+0)= \lim\limits_{t\to+0} f(t),~ f(+\infty)= \lim\limits_{t\to+\infty} f(t)[/math], если пределы существуют и конечны.

3. Совокупность всех оригиналов называется пространством оригиналов.

4. В точке [math]t_0[/math] разрыва первого рода функция имеет конечные односторонние пределы: [math]\lim\limits_{t\to t_0+0} f(t),~ \lim\limits_{t\to t_0-0} f(t)[/math].


▼ Пример 5.1

Проверить, являются ли функции оригиналами:


[math]f_1(t)= \begin{cases}2e^{5t},& t \geqslant 0,\\ 0,&t<0;\end{cases} \qquad f_2(t)= \begin{cases}\dfrac{1}{t-2},& t \geqslant 0,\\ 0,&t<0;\end{cases} \qquad f_3(t)= \begin{cases}3^{4^t},& t \geqslant 0,\\ 0,&t<0;\end{cases}[/math]

Решение. Функция [math]f_1(t)[/math] является оригиналом, так как условия пп. "а"-"в" выполнены: [math]M=2,~ \sigma_0=5[/math]; функция [math]f_2(t)[/math] не является оригиналом, так как в точке [math]t=2[/math] имеет разрыв второго рода (не выполняется условие "б"); функция [math]f_3(t)[/math] не является оригиналом, так как растет быстрее показательной функции (не выполняется условие "в", поскольку [math]3^{4^t}>M\,e^{\sigma t}[/math] для любых [math]M[/math] и [math]\sigma,~t>0[/math]).


2. Изображение функции [math]f(t)[/math] — функция [math]F(p)[/math] комплексного переменного [math]p[/math], определяемая равенством


[math]F(p)= \int\limits_{0}^{+\infty} e^{-pt}f(t)\,dxt\,.[/math]
(5.1)

Область существования этой функции определяется областью сходимости интеграла Лапласа, стоящего в правой части равенства (5.1). Исследование интeгpaлa позволяет определить эту область и установить свойства функции [math]F(p)[/math]. Имеет место следующее утверждение.


Утверждение 5.1. Если функция [math]f(t)[/math], является оригиналом, то интеграл Лапласа сходится абсолютно в области [math]\operatorname{Re}p=\sigma> \sigma_0[/math] (рис. 5.1) , где [math]\sigma_0[/math] — показатель роста оригинала. Внутри этой области, т.е. на любом замкнутом подмножестве [math]\operatorname{Re}p=\sigma \geqslant a>\sigma_0[/math], интеграл сходится равномерно и определяет аналитическую функцию [math]F(p)[/math].


Замечания 5.2


1. Утверждение 5.1 аналогично свойствам степенных рядов, сходящихся в круге и равномерно сходящихся внутри этого круга, где сумма ряда является аналитической функцией.

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

3. Совокупность всех изображений [math]F(p)[/math] называется пространством изображений.

4. Переход, определяющий изображение [math]F(p)[/math] по оригиналу [math]f(t)[/math], называется преобразованием Лапласа:

[math]F(p)= L\bigl[f(t)\bigr]= \int\limits_{0}^{+\infty} e^{-pt}f(t)\,dt\,.[/math]
(5.2)

Запись [math]F(p)=L[f(t)][/math] означает, что оригиналу [math]f(t)[/math] соответствует изображение [math]F(p)[/math].

5. Оригинал по изображению находится с помощью обратного преобразования Лапласа по формуле обращения

[math]f(t)=L^{-1}\bigl[F(p)\bigr]= \frac{1}{2\pi i} \int\limits_{\sigma-i\infty}^{\sigma+i\infty} e^{pt}F(p)\,dp\,,[/math]
(5.3)

где путь интегрирования — любая прямая [math]\operatorname{Re}p=\sigma[/math], параллельная мнимой оси и лежащая правее прямой [math]\operatorname{Re}p=\sigma_0[/math] (рис. 5.1).

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


Замечания 5.3


1. Для преобразования Лапласа используются различные обозначения, на пример [math]f(t)\risingdotseq F(p)[/math] и [math]F(p)\fallingdotseq f(t)[/math], что означает: оригиналу [math]f(t)[/math] соответствует изображение [math]F(p)[/math] и изображению [math]F(p)[/math] соответствует оригинал [math]f(t)[/math]. В некоторых учебниках вместо аргумента [math]p[/math] применяется [math]s[/math], то есть [math]F(s)=L[f(t)][/math] и [math]L^{-1}[F(s)]=f(t)[/math].


2. Для компактной записи оригиналов используется единичная ступенчатая функция [math]\boldsymbol{1}(t-\tau)\colon[/math]


[math]\boldsymbol{1}(t-\tau)= \begin{cases}1,& t>\tau,\\ 0,& t \leqslant \tau,\end{cases}[/math]
(5.4)

где [math]\tau[/math] — точка приложения (рис. 5.2). Так как во многих практических задачах аргумент [math]t[/math] имеет смысл текущего времени, то [math]\tau[/math] также называется моментом приложения единичной ступенчатой функции. В системах автоматического регулирования и управления функция [math]\boldsymbol{1}(t-\tau)[/math] рассматривается как типовой входной сигнал.


При [math]\tau=0[/math] функция [math]\boldsymbol{1}(t-\tau)[/math] является функцией Хевисайда:


[math]\boldsymbol{1}(t)= \begin{cases}1,& t>0,\\ 0,& t \leqslant 0.\end{cases}[/math]
(5.5)

Тогда, если функция [math]f(t)[/math] удовлетворяет условиям "б", "в" в определении оригинала (п. 1), но не удовлетворяет условию "а", то функция [math]f(t)\cdot \boldsymbol{1}(t)[/math] будет оригиналом, так как


[math]f(t)\cdot\boldsymbol{1}(t)= \begin{cases}f(t),& t>0,\\ 0,& t\leqslant0. \end{cases}[/math]

Далее под заданной с помощью аналитической формулы функцией [math]f(t)[/math], там, где это не вызывает недоразумений, будем понимать произведение этой функции на функцию Хевисайда, а множитель [math]\boldsymbol{1}(t)[/math] опускать.


3. Функции [math]F(p)[/math], являющиеся изображениями, удовлетворяют необходимому условию: если [math]F(p)[/math] есть изображение, то [math]F(p)\to0[/math] при [math]\operatorname{Re}p=\sigma\to+\infty[/math]. Поэтому функции [math]F_1(p)=1,~ F_2(p)=p,~ F_3(p)=\sin p,~ F_4(p)= \frac{p}{p-1}[/math] не являются изображениями. Однако в практических задачах функции типа [math]F_1(p)=1,~ F_2(p)=p[/math] и другие встречаются. Это требует расширения понятий оригинала и изображения.


Класс оригиналов можно расширить, включив в него функции, которые могут быть не ограничены в окрестности некоторых конечных точек, но такие, что интеграл Лапласа от них, тем не менее, сходится абсолютно в некоторой полуплоскости [math]\operatorname{Re}p>\sigma_0[/math]. К числу таких обобщенных оригиналов относятся степенная функция [math]f(t)=t^{\mu}[/math] при [math]\mu>-1,~\ln t[/math] и некоторые другие.


4. Во всякой точке [math]t_0[/math], являющейся точкой разрыва функции [math]f(t)[/math], правая часть формулы (5.3) равна [math]\frac{1}{2}\bigl[f(t_0-0)+f(t_0+0)\bigr][/math].


▼ Примеры 5.2-5.3

Пример 5.2. Найти изображение единичной функции Хевисайда [math]f(t)= \boldsymbol{1}(t)[/math].


Решение. Так как функция [math]\boldsymbol{1}(t)[/math] ограничена, то в качестве показателя роста можно положить [math]\sigma_0=0[/math]. По формуле (5.2) имеем


[math]F(p)= L[\boldsymbol{1}(t)]= \int\limits_{0}^{+\infty} \boldsymbol{1}(t) e^{-pt}\,dt=\left.{-\frac{1}{p}\,e^{-pt}}\right|_{0}^{+\infty}= \frac{1}{p}\,.[/math]

так как из равенства [math]|e^{-pt}|= e^{\operatorname{Re}(-pt)}= e^{-t \operatorname{Re}p}= e^{-\sigma t}[/math] следует, что [math]\lim\limits_{t\to+\infty} e^{-pt}=0[/math] при [math]\operatorname{Re}p= \sigma>\sigma_0=0[/math].


Пример 5.3. Найти изображение функции [math]f(t)=e^{at}[/math], где [math]a[/math] — действительное число.


Решение. Показателем роста можно считать [math]\sigma_0=a[/math]. По формуле (5.2)


[math]F(p)= \int\limits_{0}^{+\infty} e^{at}e^{-pt}\,dt= \int\limits_{0}^{+\infty} e^{(a-p)t}\,dt= \left.{\frac{e^{(a-p)t}}{a-p}}\right|_{0}^{+\infty}=-\frac{1}{a-p}=\frac{1}{p-a}\,.[/math]

так как из равенства [math]|e^{(a-p)t}|= e^{-t(\operatorname{Re}p-a)}= e^{-t(\sigma-a)}[/math] следует, что [math]\lim\limits_{t\to+\infty} e^{(a-p)t}=0[/math] при [math]\operatorname{Re}p=\sigma>a[/math].




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


Будем предполагать, что рассматриваемые далее функции [math]f(t),f_1(t),\ldots,f_n(t)[/math] являются оригиналами. Соответствующие им изображения (при [math]\operatorname{Re}p> \sigma_i,~ i=0,1,\ldots,n[/math]) обозначим [math]F(p),F_1(p),\ldots,F_n(p)[/math].


1. Линейность. Если [math]f_1(t),\ldots,f_n(t)[/math] — оригиналы, то для любых комплексных чисел [math]c_i,~ i=1,2,\ldots,n[/math], функция [math]\textstyle{f(t)= \sum\limits_{k=1}^{n} c_kf_k(t)}[/math] также является оригиналом и справедливо равенство [math]\textstyle{L\!\left[ \sum\limits_{k=1}^{n} c_kf_k(t)\right]= \sum\limits_{k=1}^{n}c_k L[f_k(t)]}[/math]


[math]L \bigl[c_1f_1(t)+\ldots+ c_nf_n(t)\bigr]= c_1F_1(p)+\ldots+ c_nF_n(p),\quad \operatorname{Re}p>\max\{\sigma_1,\ldots,\sigma_n\}.[/math]
(5.6)

Заметим, что для функции [math]\textstyle{f(t)= \sum\limits_{k=1}^{n} c_kf_k(t)}[/math] существенно, что все слагаемыс являются оригиналами, так как, например, функция [math]f(t)=\frac{e^{t}-1}{t}[/math] является оригиналом, а слагаемые [math]f_1(t)=\frac{e^t}{t}[/math] и [math]f_2(t)=-\frac{1}{t}[/math] не являются.


Справедливо и обратное утверждение: если [math]F_1(p),\ldots,F_n(p)[/math] — изображения, то


[math]L^{-1} \bigl[c_1F_1(p)+\ldots+ c_nF_n(p)\bigr]= c_1f_1(t)+\ldots+ c_nf_n(t).[/math]

Здесь также важно, что слагаемые функции [math]c_1F_1(p)+\ldots+ c_nF_n(p)[/math] являются изображениями, поскольку из того, что [math]F(p)= c_1F_1(p)+\ldots+ c_nF_n(p)[/math] — изображение, не следует, что [math]F_1(p)+\ldots+ F_n(p)[/math] — изображения. Например, функция [math]F(p)= \ln \frac{p-1}{p}[/math] является изображением, а слагаемые [math]F_1(p)= \ln(p-1)[/math] и [math]F_2(p)=-\ln p[/math] не являются.


▼ Примеры 5.4-5.6

Пример 5.4. Найти изображение функции [math]f(t)=3+2e^{-t}[/math].


Решение. Из примера 5.2 имеем [math]L[\boldsymbol{1}(t)]= \frac{1}{p}[/math], а из примера 5.3 при [math]a=-1[/math] имеем [math]L[e^{-t}]= \frac{1}{p+1}[/math]. Тогда согласно свойству линейности для оригинала [math]f(t)=3\cdot\boldsymbol{1}(t)+2\cdot e^{-t}[/math] получаем [math]F(p)= 3\cdot \frac{1}{p}+2\cdot \frac{1}{p+1}[/math].


Пример 5.5. Найти изображение функции [math]f(t)=\cos t[/math].


Решение. Используя формулу Эйлера (2.11), получаем [math]f(t)=\cos t=\frac{e^{it}+e^{-it}}{2}= \frac{1}{2}\,e^{it}+\frac{1}{2}\,e^{-it}[/math].


Из примера 5.3 при [math]a=i[/math] и [math]a=-i[/math] следует: [math]L[e^{it}]= \frac{1}{p-i}\,,~ L[e^{-it}]=\frac{1}{p+i}[/math]. Тогда по свойству линейности


[math]L[\cos t]= \frac{1}{2}L[e^{it}]+\frac{1}{2}L[e^{-it}]= \frac{1}{2}\cdot \frac{1}{p-i}+ \frac{1}{2}\cdot \frac{1}{p+i}= \frac{p+i+p-i}{2(p^2+1)}= \frac{p}{p^2+1}\,.[/math]

2. Подобие (теорема подобия). Для любого [math]a>0[/math] из [math]F(p)=L[f(t)][/math] следует


[math]L \bigl[f(at)\bigr]= \frac{1}{a}\,F\! \left(\frac{p}{a}\right)\!,\quad \operatorname{Re}p>a\sigma_0,[/math]
(5.7)

и обратно: [math]L^{-1} \bigl[F(ap)\bigr]= \frac{1}{a}\,f\! \left(\frac{t}{a}\right)[/math].


▼ Пример 5.6

Найти изображение функции [math]f(t)=\cos at[/math].


Решение. Из примера 5.5 следует, что [math]L[\cos t]= \frac{p}{p^2+1}[/math]. Тогда по теореме подобия


[math]L[\cos at]= \frac{1}{a}\cdot \frac{p\!\!\not{\phantom{|}}\,a}{(p\!\!\not{\phantom{|}}\, a)^2+1}= \frac{p}{p^2+a^2}\,.[/math]

3. Смещение (теорема смещения). При любом комплексном [math]a[/math] из [math]F(p)= L[f(t)][/math] следует


[math]L \bigl[e^{at}f(t)\bigr]= F(p-a),\quad \operatorname{Re}(p-a)>\sigma_0,[/math]
(5.8)

то есть умножению оригинала на [math]e^{at}[/math] соответствует смещение изображения на [math]a[/math].


▼ Пример 5.7

Найти изображение функции [math]f(t)=e^{at}\cos bt[/math].


Решение. Из примера 5.6 следует [math]L[\cos bt]=\frac{p}{p^2+a^2}[/math]. Тогда по теореме смещения [math]L[e^{at}\cos bt]= \frac{p-a}{(p-a)^2+b^2}[/math].




Запаздывание оригинала


4. Запаздывание (теорема запаздывания). Для любого [math]\tau>0[/math] из [math]F(p)=L[f(t)][/math] следует


[math]L\bigl[f(t-\tau)\bigr]= e^{-p\tau}\cdot F(p),\quad \operatorname{Re}p>\sigma_0,[/math]
(5.9)

где [math]f(t-\tau)=f(t-\tau)\cdot \boldsymbol{1}(t-\tau)[/math] (рис. 5.3), т.е. запаздыванию оригинала на [math]\tau>0[/math] соответствует умножение изображения на [math]e^{-p\tau}[/math].


▼ Примеры 5.8-5.10

Пример 5.8. Найти изображение функции [math]f(t)=\cos(t-3)\cdot \boldsymbol{1}(t-3)= \cos(t-3)[/math].


Решение. В примере 5.5 получено [math]L[\cos t]=\frac{p}{p^2+1}[/math]. По теореме запаздывания при [math]\tau=3[/math] имеем [math]L[\cos(t-3)]= \frac{e^{-3p}\cdot p}{p^2+1}[/math].


Пример 5.9. Найти оригиналы по изображениям: a) [math]F(p)=\frac{(p-1)e^{-4p}}{(p-1)^2+4}[/math]; б) [math]F(p)=\frac{e^{-p}}{p^2}[/math].


▼ Решение

а) Из примера 5.7 следует, что при [math]a=1,~b=2[/math] изображению [math]\frac{p-1}{(p-1)^2+4}[/math] соответствует оригинал [math]e^t\cos2t[/math]. Тогда по теореме запаздывания при [math]\tau=4[/math] имеем


[math]f(t)= L^{-1}[F(p)]= e^{t-4}\cos(t-4)\cdot \boldsymbol{1}(t-4).[/math]

б) По формуле 3 из табл. 5.1 [math]L^{-1}\!\left[\frac{1}{p^2}\right]=t[/math].По теореме запаздывания при [math]\tau=1[/math] получаем [math]f(t)=L^{-1}[F(p)]= (t-1)\cdot \boldsymbol{1} (t-1)[/math]. Заметим, что для похожего, но отличного от полученного, оригинала [math]f(t)=t-1[/math] (его можно записать в виде [math](t-1)\cdot \boldsymbol{1}(t)[/math]) изображение имеет вид [math]F(p)=\frac{1}{p^2}-\frac{1}{p}\ne \frac{e^{-p}}{p^2}[/math].


Пример 5.10. Найти изображение функции [math]\delta_h(t)= \begin{cases}\frac{1}{h},& 0<t\leqslant h,\\ 0,& t<0,\,t>h,\end{cases}[/math], график которой представлен на рис. 5.4.


▼ Решение

С учетом (5.4) представим функцию [math]\delta_h(t)[/math] в виде [math]\delta_h(t)= \frac{\boldsymbol{1}(t)-\boldsymbol{1}(t-h)}{h}[/math].


Из примера 5.2 имеем [math]L[\boldsymbol{1}(t)]=\frac{1}{p}[/math]. Применяя свойства линейности и запаздывания, получаем


[math]L\bigl[\delta_h(t)\bigr]= \frac{1}{h}L\bigl[\boldsymbol{1}(t)\bigr]-\frac{1}{h}L\bigl[\boldsymbol{1}(t-h)\bigr]= \frac{1}{h}\cdot \frac{1}{p}-\frac{1}{h}\cdot \frac{1}{p}\,e^{-ph}= \frac{1-e^{-ph}}{ph}\,.[/math]

Заметим, что, находя предел при [math]h\to0[/math] в последнем выражении, можно получить изображение δ-функции [math]\delta(t)= \lim\limits_{h\to0} \delta_h(t)\colon[/math]


[math]L\bigl[\delta(t)\bigr]= \lim\limits_{h\to0} \frac{1-e^{-ph}}{ph}= \lim\limits_{h\to0} \frac{p\,e^{-ph}}{p}=1.[/math]

Замечание 5.4. Дельта-функция часто встречается в инженерных приложениях как идеализация импульса конечной длительности. В теории автоматического регулирования и управления δ-функция вместе с единичной ступенчатой являются типовыми входными воздействиями.


Очевидно, изображение дельта-функции не удовлетворяет необходимому условию (п.2 замечаний 5.3). Этот факт свидетельствует о практическом требовании расширения понятия оригинала. Дельта-функция относится к обобщенным функциям и задается соотношением


[math]\int\limits_{a}^{b} f(t)\delta(t-\tau)\,dt= \begin{cases}f(\tau+0),& a \leqslant \tau<b,\\ 0,& \tau<a,\, \tau\geqslant b.\end{cases}[/math]
(5.10)



Дифференцирование оригинала


5. Если функции [math]f(t),f'(t),\ldots,f^{(n)}(t)[/math] являются оригиналами и [math]F(p)=L[f(t)][/math], то


[math]\begin{aligned}& L[f'(t)]= pF(p)-f(+0),\\ & L[f''(t)]= p^2F(p)-pf(+0)-f'(+0),\\ & \quad\vdots\\ & L[f^{(n)}(t)]= p^nF(p)-p^{n-1}f(+0)-\ldots-f^{(n-1)}(+0), \end{aligned}[/math]
(5.11)

где [math]f^{i}(+0)= \lim\limits_{t\to+0} f^{(i)}(t),~ i=0,1,2,\ldots,n-1[/math].


▼ Примеры 5.11-5.12

Пример 5.11. Найти изображение [math]f'(t)[/math], если [math]f(t)=e^{-t}\cos3t[/math].


Решение. Из примера 5.7 следует, что при [math]a=-1,~b=3[/math] имеем [math]L[e^{-t}\cos3t]= \frac{p+1}{(p+1)^2+9}[/math].


Найдем [math]f(+0)=\lim\limits_{t\to+0}e^{-t}\cos3t=1[/math]. Согласно (5.11) [math]L[f'(t)]= p\cdot \frac{p+1}{(p+1)^2+9}-1[/math].


Пример 5.12. Найти изображение выражения [math]x''+3x'+2x+1[/math] с начальными условиями [math]x(+0)=1,~x'(+0)=4[/math].


Решение. Пусть [math]X(p)=L[x(t)][/math], тогда [math]L[x'(t)]= pX(p)-1;~ L[x''(t)]= p^2X(p)-p\cdot1-4[/math]. В примере 5.2 получено [math]L[\boldsymbol{1}(t)]=1\!\!\not{\phantom{|}}\,p[/math]. Используя свойство линейности, имеем


[math]\begin{aligned}L[x''+3x'+2x+1]&= L[x'']+3L[x']+2L[x]+L[1]=\\ &=p^2X(p)-p-4+3pX(p)-3+2X(p)+\frac{1}{p}=\\ &=(p^2+3p+2)X(p)-p-7+\frac{1}{p}\,.\end{aligned}[/math]

Интегрирование оригинала


Если функция [math]f(t)[/math] является оригиналом и [math]F(p)=L[f(t)][/math], то


[math]\textstyle{L\! \left[\int\limits_{0}^{t} f(\tau)\,d\tau\right]= \dfrac{F(p)}{p},\quad \operatorname{Re}p>\sigma_0,}[/math]
(5.12)

т.е. интегрированию оригинала соответствует деление изображения на [math]p[/math].


▼ Пример 5.13

Найти изображение интеграла [math]\textstyle{\int\limits_{0}^{t} f(\tau)\,d\tau}[/math] от функции [math]f(t)=\cos t[/math].


Решение. Из примера 5.6 следует, что [math]L[\cos t]=\frac{p}{p^2+1}=F(p)[/math]. Тогда


[math]L\Biggl[\int\limits_{0}^{t} \cos\tau\,d\tau\Biggr]= L[\sin t]= \frac{p}{p(p^2+1)}= \frac{1}{p^2+1}[/math], то есть [math]L[\sin t]= \frac{1}{p^2+1}[/math].



Дифференцирование изображения


Если функция [math]f(t)[/math] является оригиналом и [math]F(p)=L[f(t)][/math]. то


[math]L\bigl[(-1)^nt^nf(t)\bigr]= F^{(n)}(p).[/math]
(5.13)

▼ Примеры 5.14-5.15

Пример 5.14. Найти изображения функций [math]f_1(t)=t, f_2(t)=t^2,\ldots,f_n(t)=t^n[/math].


Решение. Из примера 5.2 следует, что [math]L[\boldsymbol{1}(t)]=\frac{1}{p}=F(p)[/math]. Согласно (5.13) при [math]n=1[/math] получаем [math]L[-t]=F'(p)=-\frac{1}{p^2}[/math] или по свойству линейности [math]L[t]=\frac{1}{p^2}[/math]. При [math]n=2[/math] имеем [math]L[t^2]=F''(p)= \frac{2}{p^3}[/math]. Аналогично [math]L[t^n]=\frac{n!}{p^{n+1}}[/math].


Пример 5.15. Найти изображение функции [math]f(t)=t\cos3t[/math].


Решение. Из примера 5.6 имеем [math]L[\cos3t]= \frac{p}{p^2+9}=F(p)[/math]. При [math]n=1[/math] из (5.13) вытекает [math]L[(-t)f(t)]=F''(p)[/math] или


[math]L[-t\cos3t]= \left(\frac{p}{p^2+9}\right)'= \frac{p^2+9-2p^3}{(p^2+9)^2}=-\frac{p^2-9}{(p^2+9)^2}\,.[/math]

Окончательно получаем [math]L[t\cos3t]=\frac{p^2-9}{(p^2+9)^2}[/math].


Интегрирование изображения


Если функция [math]\frac{f(t)}{t}[/math] является оригиналом, то из [math]F(p)=L[f(t)][/math] следует


[math]L\! \left[\frac{f(t)}{t}\right]= \int\limits_{p}^{\infty} F(z)\,dz\,.[/math]
(5.14)

▼ Пример 5.16

Найти изображение функции [math]\frac{\sin t}{t}[/math].


Решение. Функция [math]\frac{\sin t}{t}[/math] является оригиналом, так как [math]\left|\frac{\sin t}{t}\right|<1[/math] (условие "в") и точка [math]t=0[/math] является точкой разрыва первого рода (условие "б"). Из примера 5.13 следует [math]L[\sin t]=\frac{1}{p^2+1}=F(p)[/math].


Отсюда [math]L\! \left[\frac{\sin t}{t}\right]= \int\limits_{p}^{+\infty} \frac{dz}{z^2+1}= \Bigl.{\operatorname{arctg}z}\Bigr|_{p}^{+\infty}=\frac{\pi}{2}-\operatorname{arctg}p.[/math].




Умножение изображений (теорема Бореля)


Из [math]F_1(p)=L[f_1(t)][/math] и [math]F_2(p)=L[f_2(t)][/math] следует


[math]L\bigl[f_1(t)\ast f_2(t)\bigr]= F_1(p)\cdot F_2(p),[/math]
(5.15)

т.е. свертке оригиналов соответствует произведение изображений. Функция [math]f_1(t)\ast f_2(t)[/math] определяется формулой


[math]f_1(t)\ast f_2(t)= \int\limits_{0}^{t} f_1(\tau)f_2(t-\tau)\,d\tau= \int\limits_{0}^{t} f_1(t-\tau)f_2(\tau)\,d\tau[/math]
(5.16)

и называется сверткой оригиналов [math]f_1(t)[/math] и [math]f_2(t)[/math].


▼ Пример 5.17

Найти оригинал, соответствующий изображению [math]F(p)= \frac{p}{(p^2+1)^2}[/math].


Решение. Представим [math]F(p)[/math] в виде произведения изображений: [math]F(p)= F_1(p)\cdot F_2(p)= \frac{1}{p^2+2}\cdot \frac{p}{p^2+1}[/math].


Из примеров 5.6 и 5.13 следует [math]f_1(t)=L^{-1}\! \left[\frac{1}{p^2+2}\right]=\sin t,~ f_2(t)=L^{-1}\! \left[\frac{p}{p^2+2}\right]=\cos t[/math].


Согласно (5.15),(5.1б) получаем искомый оригинал:


[math]\begin{aligned}L^{-1}\! \left[\frac{p}{(p^2+1)^2}\right]&= f_1(t)\ast f_2(t)= \int\limits_{0}^{t} \sin\tau\cos(t-\tau)\,\tau= \frac{1}{2} \int\limits_{0}^{t} \bigl[\sin t+\sin(2\tau-t)\bigr]d\tau=\\ \left.{\frac{1}{2}\! \left(\sin t\cdot\tau-\frac{1}{2}\cos(2\tau-t)\right) }\right|_{0}^{t}= \frac{1}{2}\! \left(t\sin t-\frac{1}{2}\cos t+\frac{1}{2}\cos(-t)\right)= \frac{t}{2}\sin t.\end{aligned}[/math]



Дифференцирование свертки (интеграл Дюамеля)


10. Согласно свойствам 9 и 5 найдем преобразование Лапласа от производной свертки двух функций:


[math]L\left\{\bigl[f_1(t)\ast f_2(t)\bigr]'\right\}= pF_1(p)\cdot F_2(p).[/math]

С другой стороны,

[math]\bigl[f_1(t)\ast f_2(t)\bigr]'= \frac{d}{dt} \int\limits_{0}^{t} f_1(\tau) f_2(t-\tau)\,d\tau= f_2(0)f_1(t)+ \int\limits_{0}^{t} f'_2(\tau) f_1(t-\tau)\,d\tau[/math]

или, применяя правило дифференцирования интеграла, зависящего от параметра, имеем


[math]\bigl[f_1(t)\ast f_2(t)\bigr]'= \frac{d}{dt} \int\limits_{0}^{t} f_1(t-\tau)f_2(\tau)\,d\tau= f_1(0)f_2(t)+\int\limits_{0}^{t} f'_1(\tau)f_2(t-\tau)\,d\tau\,.[/math]

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


[math]\frac{d}{d\lambda} \int\limits_{u(\lambda)}^{v(\lambda)}f(x,\lambda)\,dx= f \bigl(v(\lambda),\lambda\bigr)\frac{dv}{d\lambda}-f \bigl(u(\lambda),\lambda\bigr)\frac{du}{d\lambda}+ \int\limits_{u(\lambda)}^{v(\lambda)} \frac{\partial}{\partial\lambda}f(x,\lambda)\,dx\,.[/math]

Объединяя полученные результаты, можно записать:


[math]\begin{aligned}L^{-1} \bigl[pF_1(p)F_2(p)\bigr]&= f_1(0)f_2(0)+ f'_1\ast f_2= f_2(0)f_1(t)+ f'_2\ast f_1=\\ &=f_1(0)f_2(0)+ \int\limits_{0}^{t}f'_1(\tau)f_2(t-\tau)\,d\tau= f_2(0)f_1(t)+ \int\limits_{0}^{t} f'_2(\tau) f_1(t-\tau)\,d\tau\,.\end{aligned}[/math]
(5.17)

Формула (5.17) называется интегралом Дюамеля. Интеграл Дюамеля применяется для решения дифференциальных уравнений.


▼ Пример 5.18

Найти оригиналы, соответствующие изображениям: a) [math]F(p)= \frac{p}{(p-1)(p-2)}[/math]; б) [math]F(p)= \frac{p^3}{(p^2+1)(p^2+4)}[/math].


Решение. а) Заметим, что здесь нельзя непосредственно воспользоваться теоре мой Бореля, так как в произведении [math]F(p)=\frac{p}{p-1}\cdot \frac{1}{p-2}[/math] множитель [math]F_1(p)= \frac{p}{p-1}[/math] не является изображением (не выполняется необходимое условие).


Представим изображение в виде произведения [math]F(p)=p\cdot\frac{1}{p-1}\cdot \frac{1}{p-2}[/math]. Из примера 5.3 вытекает


[math]f_1(t)=L^{-1}\bigl[F_1(p)\bigr]= L^{-1}\!\left[\frac{1}{p-1}\right]= e^t;\qquad f_2(t)= L^{-1} \bigl[F_2(p)\bigr]= L^{-1}\! \left[\frac{1}{p-2}\right]=e^{2t}.[/math]

Тогда по формуле (5.17) имеем


[math]L^{-1} \bigl[F(p)\bigr]= f_1(0)f_2(t)+ f'_1\ast f_2= 1\cdot e^{2t}+ \int\limits_{0}^{t} e^{\tau}\cdot e^{2(t-\tau)}\,d\tau= e^{2t}+ e^{2t} \int\limits_{0}^{t} e^{-\tau}\,d\tau=-e^{t}+2e^{2t}.[/math]

Можно решить этот пример с помощью теоремы Бореля, представив изображение в виде


[math]F(p)= \frac{p}{(p-1)(p-2)}= \frac{p-1+1}{(p-1)(p-2)}= \frac{1}{p-2}+ \frac{1}{p-1}\cdot \frac{1}{p-2}\,.[/math]

Тогда, используя свойство линейности и теорему Бореля, получаем


[math]L^{-1}\bigl[F(p)\bigr]&= e^{2t}+e^{t}\ast e^{2t}= e^{2t}+ \int\limits_{0}^{t} e^{2\tau}\cdot e^{t-\tau}\,d\tau= e^{2t}+ e^{t} \int\limits_{0}^{t} e^{\tau}\,d\tau= e^{2t}+e^{t}(e^{t}-1)= 2e^{2t}-e^{t}.[/math]

б) Представим изображение в виде произведения:


[math]F(p)= \frac{p^3}{(p^2+1)(p^2+4)}= p\cdot \frac{p}{p^2+1}\cdot \frac{p}{p^2+1}= p\cdot F_1(p)\cdot F_2(p).[/math]

Из примера 5.6 при [math]a=1[/math] и [math]a=2[/math] следует


[math]L^{-1}\! \left[\frac{p}{p^2+1}\right]=\cos t=f_1(t),\qquad L^{-1}\! \left[\frac{p}{p^2+ 4}\right]= \cos2t= f_2(t).[/math]

Тогда по формуле (5.17) получаем


[math]\begin{aligned}f(t)&= L^{-1} \bigl[F(p)\bigr]= f_1(0)f_2(t)+ f'_1\ast f_2= \cos2t-\int\limits_{0}^{t} \sin\tau\cos[2(t-\tau)]\,d\tau=\\ &=\cos2t-\frac{1}{2} \int\limits_{0}^{t} \bigl[\sin(3\tau-2t)+\sin(2t-\tau)\bigr]d\tau=\\ &=\cos2t+ \left.{\frac{1}{6}\cos(3\tau-2t)}\right|_{0}^{t}-\left.{\frac{1}{2} \cos(2t-\tau)}\right|_{0}^{t}= \frac{4}{3}\cos2t-\frac{1}{3}\cos t\,.\end{aligned}[/math]



Теорема о связи "начальных" и "конечных" значений оригинала и изображения


Начальное значение оригинала находится по формуле


[math]f(+0)= \lim\limits_{p\to\infty} pF(p).[/math]
(5.18)

Если существует конечный предел [math]\lim\limits_{t\to+\infty} f(t)= f(+\infty)[/math], то


[math]f(+\infty)= \lim\limits_{p\to0} pF(p).[/math]
(5.19)

Из соотношений (5.18),(5.19) следует, что для нахождения начальных и конечных значений оригинала не требуется знания оригинала, а достаточно иметь соответствующее изображение. На практике соотношение (5.19) применяется, например, для нахождения установившегося значения выходного сигнала в системах автоматического регулирования.


▼ Пример 5.19

Найти начальное и конечное значения оригинала, которому соответствует изображение [math]F(p)= \frac{p+1}{(p+1)^2+9}[/math].


Решение. Согласно (5.18) и (5.19) имеем


[math]f(+0)= \lim\limits_{p\to\infty} \frac{p(p+1)}{(p+1)^2+9}=1;\qquad f(+\infty)= \lim\limits_{p\to0} \frac{p(p+1)}{(p+1)^2+9}=0.[/math]

С другой стороны, из примера 5.7 следует, что


[math]f(t)= L^{-1}\! \left[\frac{p+1}{(p+1)^2+9}\right]=e^{-t}\cos3t\,,[/math]

поэтому легко убедиться в правильности полученного результата.


Полученные решения примеров 5.2–5.17 позволяют сформировать таблицу преобразования Лапласа. Табл. 5.1 является фрагментом более полных таблиц, используемых далее при решении примеров и задач.




Нахождение изображения по оригиналу


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


Таблица 5.1. Таблица основных преобразований Лапласа


▼ Посмотреть
[math]f(t)[/math][math]F(p)[/math][math]f(t)[/math][math]F(p)[/math]
1[math]\boldsymbol{1}(t)[/math][math]\frac{1}{p}[/math]17[math]\frac{1}{a^2}(1-\cos at)[/math][math]\frac{1}{p(p^2+a^2)}[/math]
2[math]C~(\text{const})[/math][math]\frac{C}{p}[/math]18[math]\frac{1}{a^2}(e^{at}-1-a)[/math][math]\frac{1}{p^2(p-a)}[/math]
3[math]t[/math][math]\frac{1}{p^2}[/math]19[math]\operatorname{sh}at[/math][math]\frac{a}{p^2-a^2}[/math]
4[math]t^n[/math][math]\frac{n!}{p^{n+1}}[/math]20[math]\operatorname{ch}at[/math][math]\frac{p}{p^2-a^2}[/math]
5[math]\delta(t)[/math][math]1[/math]21[math]\left(t+\frac{1}{2}at^2\right)\!e^{at}[/math][math]\frac{p}{(p-a)^3}[/math]
6[math]e^{at}[/math][math]\frac{1}{p-a}[/math]22[math]\left(1+2at+\frac{1}{2}a^2t^2\right)\!e^{at}[/math][math]\frac{p^2}{(p-a)^3}[/math]
7[math]t^ne^{at}[/math][math]\frac{n!}{(p-a)^{n+1}}[/math]23[math](1+at)e^{at}[/math][math]\frac{p}{(p-a)^2}[/math]
8[math]\sin at[/math][math]\frac{a}{p^2+a^2}[/math]24[math]\cos^2at[/math][math]\frac{p^2+2a^2}{p(p^2+4a^2)}[/math]
9[math]\cos at[/math][math]\frac{p}{p^2+a^2}[/math]25[math]\sin^2at[/math][math]\frac{2a^2}{p(p^2+4a^2)}[/math]
10[math]t\sin at[/math][math]\frac{2pa}{(p^2+a^2)^2}[/math]26[math]\sin \frac{at}{\sqrt{2}} \operatorname{sh}\frac{at}{\sqrt{2}}[/math][math]\frac{a^2p}{p^4+a^4}[/math]
11[math]t\cos at[/math][math]\frac{p^2-a^2}{(p^2+a^2)^2}[/math]27[math]\cos\frac{at}{\sqrt{2}} \operatorname{ch}\frac{at}{\sqrt{2}}[/math][math]\frac{p^3}{p^4+a^4}[/math]
12[math]e^{at}\sin bt[/math][math]\frac{b}{(p-a)^2+b^2}[/math]28[math]\frac{1}{2}(\operatorname{sh}at-\sin at)[/math][math]\frac{a^3}{p^4-a^4}[/math]
13[math]e^{at}\cos bt[/math][math]\frac{p-a}{(p-a)^2+b^2}[/math]29[math]\frac{1}{2}(\operatorname{ch}at-\cos at)[/math][math]\frac{a^2p}{p^4-a^4}[/math]
14[math]\frac{1}{a}\,e^{-t\!\not{\phantom{|}}\,\,a}[/math][math]\frac{1}{1+ap}[/math]30[math]\frac{1}{2}(\operatorname{sh}at+ \sin at)[/math][math]\frac{ap^2}{p^4-a^4}[/math]
15[math]\frac{1}{a}(e^{at}-1)[/math][math]\frac{1}{p(p-a)}[/math]31[math]\frac{1}{2}(\operatorname{ch}at+ \cos at)[/math][math]\frac{p^3}{p^4-a^4}[/math]
16[math]\frac{e^{at}-e^{bt}}{a-b}[/math][math]\frac{1}{(p-a)(p-b)}[/math]32[math]\frac{a\,e^{at}-b\,e^{bt}}{a-b}[/math][math]\frac{p}{(p-a)(p-b)}[/math]

Продолжение таблицы 5.1


[math]f(t)[/math][math]F(p)[/math][math]f(t)[/math][math]F(p)[/math]
33[math]\frac{(c-b)e^{at}+ (a-c)e^{bt}+ (b-a)e^{ct}}{(a-b)(a-c)(c-b)}[/math][math]\frac{1}{(p-a)(p-b)(p-c)}[/math]42[math]\frac{a \operatorname{sh}at-b\operatorname{sh}bt}{a^2-b^2}[/math][math]\frac{p^2}{(p^2-a^2)(p^2-b^2)}[/math]
34[math]\frac{a(c-b)e^{at}+ b(a-c)e^{bt}+ c(b-a)e^{ct}}{(a-b)(a-c)(c-b)}[/math][math]\frac{p}{(p-a)(p-b)(p-c)}[/math]43[math]\frac{a^2 \operatorname{ch}at-b^2\operatorname{ch}bt}{a^2-b^2}[/math][math]\frac{p^3}{(p^2-a^2)(p^2-b^2)}[/math]
35[math]\frac{a^2(c-b)e^{at}+ b^2(a-c)e^{bt}+ c^2(b-a)e^{ct}}{(a-b)(a-c)(c-b)}[/math][math]\frac{p^2}{(p-a)(p-b)(p-c)}[/math]44[math]t-\frac{1}{a}\sin at[/math][math]\frac{a^2}{p^2(p^2+a^2)}[/math]
36[math]\frac{a\sin bt-b\sin at}{a^2-b^2}[/math][math]\frac{ab}{(p^2+a^2)(p^2+b^2)}[/math]45[math]\frac{1}{a}\operatorname{sh}at-t[/math][math]\frac{a^2}{p^2(p^2-a^2)}[/math]
37[math]\frac{a\cos bt-b\cos at}{a^2-b^2}[/math][math]\frac{p}{(p^2+a^2)(p^2+b^2)}[/math]46[math]1-\cos at-\frac{at}{2}\sin at[/math][math]\frac{a^4}{p(p^2+a^2)^2}[/math]
38[math]\frac{a\sin at-b\sin bt}{a^2-b^2}[/math][math]\frac{p^2}{(p^2+a^2)(p^2+b^2)}[/math]47[math]1-\operatorname{ch}at+ \frac{at}{2}\operatorname{sh}at[/math][math]\frac{a^4}{p(p^2-a^2)^2}[/math]
39[math]\frac{a^2\cos at-b^2\cos bt}{a^2-b^2}[/math][math]\frac{p^3}{(p^2+a^2)(p^2+b^2)}[/math]48[math]1+\frac{b^2\coa at-a^2\cos bt}{a^2-b^2}[/math][math]\frac{a^2b^2}{p(p^2+a^2)(p^2+b^2)}[/math]
40[math]\frac{b \operatorname{sh}at-a \operatorname{sh}bt}{a^2-b^2}[/math][math]\frac{ab}{(p^2-a^2)(p^2-b^2)}[/math]49[math]1+\frac{b^2 \operatorname{ch}at-a^2 \operatorname{ch}bt}{a^2-b^2}[/math][math]\frac{a^2b^2}{p(p^2-a^2)(p^2-b^2)}[/math]
41[math]\frac{\operatorname{ch}at-\operatorname{ch}bt}{a^2-b^2}[/math][math]\frac{p}{(p^2-a^2)(p^2-b^2)}[/math]

▼ Пример 5.21

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


a) [math]f(t)=t\,e^{2t}\sin3t[/math]; б) [math]f(t)=\frac{2}{t}(1-\cos t)[/math];


в) [math]f(t)=\begin{cases}\cos[4(t-2)],&t>3,\\ 0,&t \leqslant 2;\end{cases}[/math]; г) [math]f(t)= \begin{cases}e^{-3t}\cos[4(t-2)],& t>2,\\ 0,& t \leqslant 2;\end{cases}[/math];


д) [math]f(t)= \frac{1}{2}(t-2)^2e^{-(t-2)}\boldsymbol{1}(t-2)[/math]; е) [math]f(t)= e^{2t}+ \boldsymbol{1}(t-1)+ \boldsymbol{1}(t-4)\sin[3(t-4)][/math];


ж) [math]\int\limits_{0}^{t} \frac{\operatorname{sh}\tau}{\tau}\,d\tau[/math]; з) [math]f(t)= \cos(4t-8)[/math].


▼ Решение

а) Согласно формуле 12 из табл. 5.1 [math]L \bigl[e^{2t}\sin3t\bigr]= \frac{3}{(p-2)^2+9}[/math]. По свойству дифференцирования изображения (формула (5.13) при [math]n=1[/math]):


[math]L \bigl[(-1)tf(t)\bigr]=F'(p)[/math] или [math]L \bigl[tf(t)\bigr]=-F'(p)[/math].

Поэтому [math]L\bigl[t\,e^{2t}\sin3t\bigr]=-\left[\frac{3}{(p-2)^2+9}\right]'=-\frac{-2(p-2)\cdot3}{\bigl[(p-2)^2+9\bigr]^2}= \frac{6(p-2)}{\bigl[(p-2)^2+9\bigr]^2}[/math].


Можно решать иначе, используя формулу 10 из табл. 5.1 при а = 3 и свойство смещения при [math]a=2\colon[/math]


[math]L \bigl[t\sin3t\bigr]= \frac{6p}{(p^2+9)^2}[/math] и [math]L \bigl[t\,e^{2t}\sin3t\bigr]= \frac{6(p-2)}{\bigl[(p-2)^2+9\bigr]^2}[/math].

б) Применим свойства линейности и интегрирования изображения. Так как [math]L[1-\cos t]=\frac{1}{p}-\frac{p}{p^2+1}[/math] в силу формул 1 и 9 из табл. 5.1 , то согласно (5.14)


[math]L\! \left[\frac{2}{t}(1-\cos t)\right]= 2 \int\limits_{p}^{+\infty}\! \left(\frac{1}{z}-\frac{z}{z^2+1}\right)\!dz= \ldots= \ln\frac{p^2+1}{p^2}\,.[/math]

в) По формуле 9 из табл. 5.1 [math]L[\cos4t]= \frac{p}{p^2+16}[/math]. Функцию fit) можно записать, используя единичную ступенчатую функцию: [math]f(t)= \cos[4(t-2)]\cdot \boldsymbol{1}(t-2)[/math]. Поэтому для нахождения изображения следует применить теорему запаздывания (5.9) к оригиналу [math]\cos4t[/math] при [math]\tau=2\colon[/math]


[math]L \bigl[\cos[4(t-2)]\bigr]= e^{-2p}\cdot \frac{p}{p^2+16}\,.[/math]

г) Согласно теореме смещения (5.8) и с учетом результата п. "в" имеем


[math]L \bigl[e^{-3t}\cos[4(t-2)]\bigr]= e^{-2(p+3)}\cdot \frac{p+3}{(p+3)^2+16}[/math], так как [math]a=-3[/math].

д) По теореме запаздывания (5.9) при [math]\tau=2[/math] и по формуле 7 из табл. 5.1 при [math]a=2,~ b=-1[/math] получаем


[math]L\! \left[\frac{1}{2}(t-2)^2e^{-(t-2)} \boldsymbol{1}(t-2)\right]= e^{-2p}\cdot L\! \left[\frac{1}{2}\,t^2e^{-t}\right]= e^{-2p}\cdot \frac{1}{2}\cdot \frac{2!}{(p+1)^3}= \frac{e^{-2p}}{(p+1)^3}\,.[/math]

е) Используя свойства линейности, запаздывания и формулы 6,1,8 из табл. 5.1, получаем


[math]F(p)= \frac{1}{p-2}+ \frac{1}{p}\,e^{-p}+ e^{-4p}\cdot \frac{3}{p^2+9}\,.[/math]

ж) По формуле 19 из табл. 5.1 находим [math]L[\operatorname{sh}t]=\frac{1}{p^2-1}[/math]. По свойству интегрирования изображения (формула (5.14)):


[math]L\! \left[\frac{\operatorname{sh}t}{t}\right]= \int\limits_{p}^{+\infty} \frac{dz}{z^2-1}= \left.{\frac{1}{2}\ln \frac{z-1}{z+1}}\right|_{p}^{+\infty}= \ldots=-\frac{1}{2} \ln \frac{p-1}{p+1}= \frac{1}{2}\ln \frac{p+1}{p-1}\,.[/math]

По свойству интегрирования оригинала (5.12): [math]L \Biggl[\int\limits_{0}^{t} \frac{\operatorname{sh}\tau}{\tau}\,d\tau\Biggr]= \frac{1}{2p}\ln \frac{p+1}{p-1}[/math].


з) Используем формулу косинуса разности и запишем оригинал в виде суммы:


[math]f(t)= \cos(4t-8)= \cos4t\cos8+\sin4t\sin8[/math]

По свойству линейности получаем:


[math]L \bigl[\cos(4t-8)\bigr]= \cos8\cdot \frac{p}{p^2+16}+ \sin8\cdot \frac{4}{p^2+16}= \frac{p\cos8+4\sin8}{p^2+16}\,.[/math]

Заметим, что здесь [math]f(t)=f(t)\cdot \boldsymbol{1}(t)[/math] и результаты пп. "в" и "з" различны так как оригиналами являются разные функции.




Нахождение изображений функций, заданных графиком


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


Пример 5.22. Найти изображения функций, заданных графиками на рис 5.5.


▼ Решение

а) Представим функцию в виде [math]f(t)=(1-t)\cdot \boldsymbol{1}(t-1)=-(t-1)\cdot \boldsymbol{1}(t-1)[/math].


По формуле 3 из табл.5.1 и теореме запаздывания (формула (5.9) при [math]\tau=1[/math])


[math]F(p)=-\frac{1}{p^2}\cdot e^{-p}.[/math]

б) Запишем функцию в виде [math]f(t)=(1-t) \bigl[\boldsymbol{1}(t)+\boldsymbol{1}(t-1)\bigr]=(1-t)\cdot \boldsymbol{1}(t)+ (t-1)\cdot \boldsymbol{1}(t-1)[/math].


По формулам 1,3 из табл. 5.1 и (5.6), (5.9) имеем [math]F(p)=\frac{1}{p}-\frac{1}{p^2}+\frac{e^{-p}}{p^2}[/math].


в) Запишем изображенную функцию в виде [math]f(t)=-t\cdot \boldsymbol{1}(t-1)= (-t+1-1)\cdot \boldsymbol{1}(t-1)=-(t-1)\cdot \boldsymbol{1}(t-1)-\boldsymbol{1}(t-1)[/math].


По формулам 1,3 из табл. 5.1 и (5.6), (5.9) имеем [math]F(p)=-\frac{e^{-p}}{p^2}-\frac{e^{-p}}{p}[/math].


г) Представим функцию в виде


[math]\begin{aligned}f(t)&= (t-1)\cdot \bigl[\boldsymbol{1}(t-1)-\boldsymbol{1}(t-2)\bigr]+ (3-t)\cdot \bigl[\boldsymbol{1}(t-2)-\boldsymbol{1}(t-3)\bigr]=\\ &=(t-1)\cdot \boldsymbol{1}(t-1)-2(t-2)\cdot \boldsymbol{1}(t-2)+ (t-3)\cdot \boldsymbol{1}(t-3). \end{aligned}[/math]

По формулам 3 из табл. 5.1 и (5.6), (5.9) имеем [math]F(p)= \frac{e^{-p}}{p^2}-\frac{2e^{-2p}}{p^2}+\frac{e^{-3p}}{p^2}[/math].


д) Запишем функцию в форме


[math]\begin{aligned}f(t)&= t \bigl[\boldsymbol{1}(t)-\boldsymbol{1}(t-1)\bigr]+ \bigl[\boldsymbol{1}(t-1)-\boldsymbol{1}(t-2)\bigr]+ (3-t)\cdot \bigl[\boldsymbol{1}(t-2)-\boldsymbol{1}(t-3)\bigr]=\\[2pt] &=t\cdot\boldsymbol{1}(t)-(t-1)\cdot \boldsymbol{1}(t-1)-\boldsymbol{1}(t-1)+\boldsymbol{1}(t-1)-\boldsymbol{1}(t-2)\,+\\ &\quad+(2-t)\cdot \boldsymbol{1}(t-2)+ \boldsymbol{1}(t-2)-(3-t)\cdot \boldsymbol{1}(t-3)=\\[2pt] &=t\cdot \boldsymbol{1}(t)-(t-1)\cdot \boldsymbol{1}(t-1)-(t-2)\cdot \boldsymbol{1}(t-2)+ (t-3)\cdot \boldsymbol{1}(t-3). \end{aligned}[/math]

По формулам З из табл. 5.1 и (5.6), (5.9) [math]F(p)=\frac{1}{p^2}-\frac{e^{-p}}{p^2}-\frac{e^{-2p}}{p^2}+\frac{e^{-3p}}{p^2}[/math].


е) Представим изображенную функцию в виде


[math]\begin{aligned}f(t)&= (1-t) \bigl[\boldsymbol{1}(t)-\boldsymbol{1}(t-1)\bigr]-1\cdot \bigl[\boldsymbol{1}(t-1)-\boldsymbol{1}(t-2)\bigr]+ (t-3)\cdot \boldsymbol{1}(t-2)=\\[2pt] &=(1-t)\cdot\boldsymbol{1}(t)+ (t-1)\cdot\boldsymbol{1}(t-1)-\boldsymbol{1}(t-1)+ \boldsymbol{1}(t-2)+ (t-2)\cdot \boldsymbol{1}(t-2)-\boldsymbol{1}(t-2)=\\[2pt] &=(1-t)\cdot\boldsymbol{1}(t)+ (t-1)\cdot\boldsymbol{1}(t-1)-\boldsymbol{1}(t-1)+ (t-2)\cdot \boldsymbol{1}(t-2).\end{aligned}[/math]

По формулам 1,3 из табл. 5.1 и (5.6), (5.9) [math]F(p)=\frac{1}{p}-\frac{1}{p^2}+\frac{1}{p^2}\,e^{-p}-\frac{1}{p}\,e^{-p}+\frac{1}{p^2}\,e^{-2p}[/math].


ж) Запишем функцию в форме [math]f(t)= \sin t\cdot \bigl[\boldsymbol{1}(t)-\boldsymbol{1}(t-\pi)\bigr]= \sin t\cdot\boldsymbol{1}(t)+ \sin(t-\pi)\cdot \boldsymbol{1}(t-\pi)[/math].


По формуле 8 из табл. 5.1 и по теореме запаздывания [math]F(p)= \frac{1}{p^2+1}+ \frac{e^{\pi p}}{p^2+1}[/math].


з) Представим функцию в виде [math]f(t)= e^{-t} \bigl[\boldsymbol{1}(t)-\boldsymbol{1}(t-1)\bigr]= e^{-t}\cdot\boldsymbol{1}(t)-\frac{e^{-(t-1)}}{e}\cdot \boldsymbol{1}(t-1).[/math].


По формуле 6 из табл. 5.1 при [math]a=-1[/math] и (5.9) при [math]\tau=1[/math] имеем


[math]F(p)= \frac{1}{p+1}-\frac{1}{e}\cdot \frac{e^{-p}}{p+1}= \frac{1}{p+1}\bigl(1-e^{-p-1}\bigr).[/math]

и) Представим функцию в виде [math]f(t)=1\cdot \bigl[\boldsymbol{1}(t)-\boldsymbol{1}(t-1)\bigr]-1\cdot \bigl[\boldsymbol{1}(t-1)-\boldsymbol{1}(t-2)\bigr]= \boldsymbol{1}(t)-2\cdot \boldsymbol{1}(t-1)+\boldsymbol{1}(t-2).[/math].


Используя формулы 2 из табл. 5.1 и (5.6), (5.9), получаем [math]F(p)=\frac{1}{p}-\frac{2}{p}e^{-p}+ \frac{1}{p}\,e^{-2p}[/math].




Нахождение изображений периодических функций


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


Пусть [math]f(t)[/math] — оригинал с периодом [math]T[/math] (рис. 5.6,в), образованный повторением функции [math]f_0(t)[/math] (рис. 5.6,б):


[math]f_0(t)= \begin{cases}0,& t\leqslant 0,\\ f(t),& 0<t\leqslant T,\\ 0,& t>T.\end{cases}[/math]

Для нахождения изображения [math]F(p)[/math] периодической функции [math]f(t)[/math] следует:


1. Найти изображение функции [math]f_0(t)\colon\, F_0(p)=L\bigl[f_0(t)\bigr][/math].
2. Найти изображение [math]F(p)[/math] по формуле

[math]F(p)= \frac{F_0(p)}{1-e^{-Tp}}[/math]
(5.20)

Пример 5.23. Найти изображения функций, представленных на рис. 5.7.


▼ Решение

а) По графику (рис. 5.7,в) получаем


[math]f_0(t)= t\cdot \bigl[\boldsymbol{1}(t)-\boldsymbol{1}(t-1)\bigr]= 1\cdot \boldsymbol{1}(t)-(t-1)\cdot \boldsymbol{1}(t-1)-\boldsymbol{1}(t-1).[/math]

Поэтому [math]F_0(p)= \frac{1}{p^2}-\frac{1}{p^2}\,e^{-p}-\frac{1}{p}\,e^{-p}[/math].


Поскольку [math]T=1[/math], по формуле (5.20) находим


[math]F(p)= \frac{\dfrac{1}{p^2}(1-e^{-p}-pe^{-p})}{1-e^{-p}}= \frac{e^p(1-e^{-p}-pe^{-p})}{p^2(e^p-1)}= \frac{e^p-1-p}{p^2(e^p-1)}\,.[/math]

б) По графику (рис. 5.7,б) имеем [math]f_0(t)= \boldsymbol{1}(t)-\boldsymbol{1}(t-\tau)[/math], тогда [math]F_0(p)=\frac{1}{p}-\frac{1}{p}\,e^{-p\tau}[/math]. По формуле (5.20) при [math]T=2\tau[/math] имеем


[math]F(p)= \frac{1}{p}(1-e^{-p\tau})\frac{1}{1-e^{-2\tau p}}= \frac{1}{p(1+e^{-p\tau})}\,.[/math]

в) Функция, изображенная на рис. 5.7,в , имеет период [math]T=2c[/math]. Запишем аналитическое выражение для [math]f_0(t)[/math] и соответствующее изображение [math]F_0(p)\colon[/math]


[math]\begin{aligned}&f_0(t)= h \bigl[\boldsymbol{1}(t)-\boldsymbol{1}(t-c)\bigr]-h \bigl[\boldsymbol{1}(t-c)-\boldsymbol{1}(t-2c)\bigr]= h\cdot\boldsymbol{1}(t)-2h\cdot \boldsymbol{1}(t-c)+h\cdot \boldsymbol{1}(t-2c),\\ &F_0(p)= \frac{h}{p}-\frac{2h}{p}\,e^{-pc}+\frac{h}{p}\,e^{-2pc}.\end{aligned}[/math]

По формуле (5.20) получаем [math]F(p)=\frac{h(1+e^{-2pc}-2e^{-pc})}{p(1-e^{2pc})}[/math].


г) Для функции, изображенной на рис. 5.7,г, изображением для [math]f_0(t)[/math] является [math]F_0(p)= \frac{1}{p^2+1}(1+e^{-\pi p})[/math] (см. пример 5.22 п."ж"). Тогда по формуле (5.20) при [math]T=\pi[/math] получаем [math]F(p)= \frac{1+e^{-\pi p}}{(p^2+1)(1-e^{-\pi p})}[/math].




Нахождение оригинала по изображению


Непосредственное применение формулы обращения (5.3) затруднительно, поэтому для нахождения оригинала применяются теоремы разложения и правила преобразования изображения к виду, представленному в табл. 5.1.


Применение теорем разложения


Теорема 5.1 (первая теорема разложения). Если функция [math]F(p)[/math] аналитична в некоторой окрестности бесконечно удаленной точки и ее разложение в ряд по степеням [math]\frac{1}{p}[/math] имеет вид [math]\textstyle{F(p)= \sum\limits_{n=0}^{\infty} \dfrac{a_n}{p^{n+1}}}[/math], то функция


[math]f(t)= \sum\limits_{n=0}^{\infty} a_n \frac{t^n}{n!},\quad t\geqslant0[/math]
(5.21)

является оригиналом, соответствующим изображению [math]F(p)[/math].


Теорема 5.2 (вторая теорема разложения). Если изображение [math]F(p)[/math] является однозначной функцией и имеет лишь конечное число особых точек [math]p_1,p_2,\ldots,p_n[/math] лежащих в конечной части плоскости, то


[math]f(t)= \sum\limits_{k=1}^{n} \mathop{\operatorname{res}}\limits_{p=p_k} \bigl[e^{pt}F(p)\bigr].[/math]
(5.22)

Замечания 5.5

1. Формула (5.21) может быть записана в виде [math]\textstyle{L^{-1}\! \left[\sum\limits_{n=0}^{\infty} \dfrac{a_n}{p^{n+1}}\right]= \sum\limits_{n=0}^{\infty}a_nL^{-1}\! \left[\dfrac{1}{p^{n+1}}\right]}[/math]. Задача нахождения оригинала при выполнении условий теоремы сводится к нахождению коэффициентов разложения функции в ряд Лорана в окрестности бесконечно удаленной точки.


2. Формула (5.22) принимает наиболее простой вид в случае [math]F(p)=R(p)[/math] — рационального изображения, т.е. [math]F(p)= R(p)= \frac{P_m(p)}{Q_n(p)}[/math], где [math]P_m(p),\,Q_n(p)[/math] — многочлены степеней /пил соответственно, не имеющие общих корней. Если все полюсы [math]p_1,p_2,\ldots,p_n[/math] функции [math]F(p)[/math] простые, то по формуле (4.24) получаем [math]\mathop{\operatorname{res}}\limits_{p=p_k} \frac{P_m(p)}{Q_n(p)}\,e^{pt}= \frac{P_m(p_k)}{Q'_n(p_k)}\,e^{p_kt}[/math], а формула (5.22) принимает вид


[math]f(t)= \sum\limits_{k=1}^{n} \frac{P_m(p_k)}{Q'_n(p_k)}\,e^{p_kt}.[/math]
(5.23)

3. Если при выполнении условий п.2 коэффициенты многочлена [math]Q_n(p)[/math] — лействительные числа, то его комплексные корни, как известно, являются по парно сопряженными. Нахождение суммы вычетов в таких точках можно заме нить нахождением действительной части вычета в одной из них. Действительио, вычет в точке [math]\overline{p}_k[/math], используя свойства сопряженных чисел, можно записать следующим образом:


[math]\mathop{\operatorname{res}}\limits_{p=\overline{p}_k} \frac{P_m(p)}{Q_n(p)}\,e^{pt}= \frac{P_m(\overline{p}_k)}{Q'_n(\overline{p}_k)}\,e^{\overline{p}_kt}= \frac{\overline{P_m(p_k)}}{\overline{Q'_n(p_k)}}\,e^{\overline{p_kt}}= \overline{\frac{P_m(p_k)}{Q'_n(p_k)}\,e^{p_kt}}.[/math]

Это означает, что вычет в точке [math]\overline{p}_k[/math] есть число, сопряженное вычету в точке [math]p_k[/math], а сумма таких чисел равна их удвоенной действительной части:


[math]\mathop{\operatorname{res}}\limits_{p=p_k} \frac{P_m(p)}{Q_n(p)}\,e^{pt}+ \mathop{\operatorname{res}}\limits_{p=\overline{p}_k} \frac{P_m(p)}{Q_n(p)}\,e^{pt}= 2 \operatorname{Re} \mathop{\operatorname{res}}\limits_{p=p_k} \frac{P_m(p)}{Q_n(p)}\,e^{pt}.[/math]

▼ Пример 5.24

Найти оригиналы для функций:


a) [math]F_1(p)=\frac{1}{p}\exp \frac{1}{p^2},\quad F_2(p)=\frac{1}{p}\cos \frac{1}{p},\quad F_2(p)= \frac{1}{\sqrt{p}}\sin \frac{1}{\sqrt{p}}[/math];


б) [math]F_1(p)=\frac{p}{p^2+4p+5},\quad F_2(p)=\frac{p+2}{(p+1)(p-2)(p^2+4)},\quad F_3(p)=\frac{p^2+p+1}{(p-1)(p+1)^2}[/math].


Решение. В случае "а" для решения задачи используем теорему 5.1, а в случае "б" — теорему 5.2.


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


[math]e^z=\sum\limits_{n=0}^{\infty} \frac{z^n}{n!},\qquad \cos z=\sum\limits_{n=0}^{\infty} \frac{(-1)^nz^{2n}}{(2n)!},\qquad \sin z=\sum\limits_{n=1}^{\infty} \frac{(-1)^{n+1}z^{2n-1}}{(2n-1)!}\,.[/math]

Для заданных изображений получаем:


[math]\begin{aligned}F_1(p)&= \frac{1}{p} \sum\limits_{n=0}^{\infty} \frac{1}{p^{2n}n!}= \sum\limits_{n=0}^{\infty} \frac{1}{p^{2n+1}n!},\quad a_{2n}=\frac{1}{n!}\,;\\[2pt] F_2(p)&= \frac{1}{p} \sum\limits_{n=0}^{\infty} \frac{(-1)^n}{p^{2n}(2n)!}= \sum\limits_{n=0}^{\infty} \frac{(-1)^n}{p^{2n+1}(2n)!},\quad a_{2n}= \frac{(-1)^n}{(2n)!}\,;\\[2pt] F_3(p)&= \frac{1}{\sqrt{p}} \sum\limits_{n=1}^{\infty} \frac{(-1)^n}{\sqrt{p^{2n+1}}(2n+1)!}= \sum\limits_{n=1}^{\infty} \frac{(-1)^n}{(2n+1)!p^{n+1}},\quad a_n= \frac{(-1)^n}{(2n+1)!}\,.\end{aligned}[/math]

Согласно первой теореме разложения


[math]f_1(t)= \sum\limits_{n=0}^{\infty} \frac{1}{n!}\frac{t^{2n}}{(2n)!}\,,\qquad f_2(t)= \sum\limits_{n=0}^{\infty} \frac{(-1)^nt^{2n}}{[(2n)!]^2}\,,\qquad f_3(t)= \sum\limits_{n=0}^{\infty} \frac{(-1)^nt^n}{(2n+1)!n!}\,.[/math]

б) Представим [math]F_1(p)[/math] в виде


[math]F_1(p)= \frac{p}{p^2+4p+5}= \frac{P_1(p)}{Q_2(p)}= \frac{p}{\bigl[p-(-2+i)\bigr] \bigl[p-(-2-i)\bigr]}\,.[/math]

где [math]p_1=-2+i,~ p_2=-2-i[/math] — простые полюсы функции [math]F(p)[/math]. По второй теореме разложения


[math]\begin{aligned}f_1(t)&= \sum\limits_{k=1}^{2}\frac{P_1(p_k)}{Q'_2(p_k)}\,e^{p_kt}= \sum\limits_{k=1}^{2} \frac{p_k}{2p_k+4}\,e^{p_kt}= \frac{-2+i}{-4+2i+4}\,e^{(-2+i)t}+ \frac{-2-i}{-4-2i+4}\,e^{(-2-i)t}=\\ &=e^{-2t}\! \left(\frac{-2+i}{2i}\,e^{it}+\frac{2+i}{2i}\,e^{-it}\right)= e^{-2t}\! \left(\frac{e^{it}+e^{-it}}{2}-2\cdot\frac{e^{it}-e^{-it}}{2i}\right)= e^{-2t} \bigl[\cos t-2\sin t\bigr]. \end{aligned}[/math]

Тот же результат можно получить, пользуясь пп. 2 и 3 замечаний 5.5:


[math]\begin{aligned}f_1(t)&= 2 \operatorname{Re} \mathop{\operatorname{res}}\limits_{p=-2+i} \frac{p\,e^{pt}}{p^2+4p+5}= \left.{2 \operatorname{Re} \frac{p\,e^{pt}}{2p+4}}\right|_{-2+i}= 2 \operatorname{Re} \frac{-2+i}{2i}\,e^{(-2+i)t}=\\ &=2 \operatorname{Re}\! \left[e^{-2t}\, \frac{1+2i}{2}\,e^{it}\right]= \operatorname{Re} \bigl[e^{-2t}(1+2i)(\cos t+i\sin t)\bigr]= e^{-2t}(\cos t-2\sin t). \end{aligned}[/math]

Функция [math]F_2(p)[/math] имеет четыре простых полюса: [math]p_1=-1,~ p_2=2,~ p_3=2i[/math].


Так как вычет в простом полюсе находится по формуле [math]\mathop{\operatorname{res}}\limits_{p=p_k} F(p)= \lim\limits_{p\to p_k} F(p)(p-p_k)[/math], то по второй теореме разложения


[math]\begin{aligned}f_2(t)&= \sum\limits_{k=1}^{4} \mathop{\operatorname{res}}\limits_{p_k} \bigl[F(p)e^{pt}\bigr]= \lim\limits_{p\to-1} \frac{(p+2)e^{pt}}{(p-2)(p^2+4)}+ \lim\limits_{p\to2} \frac{(p+2)e^{pt}}{(p+1)(p^2+4)}\,+\\ &\qquad + \lim\limits_{p\to2i} \frac{(p+2)e^{pt}}{(p-1)(p+1)(p+2i)}+ \lim\limits_{p\to-2i} \frac{(p+2)e^{pt}}{(p-1)(p+1)(p-2i)}=\\[2pt] &=-\frac{1}{15}\,e^{-t}+ \frac{1}{6}\,e^{2t}+\frac{2+2i}{8-24i}\,e^{2it}+ \frac{2-2i}{8+24i}\,e^{-t}=\\[2pt] &=-\frac{1}{15}\,e^{-t}+ \frac{1}{6}\,e^{2t}+\frac{-32+64i}{640}\,e^{2it}+ \frac{-32-64i}{640}\,e^{-t}=\\[2pt] &=-\frac{1}{15}\,e^{-t}+ \frac{1}{6}\,e^{2t}-\frac{64}{640}\cdot \frac{e^{2it}+e^{-2it}}{2}-\frac{128}{640}\cdot \frac{e^{2it}-e^{-2it}}{2i}=\\[2pt] &=-\frac{1}{15}\,e^{-t}+ \frac{1}{6}\,e^{2t}-\frac{1}{10}\cos2t-\frac{1}{5}\sin2t\,.\end{aligned}[/math]

Функция [math]F_3(p)[/math] имеет два полюса: простой [math]p_1=1[/math] и полюс второго порядка [math]p_2=-1[/math].


По второй теореме разложения [math]f_3(t)= \mathop{\operatorname{res}}\limits_{p_1=1} \frac{(p^2+p+1)e^{pt}}{(p-1)(p+1)^2}+ \mathop{\operatorname{res}}\limits_{p_2=-1} \frac{(p^2+p+1)e^{pt}}{(p-1)(p+1)^2}[/math]. Находим вычеты


[math]\begin{aligned}\mathop{\operatorname{res}}\limits_{p_1=1} \frac{(p^2+p+1)e^{pt}}{(p-1)(p+1)^2}&= \lim\limits_{p\to1} \frac{(p^2+p+1)(p-1)e^{pt}}{(p-1)(p+1)^2}= \lim\limits_{p\to1} \frac{(p^2+p+1)e^{pt}}{(p+1)^2}= \frac{3}{4}\,e^{t},\\[2pt] \mathop{\operatorname{res}}\limits_{p_2=-1} \frac{(p^2+p+1)e^{pt}}{(p-1)(p+1)^2}&= \lim\limits_{p\to-1} \frac{d}{dp}\! \left[\frac{(p^2+p+1)e^{pt}(p+1)^2}{(p-1)(p+1)^2}\right]= \ldots= \frac{1}{4}\,e^{-t}-\frac{t}{2}\,e^{-t}, \end{aligned}[/math]

получаем окончательный ответ [math]f_3(t)= \frac{3}{4}\,e^{t}+\frac{1}{4}\,e^{-t}-\frac{t}{2}\,e^{-t}[/math].




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


Приведем ряд известных приемов нахождения оригинала.


1. Если изображение отличается от табличного на постоянный множитель, то его следует умножить и одновременно поделить на этот множитель, а затем воспользоваться свойством линейности.


▼ Пример 5.25

Найти оригиналы для функций a) [math]F(p)= \frac{3}{(p+2)^3}[/math]; б) [math]F(p)= \frac{5}{p^2-7}[/math].


Решение. а) Запишем изображение в виде [math]F(p)= \frac{3}{2!}\cdot \frac{2!}{(p+2)^3}[/math] и при [math]a--2,~ n=2[/math] по формуле 7 из табл. 5.1 получаем [math]f(t)= \frac{3}{2}\,t^2e^{-2t}[/math].


б) Представим изображение в форме [math]F(p)= \frac{5}{p^2-7}= \frac{5}{\sqrt{7}}\cdot \frac{\sqrt{7}}{p^2-(\sqrt{7})^2}[/math] и при [math]a=\sqrt{7}[/math] в формуле 19 из табл. 5.1 получаем [math]f(t)= \frac{5}{\sqrt{7}}\operatorname{sh}(t\sqrt{7})[/math].


2. Изображение, заданное в виде дроби [math]\frac{a\pm b}{c}[/math], разлагается на сумму дробей.


▼ Пример 5.26

Найти оригинал для функций: а) [math]F(p)= \frac{3p}{(p+5)^2}[/math]; б) [math]F(p)= \frac{3p-2}{(p+5)^2}[/math]; в) [math]F(p)= \frac{p^3+2p+2}{p^3(p+1)}[/math].


Решение. Представим дроби в виде суммы двух слагаемых, а затем воспользуемся свойством линейности и формулами из табл. 5.1:


а) [math]F(p)= \frac{3(p+5)-15}{(p+5)^2}= 3\cdot\frac{1}{p+5}-15\cdot \frac{1}{(p+5)^2}\quad \Rightarrow\quad f(t)=3e^{-5t}-15te^{-5t}[/math];


б) [math]F(p)= 3\cdot \frac{p}{(p+5)^2}-\frac{2}{1!}\cdot \frac{1!}{(p+5)^2}\quad \Rightarrow\quad f(t)= 3(1-5t)e^{-5t}-2te^{-5t}= e^{-5t}(3-17t)[/math].


в) представим [math]F(p)[/math] в виде [math]F(p)= \frac{p^3+2p+2}{p^3(p+1)}= \frac{p^3+2(p+ 1)}{p^3(p+1)}= \frac{2}{p^3}+ \frac{1}{p+1}[/math]. По формулам 4,6 из табл. 5.1 находим [math]f(t)= t^2+e^{-t}[/math].


3. Если знаменатель дроби содержит квадратный трехчлен, то в нем выде ляется полный квадрат: [math]ap^2+bp+c= a(p\pm\alpha)^2\pm\omega^2[/math]. При этом числитель дроби представляется в виде многочлена от [math](p\pm\alpha)[/math].


▼ Пример 5.27

Найти оригиналы для функций: a) [math]F(p)=\frac{3}{p^2+4p+7}[/math]; б) [math]F(p)= \frac{3p+2}{2p^2-8p+6}[/math].


Решение. а) Выделим полный квадрат в знаменателе дроби и воспользуемся табл. 5.1 (по формуле 12 из табл. 5.1 при [math]a=-2,~ b=\sqrt{3}[/math]):


[math]F(p)= \frac{3}{(p+2)^2-4+7}= \frac{3}{(p+2)^2+3}= \sqrt{3}\cdot \frac{\sqrt{3}}{(p+2)^2+(\sqrt{3})^2}\quad \Rightarrow\quad f(t)=\sqrt{3}\,e^{-2t}\sin(t \sqrt{3}).[/math]

б) Используем представление


[math]F(p)= \frac{3p+2}{2(p^2-4p+3)}= \frac{1,\!5p+1}{(p-2)^2-1}= \frac{1,\!5(p-2)+3}{(p-2)^2-1}+ \frac{1}{(p-2)^2-1}= 1,\!5\cdot\frac{p-2}{(p-2)^2-1}+4\cdot \frac{1}{(p-2)^2-1}\,.[/math]

По формулам 19,20 из табл. 5.1 и по теореме смещения (формула (5.8))


[math]f(t)=\frac{3}{2}\,e^{2t}\operatorname{ch}t+ 4e^{2t}\operatorname{sh}t= \frac{3}{2}\,e^{2t}\,\frac{e^t+e^{-t}}{2}+ 4e^{2t}\,\frac{e^t-e^{-t}}{2}= \frac{11}{4}\,e^{3t}-\frac{5}{4}\,e^{t}.[/math]

Можно решить эту задачу иначе, используя вторую теорему разложения (см. п. "б" примера 5.24):


[math]\begin{gathered}F(p)= \frac{3p+2}{2(p^2-4p+3)}= \frac{3p+2}{2(p-1)(p-3)}\\ \Downarrow\\ f(t)= \mathop{\operatorname{res}}\limits_{p=1} F(p)e^{pt}+ \mathop{\operatorname{res}}\limits_{p=3} F(p)e^{pt}= \frac{5e^t}{2\cdot(-2)}+ \frac{11e^{3t}}{2\cdot2}= \frac{11}{4}\,e^{3t}-\frac{5}{4}\,e^{t}.\end{gathered}[/math]

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


▼ Примеры 5.28-5.29

Пример 5.28. Найти оригиналы для функций:


а) [math]F(p)=\frac{3p^2+3p-13}{p(p^4+4p+13)}[/math]; б) [math]F(p)=\frac{p}{(p+2)^2(p-1)}[/math]; в) [math]F(p)=\frac{p^2+4}{(p-1)(p+2)(p-3)}[/math]; г) [math]F(p)=\frac{3p}{2p^2-2p-4}[/math].


▼ Решение

а) Представим [math]F(p)[/math] в виде [math]F(p)=\frac{3p^2+3p-13}{p(p^4+4p+13)} = \frac{A}{p}+ \frac{Bp+C}{p^2+4p+13}[/math], где [math]A,\,B,\,C[/math] — неопределенные коэффициенты.


Отсюда следует равенство [math]3p^2+3p-13= Ap^2+4Ap+13A+Bp^2+Cp[/math].


Приравнивая коэффициенты при одинаковых степенях [math]p[/math], получаем систему уравнений для нахождения неопределенных коэффициентов:


[math]\begin{cases}A+B=3,\\ 4A+C=3,\\ 13A=-13.\end{cases}[/math] Решая ее, получаем [math]A=-1,~B=4,~C=7[/math] и

[math]F(p)=-\frac{1}{p}+\frac{4p+7}{p^2+4p+13}=-\frac{1}{p}+ \frac{4(p+2)-1}{(p+2)^2+9}=-\frac{1}{p}+4\cdot \frac{p+2}{(p+2)^2+3^2}-\frac{1}{3}\cdot \frac{3}{(p+2)^2+3^2}\,.[/math]

По формулам 1,12,13 из табл. 5.1 [math]f(t)=-1+4e^{-2t}\cos3t-\frac{1}{3}\,e^{-2t}\sin3t[/math].


б) Представим [math]F(p)[/math] в виде [math]F(p)= \frac{p}{(p+2)^2(p-1)}= \frac{A}{p-1}+ \frac{B}{p+2}+ \frac{C}{(p+2)^2}[/math], где [math]A,\,B,\,C[/math] — неопределенные коэффициенты.


Отсюда [math]A(p+2)^2+ B(p-1)(p+2)+ C(p-1)=p[/math].


Подставляя последовательно [math]p=1,~p=-2,~ p=0[/math], получаем [math]A=\frac{1}{9},~ B=\frac{2}{3},~ C=-\frac{1}{9}[/math] и поэтому


[math]F(p)=\frac{1}{9}\cdot \frac{1}{p-1}-\frac{1}{9}\cdot \frac{1}{p+2}+ \frac{2}{3}\cdot \frac{1}{(p+2)^2}\,.[/math]

По формулам 6,7 из табл. 5.1 находим [math]\frac{1}{9}\,e^{t}-\frac{1}{9}\,e^{-2t}+ \frac{2}{3}\,t\,e^{-2t}.[/math].


в) Представим изображение в виде [math]F(p)=\frac{p^2+4}{(p-1)(p+2)(p-3)}= \frac{A}{p-1}+ \frac{B}{p+2}+ \frac{C}{p-3}[/math]. Отсюда


[math]A(p+2)(p-3)+ B(p-1)(p-3)+ C(p-1)(p+2)= p^2+4.[/math]

При [math]p=1,~p=-2,~p=3[/math] получаем [math]A=-\frac{5}{6},~ B=\frac{8}{15},~ C=\frac{13}{10}[/math], поэтому


[math]F(p)=-\frac{5}{6}\cdot \frac{1}{p-1}+ \frac{8}{15}\cdot \frac{1}{p+2}+ \frac{13}{10}\cdot \frac{1}{p-3}\,.[/math]

По свойству линейности и по формуле 6 из табл. 5.1 получаем


[math]f(t)=-\frac{5}{6}\,e^{t}+ \frac{8}{15}\,e^{-2t}+ \frac{13}{10}\,e^{3t}.[/math]

Можно решить эту задачу иначе, используя вторую теорему разложения (см. п. "б" примера 5.24):


[math]\begin{gathered}F(p)= \frac{p^2+4}{(p-1)(p+2)(p-3)}\\ \Downarrow\\ \begin{aligned}f(t)&= \mathop{\operatorname{res}}\limits_{p=1} F(p)e^{pt}+ \mathop{\operatorname{res}}\limits_{p=-2} F(p)e^{pt}+ \mathop{\operatorname{res}}\limits_{p=3} F(p)e^{pt}=\\ &= \frac{5e^t}{3\cdot(-2)}+ \frac{8e^{-2t}}{(-3)\cdot(-5)}+ \frac{13e^{3t}}{2\cdot5}=-\frac{5}{6}\,e^t+ \frac{8}{15}\,e^{-2t}+ \frac{13}{10}\,e^{3t}.\end{aligned}\end{gathered}[/math]

г) Представим [math]F(p)[/math] в виде [math]F(p)= \frac{3p}{2p^2-2p-4}= \frac{3p}{2(p-2)(p+1)}= \frac{A}{p-2}+ \frac{B}{p+1}[/math], где [math]A,\,B[/math] — неопределенные коэффициенты.


Из равенства [math]A(p+1)+B(p-2)=\frac{3p}{2}[/math] при [math]p=-1,~p=2[/math] получаем [math]A=1,~ B=\frac{1}{2}[/math], поэтому [math]F(p)=\frac{1}{p-2}+ \frac{1}{2}\cdot \frac{1}{p+1}[/math].


По формуле 6 из табл. 5.1 имеем [math]f(t)=e^{2t}+\frac{1}{2}\,e^{-t}[/math].


Можно также решить эту задачу иначе, используя вторую теорему разложения (см. п. "б" примера 5.24 и п. "в" данного примера):


[math]\begin{gathered}F(p)= \frac{3p}{2(p^2-p-2)}= \frac{3p}{2(p+1)(p-2)}\\ \Downarrow\\ f(t)= \mathop{\operatorname{res}}\limits_{p=-1} F(p)e^{pt}+ \mathop{\operatorname{res}}\limits_{p=2} F(p)e^{pt}= \frac{-3e^{-t}}{2\cdot(-3)}+ \frac{6e^{2t}}{2\cdot3}= e^{2t}+\frac{1}{2}\,e^{-t}. \end{gathered}[/math]

Пример 5.29. Найти оригиналы для функций: a) [math]F(p)= \frac{1}{p^2(p-1)^2}[/math]; б) [math]F(p)= \frac{p}{(p^2+1)(p^2+4)}[/math].


▼ Решение

а) Решим пример различными способами.


Первый способ. Воспользуемся разложением дроби на элементарные:


[math]F(p)= \frac{1}{p^2(p-1)^2}= \frac{2}{p}+ \frac{1}{p^2}-\frac{2}{p-1}+ \frac{1}{(p-1)^2}\,.[/math]

По формулам 2,3,6,7 из табл. 5.1 получаем [math]f(t)=2+t-2e^t+te^t[/math].


Второй способ. Применим вторую теорему разложения, учитывая, что [math]p_1=0[/math] и [math]p_2=1[/math] — полюсы второго порядка функции [math]F(p)\colon[/math]


[math]\begin{aligned} f(t)= \mathop{\operatorname{res}}\limits_{p=0}F(p)e^{pt}+ \mathop{\operatorname{res}}\limits_{p=1}F(p)e^{pt}= \lim\limits_{p\to0} \frac{d}{dp}\! \left[\frac{p^2e^{pt}}{p^2(p-1)^2}\right]+ \lim\limits_{p\to1} \frac{d}{dp}\! \left[\frac{(p-1)^2e^{pt}}{p^2(p-1)^2}\right]= \ldots= t+2+te^t-2e^t.\end{aligned}[/math]

Третий способ. Обозначим [math]F_1(p)= \frac{1}{(p-1)^2}[/math]. Тогда [math]f_1(t)=t\,e^t[/math]. Рассмотрим функцию [math]F_2(p)= \frac{1}{p}F_1(p)[/math]. По свойству интефирования оригинала (формула (5.12)) получаем


[math]f_2(t)= \int\limits_{0}^{t}\tau\,e^{\tau}\,d\tau= \Bigl.{\tau\,e^{\tau}}\Bigr|_{0}^{t}-\Bigl.{e^{\tau}}\Bigr|_{0}^{t}= t\,e^t-e^t+1.[/math]

Заметим, что [math]F(p)= \frac{1}{p}F_2(p)[/math]. Применяя еще раз свойство интегрирования оригинала, имеем


[math]f(t)= \int\limits_{0}^{t}(\tau\,e^{\tau}-e^{\tau}+1)d\tau= \ldots= t\,e^t-2e^t+2+t.[/math]

Четвертый способ. Представим изображение в виде произведения


[math]F(p)= \frac{1}{p^2(p-1)^2}= \frac{1}{p^2}\cdot \frac{1}{(p-1)^2}= F_1(p)\cdot F_2(p)[/math], где [math]F_1(p)=\frac{1}{p^2},~~ F_2(p)=\frac{1}{(p-1)^2}[/math].

По формулам 3 и 7 из табл. 5.1 [math]f_1(t)=t,~ f_2(t)=t\,e^t[/math]. Далее по теореме Бореля (формула (5.15))


[math]f(t)= f_1(t)\ast f_2(t)= \int\limits_{0}^{t}\tau\,e^{\tau}(t-\tau)d\tau= \ldots= t\,e^t-2e^t+2+t.[/math]

б) Решим пример также несколькими способами.


Первый способ. Воспользуемся разложением дроби на элементарные:


[math]F(p)= \frac{p}{(p^2+1)(p^2+4)}= \frac{p[p^2+4-(p^2+1)]}{3(p^2+1)(p^2+4)}= \frac{1}{3}\cdot \frac{p}{p^2+1}-\frac{1}{3}\cdot \frac{p}{p^2+4}\,.[/math]

По формуле 9 из табл. 5.1 получаем [math]f(t)= \frac{1}{3}\cos t-\frac{1}{3}\cos2t[/math].


Второй способ. Применим вторую теорему разложения с учетом пп. 2,3 замечаний 5.5:


[math]f(t)= 2\operatorname{Re} \mathop{\operatorname{res}}\limits_{p=i} \frac{p\,e^{pt}}{(p^2+1)(p^2+4)}+ 2\operatorname{Re} \mathop{\operatorname{res}}\limits_{p=2i} \frac{p\,e^{pt}}{(p^2+1)(p^2+4)}= \ldots= \frac{1}{3}\cos t-\frac{1}{3}\cos2t.[/math]

Третий способ. Представим изображение в виде произведения:


[math]F(p)= \frac{p}{(p^2+1)(p^2+4)}= \frac{p}{p^2+4}\cdot \frac{1}{p^2+1}= F_1(p)\cdot F_2(p).[/math]

Отсюда [math]f_1(t)= L^{-1}\! \left[\frac{p}{p^2+4}]=\cos2t,~ f_2(t)= L^{-1}\! \left[\frac{1}{p^2+1}]= \sin t[/math]. По теореме Бореля


[math]f(t)=f_1(t)\ast f_2(t)= \int\limits_{0}^{t} \cos2\tau\sin(t-\tau)d\tau= \ldots= \frac{1}{3}\cos t-\frac{1}{3}\cos2t.[/math]

Четвертый способ. Используем формулу 37 из табл. 5.1. При [math]a=1,~b=2[/math] получаем


[math]f(t)=\frac{\cos2t-\cos t}{1-4}= \frac{1}{3}\cos t-\frac{1}{3}\cos2t.[/math]

Часовой пояс: UTC + 4 часа [ Летнее время ]


Яндекс.Метрика

Copyright © 2010-2016 MathHelpPlanet.com. All rights reserved