Метод аннигилятора

Метод решения неоднородных обыкновенных дифференциальных уравнений

В математике метод аннигиляторов — это процедура, используемая для нахождения частного решения для определенных типов неоднородных обыкновенных дифференциальных уравнений (ОДУ). [1] Он похож на метод неопределенных коэффициентов , но вместо угадывания частного решения в методе неопределенных коэффициентов , в этой технике частное решение определяется систематически. Фраза неопределенные коэффициенты также может использоваться для обозначения шага в методе аннигиляторов, на котором вычисляются коэффициенты.

Метод аннулятора используется следующим образом. Дано ОДУ , найдите другой дифференциальный оператор такой, что . Этот оператор называется аннулятором , отсюда и название метода. Применение к обеим сторонам ОДУ дает однородное ОДУ , для которого мы находим базис решения, как и прежде. Затем исходное неоднородное ОДУ используется для построения системы уравнений, ограничивающих коэффициенты линейной комбинации для удовлетворения ОДУ. П ( Д ) у = ф ( х ) {\displaystyle P(D)y=f(x)} А ( Д ) {\displaystyle A(D)} А ( Д ) ф ( х ) = 0 {\displaystyle A(D)f(x)=0} А ( Д ) {\displaystyle A(D)} ( А ( Д ) П ( Д ) ) у = 0 {\displaystyle {\big (}A(D)P(D){\big )}y=0} { у 1 , , у н } {\displaystyle \{y_{1},\ldots ,y_{n}\}}

Этот метод не столь универсален, как вариация параметров, в том смысле, что аннигилятор не всегда существует.

Таблица аннигиляторов

ж ( х )А ( Д )
а н х н + а н 1 х н 1 + + а 1 х + а 0 {\displaystyle a_{n}x^{n}+a_{n-1}x^{n-1}+\cdots +a_{1}x+a_{0}\!} Д н + 1 {\displaystyle D^{n+1}\!}
е а х {\displaystyle e^{топор}\!} Д а {\displaystyle Да\!}
х к . е а х {\displaystyle x^{k}.e^{ax}\!} ( Д а ) к + 1 {\displaystyle (Да)^{k+1}\!}
потому что ( б х ) о г грех ( б х ) {\displaystyle \cos(bx)\;\;\mathrm {или} \;\;\sin(bx)\!} Д 2 + б 2 {\displaystyle D^{2}+b^{2}\!}
х к потому что ( б х ) о г х к грех ( б х ) {\displaystyle x^{k}\cos(bx)\;\;\mathrm {или} \;\;x^{k}\sin(bx)\!} ( Д 2 + б 2 ) к + 1 {\displaystyle (D^{2}+b^{2})^{k+1}\!}
е а х потому что ( б х ) о г е а х грех ( б х ) {\displaystyle е^{ax}\cos(bx)\;\;\mathrm {or} \;\;e^{ax}\sin(bx)\!} ( Д а ) 2 + б 2 = Д 2 2 а Д + а 2 + б 2 {\displaystyle (Da)^{2}+b^{2}=D^{2}-2aD+a^{2}+b^{2}\!}
х к е а х потому что ( б х ) о г х к е а х грех ( б х ) {\displaystyle x^{k}e^{ax}\cos(bx)\;\;\mathrm {or} \;\;x^{k}e^{ax}\sin(bx)\!} [ ( Д а ) 2 + б 2 ] к + 1 = [ Д 2 2 а Д + а 2 + б 2 ] к + 1 {\displaystyle \left[(Da)^{2}+b^{2}\right]^{k+1}=\left[D^{2}-2aD+a^{2}+b^{2}\right]^{k+1}\!}
а н х н + + а 1 х + а 0 + б 1 е ± с 1 х + + б к е с к х {\displaystyle a_{n}x^{n}+\cdots +a_{1}x+a_{0}+b_{1}e^{\pm c_{1}x}+\cdots +b_{k}e^{\mp c_{k}x}\!} Д н + 1 ( Д с 1 ) . . ( Д ± с к ) {\displaystyle D^{n+1}(D\mp c_{1}).\cdots .(D\pm c_{k})\!}

