Исследуйте поведение решений обыкновенного дифференциального уравнения y' = y^2 при разных начальных условиях; обсудите существование глобальных решений, конечное время взрыва и методы оценки времени разрыва
Решение уравнения y′=y2y' = y^2y′=y2. 1) Явное решение (сепарация переменных): dyy2=dt⇒−1y=t+C.
\frac{dy}{y^2}=dt \quad\Rightarrow\quad -\frac{1}{y}=t+C. y2dy=dt⇒−y1=t+C.
При начальном условии y(t0)=y0y(t_0)=y_0y(t0)=y0 получаем −1y(t)=t−t0−1y0⇒y(t)=11y0−(t−t0).
-\frac{1}{y(t)}=t-t_0-\frac{1}{y_0}\quad\Rightarrow\quad y(t)=\frac{1}{\dfrac{1}{y_0}-(t-t_0)}. −y(t)1=t−t0−y01⇒y(t)=y01−(t−t0)1. 2) Существование и конечное время взрыва (blow-up): - Если y0>0y_0>0y0>0, знаменатель обнуляется при T+=t0+1y0,
T^+=t_0+\frac{1}{y_0}, T+=t0+y01,
поэтому решение стремится к +∞+\infty+∞ при t→T+−t\to T^+-t→T+−. Решение определено на промежутке (−∞,T+)(-\infty,T^+)(−∞,T+) (вперёд до T+T^+T+ — конечный взрыв, назад бесконечно). - Если y0=0y_0=0y0=0, то y(t)≡0y(t)\equiv0y(t)≡0 глобально для всех ttt. - Если y0<0y_0<0y0<0, то для t>t0t>t_0t>t0 знаменатель никогда не равен нулю (обнулится при t=t0+1/y0<t0t=t_0+1/y_0<t_0t=t0+1/y0<t0), поэтому решение существует глобально вперёд: y(t)→0−y(t)\to0^-y(t)→0− при t→+∞t\to+\inftyt→+∞. Влево по времени оно взрывается в конечный момент T−=t0+1y0<t0T^-=t_0+\dfrac{1}{y_0}<t_0T−=t0+y01<t0. 3) Оценки и методы определения времени разрыва: - Интегральный тест (общий приём для y′=f(y)y'=f(y)y′=f(y)): время до взрыва определяется интегралом T+−t0=∫y0+∞duf(u).
T^+-t_0=\int_{y_0}^{+\infty}\frac{du}{f(u)}. T+−t0=∫y0+∞f(u)du.
Для f(u)=u2f(u)=u^2f(u)=u2 это даёт ∫y0∞u−2du=1y0,
\int_{y_0}^{\infty}u^{-2}du=\frac{1}{y_0}, ∫y0∞u−2du=y01,
что совпадает с явной формулой. - Сравнение/неравенства Гронваля: если известны оценки a(y)≤f(y)≤b(y)a(y)\le f(y)\le b(y)a(y)≤f(y)≤b(y), то времена разрыва можно оценить интегралами ∫dy/a(y)\int dy/a(y)∫dy/a(y) и ∫dy/b(y)\int dy/b(y)∫dy/b(y). - Поведение у границы взрыва: при t→T−t\to T^-t→T−y(t)∼1T−t,
y(t)\sim\frac{1}{T-t}, y(t)∼T−t1,
то есть типичный обратнопропорциональный профиль взрыва. 4) Замечания: - Поскольку f(y)=y2f(y)=y^2f(y)=y2 локально липшицева, решение единственно до момента возможного взрыва; единственная причина конечного времени существования — рост решения до бесконечности. - Аналогичные рассуждения применимы к нелинейностям y′=ypy'=y^py′=yp (p>1p>1p>1) с тем, что интеграл ∫y0∞u−pdu\int_{y_0}^{\infty}u^{-p}du∫y0∞u−pdu даёт конечное или бесконечное время взрыва в зависимости от ppp. Если нужно, могу привести краткие численные оценки для приближённого определения TTT при заданном y0y_0y0 или показать графики поведения решений.
1) Явное решение (сепарация переменных):
dyy2=dt⇒−1y=t+C. \frac{dy}{y^2}=dt \quad\Rightarrow\quad -\frac{1}{y}=t+C.
y2dy =dt⇒−y1 =t+C. При начальном условии y(t0)=y0y(t_0)=y_0y(t0 )=y0 получаем
−1y(t)=t−t0−1y0⇒y(t)=11y0−(t−t0). -\frac{1}{y(t)}=t-t_0-\frac{1}{y_0}\quad\Rightarrow\quad
y(t)=\frac{1}{\dfrac{1}{y_0}-(t-t_0)}.
−y(t)1 =t−t0 −y0 1 ⇒y(t)=y0 1 −(t−t0 )1 .
2) Существование и конечное время взрыва (blow-up):
- Если y0>0y_0>0y0 >0, знаменатель обнуляется при
T+=t0+1y0, T^+=t_0+\frac{1}{y_0},
T+=t0 +y0 1 , поэтому решение стремится к +∞+\infty+∞ при t→T+−t\to T^+-t→T+−. Решение определено на промежутке (−∞,T+)(-\infty,T^+)(−∞,T+) (вперёд до T+T^+T+ — конечный взрыв, назад бесконечно).
- Если y0=0y_0=0y0 =0, то y(t)≡0y(t)\equiv0y(t)≡0 глобально для всех ttt.
- Если y0<0y_0<0y0 <0, то для t>t0t>t_0t>t0 знаменатель никогда не равен нулю (обнулится при t=t0+1/y0<t0t=t_0+1/y_0<t_0t=t0 +1/y0 <t0 ), поэтому решение существует глобально вперёд: y(t)→0−y(t)\to0^-y(t)→0− при t→+∞t\to+\inftyt→+∞. Влево по времени оно взрывается в конечный момент T−=t0+1y0<t0T^-=t_0+\dfrac{1}{y_0}<t_0T−=t0 +y0 1 <t0 .
3) Оценки и методы определения времени разрыва:
- Интегральный тест (общий приём для y′=f(y)y'=f(y)y′=f(y)): время до взрыва определяется интегралом
T+−t0=∫y0+∞duf(u). T^+-t_0=\int_{y_0}^{+\infty}\frac{du}{f(u)}.
T+−t0 =∫y0 +∞ f(u)du . Для f(u)=u2f(u)=u^2f(u)=u2 это даёт
∫y0∞u−2du=1y0, \int_{y_0}^{\infty}u^{-2}du=\frac{1}{y_0},
∫y0 ∞ u−2du=y0 1 , что совпадает с явной формулой.
- Сравнение/неравенства Гронваля: если известны оценки a(y)≤f(y)≤b(y)a(y)\le f(y)\le b(y)a(y)≤f(y)≤b(y), то времена разрыва можно оценить интегралами ∫dy/a(y)\int dy/a(y)∫dy/a(y) и ∫dy/b(y)\int dy/b(y)∫dy/b(y).
- Поведение у границы взрыва: при t→T−t\to T^-t→T− y(t)∼1T−t, y(t)\sim\frac{1}{T-t},
y(t)∼T−t1 , то есть типичный обратнопропорциональный профиль взрыва.
4) Замечания:
- Поскольку f(y)=y2f(y)=y^2f(y)=y2 локально липшицева, решение единственно до момента возможного взрыва; единственная причина конечного времени существования — рост решения до бесконечности.
- Аналогичные рассуждения применимы к нелинейностям y′=ypy'=y^py′=yp (p>1p>1p>1) с тем, что интеграл ∫y0∞u−pdu\int_{y_0}^{\infty}u^{-p}du∫y0 ∞ u−pdu даёт конечное или бесконечное время взрыва в зависимости от ppp.
Если нужно, могу привести краткие численные оценки для приближённого определения TTT при заданном y0y_0y0 или показать графики поведения решений.