Разложение Лоуи

редактировать

При изучении дифференциальных уравнений разложение Лоуи разбивается на каждые линейные обыкновенное дифференциальное уравнение (ОДУ) на так называемые наибольшие полностью приводимые компоненты. Он был введен Альфредом Лоуи.

Решение дифференциальных уравнений - одно из самых важных подполей в математике. Особый интерес представляют решения в закрытой форме. Разбиение ОДУ на самые большие неприводимые компоненты сводит процесс решения исходного уравнения к решению неприводимых уравнений самого низкого возможного порядка. Эта процедура является алгоритмической, поэтому гарантируется наилучший ответ для решения приводимого уравнения. Подробное обсуждение можно найти в.

Результаты Лоуи были расширены до линейных дифференциальных уравнений в частных (PDE) с двумя независимыми переменными. Таким образом, стали доступны алгоритмические методы решения больших классов линейных pde.

Содержание
  • 1 Разложение линейных обыкновенных дифференциальных уравнений
  • 2 Основные факты из дифференциальной алгебры
  • 3 Разложение линейных дифференциальных уравнений в частных производных второго порядка на плоскости
  • 4 Разложение линейных pde порядка выше 2
  • 5 Ссылки
Разложение линейных обыкновенных дифференциальных уравнений

Пусть D ≡ ddx {\ displaystyle D \ Equiv {\ frac {d} {dx}}}{\ displaystyle D \ Equiv {\ frac {d} {dx}}} обозначает производная по переменная x {\ displaystyle x}x . Дифференциальный оператор порядка n {\ displaystyle n}n- это многочлен вида

L ≡ D n + a 1 D n - 1 + ⋯ + an - 1 D + an { \ Displaystyle L \ Equiv D ^ {n} + a_ {1} D ^ {n-1} + \ cdots + a_ {n-1} D + a_ {n}}{\ displaystyle L \ Equ D ^ {n} + a_ {1} D ^ {n-1} + \ cdots + a_ {n-1} D + a_ {n}}

где коэффициенты ai { \ displaystyle a_ {i}}a_ {i} , i = 1,…, n {\ displaystyle i = 1, \ ldots, n}i = 1, \ ldots, n взяты из некоторого функционального поля, базового поля L { \ Displaystyle L}L . Обычно это поле рациональных функций в переменной x {\ displaystyle x}x , т.е. ai ∈ Q (x) {\ displaystyle a_ {i} \ in {\ mathbb { Q}} (x)}{\ displaystyle a_ {i} \ in {\ mathbb {Q}} (x)} . Если y {\ displaystyle y}y является неопределенным с dydx ≠ 0 {\ displaystyle {\ frac {dy} {dx}} \ neq 0}{\ displaystyle {\ frac {dy} {dx}} \ neq 0} , L y {\ displaystyle Ly}{\ displaystyle Ly} становится дифференциальным многочленом, а L y = 0 {\ displaystyle Ly = 0}Ly = 0 является дифференциальным уравнением, соответствующим L {\ displaystyle L}L .

Оператор L {\ displaystyle L}L порядка n {\ displaystyle n}nназывается сокращаемым, если он может быть представлен как произведение двух операторы L 1 {\ displaystyle L_ {1}}L_{1}и L 2 {\ displaystyle L_ {2}}L_ {2} , оба порядка ниже n {\ displaystyle n}n. Затем пишут L = L 1 L 2 {\ displaystyle L = L_ {1} L_ {2}}{\ displaystyle L = L_ {1} L_ {2}} , т.е. сопоставление означает произведение оператора, оно определяется правилом D ai = ai D + ai ′ {\ displaystyle Da_ {i} = a_ {i} D + a_ {i} '}{\displaystyle Da_{i}=a_{i}D+a_{i}'}; L 1 {\ displaystyle L_ {1}}L_{1}называется левым множитель L {\ displaystyle L}L , L 2 {\ displaystyle L_ {2}}L_ {2} правый множитель. По умолчанию предполагается, что область коэффициентов факторов является базовым полем L {\ displaystyle L}L , возможно, расширенным некоторыми алгебраическими числами, например Q ¯ (x) { \ displaystyle {\ bar {\ mathbb {Q}}} (x)}{\ displaystyle {\ bar {\ mathbb {Q}}} (x)} разрешено. Если оператор не допускает правого множителя, он называется неприводимым.

Для любых двух операторов L 1 {\ displaystyle L_ {1}}L_{1}и L 2 {\ displaystyle L_ {2}}L_ {2} наименьшее общее левое кратное L clm (L 1, L 2) {\ displaystyle Lclm (L_ {1}, L_ {2})}{\ displaystyle Lclm (L_ {1}, L_ {2})} - оператор самого низкого порядка, такой, что оба L 1 {\ displaystyle L_ {1}}L_{1}и L 2 {\ displaystyle L_ {2}}L_ {2} делят его справа. Наибольший общий правый делитель G crd (L 1, L 2) {\ displaystyle Gcrd (L_ {1}, L_ {2})}{\ displaystyle Gcrd (L_ {1}, L_ {2})} - оператор высшего порядка, который делит оба L 1 {\ displaystyle L_ {1}}L_{1}и L 2 {\ displaystyle L_ {2}}L_ {2} справа. Если оператор может быть представлен как L c l m {\ displaystyle Lclm}{\ displaystyle Lclm} неприводимых операторов, он называется полностью приводимым. По определению неприводимый оператор называется вполне приводимым.

Если оператор не является полностью сводимым, L c l m {\ displaystyle Lclm}{\ displaystyle Lclm} его неприводимых правых множителей разделяется, и та же процедура повторяется с частным. Из-за понижения порядка на каждом шаге эта процедура прекращается после конечного числа итераций и достигается желаемое разложение. На основании этих соображений Леви получил следующий фундаментальный результат.

Теорема 1 (Loewy 1906). Пусть D = ddx {\ displaystyle D = {\ frac {d} {dx}}}{\ displaystyle D = {\ frac {d} {dx}}} является производной и ai ∈ Q (x) {\ displaystyle a_ {i} \ in {\ mathbb {Q} (x)}}{ \ displaystyle a_ {i} \ in {\ mathbb {Q} (x)}} . Дифференциальный оператор

L ≡ D n + a 1 D n - 1 + ⋯ + an - 1 D + an {\ displaystyle L \ Equiv D ^ {n} + a_ {1} D ^ {n-1} + \ cdots + a_ {n-1} D + a_ {n}}{\ displaystyle L \ Equ D ^ {n} + a_ {1} D ^ {n-1} + \ cdots + a_ {n-1} D + a_ {n}}

порядка n {\ displaystyle n}nможно записать однозначно как произведение полностью сокращаемых множителей L k (dk) {\ displaystyle L_ {k} ^ {(d_ {k})}}{\ displaystyle L_ {k} ^ {(d_ {k})}} максимального порядка dk {\ displaystyle d_ {k}}d_ { k} сверх Q (x) {\ displaystyle {\ mathbb {Q}} (x)}{\ displaystyle {\ mathbb {Q}} (x)} в форме

L = L m (dm) L m - 1 (dm - 1)… L 1 (d 1) {\ displaystyle L = L_ {m} ^ {(d_ {m})} L_ {m-1} ^ {(d_ {m-1})} \ ldots L_ {1} ^ {( d_ {1})}}{\ displaystyle L = L_ {m} ^ {(d_ {m})} L_ {m-1} ^ {(d_ {m-1})} \ ldots L_ {1} ^ {(d_ {1})}}

с d 1 +… + dm = n {\ displaystyle d_ {1} + \ ldots + d_ {m} = n}{\ displaystyle d_ {1} + \ ldots + d_ {m} = n} . Факторы L k (d k) {\ displaystyle L_ {k} ^ {(d_ {k})}}{\ displaystyle L_ {k} ^ {(d_ {k})}} уникальны. Любой множитель L k (dk) {\ displaystyle L_ {k} ^ {(d_ {k})}}{\ displaystyle L_ {k} ^ {(d_ {k})}} , k = 1,…, m {\ displaystyle k = 1, \ ldots, m}k = 1, \ ldots, m можно записать как

L k (dk) = L clm (lj 1 (e 1), lj 2 (e 2),…, ljk (ek)) {\ displaystyle L_ {k} ^ {(d_ {k})} = Lclm (l_ {j_ {1}} ^ {(e_ {1})}, l_ {j_ {2}} ^ {(e_ {2})}, \ ldots, l_ { j_ {k}} ^ {(e_ {k})})}{\ displaystyle L_ {k} ^ {(d_ {k})} = Lclm (l_ {j_ {1}} ^ {(e_ {1})}, l_ {j_ {2}} ^ {(е_ {2})}, \ ldots, l_ {j_ {k}} ^ {(e_ {k})})}

с e 1 + e 2 +… + ek = dk {\ displaystyle e_ {1} + e_ {2} + \ ldots + e_ {k} = d_ {k}}{\ displaystyle e_ {1} + e_ {2} + \ ldots + e_ {k} = d_ {k}} ; lji (ei) {\ displaystyle l_ {j_ {i}} ^ {(e_ {i})}}{\ displaystyle l_ {j_ {i}} ^ {(e_ {i})}} для i = 1,…, K {\ displaystyle i = 1, \ ldots, k}я = 1, \ ldots, k , обозначает неприводимый оператор порядка ei {\ displaystyle e_ {i}}e_ {i} over Q (x) {\ displaystyle {\ mathbb {Q}} (x)}{\ displaystyle {\ mathbb {Q}} (x)} .

Разложение, определенное в этой теореме, называется разложением Лоуи L {\ displaystyle L}L . Он предоставляет подробное описание функционального пространства, содержащего решение приводимого линейного дифференциального уравнения L y = 0 {\ displaystyle Ly = 0}Ly = 0 .

Для операторов фиксированного порядка возможные разложения Лоуи, различающиеся числом и порядок факторов может быть указан явно; некоторые факторы могут содержать параметры. Каждая альтернатива называется разновидностью разложения Лоуи. Полный ответ для n = 2 {\ displaystyle n = 2}n = 2 подробно описан в следующем следствии приведенной выше теоремы.