Где — в натуральных числах , а — в действительных числах . k {\displaystyle k} k , b , a , c 1 , , c k {\displaystyle k,b,a,c_{1},\cdots ,c_{k}}

Если состоит из суммы выражений, приведенных в таблице, то аннулятор является произведением соответствующих аннуляторов. f ( x ) {\displaystyle f(x)}

Пример

Дано , . Простейший аннигилятор — это . Нули — это , поэтому базис решения — это y 4 y + 5 y = sin ( k x ) {\displaystyle y''-4y'+5y=\sin(kx)} P ( D ) = D 2 4 D + 5 {\displaystyle P(D)=D^{2}-4D+5} sin ( k x ) {\displaystyle \sin(kx)} A ( D ) = D 2 + k 2 {\displaystyle A(D)=D^{2}+k^{2}} A ( z ) P ( z ) {\displaystyle A(z)P(z)} { 2 + i , 2 i , i k , i k } {\displaystyle \{2+i,2-i,ik,-ik\}} A ( D ) P ( D ) {\displaystyle A(D)P(D)} { y 1 , y 2 , y 3 , y 4 } = { e ( 2 + i ) x , e ( 2 i ) x , e i k x , e i k x } . {\displaystyle \{y_{1},y_{2},y_{3},y_{4}\}=\{e^{(2+i)x},e^{(2-i)x},e^{ikx},e^{-ikx}\}.}

Установку мы находим y = c 1 y 1 + c 2 y 2 + c 3 y 3 + c 4 y 4 {\displaystyle y=c_{1}y_{1}+c_{2}y_{2}+c_{3}y_{3}+c_{4}y_{4}}

sin ( k x ) = P ( D ) y = P ( D ) ( c 1 y 1 + c 2 y 2 + c 3 y 3 + c 4 y 4 ) = c 1 P ( D ) y 1 + c 2 P ( D ) y 2 + c 3 P ( D ) y 3 + c 4 P ( D ) y 4 = 0 + 0 + c 3 ( k 2 4 i k + 5 ) y 3 + c 4 ( k 2 + 4 i k + 5 ) y 4 = c 3 ( k 2 4 i k + 5 ) ( cos ( k x ) + i sin ( k x ) ) + c 4 ( k 2 + 4 i k + 5 ) ( cos ( k x ) i sin ( k x ) ) {\displaystyle {\begin{aligned}\sin(kx)&=P(D)y\\[8pt]&=P(D)(c_{1}y_{1}+c_{2}y_{2}+c_{3}y_{3}+c_{4}y_{4})\\[8pt]&=c_{1}P(D)y_{1}+c_{2}P(D)y_{2}+c_{3}P(D)y_{3}+c_{4}P(D)y_{4}\\[8pt]&=0+0+c_{3}(-k^{2}-4ik+5)y_{3}+c_{4}(-k^{2}+4ik+5)y_{4}\\[8pt]&=c_{3}(-k^{2}-4ik+5)(\cos(kx)+i\sin(kx))+c_{4}(-k^{2}+4ik+5)(\cos(kx)-i\sin(kx))\end{aligned}}}

давая системе

i = ( k 2 + 4 i k 5 ) c 3 + ( k 2 + 4 i k + 5 ) c 4 {\displaystyle i=(k^{2}+4ik-5)c_{3}+(-k^{2}+4ik+5)c_{4}}
0 = ( k 2 + 4 i k 5 ) c 3 + ( k 2 4 i k 5 ) c 4 {\displaystyle 0=(k^{2}+4ik-5)c_{3}+(k^{2}-4ik-5)c_{4}}

который имеет решения

c 3 = i 2 ( k 2 + 4 i k 5 ) {\displaystyle c_{3}={\frac {i}{2(k^{2}+4ik-5)}}} , c 4 = i 2 ( k 2 + 4 i k + 5 ) {\displaystyle c_{4}={\frac {i}{2(-k^{2}+4ik+5)}}}

давая набор решений

y = c 1 y 1 + c 2 y 2 + i 2 ( k 2 + 4 i k 5 ) y 3 + i 2 ( k 2 + 4 i k + 5 ) y 4 = c 1 y 1 + c 2 y 2 + 4 k cos ( k x ) ( k 2 5 ) sin ( k x ) ( k 2 + 4 i k 5 ) ( k 2 4 i k 5 ) = c 1 y 1 + c 2 y 2 + 4 k cos ( k x ) + ( 5 k 2 ) sin ( k x ) k 4 + 6 k 2 + 25 . {\displaystyle {\begin{aligned}y&=c_{1}y_{1}+c_{2}y_{2}+{\frac {i}{2(k^{2}+4ik-5)}}y_{3}+{\frac {i}{2(-k^{2}+4ik+5)}}y_{4}\\[8pt]&=c_{1}y_{1}+c_{2}y_{2}+{\frac {4k\cos(kx)-(k^{2}-5)\sin(kx)}{(k^{2}+4ik-5)(k^{2}-4ik-5)}}\\[8pt]&=c_{1}y_{1}+c_{2}y_{2}+{\frac {4k\cos(kx)+(5-k^{2})\sin(kx)}{k^{4}+6k^{2}+25}}.\end{aligned}}}

Это решение можно разбить на однородную и неоднородную части. В частности, является частным интегралом для неоднородного дифференциального уравнения, а является дополнительным решением соответствующего однородного уравнения. Значения и определяются обычно через набор начальных условий. Поскольку это уравнение второго порядка, для определения этих значений необходимы два таких условия. y p = 4 k cos ( k x ) + ( 5 k 2 ) sin ( k x ) k 4 + 6 k 2 + 25 {\displaystyle y_{p}={\frac {4k\cos(kx)+(5-k^{2})\sin(kx)}{k^{4}+6k^{2}+25}}} y c = c 1 y 1 + c 2 y 2 {\displaystyle y_{c}=c_{1}y_{1}+c_{2}y_{2}} c 1 {\displaystyle c_{1}} c 2 {\displaystyle c_{2}}

Фундаментальные решения и можно далее переписать с использованием формулы Эйлера : y 1 = e ( 2 + i ) x {\displaystyle y_{1}=e^{(2+i)x}} y 2 = e ( 2 i ) x {\displaystyle y_{2}=e^{(2-i)x}}

e ( 2 + i ) x = e 2 x e i x = e 2 x ( cos x + i sin x ) {\displaystyle e^{(2+i)x}=e^{2x}e^{ix}=e^{2x}(\cos x+i\sin x)}
e ( 2 i ) x = e 2 x e i x = e 2 x ( cos x i sin x ) {\displaystyle e^{(2-i)x}=e^{2x}e^{-ix}=e^{2x}(\cos x-i\sin x)}

Тогда , и соответствующее переназначение констант дает более простую и понятную форму дополнительного решения, . c 1 y 1 + c 2 y 2 = c 1 e 2 x ( cos x + i sin x ) + c 2 e 2 x ( cos x i sin x ) = ( c 1 + c 2 ) e 2 x cos x + i ( c 1 c 2 ) e 2 x sin x {\displaystyle c_{1}y_{1}+c_{2}y_{2}=c_{1}e^{2x}(\cos x+i\sin x)+c_{2}e^{2x}(\cos x-i\sin x)=(c_{1}+c_{2})e^{2x}\cos x+i(c_{1}-c_{2})e^{2x}\sin x} y c = e 2 x ( c 1 cos x + c 2 sin x ) {\displaystyle y_{c}=e^{2x}(c_{1}\cos x+c_{2}\sin x)}

Ссылки

  1. ^ Golberg, MA; Chen, CS; Rashed, YF (1 марта 1999 г.). «Метод аннигилятора для вычисления частных решений уравнений с частными производными». Инженерный анализ с граничными элементами . 23 (3): 275–279. doi :10.1016/S0955-7997(98)00081-2. ISSN  0955-7997.
Retrieved from "https://en.wikipedia.org/w/index.php?title=Annihilator_method&oldid=1256555162"