Главная страница
Навигация по странице:

  • Определение. Функцией Хевисайда

  • Определение. Выражение называется сверткой

  • ду. ДУ много решенных. Обыкновенные дифференциальные уравнения. Основные определения


    Скачать 1.42 Mb.
    НазваниеОбыкновенные дифференциальные уравнения. Основные определения
    Дата19.07.2022
    Размер1.42 Mb.
    Формат файлаpdf
    Имя файлаДУ много решенных.pdf
    ТипРешение
    #633601
    страница16 из 19
    1   ...   11   12   13   14   15   16   17   18   19
    Определение.
    Функция F(p) называется изображением Лапласа функции f(t).
    Также функцию F(p) называют L – изображением или преобразованием Лапласа.
    Обозначается
    );
    (
    )
    (
    );
    (
    )
    (
    )};
    (
    {
    )
    (
    t
    f
    p
    F
    t
    f
    p
    F
    t
    f
    L
    p
    F




    =

    =
    При этом функция f(t) называется начальной функцией или оригиналом, а процесс нахождения оригинала по известному изображению называется операционным
    исчислением.
    Теорема.
    (Теорема единственности) Если две непрерывнные функции f(x) и g(x)
    имеют одно и то же L – изображение F(p), то они тождественно равны.
    103

    Ларин Александр Александрович “Курс высшей математики. Часть 3.”
    Определение.
    Функцией Хевисайда (Оливер Хевисайд (1850 – 1925) – английский физик) называется функция



    <

    =
    σ
    0
    ,
    0 0
    ,
    1
    )
    (
    0
    t
    t
    t
    Свойства изображений.
    Если
    ) , то справедливы следующие свойства:
    (
    )
    (
    t
    f
    p
    F


    =
    1)
    Свойство подобия.
    ;
    0
    ;
    1
    )
    (
    >
    α






    α
    α
    =
    α


    p
    F
    t
    f
    2)
    Свойство линейности.
    )].
    (
    [
    )]
    (
    [
    )]
    (
    )
    (
    [
    t
    g
    BL
    t
    f
    AL
    t
    Bg
    t
    Af
    L
    +
    =
    +
    3)
    Смещение изображения.
    )
    (
    )
    (
    α
    +
    =


    α

    p
    F
    e
    t
    f
    t
    4)
    Дифференцирование изображения.
    )
    (
    )
    (
    )
    1
    (
    t
    f
    t
    p
    F
    dp
    d
    n
    n
    n
    n


    =

    5)
    Дифференцирование оригинала.
    )
    (
    )
    0
    (
    )
    (
    t
    f
    f
    p
    pF

    =



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




    =
    p
    dq
    q
    F
    t
    t
    f
    )
    (
    )
    (
    (Справедливо при условии, что интеграл сходится)
    7)
    Интегрирование оригинала.
    p
    p
    F
    d
    f
    t
    )
    (
    )
    (
    0


    =
    τ
    τ

    Таблица изображений некоторых функций.
    Для большинства функций изображение находится непосредственным интегрированием.
    Пример. Найти изображение функции f(t) = sint.
    104

    Ларин Александр Александрович “Курс высшей математики. Часть 3.”
    =


    =
    ⎪⎭



    ⎪⎩




    =

    =
    =
    =
    =
    =










    dt
    t
    pe
    t
    e
    t
    v
    dt
    pe
    du
    tdt
    dv
    e
    u
    tdt
    e
    p
    F
    pt
    pt
    pt
    pt
    pt
    0 0
    0
    cos cos
    ;
    cos
    ;
    ;
    sin
    ;
    sin
    )
    (


    =
    ⎪⎭



    ⎪⎩



    =

    =
    =
    =
    =

    =







    0 0
    sin
    1
    ;
    sin
    ;
    ;
    cos
    ;
    cos
    1
    t
    pe
    t
    v
    dt
    pe
    du
    tdt
    dv
    e
    u
    tdt
    e
    p
    pt
    pt
    pt
    pt



    0 2
    sin tdt
    e
    p
    pt
    1
    sin
    )
    1
    (
    0 2
    =
    +



    tdt
    e
    p
    pt
    ;
    1 1
    sin
    2 0
    p
    tdt
    e
    pt
    +
    =



    ;
    1 1
    sin
    2
    p
    t
    +
    =


    Для многих функций изображения посчитаны и приведены в соответствующих таблицах.
    № f(t)
    F(p) № f(t)
    F(p)
    1 1
    p
    1 9
    n
    t
    1
    !
    +
    n
    p
    n
    2 sin
    αt
    2 2
    α
    +
    α
    p
    10
    at
    t sin
    2 2
    2
    )
    (
    2
    a
    p
    pa
    +
    3 cos
    αt
    2 2
    α
    +
    p
    p
    11
    at
    t cos
    2 2
    2 2
    2
    )
    (
    a
    p
    p
    a
    +


    4
    e
    -
    α
    t
    α
    +
    p
    1 12
    t
    te
    α

    2
    )
    (
    1
    α
    +
    p
    5 sh
    αt
    2 2
    α

    α
    p
    13
    )
    cos
    (sin
    2 1
    3
    at
    at
    at
    a


    2 2
    2
    )
    (
    1
    a
    p
    +
    6 ch
    αt
    2 2
    α

    p
    p
    14
    )
    (t
    f
    t
    n
    )
    (
    )
    1
    (
    p
    F
    dp
    d
    n
    n
    n

    7
    at
    e
    t
    sin
    α

    2 2
    )
    (
    a
    p
    a
    +
    α
    +
    15

    τ
    τ

    τ
    t
    d
    t
    f
    f
    0 2
    1
    )
    (
    )
    (
    )
    (
    )
    (
    2 1
    p
    F
    p
    F
    8
    at
    e
    t
    cos
    α

    2 2
    )
    (
    a
    p
    p
    +
    α
    +
    α
    +
    16
    )
    (
    )
    (
    t
    f
    n
    )
    ( p
    F
    p
    n
    *
    * - при условии, что
    0
    )
    0
    (
    )
    0
    (
    )
    0
    (
    )
    1
    (
    =
    =
    =

    =

    n
    f
    f
    f
    Теоремы свертки и запаздывания.
    Теорема.
    (теорема запаздывания) Если f(t) = 0 при t < 0, то справедлива
    формула
    )]
    (
    [
    )]
    (
    [
    0 0
    t
    f
    L
    e
    t
    t
    f
    L
    pt

    =

    где t
    0
    – некоторая точка.
    105

    Ларин Александр Александрович “Курс высшей математики. Часть 3.”
    Определение.
    Выражение называется сверткой функций f
    1
    (t)
    f
    2
    (t) и обозначается f
    1

    f
    2

    τ
    τ

    τ
    t
    d
    t
    f
    f
    0 2
    1
    )
    (
    )
    (
    и
    Теорема.
    (теорема свертки) Преобразование Лапласа от свертки равно
    оизведению преобразований Лапласа от функций f
    1
    (t) и f
    2
    (t) .
    пр

    τ
    τ

    τ
    =


    t
    d
    t
    f
    f
    p
    F
    p
    F
    0 2
    1 2
    1
    )
    (
    )
    (
    )
    (
    )
    (
    Теорема.
    (
    (1797 – 1872) –
    Интеграл Дюамеля (Дюамель французский атематик)). Если то верно равенство осредственным нтегрированием применяются приведенные выще теоремы и свойства.
    Пример.
    )
    (
    )
    (
    );
    (
    )
    (
    t
    g
    p
    G
    t
    f
    p
    F




    =
    =
    , м

    τ
    τ


    τ
    +
    =


    t
    d
    t
    g
    f
    g
    t
    f
    p
    G
    p
    pF
    0
    )
    (
    )
    (
    )
    0
    (
    )
    (
    )
    (
    )
    (
    Для нахождения изображений различных функций наряду с неп и
    t
    t
    sin
    Найти изображение функции
    1 1
    sin
    2
    +
    =


    p
    t
    Из таблицы изображений получаем:
    По свойству интегрирования изображения получаем:




    =
    p
    dq
    q
    F
    t
    t
    f
    )
    (
    )
    (
    ;
    2 1
    1
    sin
    2
    arctgp
    arctgq
    dq
    q
    t
    t
    p
    p

    π
    =
    =
    +
    =





    t
    2
    sin
    Пример.
    Найти изображение функции известна формула
    2 2
    cos
    1
    sin
    2
    t
    t

    =
    Из тригонометрии
    Тогда
    )
    4
    (
    2 2
    1
    ]
    2
    [cos
    2 1
    ]
    1
    [
    2 1
    ]
    2
    cos
    1
    [
    2 1
    sin
    2 2
    +

    =

    =

    =


    p
    p
    p
    t
    L
    L
    t
    L
    t
    =
    )
    4
    (
    2
    )
    4
    (
    2 4
    2 2
    2 2
    +
    =
    +

    +
    p
    p
    p
    p
    p
    p
    Операционное исчисление используется как для нахождения значений нтегралов, так и для решение дифференциальных уравнений. о линейное дифференциальное уравнение с постоянными коэффициентами. и
    Пусть дан
    )
    (
    )
    (
    )
    (
    )
    (
    0 1
    )
    (
    t
    f
    t
    x
    a
    t
    x
    a
    t
    x
    a
    n
    n
    =
    +

    +
    +
    106

    Ларин Александр Александрович “Курс высшей математики. Часть 3.”
    Требуется найти решение этого дифференциального уравнения, удовлетворяющее начальным условиям: и функция f(t) имеет (по теореме единственности) одно и то же изображение
    Лапласа.
    )
    0
    (
    ;
    )
    0
    (
    ;
    )
    0
    (
    )
    1
    (
    0
    )
    1
    (
    0 0


    =

    =

    =
    n
    n
    x
    x
    x
    x
    x
    x
    Если функция x(t) является решением этого дифференциального уравнения, то оно обращает исходное уравнение в тождество, значит функция, стоящая в левой части уравнения
    )]
    (
    [
    0
    k

    =
    t
    f
    L
    dt
    x
    d
    a
    L
    n
    k
    k
    k
    =






    Из теоремы о дифференцировании оригинала {
    } можно сделать вывод, что
    )
    (
    )
    0
    (
    )
    (
    t
    f
    f
    p
    pF

    =

    ).
    0
    (
    )
    0
    (
    )
    1
    (


    k
    x
    )
    0
    (
    ]
    [
    )
    2
    (
    1





    =






    k
    k
    k
    k
    k
    px
    x
    p
    x
    L
    p
    dt
    x
    d
    L
    Тогда
    ].
    [
    ]
    [
    0
    f
    L
    x
    L
    a
    dt
    n

    x
    d
    L
    a
    n
    n
    =
    +
    +





    Обозначим
    ).
    (
    ]
    [
    ),
    (
    ]
    [
    p
    F
    f
    L
    p
    x
    x
    L
    =
    =
    +
    +
    +

    +
    =
    +
    +
    +
    +





    ]
    [
    ]
    )[
    (
    )
    1
    (
    0 0
    2 0
    1 0
    1 1
    1
    n
    n
    n
    n
    n
    n
    n
    n
    x
    x
    p
    x
    p
    a
    a
    p
    a
    p
    a
    p
    a
    p
    x
    Получаем: е называется вспомогательным (изображающим) или операторным
    Отсюда получаем изображение
    ).
    (
    ]
    [
    ]
    [
    0 1
    0 0
    2
    )
    2
    (
    0 0
    3 0
    2 1
    p
    F
    x
    a
    x
    px
    a
    x
    x
    p
    x
    p
    a
    n
    n
    n
    n
    +
    +

    +
    +
    +
    +
    +

    +
    +




    Это уравнени
    уравнением.
    )
    ( p
    x
    , а по нему и искомую функцию x(t).
    Изображение получаем в виде:
    )
    (
    )
    (
    )
    (
    )
    (
    )
    (
    1
    p
    R
    p
    p
    R
    p
    F
    p
    x
    n
    n
    n

    Ψ
    +
    =
    Где
    . Если эти условия нулевые, то многочлен равен нулю, и формула принимает вид:
    ;
    )
    (
    0 1
    1 1
    a
    p
    a
    p
    a
    p
    a
    p
    R
    n
    n
    n
    n
    n
    +
    +
    +
    +
    =


    )
    )
    1
    (
    0


    n
    (
    )
    (
    )
    (
    )
    (
    )
    2
    (
    0 0
    2 0
    1 0
    0 0
    2 3
    0 0
    2 0
    1 1




    +
    +
    +
    +
    +
    +
    ′′
    +

    +
    +

    +
    +
    =
    Ψ
    n
    n
    n
    n
    n
    x
    px
    x
    p
    x
    p
    a
    x
    x
    p
    x
    p
    a
    x
    px
    a
    x
    a
    p
    Этот многочлен зависит от начальных условий
    )
    ( p
    R
    n
    )
    (
    )
    (
    p
    F
    p
    x
    =
    ассмотрим применение этого метода на примерах.
    Пример.
    Р
    0
    )
    0
    (
    )
    0
    (
    ;
    2 4
    =

    =
    =
    +
    ′′
    y
    y
    y
    y
    Решить уравнение
    Изображение искомой функции будем искать в виде:
    )
    (
    )
    (
    p
    R
    p
    F
    y
    n
    =
    107

    Ларин Александр Александрович “Курс высшей математики. Часть 3.”
    4 4
    0 1
    )
    (
    ;
    2
    ]
    2
    [
    ]
    [
    )
    (
    2 2
    +
    =
    +

    +

    =
    =
    =
    =
    p
    p
    p
    p
    R
    p
    L
    f
    L
    p
    F
    n






    +

    =
    +
    =
    4 1
    2 1
    )
    4
    (
    2 2
    2
    p
    p
    p
    p
    p
    y
    Находим оригинал, т.е. искомую функцию:
    )
    2
    cos
    1
    (
    2 1
    x
    y
    y

    =
    =


    Пример. Решить уравнение
    1
    )
    0
    (
    ;
    0 2
    =
    =


    y
    y
    y
    ;
    1
    ;
    2
    )
    (
    ;
    0
    ]
    0
    [
    ]
    [
    )
    (
    0 1
    1
    =
    =
    Ψ

    =
    =
    =
    =

    y
    a
    p
    p
    R
    L
    f
    L
    p
    F
    n
    n
    ;
    2 1

    =
    p
    y
    ;
    2x
    e
    y
    y
    =
    =


    Пример. Решить уравнение:
    ;
    0
    )
    0
    (
    ;
    1
    )
    0
    (
    ;
    0
    )
    0
    (
    ;
    0 6
    11 6
    =
    ′′
    =

    =
    =


    +
    ′′

    ′′′
    y
    y
    y
    y
    y
    y
    y
    ;
    6 11 6
    )
    (
    ;
    0
    ]
    0
    [
    ]
    [
    )
    (
    2 3

    +

    =
    =
    =
    =
    p
    p
    p
    p
    R
    L
    f
    L
    p
    F
    n
    6
    )
    (
    )
    (
    )
    (
    0 0
    0 2
    3 0
    0 2
    0 1
    1
    p
    y
    y
    p
    y
    p
    a
    y
    py
    a
    y
    a
    p
    n
    +

    =
    ′′
    +

    +
    +

    +
    +
    =
    Ψ

    Изображение искомой функции
    6 11 6
    6 2
    3

    +

    +

    =
    p
    p
    p
    p
    y
    Для нахождения оригинала необходимо разложить полученную дробь на элементарные дроби. Воспользуемся делением многочленов (знаменатель делится без остатка на p – 1): p
    3
    – 6p
    2
    + 11p – 6 p - 1 p
    3
    – p
    2
    p
    2
    – 5p + 6
    -5p
    2
    + 11p
    -5p
    2
    + 5p
    6p - 6 6p - 6 0
    В свою очередь
    )
    3
    )(
    2
    (
    6 5
    2


    =
    +

    p
    p
    p
    p
    Получаем:
    ).
    3
    )(
    2
    )(
    1
    (
    6 11 6
    2 3



    =

    +

    p
    p
    p
    p
    p
    p
    Тогда:
    ;
    3 2
    1 6
    11 6
    6 2
    3

    +

    +

    =

    +

    +

    =
    p
    C
    p
    B
    p
    A
    p
    p
    p
    p
    y
    Определим коэффициенты А, В и С.
    p
    p
    p
    C
    p
    p
    B
    p
    p
    A
    +

    =


    +


    +


    6
    )
    2
    )(
    1
    (
    )
    3
    )(
    1
    (
    )
    3
    )(
    2
    (
    p
    C
    Cp
    Cp
    B
    Bp
    Bp
    A
    Ap
    Ap
    +

    =
    +

    +
    +

    +
    +

    6 2
    3 3
    4 6
    5 2
    2 2
    p
    C
    B
    A
    C
    B
    A
    p
    C
    B
    A
    p
    +

    =
    +
    +
    +
    +
    +

    +
    +
    6 2
    3 6
    )
    3 4
    5
    (
    )
    (
    2 108

    Ларин Александр Александрович “Курс высшей математики. Часть 3.”








    =
    =

    =






    =



    =


    =






    =
    +

    =
    +


    =






    =
    +
    +

    =
    +
    +
    =
    +
    +
    2 3
    4 2
    5 6
    1 2
    2 1
    6 4
    1 2
    6 2
    3 6
    1 3
    4 5
    0
    C
    B
    A
    A
    A
    B
    B
    A
    C
    B
    A
    B
    A
    B
    A
    C
    C
    B
    A
    C
    B
    A
    C
    B
    A
    Тогда
    ;
    3 2
    3 2
    4 1
    2 5
    6 11 6
    6 11 2
    3


    +

    +


    =

    +


    =
    p
    p
    p
    p
    p
    p
    p
    y
    ;
    2 3
    4 2
    5 3
    2
    x
    x
    x
    e
    e
    e
    y
    y

    +

    =
    =


    Приемы операционного исчисления можно также использовать для решения систем дифференциальных уравнений.
    Пример. Решить систему уравнений:
    ;
    3 4
    4 3




    =

    +
    =

    y
    x
    y
    y
    x
    x
    1
    )
    0
    (
    )
    0
    (
    =
    = y
    x
    Обозначим
    )
    (
    ),
    (
    p
    y
    p
    x
    - изображения искомых функций и решим вспомогательные уравнения:




    =

    +
    =





    =

    +
    =

    )
    (
    3
    )
    (
    4
    )
    0
    (
    )
    (
    )
    (
    4
    )
    (
    3
    )
    0
    (
    )
    (
    ;
    ]
    [
    3
    ]
    [
    4
    ]
    [
    ]
    [
    4
    ]
    [
    3
    ]
    [
    p
    y
    p
    x
    y
    p
    y
    p
    p
    y
    p
    x
    x
    p
    x
    p
    y
    L
    x
    L
    y
    L
    y
    L
    x
    L
    x
    L
    Решим полученную систему алгебраических уравнений.
    ;
    )
    3
    )(
    25
    (
    21 4
    )
    (
    25 1
    )
    (
    ;
    )
    (
    3 3
    1
    )
    (
    4 4
    1
    )
    (
    3 1
    )
    (
    4
    )
    (
    2 2
    2



    ⎪⎪





    +
    =

    +
    =



    ⎪⎪




    +
    =


    +
    =
    p
    p
    p
    p
    p
    x
    p
    p
    p
    y
    p
    y
    p
    p
    y
    p
    y
    p
    p
    p
    y
    p
    x
    ;
    25 7
    25 25 7
    )
    3
    )(
    25
    (
    )
    3
    )(
    7
    (
    )
    (
    2 2
    2 2

    +

    =

    +
    =



    +
    =
    p
    p
    p
    p
    p
    p
    p
    p
    p
    p
    x
    ;
    5 5
    7 5
    )
    (
    )
    (
    t
    sh
    t
    ch
    t
    x
    p
    x
    +
    =
    =


    ;
    25 1
    25
    )
    (
    2 2

    +

    =
    p
    p
    p
    p
    y
    ;
    5 5
    1 5
    )
    (
    )
    (
    t
    sh
    t
    ch
    t
    y
    p
    y
    +
    =
    =


    Если применить к полученным результатам формулы
    109

    Ларин Александр Александрович “Курс высшей математики. Часть 3.”
    ;
    2
    ;
    2
    z
    z
    z
    z
    e
    e
    shz
    e
    e
    chz



    =
    +
    =
    то ответ можно представить в виде:
    ;
    5 2
    5 3
    5 1
    5 6
    5 5
    5 5



    ⎪⎪


    +
    =

    =


    t
    t
    t
    t
    e
    e
    y
    e
    e
    x
    Как видно, гиперболические функции в ответе могут быть легко заменены на показательные.
    Пример. Решить систему уравнений



    +
    =

    +
    =

    y
    x
    y
    y
    x
    x
    2 2
    2 5
    при x(0) = y(0) = 1
    Составим систему вспомогательных уравнений:
    ;
    )
    (
    2
    )
    (
    2
    )
    0
    (
    )
    (
    )
    (
    2
    )
    (
    5
    )
    0
    (
    )
    (
    ;
    ]
    [
    2
    ]
    [
    2
    ]
    [
    ]
    [
    2
    ]
    [
    5
    ]
    [



    +
    =

    +
    =




    +
    =

    +
    =

    p
    y
    p
    x
    y
    p
    y
    p
    p
    y
    p
    x
    x
    p
    x
    p
    y
    L
    x
    L
    y
    L
    y
    L
    x
    L
    x
    L
    ;
    )
    6
    )(
    1
    (
    )
    (
    )
    6
    )(
    1
    (
    3
    )
    (
    ;
    1
    )
    (
    2 5
    2
    )
    (
    4
    )
    (
    5 1
    )
    (
    2
    )
    (



    ⎪⎪




    =



    =



    ⎪⎪


    +
    +

    +
    =

    +
    =
    p
    p
    p
    p
    x
    p
    p
    p
    p
    y
    p
    y
    p
    p
    y
    p
    y
    p
    p
    p
    y
    p
    x
    ;
    5 3
    5 2
    6 1
    5 3
    1 1
    5 2
    6 1
    )
    (
    6
    t
    t
    e
    e
    p
    p
    p
    B
    p
    A
    p
    y
    +
    =

    +

    =

    +

    =


    ;
    5 6
    5 1
    6 1
    5 6
    1 1
    5 1
    6 1
    )
    (
    6
    t
    t
    e
    e
    p
    p
    p
    D
    p
    C
    p
    x
    +

    =

    +


    =

    +

    =


    Если обозначить
    ;
    5 3
    ;
    5 1
    2 1
    =

    =
    C
    C
    то из полученного частного решения системы можно записать и общее решение:
    ⎪⎩



    +

    =
    +
    =
    t
    t
    t
    t
    e
    C
    e
    C
    y
    e
    C
    e
    C
    x
    6 2
    1 6
    2 1
    2 2
    При рассмотрении нормальных систем дифференциальных уравнений этот пример был решен традиционным способом (См.
    Другой способ решения.
    ). Как видно, результаты совпадают.
    Отметим, что операторный способ решения систем дифференциальных уравнений применим к системам порядка выше первого, что очень важно, т.к. в этом случае применение других способов крайне затруднительно.
    110

    Ларин Александр Александрович “Курс высшей математики. Часть 3.”
    Криволинейные интегралы.
    1   ...   11   12   13   14   15   16   17   18   19


    написать администратору сайта