Следствие 1 Пусть L {\ displaystyle L}L быть оператором второго порядка. Его возможные разложения Леви обозначаются L 0 2,… L 3 2 {\ displaystyle {\ mathcal {L}} _ {0} ^ {2}, \ ldots {\ mathcal {L}} _ {3} ^ {2}}{\ displaystyle {\ mathcal {L}} _ {0 } ^ {2}, \ ldots {\ mathcal {L}} _ {3} ^ {2}} , их можно описать следующим образом; l (i) {\ displaystyle l ^ {(i)}}{\ displaystyle l ^ {(i)}} и lj (i) {\ displaystyle l_ {j} ^ {(i)}}{\ displaystyle l_ {j} ^ {(i)}} - неприводимые операторы порядка i {\ displaystyle i}i; C {\ displaystyle C}C- константа.

L 1 2: L = l 2 (1) l 1 (1); {\ displaystyle {\ mathcal {L}} _ {1} ^ {2}: L = l_ {2} ^ {(1)} l_ {1} ^ {(1)};}{\ displaystyle {\ mathcal {L}} _ {1} ^ {2}: L = l_ {2} ^ {(1)} l_ {1} ^ {(1)};}
L 2 2: L = L clm (l 2 (1), l 1 (1)); {\ displaystyle {\ mathcal {L}} _ {2} ^ {2}: L = Lclm (l_ {2} ^ {(1)}, l_ {1} ^ {(1)});}{\ displaystyle {\ mathcal {L}} _ {2} ^ {2}: L = Lclm (l_ {2} ^ {(1)}, l_ {1} ^ {(1)});}
L 3 2: L = L clm (l (1) (C)). {\ displaystyle {\ mathcal {L}} _ {3} ^ {2}: L = Lclm (l ^ {(1)} (C)).}{\ displaystyle {\ mathcal {L}} _ {3} ^ {2}: L = Lclm (l ^ {(1)} (C)).}

Тип разложения оператора - разложение L я 2 {\ displaystyle {\ mathcal {L}} _ {i} ^ {2}}{\ displaystyle {\ mathcal {L}} _ {i} ^ {2}} с наивысшим значением i {\ displaystyle i}i. Неприводимый оператор второго порядка определяется как имеющий тип разложения L 0 2 {\ displaystyle {\ mathcal {L}} _ {0} ^ {2}}{\ displaystyle {\ mathcal {L}} _ {0} ^ {2}} .

Разложения L 0 2 {\ displaystyle {\ mathcal {L}} _ {0} ^ {2}}{\ displaystyle {\ mathcal {L}} _ {0} ^ {2}} , L 2 2 {\ displaystyle {\ mathcal {L}} _ {2} ^ {2}}{\ displaystyle {\ mathcal {L }} _ {2} ^ {2}} и L 3 2 {\ displaystyle {\ mathcal {L}} _ {3} ^ {2}}{\ displaystyle {\ mathcal {L}} _ {3} ^ {2}} полностью сводимы.

Если разложение типа L i 2 {\ displaystyle {\ mathcal {L}} _ {i} ^ {2}}{\ displaystyle {\ mathcal {L}} _ {i} ^ {2}} , i = 1, 2 {\ displaystyle i = 1, 2}i = 1,2 или 3 {\ displaystyle 3}3 было получено для уравнения второго порядка L y = 0 {\ displaystyle Ly = 0}Ly = 0 , фундаментальная система может быть указана явно.

Следствие 2 Пусть L {\ displaystyle L}L будет дифференциальным оператором второго порядка, D ≡ ddx {\ displaystyle D \ Equiv {\ frac {d} {dx}}}{\ displaystyle D \ Equiv {\ frac {d} {dx}}} , y {\ displaystyle y}y неопределенный дифференциал и ai ∈ Q (x) {\ displaystyle a_ {i} \ in {\ mathbb {Q}} (x)}{\ displaystyle a_ {i} \ in {\ mathbb {Q}} (x)} . Определим ε я (Икс) ≡ ехр ⁡ (- ∫ aidx) {\ Displaystyle \ varepsilon _ {i} (x) \ Equiv \ exp {(- \ int a_ {i} dx)}}{\ displaystyle \ varepsilon _ {я} (х) \ эквив \ ехр {(- \ int a_ {i} dx)}} для я = 1, 2 {\ displaystyle i = 1,2}i = 1,2 и ε (x, C) ≡ exp ⁡ (- ∫ a (C) dx) {\ displaystyle \ varepsilon (x, C) \ Equiv \ exp {(- \ int a (C) dx)}}{\ displaystyle \ varepsilon (x, C) \ Equiv \ exp {(- \ int a (C) dx)}} , C {\ displaystyle C}C- параметр; запрещенные количества C ¯ {\ displaystyle {\ bar {C}}}{\ bar {C}} и C ¯ ¯ {\ displaystyle {\ bar {\ bar {C}}}}{\ displaystyle {\ bar {\ bar {C}}}} - произвольные числа, C ¯ ≠ C ¯ ¯ {\ displaystyle {\ bar {C}} \ neq {\ bar {\ bar {C}}}}{\ displaystyle {\ bar {C}} \ neq {\ bar {\ bar {C}}}} . Для трех нетривиальных разложений следствия 1 следующие элементы y 1 {\ displaystyle y_ {1}}y_ {1} и y 2 {\ displaystyle y_ {2}}y_{2}фундаментальной системы.

L 1 2 {\ displaystyle {\ mathcal {L}} _ {1} ^ {2}}{\ displaystyle {\ mathcal {L}} _ {1} ^ {2}} : L y = (D + a 2) (D + a 1) y = 0 {\ displaystyle Ly Знак равно (D + a_ {2}) (D + a_ {1}) y = 0}{\ displaystyle Ly = (D + a_ {2}) (D + a_ {1}) y = 0} ;

y 1 = ε 1 (x), {\ displaystyle y_ {1} = \ varepsilon _ {1} (x),}{\ displaystyle y_ {1} = \ varepsilon _ {1} (x),}

y 2 = ε 1 (x) ∫ ε 2 (x) ε 1 (x) dx. {\ displaystyle y_ {2} = \ varepsilon _ {1} (x) \ int {\ frac {\ varepsilon _ {2} (x)} {\ varepsilon _ {1} (x)}} \, dx.}{\ displaystyle y_ {2} = \ varepsilon _ {1} (x) \ int {\ frac {\ varepsilon _ {2} (x)} {\ varepsilon _ {1} ( x)}} \, dx.}
L 2 2 {\ displaystyle {\ mathcal {L}} _ {2} ^ {2}}{\ displaystyle {\ mathcal {L }} _ {2} ^ {2}} : L y = L clm (D + a 2, D + a 1) y = 0; {\ Displaystyle Ly = Lclm (D + a_ {2}, D + a_ {1}) y = 0;}{\ displaystyle Ly = Lclm (D + a_ {2}, D + a_ {1}) y = 0; }
y i = ε i (x); {\ displaystyle y_ {i} = \ varepsilon _ {i} (x);}{\ displaystyle y_ {i} = \ varepsilon _ {i} (x);}

a 1 {\ displaystyle a_ {1}}a_ {1} не эквивалентно a 2 {\ displaystyle a_ {2}}a_ {2} .

L 3 2 {\ displaystyle {\ mathcal {L}} _ {3} ^ {2}}{\ displaystyle {\ mathcal {L}} _ {3} ^ {2}} : L y = L clm (D + a (C)) y = 0; {\ displaystyle Ly = Lclm (D + a (C)) y = 0;}{\ displaystyle Ly = Lclm (D + a (C)) y = 0;}
y 1 = ε (x, C ¯) {\ displaystyle y_ {1} = \ varepsilon (x, {\ bar {C) }})}{\ displaystyle y_ {1} = \ varepsilon (x, {\ bar {C}})}
y 2 = ε (x, C ¯ ¯). {\ displaystyle y_ {2} = \ varepsilon (x, {\ bar {\ bar {C}}}).}{\ displaystyle y_ {2} = \ varepsilon (x, {\ bar {\ bar {C}}}).}

Здесь две рациональные функции p, q ∈ Q (x) {\ displaystyle p, q \ in {\ mathbb {Q}} (x)}{\ displaystyle p, q \ in {\ mathbb {Q}} (x)} называются эквивалентными, если существует другая рациональная функция r ∈ Q (x) {\ displaystyle r \ in {\ mathbb {Q} } (x)}{\ displaystyle r \ in {\ mathbb {Q}} (x)} такой, что

p - q = r ′ r {\ displaystyle pq = {\ frac {r '} {r}}}{\displaystyle p-q={\frac {r'}{r}}}.

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

Пример 1 Уравнение 2.201 из коллекции Камке. имеет L 2 2 {\ displaystyle {\ mathcal {L}} _ {2} ^ {2}}{\ displaystyle {\ mathcal {L }} _ {2} ^ {2}} разложение

y ″ + (2 + 1 x) y ′ - 4 Икс 2 Y знак равно L clm (D + 2 x - 2 x - 2 x 2 - 2 x + 3 2, D + 2 + 2 x - 1 x + 3 2) y = 0. {\ displaystyle y '' + ( 2 + {\ frac {1} {x}}) y '- {\ frac {4} {x ^ {2}}} y = Lclm {\ Big (} D + {\ frac {2} {x}} - {\ frac {2x-2} {x ^ {2} -2x + {\ frac {3} {2}}}}, D + 2 + {\ frac {2} {x}} - {\ frac {1} {x + {\ frac {3} {2}}}} {\ Big)} y = 0.}{\displaystyle y''+(2+{\frac {1}{x}})y'-{\frac {4}{x^{2}}}y=Lclm{\Big (}D+{\frac {2}{x}}-{\frac {2x-2}{x^{2}-2x+{\frac {3}{2}}}},D+2+{\frac {2}{x}}-{\frac {1}{x+{\frac {3}{2}}}}{\Big)}y=0.}

