"ref_1.gif"

"ref_2.gif"

      Псевдодифференциальный оператор.

Операторный метод нахождения  частного решения линейного дифференциального уравнения с полиномиальной правой частью в символьном виде.

  Линейный дифференциальный оператор с постоянными коэффициентами 5-го порядка.

"PseudoD2_1.gif"

"PseudoD2_2.gif"

   Полином в правой части уравнения 4-го порядка.

"PseudoD2_3.gif"

"PseudoD2_4.gif"

  Пусть требуется найти частное решение уравнения

"PseudoD2_5.gif"

"PseudoD2_6.gif"

Построение соответственного характеристического полинома.

"PseudoD2_7.gif"

"PseudoD2_8.gif"

"PseudoD2_9.gif"

"PseudoD2_10.gif"

"PseudoD2_11.gif"

Найдём решение, применяя обратный оператор "PseudoD2_12.gif" к правой части  R

"PseudoD2_13.gif"

"PseudoD2_14.gif"

"PseudoD2_15.gif"

"PseudoD2_16.gif"

Проверяем, удовлетворяет ли полученное решение нашему уравнению.

"PseudoD2_17.gif"

"PseudoD2_18.gif"

  Очевидно, выполняется тождество

"PseudoD2_19.gif"

"PseudoD2_20.gif"

  Выполняется также тождество с операторами, действующими в обратном порядке.

"PseudoD2_21.gif"

"PseudoD2_22.gif"

Такой же пример, но в численном виде.

  Линейный дифференциальный оператор с постоянными коэффициентами 3-го порядка.

"PseudoD2_23.gif"

"PseudoD2_24.gif"

"PseudoD2_25.gif"

"PseudoD2_26.gif"

  Полином в правой части уравнения 10-го порядка.

"PseudoD2_27.gif"

"PseudoD2_28.gif"

  Требуется найти частное решение уравнения

"PseudoD2_29.gif"

"PseudoD2_30.gif"

Построение соответственного характеристического полинома.

"PseudoD2_31.gif"

"PseudoD2_32.gif"

"PseudoD2_33.gif"

"PseudoD2_34.gif"

"PseudoD2_35.gif"

Найдём решение, применяя обратный оператор "PseudoD2_36.gif" к правой части  R

"PseudoD2_37.gif"

"PseudoD2_38.gif"

  Проверяем, удовлетворяет ли полученное решение нашему уравнению.

"PseudoD2_39.gif"

"PseudoD2_40.gif"

  Очевидно, выполняется тождество

"PseudoD2_41.gif"

"PseudoD2_42.gif"

Выполняется также тождество с операторами, действующими в обратном порядке.

"PseudoD2_43.gif"

"PseudoD2_44.gif"

Выясним, какой задаче Коши удовлетворяет полученное решение.

"PseudoD2_45.gif"

"PseudoD2_46.gif"

Это начальные условия Клши, отвечающее полученному решению.

Псевдодифференциальные уравнения с полиномиальной правой частью.

1. Рассмотрим псевдодифференциальное уравнение, где правая часть R(x) –полином 10-порядка.

"PseudoD2_47.gif"

"PseudoD2_48.gif"

"PseudoD2_49.gif"

Разложим функцию

"PseudoD2_50.gif"

в ряд Тейлора до 10-го  порядка включительно (остаток в форме "PseudoD2_51.gif") так как 11-ая производная R(x)равна нулю:

"PseudoD2_52.gif"

"PseudoD2_53.gif"

Заменим каждый член разложения "PseudoD2_54.gif" оператором дифференцирования "PseudoD2_55.gif" и применим полученный оператор к R.

"PseudoD2_56.gif"

"PseudoD2_57.gif"

Находим решение.

"PseudoD2_58.gif"

"PseudoD2_59.gif"

  Проверим правильность решения.

"PseudoD2_60.gif"

"PseudoD2_61.gif"

"PseudoD2_62.gif"

"PseudoD2_63.gif"

  Если вместо R взять экспоненту exp(x), то так как

"PseudoD2_64.gif"

а "PseudoD2_65.gif", то

"PseudoD2_66.gif"

2. Рассмотрим псевдодифференциальный оператор

"PseudoD2_67.gif"

Выясним, как он деёствует на произвольную функцию u(x). Положим

"PseudoD2_68.gif"

Дифференцируя функцию v(x,t)  по t, получим

"PseudoD2_69.gif"

Имеем Задачу Коши

"PseudoD2_70.gif"

Она легко решается, например, методом характермстик. Её решение

"PseudoD2_71.gif"

"PseudoD2_72.gif"

"PseudoD2_73.gif"

Таким образом, данный псевдодифференциальный оператор действует на функцию прибавляя к её аргументу t. Он называется оператором сдвига.   
Вообще, для некоторого специального класса функций, псевдодифференциальный оператор (функция от дифференциального оператора) определяется по формуле

"PseudoD2_74.gif"

  где "PseudoD2_75.gif"– преобразование Фурье функции f. Найдём решение псевдодифференциального уравнения (I), используя формулу (II)/
Находим преобразование функции  "PseudoD2_76.gif"

"PseudoD2_77.gif"

"PseudoD2_78.gif"

Применяем оператор (II) к правой части R(x). Учитывая правило, по которому деёствует оператор сдвига, получим

"PseudoD2_79.gif"

"PseudoD2_80.gif"

Сравниваем с ранее полученным результатом.

"PseudoD2_81.gif"

"PseudoD2_82.gif"

Если вместо R(x) взять "PseudoD2_83.gif", получим u(x) в качестве решения (I)

"PseudoD2_84.gif"

"PseudoD2_85.gif"

3. Рассмотрим псевдодифференциальный оператор

"PseudoD2_86.gif"

Это оператор интегрирования дробногй кратности (в данном случае "PseudoD2_87.gif"-кратное интегрирование).. Используя формулу (II) применим его к функции "PseudoD2_88.gif" . Сначала найдём Фурье-образ символа этого оператора.

"PseudoD2_89.gif"

"PseudoD2_90.gif"

  Теперь применяем орератор (II) к функции "PseudoD2_91.gif"

"PseudoD2_92.gif"

"PseudoD2_93.gif"

Это результат "PseudoD2_94.gif"-кратного интегрирования. Применим к полученной функции ещё раз оператор (II)

"PseudoD2_95.gif"

"PseudoD2_96.gif"

"PseudoD2_97.gif"

"PseudoD2_98.gif"

  Эта функция лишь на константу отличается от arctag(x) и является первообразной от  "PseudoD2_99.gif". Проверка дифференцированием:

"PseudoD2_100.gif"

"PseudoD2_101.gif"

"PseudoD2_102.gif"

"PseudoD2_103.gif"

4. Используя оператор сдвига можно найти производящую функцию для полиномов Эрмита, просуммировав ряд

"PseudoD2_104.gif"

Полиномы Эрмита можно определить по формуле

"PseudoD2_105.gif"

Подставляя это выражение в ряд, получим

"PseudoD2_106.gif"

"ref_1.gif"

"ref_2.gif"

Spikey Created with Wolfram Mathematica 7.0
Hosted by uCoz