Задача на дифференциальные уравнения: решите задачу о поведении решений при малом параметре epsilon в уравнении epsilon y'' + y' + y = 0, сравните методы регулярной и сингулярной возмущающей теории и объясните, когда требуется слой пограничного типа.

1 Дек в 02:32
2 +2
0
Ответы
1
Уравнение и характерный разбор
Уравнение: ϵy′′+y′+y=0, 0<ϵ≪1\epsilon y''+y'+y=0,\; 0<\epsilon\ll1ϵy′′+y+y=0,0<ϵ1.
Характеристическое уравнение: ϵr2+r+1=0\epsilon r^2+r+1=0ϵr2+r+1=0. Корни
r=−1±1−4ϵ2ϵ. r=\frac{-1\pm\sqrt{1-4\epsilon}}{2\epsilon}.
r=2ϵ1±14ϵ .
При ϵ→0\epsilon\to0ϵ0 имеем асимптотически
r1∼−1+O(ϵ),r2∼−1ϵ+1+O(ϵ), r_1\sim-1+O(\epsilon),\qquad r_2\sim-\frac{1}{\epsilon}+1+O(\epsilon),
r1 1+O(ϵ),r2 ϵ1 +1+O(ϵ),
т.е. общее решение комбинация e−xe^{-x}ex (медленно меняющаяся) и e−x/ϵe^{-x/\epsilon}ex/ϵ-члена (быстро затухающая на масштабе O(ϵ)O(\epsilon)O(ϵ)).
Регулярная теория возмущений
Подставляем разложение в степени ϵ\epsilonϵ: y=y0+ϵy1+…y=y_0+\epsilon y_1+\dotsy=y0 +ϵy1 +. При ϵ0\epsilon^0ϵ0 получаем редуцированное уравнение первого порядка
y0′+y0=0⇒y0=C e−x. y_0'+y_0=0\quad\Rightarrow\quad y_0=C\,e^{-x}.
y0 +y0 =0y0 =Cex.
Регулярный подход даёт только медленную компоненту e−xe^{-x}ex. Он годится, если задача задаётся, например, как начальная задача в точке x=0x=0x=0 и задано значение y(0)y(0)y(0): тогда постоянная CCC определяется и решение не требует внутренних слоёв. Но регулярная теория не в состоянии получить быстро меняющийся член e−x/ϵe^{-x/\epsilon}ex/ϵ, поэтому она ломается, когда требуется учесть пограничное условие, несовместное с y0y_0y0 .
Сингулярная (matched asymptotics) теория и слой пограничного типа
Поскольку при ϵ→0\epsilon\to0ϵ0 порядок уравнения падает (высший производный умножен на малый множитель), задача сингулярно возмущена и может требовать пограничного слоя. Рассмотрим отрезок [0,1][0,1][0,1] с граничными условиями y(0)=α, y(1)=βy(0)=\alpha,\;y(1)=\betay(0)=α,y(1)=β. Внешний (outer) разложение даёт
yout(x)=C e−x. y_{\text{out}}(x)=C\,e^{-x}.
yout (x)=Cex.
Чтобы удовлетворить условию в правой границе, выбираем C=βe1C=\beta e^{1}C=βe1, так что yout(1)=βy_{\text{out}}(1)=\betayout (1)=β. Обычно тогда yout(0)=βe1≠αy_{\text{out}}(0)=\beta e^{1}\neq\alphayout (0)=βe1=α и требуется пограничный слой у левой границы x=0x=0x=0.
Внутреннее (inner) разложение: вводим растянутую переменную X=x/ϵX=x/\epsilonX=x/ϵ. Тогда на ведущем порядке получаем
YXX+YX=0, Y_{XX}+Y_X=0,
YXX +YX =0,
общий вид решения
Y(X)=A+Be−X. Y(X)=A+B e^{-X}.
Y(X)=A+BeX.
Условие совпадения с внешним решением при X→∞X\to\inftyX даёт A=yout(0)=βe1A=y_{\text{out}}(0)=\beta e^{1}A=yout (0)=βe1. Налагая Y(0)=αY(0)=\alphaY(0)=α получаем B=α−βe1B=\alpha-\beta e^{1}B=αβe1. Возвращаясь в исходную переменную, слой имеет ширину O(ϵ)O(\epsilon)O(ϵ) и вклад вида (α−βe1)e−x/ϵ(\alpha-\beta e^{1})e^{-x/\epsilon}(αβe1)ex/ϵ.
Композиционное (составное) приближение на ведущем порядке:
ycomp(x)=βe1−x+(α−βe1)e−x/ϵ. y_{\text{comp}}(x)=\beta e^{1-x}+(\alpha-\beta e^{1})e^{-x/\epsilon}.
ycomp (x)=βe1x+(αβe1)ex/ϵ.
Это согласуется с точным решением (линейной комбинацией двух экспонент) и показывает явную роль слоя.
Когда требуется слой:
- Если задача задаётся на интервале с граничными условиями на обоих концах (BVP), то регулярное решение первого порядка может удовлетворять максимум одному граничному условию; второе обычно требует слоя.
- Положение слоя определяется знаком коэффициента при y′y'y (напр., здесь коэффициент +1>0+1>0+1>0 — слой у левого конца x=0x=0x=0); в общем для ϵy′′+ay′+…\epsilon y''+a y'+\dotsϵy′′+ay+ слой располагается на том конце, откуда «входят характеристики», т.е. где решение быстро затухает внутрь (если a>0a>0a>0 — слой у левого конца, если a<0a<0a<0 — у правого).
- Для начальной задачи (IVP) с заданным y(0)y(0)y(0) и/или y′(0)y'(0)y(0) слой может не требоваться: регулярная разложение часто достаточно.
Краткое резюме
- Уравнение сингулярно возмущено: общее решение содержит медленную компоненту e−xe^{-x}ex и быстрый слой e−x/ϵe^{-x/\epsilon}ex/ϵ.
- Регулярная теория даёт только медленную часть и годится, когда граничные/начальные условия совместимы с редуцированным (первого порядка) уравнением.
- Сингулярная (matched asymptotics) теория вводит внутреннюю переменную X=x/ϵX=x/\epsilonX=x/ϵ, строит слой толщиной O(ϵ)O(\epsilon)O(ϵ) и составное приближение; слой требуется при BVP, если внешнее решение не может удовлетворить все граничные условия.
1 Дек в 02:43
Не можешь разобраться в этой теме?
Обратись за помощью к экспертам
Гарантированные бесплатные доработки в течение 1 года
Быстрое выполнение от 2 часов
Проверка работы на плагиат
Поможем написать учебную работу
Прямой эфир