Коэффициенты a 1 = 2 + 2 x - 1 x + 3 2 {\ displaystyle a_ { 1} = 2 + {\ frac {2} {x}} - {\ frac {1} {x + {\ frac {3} {2}}}}}{\ displaystyle a_ {1} = 2 + {\ frac {2} {x}} - {\ frac {1} {x + {\ frac {3} {2}}}} } и a 2 = 2 x - 2 x - 2 x 2 - 2 x + 3 2 {\ displaystyle a_ {2} = {\ frac {2} {x}} - {\ frac {2x-2} {x ^ {2} -2x + {\ frac {3} {2}}}}}{\ displaystyle a_ {2} = {\ frac {2} {x}} - {\ frac {2x-2} {x ^ {2} -2x + {\ frac {3} {2}}}}} - рациональные решения уравнения Риккати a ′ - a 2 + (2 + 1 x) + 4 x 2 = 0 {\ displaystyle a'-a ^ {2} + {\ big (} 2 + {\ frac {1} {x}} {\ big)} + {\ frac {4} {x ^ {2}}} = 0}{\displaystyle a'-a^{2}+{\big (}2+{\frac {1}{x}}{\big)}+{\frac {4}{x^{2}}}=0}, они дают фундаментальную систему

y 1 = 2 3 - 4 3 x + 1 x 2, {\ displaystyle y_ {1} = {\ frac {2} {3}} - {\ frac {4} {3x}} + {\ frac {1} {x ^ {2}}},}{\ displaystyle y_ {1} = {\ frac {2} {3}} - {\ frac {4} {3x}} + {\ frac {1} {x ^ {2}}},}
y 2 = 2 x + 3 x 2 e - 2 х. {\ displaystyle y_ {2} = {\ frac {2} {x}} + {\ frac {3} {x ^ {2}}} e ^ {- 2x}.}{\ displaystyle y_ {2} = {\ frac {2} {x}} + {\ frac {3} {x ^ {2}}} e ^ {- 2x}.}

Пример 2 An уравнение с типом L 3 2 {\ displaystyle {\ mathcal {L}} _ {3} ^ {2}}{\ displaystyle {\ mathcal {L}} _ {3} ^ {2}} разложение:

y ″ - 6 x 2 y = L clm (D + 2 x - 5 x 4 x 5 + C) y = 0. {\ displaystyle y '' - {\ frac {6} {x ^ {2}}} y = Lclm {\ big (} D + {\ frac {2} {x}} - {\ frac {5x ^ {4}} {x ^ {5} + C}} {\ big)} y = 0.}{\displaystyle y''-{\frac {6}{x^{2}}}y=Lclm{\big (}D+{\frac {2}{x}}-{\frac {5x^{4}}{x^{5}+C}}{\big)}y=0.}

Коэффициент множителя первого порядка является рациональным решением a ′ - a 2 + 6 x 2 = 0 {\ displaystyle a'-a ^ {2} + {\ frac {6} {x ^ {2}}} = 0}{\displaystyle a'-a^{2}+{\frac {6}{x^{2}}}=0}. После интеграции фундаментальная система y 1 = x 3 {\ displaystyle y_ {1} = x ^ {3}}{\ displaystyle y_ {1} = x ^ {3}} и y 2 = 1 x 2 {\ displaystyle y_ {2} = {\ frac {1} {x ^ {2}}}}{\ displaystyle y_ {2} = {\ frac {1} {x ^ {2}}}} для C = 0 {\ displaystyle C = 0}C = 0 и C → ∞ { \ displaystyle C \ rightarrow \ infty}{\ displaystyle C \ rightarrow \ infty} соответственно получается.

Эти результаты показывают, что факторизация обеспечивает алгоритмическую схему для решения приводимых линейных од. Всякий раз, когда уравнение порядка 2 факторизуется в соответствии с одним из типов, определенных выше, элементы фундаментальной системы явно известны, то есть факторизация эквивалентна ее решению.

Подобная схема может быть создана для линейных од любого порядка, хотя количество альтернатив значительно увеличивается с порядком; для порядка n = 3 {\ displaystyle n = 3}n = 3 подробный ответ дается в.

Если уравнение неприводимо, может оказаться, что его группа Галуа нетривиальна, то могут существовать алгебраические решения. Если группа Галуа тривиальна, можно выразить решения в терминах специальной функции, например, Функции Бесселя или Лежандра см. Или.

Основные факты из дифференциальной алгебры

Чтобы обобщить результат Лоуи на линейные pde, необходимо применить более общую установку дифференциальной алгебры. Поэтому ниже приведены несколько основных понятий, которые необходимы для этой цели.

Поле F {\ displaystyle {\ mathcal {F}}}{\ mathcal F} называется дифференциальным полем, если оно оснащено оператором вывода. Оператор δ {\ displaystyle \ delta}\ delta в поле F {\ displaystyle {\ mathcal {F}}}{\ mathcal F} называется оператором вывода, если δ (a + b) знак равно δ (a) + δ (b) {\ displaystyle \ delta (a + b) = \ delta (a) + \ delta (b)}{\ displaystyle \ delta (a + b) = \ delta (a) + \ delta (b)} и δ (ab) знак равно δ (a) b + a δ (b) {\ displaystyle \ delta (ab) = \ delta (a) b + a \ delta (b)}{\ displaystyle \ delta (ab) = \ delta (a) b + a \ delta (b)} для всех элементов a, b ∈ F {\ displaystyle a, b \ in {\ mathcal {F}}}{\ displaystyle a, b \ in {\ mathcal {F}} } . Поле с одним оператором вывода называется обыкновенным дифференциальным полем; если существует конечное множество, содержащее несколько коммутирующих операторов вывода, поле называется полем в частных производных.

Здесь дифференциальные операторы с производными ∂ x = ∂ ∂ x {\ displaystyle \ partial _ {x} = {\ frac {\ partial} {\ partial x}}}{\ displaystyle \ partial _ {x} = {\ frac {\ partial} {\ partial x}}} и ∂ y = ∂ ∂ y {\ displaystyle \ partial _ {y} = {\ frac {\ partial} {\ partial y}}}{\ displaystyle \ partial _ {y} = {\ frac {\ partial} {\ partial y}}} с коэффициентами из некоторого дифференциального поля. Его элементы имеют вид ∑ i, jri, j (x, y) ∂ xi ∂ yj {\ displaystyle \ sum _ {i, j} r_ {i, j} (x, y) \ partial _ {x } ^ {i} \ partial _ {y} ^ {j}}{\ displaystyle \ sum _ {i, j} r_ {i, j} (x, y) \ partial _ {x} ^ {i} \ partial _ {y} ^ {j}} ; почти все коэффициенты r i, j {\ displaystyle r_ {i, j}}{\ displaystyle r_ {i, j}} равны нулю. Поле коэффициентов называется базовым полем. Если конструктивные и алгоритмические методы являются основной проблемой, то это Q (x, y) {\ displaystyle {\ mathbb {Q}} (x, y)}{\ displaystyle {\ mathbb {Q}} (x, y)} . Соответствующее кольцо дифференциальных операторов обозначается D = Q (x, y) [∂ x, ∂ y] {\ displaystyle {\ mathcal {D}} = {\ mathbb {Q}} (x, y) [\ partial _ {x}, \ partial _ {y}]}{\ displaystyle {\ mathcal {D}} = {\ mathbb {Q}} (x, y) [\ partial _ {x}, \ partial _ {y}]} или D = F [∂ x, ∂ y] {\ displaystyle {\ mathcal {D}} = {\ mathcal {F}} [\ partial _ {x}, \ partial _ {y}]}{\ displaystyle {\ mathcal {D}} = {\ mathcal {F}} [\ partial _ {x}, \ partial _ {y}]} . Кольцо D {\ displaystyle {\ mathcal {D}}}{\ displaystyle {\ mathcal {D} }} некоммутативно, ∂ xa = a ∂ x + ∂ a ∂ x {\ displaystyle \ partial _ {x } a = a \ partial _ {x} + {\ frac {\ partial a} {\ partial x}}}{\ displaystyle \ partial _ {x} a = a \ partial _ {x} + {\ frac {\ partial a} {\ p artial x}}} и аналогично для других переменных; a {\ displaystyle a}a взято из базового поля.

Для оператора L = ∑ я + j ≤ nri, j (x, y) ∂ xi ∂ yj {\ displaystyle L = \ sum _ {i + j \ leq n} r_ {i, j} (x, y) \ partial _ {x} ^ {i} \ partial _ {y} ^ {j}}{\ displaystyle L = \ sum _ {i + j \ leq n} r_ {i, j} (x, y) \ partial _ {x} ^ {i } \ partial _ {y} ^ {j}} порядка n {\ displaystyle n}nсимвол L - однородный алгебраический многочлен symb (L) ≡ ∑ i + j = nri, j (x, y) X i Y j {\ displaystyle symb (L) \ Equiv \ sum _ {i + j = n} r_ {i, j} (x, y) X ^ {i} Y ^ {j}}{\ displaystyle symb (L) \ Equiv \ sum _ {i + j = n} r_ {i, j} (x, y) X ^ {i} Y ^ {j}} где X {\ displaystyle X}X и Y {\ displaystyle Y}Y алгебраические неопределенности.

Пусть I {\ displaystyle I}I будет левым идеалом, который порождается li ∈ D {\ displaystyle l_ {i} \ in {\ mathcal {D }}}{\ displaystyle l_ {i} \ in {\ mathcal {D}}} , i = 1,…, p {\ displaystyle i = 1, \ ldots, p}{\ displaystyle i = 1, \ ldots, p} . Затем пишут I = ⟨l 1,…, l p⟩ {\ displaystyle I = \ langle l_ {1}, \ ldots, l_ {p} \ rangle}{\ displaystyle I = \ langle l_ {1}, \ ldots, l_ {p} \ rangle } . Поскольку правильные идеалы здесь не рассматриваются, иногда I {\ displaystyle I}I просто называют идеалом.

Связь между левыми идеалами в D {\ displaystyle {\ mathcal {D}}}{\ displaystyle {\ mathcal {D} }} и системами линейных pde устанавливается следующим образом. Элементы li ∈ D {\ displaystyle l_ {i} \ in {\ mathcal {D}}}{\ displaystyle l_ {i} \ in {\ mathcal {D}}} применяются к одному дифференциальному неопределенному z {\ displaystyle z}z . Таким образом, идеал I = ⟨l 1, l 2,…⟩ {\ displaystyle I = \ langle l_ {1}, l_ {2}, \ ldots \ rangle}{\ displaystyle I = \ langle l_ {1}, l_ {2}, \ ldots \ rangle} соответствует система pde l 1 z = 0 {\ displaystyle l_ {1} z = 0}{\ displaystyle l_ {1} z = 0} , l 2 z = 0,… {\ displaystyle l_ {2} z = 0, \ ldots}{\ displaystyle l_ {2} z = 0, \ ldots} для единственной функции z {\ displaystyle z}z .

Генераторы идеала в высшей степени неуникальны; его члены могут быть преобразованы бесконечно многими способами, взяв их линейные комбинации или их производные без изменения идеала. Поэтому М. Джанет ввел нормальную форму для систем линейных pde (см. базис Жане ). Они являются дифференциальным аналогом базисов Грёбнера в коммутативной алгебре (которые первоначально были введены Бруно Бухбергером ); поэтому их также иногда называют дифференциальным базисом Грёбнера.

Для создания основы Джанет необходимо определить ранжирование производных финансовых инструментов. Это полный порядок, такой, что для любых производных δ {\ displaystyle \ delta}\ delta , δ 1 {\ displaystyle \ delta _ {1}}\ delta _ {1} и δ 2 {\ displaystyle \ delta _ {2}}\ delta _ {2} , и любой оператор вывода θ {\ displaystyle \ theta}\ theta отношения δ ⪯ θ δ {\ displaystyle \ delta \ prevq \ theta \ delta}{\ displaystyle \ delta \ prevq \ theta \ delta} и δ 1 ⪯ δ 2 → δ δ 1 ⪯ δ δ 2 {\ displaystyle \ delta _ {1} \ prevq \ delta _ {2} \ rightarrow \ delta \ delta _ {1} \ prevq \ delta \ delta _ {2}}{ \ displaystyle \ delta _ {1} \ prevq \ delta _ {2} \ rightarrow \ delta \ delta _ {1} \ prevq \ delta \ delta _ {2}} действительны. Здесь применяется градуированный лексикографический порядок терминов g r l e x {\ displaystyle grlex}{\ displaystyle grlex} . Для частных производных одной функции их определение аналогично мономиальным порядкам в коммутативной алгебре. S-пары в коммутативной алгебре соответствуют условиям интегрируемости.

Если есть уверенность, что образующие l 1,…, lp {\ displaystyle l_ {1}, \ ldots, l_ {p}}{\ displaystyle l_ {1}, \ ldots, l_ { p}} идеального I {\ displaystyle I}I образуют основу Джанет: I = ⟨⟨l 1,…, lp⟩⟩ {\ displaystyle I = {{\ big \ langle} {\ big \ langle} }} l_ {1}, \ ldots, l_ {p} {{\ big \ rangle} {\ big \ rangle}}}{\ displaystyle I = {{\ big \ langle} {\ big \ langle}} l_ {1}, \ ldots, l_ {p} {{\ big \ rangle} {\ big \ rangle}}} применяется.

Пример 3 Рассмотрим идеал

I = ⟨l 1 ≡ ∂ xx - 1 x ∂ x - yx (x + y) ∂ y, {\ displaystyle I = {\ Big \ langle} l_ { 1} \ Equiv \ partial _ {xx} - {\ frac {1} {x}} \ partial _ {x} - {\ frac {y} {x (x + y)}} \ partial _ {y}, }{\ displaystyle I = {\ Big \ langle} l_ {1} \ Equiv \ partial _ {xx} - {\ frac {1} {x}} \ partial _ {x} - {\ frac {y} {x (x + y)}} \ partial _ {y},} 

l 2 ≡ ∂ xy + 1 x + y ∂ y, {\ displaystyle l_ {2} \ Equiv \ partial _ {xy} + {\ frac {1} {x + y}} \ partial _ {y },}{\ displaystyle l_ {2} \ Equiv \ partial _ {xy} + {\ frac {1} {x + y}} \ partial _ {y},} l 3 ≡ ∂ yy + 1 x + y ∂ y⟩ {\ displaystyle l_ {3} \ Equiv \ partial _ {yy} + {\ frac {1} {x + y}} \ partial _ {y} {\ Big \ rangle}}{\ displaystyle l_ {3} \ Equiv \ partial _ {yy} + {\ frac {1} {x + y}} \ partial _ {y} {\ Big \ rangle}}

в grlex {\ displaystyle grlex}{\ displaystyle grlex} временной порядок с x ≻ y {\ displaystyle x \ succ y}x \ succ y . Его генераторы автовосстановлены. Если условие интегрируемости

l 1, y = l 2, x - l 2, y = y + 2 xx (x + y) ∂ xy + yx (x + y) ∂ yy {\ displaystyle l_ {1, y } = l_ {2, x} -l_ {2, y} = {\ frac {y + 2x} {x (x + y)}} \ partial _ {xy} + {\ frac {y} {x (x + y)}} \ partial _ {yy}}{\ displaystyle l_ {1, y} = l_ {2, x} -l_ {2, y} = {\ frac {y + 2x} {x (x + y)}} \ partial _ {xy} + {\ frac {y} {x (x + y)}} \ partial _ {y y}} 

сокращается относительно на I {\ displaystyle I}I получается новый генератор ∂ y {\ displaystyle \ partial _ {y}}\ partial _ { y} . Добавляя его к генераторам и выполняя все возможные редукции, данный идеал представляется как I = ⟨⟨∂ xx - 1 x ∂ x, ∂ y⟩⟩ {\ displaystyle I = {{\ Big \ langle} {\ Большой \ langle}} \ partial _ {xx} - {\ frac {1} {x}} \ partial _ {x}, \ partial _ {y} {{\ Big \ rangle} {\ Big \ rangle}}}{\ displaystyle I = {{\ Big \ langle} {\ Big \ langle} }} \ partial _ {xx} - {\ frac {1} {x}} \ partial _ {x}, \ partial _ {y} {{\ Big \ rangle} {\ Big \ rangle}}} . Его генераторы авторедуцируются и выполняется единственное условие интегрируемости, т. Е. Они образуют базис Жане.

Для любого идеала I {\ displaystyle I}I может оказаться, что он должным образом содержится в некотором более широком идеале J {\ displaystyle J}J с коэффициентами в базовом поле I {\ displaystyle I}I ; тогда J {\ displaystyle J}J называется делителем I {\ displaystyle I}I . Вообще говоря, дивизор в кольце дифференциальных операторов в частных производных может не быть главным.

Наибольший общий правый делитель (Gcrd) или сумма двух идеалов I {\ displaystyle I}I и J {\ displaystyle J}J - это наименьший идеал, обладающий тем свойством, что в нем содержатся как I {\ displaystyle I}I , так и J {\ displaystyle J}J . Если у них есть представление I ≡ ⟨f 1,…, fp⟩ {\ displaystyle I \ Equiv \ langle f_ {1}, \ ldots, f_ {p} \ rangle}{\ displaystyle I \ Equiv \ langle f_ {1}, \ ldots, f_ {p} \ rangle} и J ≡ ⟨g 1,…, gq⟩, {\ displaystyle J \ Equiv \ langle g_ {1}, \ ldots, g_ {q} \ rangle,}{\ displaystyle J \ Equiv \ langle g_ {1}, \ ldots, g_ {q} \ rangle,} fi {\ displaystyle f_ {i}}f_ {i} , gj ∈ D {\ displaystyle g_ {j} \ in {\ mathcal {D}}}{ \ displaystyle g_ {j} \ in {\ mathcal {D}}} для всех i {\ displaystyle i}iи j { \ displaystyle j}j , сумма генерируется объединением генераторов I {\ displaystyle I}I и J {\ displaystyle J}J . Пространство решений уравнений, соответствующих G c r d (I, J) {\ displaystyle Gcrd (I, J)}{\ displaystyle Gcrd (I, J)} , является пересечением пространств решений его аргументов.

Наименьшее общее левое кратное (Lclm) или левое пересечение двух идеалов I {\ displaystyle I}I и J {\ displaystyle J}J - наибольший идеал со свойством, что он содержится как в I {\ displaystyle I}I , так и в J {\ displaystyle J}J . Пространство решений L c l m (I, J) z = 0 {\ displaystyle Lclm (I, J) z = 0}{\ displaystyle Lclm (I, J) z = 0} - это наименьшее пространство, содержащее пространства решений его аргументов.

Особым видом делителя является так называемый делитель Лапласа заданного оператора L {\ displaystyle L}L , стр. 34. Он определяется следующим образом.

Определение Пусть L {\ displaystyle L}L будет оператором в частных производных на плоскости; определить

лм ≡ ∂ xm + am - 1 ∂ xm - 1 +… + a 1 ∂ x + a 0 {\ displaystyle {\ mathfrak {l}} _ {m} \ Equiv \ partial _ {x ^ {m }} + a_ {m-1} \ partial _ {x ^ {m-1}} + \ ldots + a_ {1} \ partial _ {x} + a_ {0}}{\ displaystyle {\ mathfrak { l}} _ {m} \ Equiv \ partial _ {x ^ {m}} + a_ {m-1} \ partial _ {x ^ {m-1}} + \ ldots + a_ {1} \ partial _ { х} + a_ {0}} и
kn ≡ ∂ yn + bn - 1 ∂ yn - 1 +… + b 1 ∂ y + b 0 {\ displaystyle {\ mathfrak {k}} _ {n} \ Equiv \ partial _ {y ^ {n}} + b_ {n-1} \ partial _ {y ^ {n-1}} + \ ldots + b_ {1} \ partial _ {y} + b_ {0}}{\ displaystyle {\ mathfrak {k}} _ {n} \ Equiv \ partial _ {y ^ {n}} + b_ {n-1} \ partial _ { y ^ {n-1}} + \ ldots + b_ {1} \ partial _ {y} + b_ {0}} 

быть обыкновенными дифференциальными операторами относительно x {\ displaystyle x}x или y {\ displaystyle y}y ; ai, bi ∈ Q (x, y) {\ displaystyle a_ {i}, b_ {i} \ in {\ mathbb {Q}} (x, y)}{\ displaystyle a_ {i}, b_ {i} \ in {\ mathbb {Q}} (x, y)} для всех i; m {\ displaystyle m}м и n {\ displaystyle n}n- натуральные числа не менее 2. Предположим, что коэффициенты ai {\ displaystyle a_ {i}}a_ {i} , i = 0,…, m - 1 {\ displaystyle i = 0, \ ldots, m-1}{\ displaystyle i = 0, \ ldots, m-1} таковы, что L {\ displaystyle L}L и lm {\ displaystyle {\ mathfrak {l}} _ {m}}{\ displaystyl е {\ mathfrak {l}} _ {m}} образуют основу Джанет. Если m {\ displaystyle m}м - наименьшее целое число с этим свойством, то L xm (L) ≡ ⟨⟨L, lm⟩⟩ {\ displaystyle {\ mathbb {L}} _ {x ^ {m}} (L) \ Equiv {\ langle \ langle} L, {\ mathfrak {l}} _ {m} {\ rangle \ rangle}}{\ displaystyle {\ mathbb {L}} _ {x ^ {m}} (L) \ Equiv {\ langle \ langle} L, {\ mathfrak {l}} _ {m} {\ rangle \ rangle}} называется делителем Лапласа из L {\ displaystyle L}L . Аналогично, если bj {\ displaystyle b_ {j}}b_ {j} , j = 0,…, n - 1 {\ displaystyle j = 0, \ ldots, n-1}{\ displaystyle j = 0, \ ldots, n -1} таковы, что L {\ displaystyle L}L и kn {\ displaystyle {\ mathfrak {k}} _ {n}}{\ displaystyle {\ mathfrak {k}} _ {n}} образуют базис Джанет и n {\ displaystyle n}nминимально, тогда L yn (L) ≡ ⟨⟨L, kn⟩⟩ {\ displaystyle {\ mathbb {L}} _ {y ^ {n}} ( L) \ Equiv {\ langle \ langle} L, {\ mathfrak {k}} _ {n} {\ rangle \ rangle}}{\ displaystyle {\ mathbb {L}} _ {y ^ {n}} (L) \ Equiv {\ langle \ langle} L, {\ mathfrak {k}} _ {n} {\ rangle \ rangle} } также называется делителем Лапласа L {\ displaystyle L}L .

Для существования делителя Лапласа коэффициенты оператора L {\ displaystyle L}L должны подчиняться определенным ограничениям. Алгоритм определения верхней границы для делителя Лапласа в настоящее время неизвестен, поэтому в общем случае существование делителя Лапласа может быть неразрешимым

Разложение линейных дифференциальных уравнений в частных производных второго порядка на плоскости

Применяя вышеупомянутые концепции, теорию Лоуи можно обобщить на линейные pde. Здесь он применяется к отдельным линейным pde второго порядка на плоскости с координатами x {\ displaystyle x}x и y {\ displaystyle y}y , а главные идеалы, порожденные соответствующими операторами.

Уравнения второго порядка широко рассматривались в литературе XIX века. Обычно выделяют уравнения со старшими производными ∂ xx {\ displaystyle \ partial _ {xx}}{\ displaystyle \ partial _ {xx}} или ∂ xy {\ displaystyle \ partial _ {xy}}{ \ displaystyle \ partial _ {xy}} . Их общие решения содержат не только константы, но и неопределенные функции переменного числа аргументов; их определение является частью процедуры решения. Для уравнений с ведущей производной ∂ x x {\ displaystyle \ partial _ {xx}}{\ displaystyle \ partial _ {xx}} результаты Лоуи можно обобщить следующим образом.

Теорема 2 Пусть дифференциальный оператор L {\ displaystyle L}L определяется как

L ≡ ∂ xx + A 1 ∂ xy + A 2 ∂ yy + A 3 ∂ Икс + A 4 ∂ Y + A 5 {\ Displaystyle L \ Equiv \ partial _ {xx} + A_ {1} \ partial _ {xy} + A_ {2} \ partial _ {yy} + A_ {3} \ частичный _ {x} + A_ {4} \ partial _ {y} + A_ {5}}{\ displaystyle L \ Equiv \ partial _ {xx} + A_ {1} \ partial _ {xy} + A_ {2} \ p artial _ {yy} + A_ {3} \ partial _ {x} + A_ {4} \ partial _ {y} + A_ {5}} где A i ∈ Q (x, y) {\ displaystyle A_ {i} \ in {\ mathbb {Q}} (x, y)}{\ displaystyle A_ {i} \ in {\ mathbb {Q}} (x, y)} для всех i {\ displaystyle i}i.

Пусть li ≡ ∂ x + ai ∂ y + bi {\ displaystyle l_ {i} \ Equiv \ partial _ {x} + a_ {i} \ partial _ {y} + b_ {i}}{\ displaystyle l_ {i} \ Equiv \ partial _ {x} + a_ {i} \ partial _ {y} + b_ {i}} для i = 1 {\ displaystyle i = 1}i = 1 и я = 2 {\ displaystyle i = 2}i = 2 и l (Φ) ≡ ∂ x + a ∂ y + b (Φ) {\ displaystyle l (\ Phi) \ Equiv \ partial _ {x} + a \ partial _ {y} + b (\ Phi)}{\ displaystyle l (\ Phi) \ Equiv \ partial _ {x} + a \ partial _ {y} + b (\ Phi)} - операторы первого порядка с ai, bi, a ∈ Q ( x, y) {\ displaystyle a_ {i}, b_ {i}, a \ in {\ mathbb {Q}} (x, y)}{\ displaystyle a_ {i}, b_ {i}, a \ in {\ mathbb {Q}} (x, y)} ; Φ {\ displaystyle \ Phi}\ Phi равно неопределенная функция одного аргумента. Тогда L {\ displaystyle L}L имеет разложение Леви согласно одному из следующих типов.

L x x 1: L = l 2 l 1; {\ displaystyle {\ mathcal {L}} _ {xx} ^ {1}: L = l_ {2} l_ {1};}{\ displaystyle {\ mathcal {L}} _ {xx} ^ {1}: L = l_ { 2} l_ {1};} L x x 2: L = L c l m (l 2, l 1); {\ displaystyle {\ mathcal {L}} _ {xx} ^ {2}: L = Lclm (l_ {2}, l_ {1});}{\ displaystyle {\ mathcal {L}} _ {xx} ^ {2} : L = Lclm (l_ {2}, l_ {1});} L xx 3: L = L clm (l (Φ)). {\ displaystyle {\ mathcal {L}} _ {xx} ^ {3}: L = Lclm (l (\ Phi)).}{\ displaystyle {\ mathcal {L}} _ {xx} ^ {3}: L = Lclm (l (\ Phi)).}

Тип разложения оператора L {\ displaystyle L}L - это разложение L xxi {\ displaystyle {\ mathcal {L}} _ {xx} ^ {i}}{\ displaystyle {\ mathcal {L}} _ {xx} ^ {i}} с наивысшим значением i {\ displaystyle i}i. Если L {\ displaystyle L}L не имеет множителей первого порядка в базовом поле, его тип разложения определяется как L xx 0 {\ displaystyle {\ mathcal {L }} _ {xx} ^ {0}}{\ displaystyle {\ mathcal {L}} _ {xx} ^ {0}} . Разложения L xx 0 {\ displaystyle {\ mathcal {L}} _ {xx} ^ {0}}{\ displaystyle {\ mathcal {L}} _ {xx} ^ {0}} , L xx 2 {\ displaystyle {\ mathcal {L}} _ {xx} ^ {2} }{\ displaystyle {\ mathcal {L}} _ {xx} ^ {2}} и L xx 3 {\ displaystyle {\ mathcal {L}} _ {xx} ^ {3}}{\ displaystyle {\ mathcal {L}} _ {xx} ^ {3}} полностью сокращаются.

Чтобы применить этот результат для решения любого данного дифференциального уравнения, включающего оператор L {\ displaystyle L}L , возникает вопрос, можно ли алгоритмически определить его факторы первого порядка. Последующее следствие дает ответ для факторов с коэффициентами либо в базовом поле, либо в универсальном расширении поля.

Следствие 3 В общем случае правые множители первого порядка линейного pde в базовом поле не могут быть определены алгоритмически. Если символьный полином является отделимым, может быть определен любой коэффициент. Если он имеет двойной корень, в общем случае невозможно определить правильные коэффициенты в базовом поле. Всегда можно решить вопрос о существовании факторов в универсальном поле, т.е. об абсолютной несводимости.

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

Утверждение 1 Пусть приводимое уравнение второго порядка

L z ≡ zxx + A 1 zxy + A 2 zyy + A 3 zx + A 4 zy + A 5 z = 0 {\ displaystyle Lz \ Equiv z_ {xx} + A_ {1} z_ {xy} + A_ {2} z_ {yy} + A_ {3} z_ {x} + A_ {4} z_ {y} + A_ {5} z = 0}{\ displaystyle Lz \ Equiv z_ {xx} + A_ {1} z_ {xy} + A_ {2} z_ {yy} + A_ {3} z_ {x} + A_ {4} z_ {y} + A_ {5} z = 0} где A 1,…, A 5 ∈ Q (x, y) {\ displaystyle A_ {1}, \ ldots, A_ {5} \ in {\ mathbb {Q}} (x, y)}{\ displaystyle A_ {1}, \ ldots, A_ {5} \ in {\ mathbb {Q}} (x, y)} .

Определить li ≡ ∂ x + ai ∂ y + bi {\ displaystyle l_ {i} \ Equiv \ partial _ {x} + a_ {i} \ partial _ {y} + b_ {i }}{\ displaystyle l_ {i} \ Equiv \ partial _ {x} + a_ {i} \ partial _ {y} + b_ {i}} , ai, bi ∈ Q (x, y) {\ displaystyle a_ {i}, b_ {i} \ in {\ mathbb {Q}} (x, y)}{\ displaystyle a_ {i}, b_ {i} \ in {\ mathbb {Q}} (x, y)} для я знак равно 1, 2 {\ Displaystyle я = 1,2}i = 1,2 ; φ я (х, y) = const {\ displaystyle \ varphi _ {я} (х, y) = const}{\ displaystyle \ varphi _ {i} (x, y) = const} - рациональный первый интеграл от dydx = ai (x, y) {\ displaystyle {\ frac {dy} {dx}} = a_ {i} (x, y)}{\ displaystyle {\ frac {dy} {dx}} = a_ {i} (x, y)} ; y ¯ ≡ φ i (Икс, Y) {\ Displaystyle {\ bar {y}} \ Equiv \ varphi _ {i} (x, y)}{\ displaystyle {\ bar {y}} \ Equiv \ varphi _ {i} (x, y)} и обратное y = ψ i (x, y ¯) {\ displaystyle y = \ psi _ {i} (x, {\ bar {y}})}{\ displaystyle y = \ psi _ {i} (x, {\ bar {y}})} ; Предполагается, что существуют оба φ i {\ displaystyle \ varphi _ {i}}\ varphi _ {i} и ψ i {\ displaystyle \ psi _ {i}}\psi_i. Кроме того, определим

E i (x, y) ≡ exp ⁡ (- ∫ b i (x, y) | y = ψ i (x, y ¯) d x) | Y ¯ знак равно φ я (Икс, Y) {\ Displaystyle {\ mathcal {E}} _ {я} (х, у) \ эквив \ ехр {\ Big (} - {\ Displaystyle \ int} b_ {я} ( x, y) {\ big |} _ {y = \ psi _ {i} (x, {\ bar {y}})} dx {\ Big)} {\ Big |} _ {{\ bar {y} } = \ varphi _ {i} (x, y)}}{\ displaystyle {\ mathcal {E}} _ {i} (x, y) \ Equiv \ exp { \ Big (} - {\ displaystyle \ int} b_ {i} (x, y) {\ big |} _ {y = \ psi _ {i} (x, {\ bar {y}})} dx {\ Большой)} {\ Big |} _ {{\ bar {y}} = \ varphi _ {i} (x, y)}} для i = 1, 2 {\ displaystyle i = 1,2}i = 1,2 .

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

L xx 1: z 1 (x, y) = E 1 (x, y) F 1 (φ 1) {\ displaystyle {\ mathcal {L}} _ {xx} ^ {1}: z_ {1 } (x, y) = {\ mathcal {E}} _ {1} (x, y) F_ {1} (\ varphi _ {1})}{\ displaystyle {\ mathcal {L}} _ {xx} ^ {1}: z_ {1} (x, y) = {\ mathcal {E}} _ {1} (x, y) F_ {1} (\ varphi _ {1})} , z 2 (x, y) = E 1 ( x, y) ∫ E 2 (x, y) E 1 (x, y) F 2 (φ 2 (x, y)) | y = ψ 1 (x, y ¯) d x | y ¯ = φ 1 (x, y); {\ displaystyle z_ {2} (x, y) = {\ mathcal {E}} _ {1} (x, y) {\ displaystyle \ int} {\ frac {{\ mathcal {E}} _ {2} (x, y)} {{\ mathcal {E}} _ {1} (x, y)}} F_ {2} {\ big (} \ varphi _ {2} (x, y) {\ big)} {\ big |} _ {y = \ psi _ {1} (x, {\ bar {y}})} dx {\ Big |} _ {{\ bar {y}} = \ varphi _ {1} ( x, y)};}{\ displaystyle z_ {2} (x, y) = {\ mathcal {E}} _ {1} (x, y) {\ displaystyle \ int} {\ frac {{\ mathcal {E}} _ {2} (x, y)} {{\ mathcal {E}} _ {1} (x, y)}} F_ {2} {\ big (} \ varphi _ {2} (x, y) {\ big)} {\ big |} _ {y = \ psi _ {1} (x, {\ bar {y}})} dx {\ Big |} _ {{\ bar {y}} = \ varphi _ { 1} (x, y)};}

L xx 2: zi (x, y) = E i (x, y) F i (φ i (x, y)), i = 1, 2; {\ displaystyle {\ mathcal {L}} _ {xx} ^ {2}: z_ {i} (x, y) = {\ mathcal {E}} _ {i} (x, y) F_ {i} { \ big (} \ varphi _ {i} (x, y) {\ big)}, i = 1,2;}{\ displaystyle {\ mathcal {L}} _ {xx} ^ {2}: z_ {i} (x, y) = {\ mathcal {E}} _ {i} (x, y) F_ {i} {\ big ( } \ varphi _ {i} (x, y) {\ big)}, i = 1,2;}

L xx 3: zi (x, y) = E i (x, y) F я (φ (х, y)), я знак равно 1, 2. {\ displaystyle {\ mathcal {L}} _ {xx} ^ {3}: z_ {i} (x, y) = {\ mathcal {E }} _ {i} (x, y) F_ {i} {\ big (} \ varphi (x, y) {\ big)}, i = 1,2.}{\ displaystyle {\ mathcal {L}} _ {xx} ^ {3}: z_ {i} (x, y) = {\ mathcal {E}} _ {i} (x, y) F_ {i} {\ big (} \ varphi (x, y) {\ big)}, i = 1,2.}

F i { \ displaystyle F_ {i}}F_ { i} - неопределенные функции одного аргумента; φ {\ displaystyle \ varphi}\ varphi , φ 1 {\ displaystyle \ varphi _ {1}}\ varphi _ {1} и φ 2 {\ displaystyle \ varphi _ {2}}\ varphi _ {2} рациональны по всем аргументам; ψ 1 {\ displaystyle \ psi _ {1}}\ psi _ {1} предполагается. Обычно φ 1 ≠ φ 2 {\ displaystyle \ varphi _ {1} \ neq \ varphi _ {2}}{\ displaystyle \ varphi _ {1} \ neq \ varphi _ {2}} , они определяются коэффициентами A 1 {\ displaystyle A_ {1}}A_ {1} , A 2 {\ displaystyle A_ {2}}A_ {2} и A 3 {\ displaystyle A_ {3}}A_ {3} данного уравнения.

Типичным примером линейного pde, в котором применяется факторизация, является уравнение, которое обсуждалось Forsyth, vol. VI, page 16,

Пример 5 (Forsyth 1906)} Рассмотрим дифференциальное уравнение zxx - zyy + 4 x + yzx = 0 {\ displaystyle z_ {xx} -z_ {yy} + {\ гидроразрыв {4} {x + y}} z_ {x} = 0}{\ displaystyle z_ {xx} -z_ {yy} + {\ frac {4} {x + y}} z_ {x} = 0} . После факторизации представление

L z ≡ l 2 l 1 z = (∂ x + ∂ y + 2 x + y) (∂ x - ∂ y + 2 x + y) z = 0 {\ displaystyle Lz \ Equiv l_ {2} l_ {1} z = {\ Big (} \ partial _ {x} + \ partial _ {y} + {\ frac {2} {x + y}} {\ Big)} {\ Big (} \ partial _ {x} - \ partial _ {y} + {\ frac {2} {x + y}} {\ Big)} z = 0}{\ displaystyle Lz \ Equiv l_ {2} l_ {1} z = {\ Big (} \ partial _ {x} + \ partial _ {y} + {\ frac {2} {x + y}} {\ Большой)} {\ Big (} \ partial _ {x} - \ partial _ {y} + {\ frac {2} {x + y}} {\ Big)} z = 0} получается. Отсюда следует

φ 1 (x, y) = x + y, ψ 1 (x, y) = y ¯ - x, E 1 (x, y) = exp ⁡ (2 yx + y) {\ displaystyle \ varphi _ {1} (x, y) = x + y, \ psi _ {1} (x, y) = {\ bar {y}} - x, {\ mathcal {E}} _ {1} (x, y) = \ exp {{\ Big (} {\ frac {2y} {x + y}} {\ Big)}}}{\ displaystyle \ varphi _ {1} (x, y) = x + y, \ psi _ {1} (x, y) = {\ bar {y}} - x, {\ mathcal { E}} _ {1} (x, y) = \ exp {{\ Big (} {\ frac {2y} {x + y}} {\ Big)}}} ,

φ 2 (x, y) = x - y, ψ 2 (x, y) = x - y ¯, E 2 (x, y) = - 1 x + y. {\ displaystyle \ varphi _ {2} (x, y) = xy, \ psi _ {2} (x, y) = x - {\ bar {y}}, {\ mathcal {E}} _ {2} (x, y) = - {\ frac {1} {x + y}}.}{\ displaystyle \ varphi _ {2} (x, y) = xy, \ psi _ {2} (x, y) = x - {\ bar {y}}, { \ mathcal {E}} _ {2} (x, y) = - {\ frac {1} {x + y}}.}

Следовательно, дифференциальная фундаментальная система имеет вид

z 1 (x, y) = exp ⁡ (2 yx + y) F (x + y), {\ displaystyle z_ {1} (x, y) = \ exp {{\ Big (} {\ frac {2y} {x + y}} {\ Big)}} F (x + y),}{\ displaystyle z_ { 1} (x, y) = \ exp {{\ Big (} {\ frac {2y} {x + y}} {\ Big)}} F (x + y),} z 2 (x, y) = 1 x + y exp ⁡ (2 yx + y) ∫ exp ⁡ (2 x - y ¯ y ¯) G (2 x - y ¯) dx | y ¯ = x + y. {\ displaystyle z_ {2} (x, y) = {\ frac {1} {x + y}} \ exp {{\ Big (} {\ frac {2y} {x + y}} {\ Big)} }{\displaystyle \int }\exp {{\Big (}{\frac {2x-{\bar {y}}}{\bar {y}}}{\Big)}}G(2x-{\bar {y}})dx{\Big |}_{{\bar {y}}=x+y}.}{\ displaystyle z_ {2} (x, y) = {\ frac {1} { x + y}} \ exp {{\ Big (} {\ frac {2y} {x + y}} {\ Big)}} {\ displaystyle \ int} \ exp {{\ Big (} {\ frac {2x) - {\ bar {y}}} {\ bar {y}}} {\ Big)}} G (2x - {\ bar {y}}) dx {\ Big |} _ {{\ bar {y}} = x + y}.}

F {\displaystyle F}F and G {\displaystyle G}G are undetermined functions.

If the only second-order derivative of an operator is ∂ x y {\displaystyle \partial _{xy}}{ \ displaystyle \ partial _ {xy}} , its possible decompositions involving only principal divisors may be described as follows.

Theorem 3Let the differential operator L {\displaystyle L}L be defined by

L ≡ ∂ x y + A 1 ∂ x + A 2 ∂ y + A 3 {\displaystyle L\equiv \partial _{xy}+A_{1}\partial _{x}+A_{2}\partial _{y}+A_{3}}{\ Displaystyle L \ Equiv \ partial _ {xy} + A_ {1} \ partial _ {x} + A_ {2} \ partial _ {y} + A_ {3}} where A i ∈ Q ( x, y) {\displaystyle A_{i}\in {\mathbb {Q} }(x,y)}{\ displaystyle A_ {i} \ in {\ mathbb {Q}} (x, y)} for all i {\displaystyle i}i.

Let l ≡ ∂ x + A 2 {\displaystyle l\equiv \partial _{x}+A_{2}}{\ displaystyle l \ Equiv \ partial _ {x} + A_ {2 }} and k ≡ ∂ y + A 1 {\displaystyle k\equiv \partial _{y}+A_{1}}{\ displaystyle k \ Equiv \ partial _ {y} + A_ {1}} are first-order operators. L {\displaystyle L}L has Loewy decompositions involving first-order principal divisors of the following form.

L x y 1 : L = k l ; {\displaystyle {\mathcal {L}}_{xy}^{1}:L=kl;}{\ displaystyle {\ mathcal {L}} _ {xy} ^ {1}: L = kl; } L x y 2 : L = l k ; {\displaystyle {\mathcal {L}}_{xy}^{2}:L=lk;}{\ displaystyle {\ mathcal {L}} _ {xy} ^ {2}: L = lk;} L x y 3 : L = L c l m ( k, l). {\displaystyle {\mathcal {L}}_{xy}^{3}:L=Lclm(k,l).}{\ displaystyle {\ mathcal {L}} _ {xy} ^ {3}: L = Lclm (k, l).}

The тип разложения оператора L {\ displaystyle L}L - это разложение L xyi {\ displaystyle {\ mathcal {L}} _ {xy} ^ {i}}{\ displaystyle {\ mathcal {L}} _ {xy} ^ {i}} с наибольшим значением i {\ displaystyle i}i. Разложение типа L xy 3 {\ displaystyle {\ mathcal {L}} _ {xy} ^ {3}}{\ displaystyle {\ mathcal {L}} _ { ху} ^ {3}} полностью сводимо

Кроме того, есть еще пять возможные типы разложения с участием неглавных делителей Лапласа, как показано ниже.

Теорема 4 Пусть дифференциальный оператор L {\ displaystyle L}L определяется как

L ≡ ∂ xy + A 1 ∂ x + A 2 ∂ y + A 3 {\ displaystyle L \ Equiv \ partial _ {xy} + A_ {1} \ partial _ {x} + A_ {2} \ partial _ {y} + A_ {3}}{\ Displaystyle L \ Equiv \ partial _ {xy} + A_ {1} \ partial _ {x} + A_ {2} \ partial _ {y} + A_ {3}} где A i ∈ Q (x, y) {\ displaystyle A_ {i} \ in {\ mathbb {Q}} (x, y)}{\ displaystyle A_ {i} \ in {\ mathbb {Q}} (x, y)} для всех i {\ displaystyle i}i.

L xm (L) {\ displaystyle \ mathbb {L} _ {x ^ {m}} (L)}{\ displaystyle \ mathbb {L} _ {x ^ {m}} (L)} и L yn (L) {\ displaystyle \ mathbb {L} _ {y ^ {n}} (L)}{\ displaystyle \ mathbb {L} _ {y ^ {n}} (L)} , а также lm {\ displaystyle {\ mathfrak {l}} _ {m}}{\ displaystyl е {\ mathfrak {l}} _ {m}} и kn {\ displaystyle {\ mathfrak {k}} _ {n}}{\ displaystyle {\ mathfrak {k}} _ {n}} определены выше; кроме того l ≡ ∂ Икс + a {\ displaystyle l \ Equiv \ partial _ {x} + a}{\ displaystyle l \ Equiv \ partial _ {x} + a} , к ≡ ∂ y + b {\ displaystyle k \ Equiv \ partial _ {y} + b}{\ displaystyle k \ Equiv \ partial _ {y } + b} , a, b ∈ Q (x, y) {\ displaystyle a, b \ in {\ mathbb {Q}} (x, y)}{\ displaystyle a, b \ in {\ mathbb {Q}} (x, y)} . L {\ displaystyle L}L имеет Loewy разложения с делителями Лапласа по одному из следующих типов; m {\ displaystyle m}м и n {\ displaystyle n}nподчиняться m, n ≥ 2 {\ displaystyle m, n \ geq 2}{\ displaystyle m, n \ geq 2} .

L xy 4: L = L clm (L xm (L), L yn (L)); {\ displaystyle {\ mathcal {L}} _ {xy} ^ {4}: L = Lclm {\ big (} \ mathbb {L} _ {x ^ {m}} (L), \ mathbb {L} _ {y ^ {n}} (L) {\ big)};}{\ Displaystyle {\ mathcal {L}} _ {ху} ^ {4}: L = Lclm {\ big (} \ mathbb {L} _ {x ^ {m}} (L), \ mathbb {L} _ {y ^ {n}} (L) {\ big) };}

L xy 5: L = E xquo (L, L xm (L)) L xm (L) = (1 0 0 ∂ y + A 1) (L lm); {\ displaystyle {\ mathcal {L}} _ {xy} ^ {5}: L = Exquo {\ big (} L, \ mathbb {L} _ {x ^ {m}} (L) {\ big)} \ mathbb {L} _ {x ^ {m}} (L) = \ left ({\ begin {array} {cc} 1 0 \\ 0 \ partial _ {y} + A_ {1} \ end {array}} \ right) \ left ({\ begin {array} {c} L \\ {\ mathfrak {l}} _ {m} \ end {array}} \ right);}{\ displaystyle {\ mathcal {L}} _ {xy } ^ {5}: L = Exquo {\ big (} L, \ mathbb {L} _ {x ^ {m}} (L) {\ big)} \ mathbb {L} _ {x ^ {m}} (L) = \ left ({\ begin {array} {cc} 1 0 \\ 0 \ partial _ {y} + A_ {1} \ end {array}} \ right) \ left ({\ begin {array} { c} L \\ {\ mathfrak {l}} _ {m} \ end {array}} \ right);}

L xy 6: L = E xquo (L, L yn (L)) L yn (L) = (1 0 0 ∂ x + A 2) (L kn); {\ displaystyle {\ mathcal {L}} _ {xy} ^ {6}: L = Exquo {\ big (} L, \ mathbb {L} _ {y ^ {n}} (L) {\ big)} \ mathbb {L} _ {y ^ {n}} (L) = \ left ({\ begin {array} {cc} 1 0 \\ 0 \ partial _ {x} + A_ {2} \ end {array}} \ right) \ left ({\ begin {array} {c} L \\ {\ mathfrak {k}} _ {n} \ end {array}} \ right);}{\ displaystyle {\ mathcal { L}} _ {xy} ^ {6}: L = Exquo {\ big (} L, \ mathbb {L} _ {y ^ {n}} (L) {\ big)} \ mathbb {L} _ { y ^ {n}} (L) = \ left ({\ begin {array} {cc} 1 0 \\ 0 \ partial _ {x} + A_ {2} \ end {array}} \ right) \ left ({ \ begin {array} {c} L \\ {\ mathfrak {k}} _ {n} \ end {array}} \ right);}

L xy 7: L = L clm (k, L xm (L)); {\ displaystyle {\ mathcal {L}} _ {xy} ^ {7}: L = Lclm {\ big (} k, \ mathbb {L} _ {x ^ {m}} (L) {\ big)} ;}{\ displaystyle {\ mathcal {L}} _ {xy} ^ {7}: L = Lclm {\ big (} k, \ mathbb {L} _ {x ^ {m}} (L) {\ big)};} L xy 8: L = L clm (l, L yn (L)). {\ displaystyle {\ mathcal {L}} _ {xy} ^ {8}: L = Lclm {\ big (} l, \ mathbb {L} _ {y ^ {n}} (L) {\ big)}.}{\ displaystyle {\ mathcal {L}} _ {xy} ^ {8}: L = Lclm {\ big (} l, \ mathbb {L} _ {y ^ {n}} (L) {\ big)}.}

Если L {\ displaystyle L}L не имеет правого множителя первого порядка и может быть показано, что делитель Лапласа не существует, его тип разложения определяется как L ху 0 {\ displaystyle {\ mathcal {L}} _ {xy} ^ {0}}{\ displaystyle {\ mathcal {L}} _ {xy} ^ {0}} . Разложения L xy 0 {\ displaystyle {\ mathcal {L}} _ {xy} ^ {0}}{\ displaystyle {\ mathcal {L}} _ {xy} ^ {0}} , L xy 4 {\ displaystyle {\ mathcal {L}} _ {xy} ^ {4 }}{\ displaystyle {\ mathcal {L}} _ {xy} ^ { 4}} , L xy 7 {\ displaystyle {\ mathcal {L}} _ {xy} ^ {7}}{\ displaystyle {\ mathcal {L}} _ {xy} ^ {7}} и L xy 8 {\ displaystyle {\ mathcal {L}} _ {xy} ^ {8}}{\ displaystyle {\ mathcal {L}} _ {xy} ^ {8}} полностью сокращаются.

Уравнение, которое не допускает разложения с участием главных делителей, но полностью сводимо относительно. неглавные делители Лапласа типа L x y 4 {\ displaystyle {\ mathcal {L}} _ {xy} ^ {4}}{\ displaystyle {\ mathcal {L}} _ {xy} ^ { 4}} были рассмотрены Форсайтом.

Пример 6 (Форсайт 1906) Определить

L ≡ ∂ xy + 2 x - y ∂ x - 2 x - y ∂ y - 4 (x - y) 2 {\ displaystyle L \ Equiv \ partial _ {xy} + {\ frac {2} {xy}} \ partial _ {x} - {\ frac {2} {xy}} \ partial _ {y} - {\ frac {4} {(xy) ^ {2}}}}{\ displaystyle L \ Equiv \ partial _ {xy} + {\ frac {2} {xy}} \ partial _ {x} - {\ frac {2} {xy}} \ partial _ {y} - {\ frac {4} {(xy) ^ {2}}}}

создание главного идеала ⟨L⟩ {\ displaystyle \ langle L \ rangle}{\ displaystyle \ langle L \ rangle} . Фактора первого порядка не существует. Однако есть делители Лапласа

L x 2 (L) ≡ ⟨⟨∂ xx - 2 x - y ∂ x + 2 (x - y) 2, L⟩⟩ {\ displaystyle {\ mathbb {L}} _ {x ^ {2}} (L) \ Equiv {{\ Big \ langle} {\ Big \ langle}} \ partial _ {xx} - {\ frac {2} {xy}} \ partial _ {x} + {\ frac {2} {(xy) ^ {2}}}, L {{\ Big \ rangle} {\ Big \ rangle}}}{\ displaystyle {\ mathbb {L}} _ {x ^ {2}} (L) \ Equiv {{\ Big \ langle} {\ Big \ langle}} \ partial _ {xx} - {\ frac {2} {xy}} \ partial _ {x} + {\ frac {2} {(xy) ^ { 2}}}, L {{\ Big \ rangle} {\ Big \ rangle}}} и L y 2 (L) ≡ ⟨⟨L, ∂ yy + 2 x - y ∂ y + 2 (x - y) 2⟩⟩. {\ displaystyle {\ mathbb {L}} _ {y ^ {2}} (L) \ Equiv {{\ Big \ langle} {\ Big \ langle}} L, \ partial _ {yy} + {\ frac { 2} {xy}} \ partial _ {y} + {\ frac {2} {(xy) ^ {2}}} {{\ Big \ rangle} {\ Big \ rangle}}.}{\ displaystyle {\ mathbb {L}} _ {y ^ {2}} (L) \ Equiv {{\ Big \ langle} {\ Big \ langle}} L, \ partial _ {yy} + {\ frac {2} {xy}} \ partial _ {y} + {\ frac {2} {(xy) ^ {2} }} {{\ Big \ rangle} {\ Big \ rangle}}.} 

Идеальный генерируется L {\ displaystyle L}L имеет представление ⟨L⟩ = L clm (L x 2 (L), L y 2 (L)) {\ displaystyle \ langle L \ rangle = Lclm {\ big (} {\ mathbb {L}} _ {x ^ {2}} (L), {\ mathbb {L}} _ {y ^ {2}} (L) {\ big) }}{\ displaystyle \ langle L \ rangle = Lclm { \ big (} {\ mathbb {L}} _ {x ^ {2}} (L), {\ mathbb {L}} _ {y ^ {2}} (L) {\ big)}} , т. Е. Полностью сводимо; его тип разложения - L x y 4 {\ displaystyle {\ mathcal {L}} _ {xy} ^ {4}}{\ displaystyle {\ mathcal {L}} _ {xy} ^ { 4}} . Следовательно, уравнение L z = 0 {\ displaystyle Lz = 0}{\ displaystyle Lz = 0} имеет дифференциальную фундаментальную систему

z 1 (x, y) = 2 (x - y) F (y) + (Икс - Y) 2 F '(Y) {\ Displaystyle Z_ {1} (x, y) = 2 (ху) F (y) + (ху) ^ {2} F' (y)}{\displaystyle z_{1}(x,y)=2(x-y)F(y)+(x-y)^{2}F'(y)}и z 2 (x, y) = 2 (y - x) G (x) + (y - x) 2 G ′ (x) {\ displaystyle z_ {2} (x, y) = 2 (yx) G (x) + (yx) ^ {2} G '(x)}{\displaystyle z_{2}(x,y)=2(y-x)G(x)+(y-x)^{2}G'(x)}.
Разложение линейных pde порядка выше 2

Оказывается, у операторов более высокого порядка сложнее разложений, и есть другие альтернативы, многие из которых в терминах неглавных делителей. Решения соответствующих уравнений усложняются. Для уравнений третьего порядка на плоскости достаточно полный ответ можно найти в. Типичный пример уравнения третьего порядка, который также представляет исторический интерес, принадлежит Блюмбергу.

Пример 7 (Blumberg 1912) В В своей диссертации Блумберг рассматривал оператор третьего порядка

L ≡ ∂ xxx + x ∂ xxy + 2 ∂ xx + 2 (x + 1) ∂ xy + ∂ x + (x + 2) ∂ y {\ displaystyle L \ Equiv \ частичный _ {xxx} + x \ partial _ {xxy} +2 \ partial _ {xx} +2 (x + 1) \ partial _ {xy} + \ partial _ {x} + (x + 2) \ partial _ {y}}{\ displaystyle L \ Equiv \ partial _ {xxx} + x \ partial _ {xxy} +2 \ partial _ {xx} +2 (x +1) \ partial _ {xy} + \ partial _ {x} + (x + 2) \ partial _ {y}} .

Он позволяет использовать два фактора первого порядка l 1 ≡ ∂ x + 1 {\ displaystyle l_ {1} \ Equiv \ partial _ {x} +1}{\ displaystyle l_ {1} \ Equiv \ partial _ {x} +1} и l 2 ≡ ∂ Икс + Икс ∂ Y {\ Displaystyle L_ {2} \ Equiv \ partial _ {x} + x \ partial _ {y}}{\ displaystyle l_ {2} \ Equiv \ partial _ {x} + x \ partial _ {y}} . Их пересечение не принципиально; определение

L 1 ≡ ∂ xxx - x 2 ∂ xyy + 3 ∂ xx + (2 x + 3) ∂ xy - x 2 ∂ yy + 2 ∂ x + (2 x + 3) ∂ y {\ displaystyle L_ { 1} \ Equiv \ partial _ {xxx} -x ^ {2} \ partial _ {xyy} +3 \ partial _ {xx} + (2x + 3) \ partial _ {xy} -x ^ {2} \ partial _ {yy} +2 \ partial _ {x} + (2x + 3) \ partial _ {y}}{\ displaystyle L_ {1} \ Equiv \ partial _ {xxx} -x ^ {2} \ partial _ {xyy} +3 \ partial _ {xx} + (2x + 3) \ partial _ {xy} -x ^ {2} \ partial _ {yy} +2 \ partial _ {x} + (2x + 3) \ partial _ {y}}

L 2 ≡ ∂ xxy + x ∂ xyy - 1 x ∂ xx - 1 x ∂ xy + x ∂ yy - 1 x ∂ x - (1 + 1 x) ∂ y⟩⟩. {\ Displaystyle L_ {2} \ Equiv \ partial _ {xxy} + x \ partial _ {xyy} - {\ frac {1} {x}} \ partial _ {xx} - {\ frac {1} {x} } \ partial _ {xy} + x \ partial _ {y} y - {\ frac {1} {x}} \ partial _ {x} - {\ big (} 1 + {\ frac {1} {x} } {\ big)} \ partial _ {y} {{\ big \ rangle} {\ big \ rangle}}.}{\ displaystyle L_ {2} \ Equiv \ partial _ {xxy} + x \ partial _ {xyy} - {\ frac {1} {x}} \ partial _ {xx} - {\ frac {1} {x}} \ partial _ {xy} + x \ partial _ {y} y - {\ frac {1} {x}} \ partial _ {x} - {\ big (} 1 + {\ frac {1} {x}} {\ big)} \ partial _ {y} {{\ big \ rangle} {\ big \ rangle}}.}

это может быть записано как L clm (l 2, l 1) = ⟨ ⟨L 1, L 2⟩⟩ {\ displaystyle Lclm (l_ {2}, l_ {1}) = {\ langle \ langle} L_ {1}, L_ {2} {\ rangle \ rangle}}{\ displaystyle Lclm (l_ {2}, l_ { 1}) = {\ langle \ langle} L_ {1}, L_ {2} {\ rangle \ rangle}} . Следовательно, разложение Лоуи оператора Блюмбергса имеет вид

L = (1 x 0 ∂ x + 1 + 1 x) (L 1 L 2). {\ displaystyle L = \ left ({\ begin {array} {cc} 1 x \\ 0 \ partial _ {x} +1 + {\ frac {1} {x}} \ end {array}} \ right) \ left ({\ begin {array} {c} L_ {1} \\ L_ {2} \ end {array}} \ right).}{\ displaystyle L = \ left ({\ begin {array} {cc} 1 x \\ 0 \ partial _ {x} +1 + {\ frac {1} {x}} \ end {array}} \ right) \ left ({\ begin {array} {c} L_ {1} \\ L_ {2} \ end {array}} \ right).}

Это дает следующую дифференциальную фундаментальную систему для дифференциального уравнения L z знак равно 0 {\ displaystyle Lz = 0}{\ displaystyle Lz = 0} .

z 1 (x, y) = F (y - 1 2 x 2) {\ displaystyle z_ {1} (x, y) = F (y - {\ frac {1} {2}} x ^ {2})}{\ displaystyle z_ {1} (x, y) = F ( y - {\ frac {1} {2}} x ^ {2})} , z 2 (x, y) = G (y) e - x {\ displaystyle z_ {2} (x, y) = G (y) e ^ {- x}}{\ displaystyle z_ {2} (x, y) = G (y) e ^ {- x}} , z 3 (x, y) = ∫ xe - x H (y ¯ + 1 2 x 2) dx | Y ¯ = Y - 1 2 Икс 2 {\ Displaystyle Z_ {3} (x, y) = {\ displaystyle \ int} xe ^ {- x} H {\ big (} {\ bar {y}} + {\ frac {1} {2}} x ^ {2} {\ big)} dx {\ Big |} _ {{\ bar {y}} = y - {\ frac {1} {2}} x ^ {2 }}}{\ displaystyle z_ {3} (x, y) = {\ displaystyle \ int} xe ^ {- x} H { \ big (} {\ bar {y}} + {\ frac {1} {2}} x ^ {2} {\ big)} dx {\ Big |} _ {{\ bar {y}} = y- {\ frac {1} {2}} x ^ {2}}}

F, G {\ displaystyle F, G}F, G и H {\ displaystyle H}H - неопределенные функции.

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

Ссылки
Последняя правка сделана 2021-05-28 05:26:43
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте