ODE 练习题
ODE 练习题：爆破估计、比较定理、Gronwall 与能量估计等。
好的，同学你好！很高兴能以常微分方程课程教师的身份为你提供帮助。
关于解的存在区间估计（Estimation of Interval of Existence），这是一个非常关键的话题。很多时候我们无法求出方程的解析解，但我们需要知道这个解能在多大的范围内保持有效（不发生爆破）。比较定理（Comparison Theorem）及其相关技巧是处理这类问题的核心工具。
下面我为你精心设计了4道题目。这些题目涵盖了从基础的直接比较、构造辅助函数到利用 Gronwall 不等式的不同侧面，解法各有千秋，难度适中。
习题一：基础非线性项的爆破估计（直接比较法）
题目： 考虑初值问题： \(\frac{dy}{dt} = t^2 + y^2, \quad y(0) = 0\) 请利用比较定理证明：该方程的解 \(y(t)\) 必定在有限时间内爆破（即解的存在区间是有界的），并给出一个爆破时间的上界估计。
提示与思路： 这个方程的右端项 \(f(t, y) = t^2 + y^2\) 既包含 \(t\) 也包含 \(y\)，很难求解析解。
- 观察： 显然 \(t^2 \ge 0\)，所以 \(\frac{dy}{dt} \ge y^2\)。但是直接解 \(z' = z^2, z(0)=0\) 会得到 \(z(t) \equiv 0\)，这对证明爆破没有帮助（因为不等号方向只能说明 \(y \ge 0\)）。
- 技巧： 寻找一个使得 \(y\) 变得很大的时刻。当 \(t\) 稍微大一点时，比如 \(t \ge 1\)，方程的性质会如何？
- 构造比较方程： 考虑在某个 \(t_0 > 0\) 之后，利用 \(y(t_0) > 0\) 的事实，构造一个形式为 \(u' = u^2\) 的方程来“从下方”顶住 \(y(t)\)，从而迫使 \(y(t)\) 爆破。
习题二：解的整体存在性（线性增长条件与比较定理）
题目： 考虑初值问题： \(\frac{dy}{dt} = \sin(t) \cdot y + \cos(y) \cdot e^{-t}, \quad y(0) = 1\) 请证明：该初值问题的解 \(y(t)\) 在区间 \([0, +\infty)\) 上整体存在。
提示与思路： 如果一个微分方程的右端项满足线性增长条件（即 \(|f(t,y)| \le A(t)|y| + B(t)\)），解通常不会在有限时间内爆破。
- 取绝对值： 考察 \(|y(t)|\) 的增长速度。我们通常计算 \(\frac{d}{dt}(y^2) = 2y y'\) 或者直接对 \(y'\) 放缩。
- 放缩： 利用 \(|\sin(t)| \le 1\) 和 \(|\cos(y)| \le 1\)。
- 应用比较原理： 你会得到一个关于 \(|y|\) 的微分不等式，形式类似于 \(|y|' \le k_1 |y| + k_2\)。
- 结论： 比较方程是一个线性方程，线性方程的解是指数增长的，虽然快，但是处处有定义（不会垂直渐近）。
习题三：利用积分不等式（Gronwall 形式的变体）
题目： 设 \(y(t)\) 是以下初值问题的非负解： \(\frac{dy}{dt} = \sqrt{y} + y, \quad y(0) = y_0 > 0\) 请估计解 \(y(t)\) 在区间 \([0, T]\) 上的上界。特别地，通过比较定理说明解是否存在爆破点？
提示与思路： 这里右端项 \(f(y) = \sqrt{y} + y\) 是非线性的。
- 直接比较的陷阱： 如果你简单放缩 \(\sqrt{y} \le y\) (当 \(y\) 很大时成立) 或 \(\sqrt{y} \le 1+y\)，你可以将问题转化为线性的。
- 另一种视角（Osgood 判别法）： 考察积分 \(\int_{y_0}^{\infty} \frac{ds}{\sqrt{s} + s}\) 的敛散性。
- 如果积分发散（趋于无穷），则解整体存在。
- 如果积分收敛（有限值），则解在有限时间爆破。
- 计算这个广义积分来看看结果。
习题四：双边比较（夹逼思想确定区间特征）
题目： 考虑初值问题： \(\frac{dy}{dt} = 1 + y^4, \quad y(0) = 0\) 令该解的最大存在区间为 \([0, T_{max})\)。请利用比较定理，通过构造两个简单的辅助方程，给出 \(T_{max}\) 的上下界估计（即找到常数 \(A, B\) 使得 \(A < T_{max} < B\)）。
提示与思路： 我们需要分别从上方和下方控制 \(y(t)\) 的增长速度。
- 求下界 \(A\)（说明解一定会爆破）： 需要一个比原方程增长得慢的方程，如果慢的方程都爆破了，快的肯定也爆破了。
- 显而易见，\(1+y^4 > y^4\)。考察 \(u' = u^4\)？注意初值是0，直接用 \(u^4\) 启动不了。
- 可以考虑 \(1+y^4 \ge 1\) (前期) 和 \(1+y^4 > y^4\) (后期)？或者直接用 \(y' \ge 1+y^2\) (利用 \(y^4\) 和 \(y^2\) 的关系，注意区间)？
- 更简单的做法： 直接积分 \(\int_0^{\infty} \frac{dy}{1+y^4} = T_{max}\)。这道题其实可以通过定积分的估计来做。
- 比较对象：
- 为了找 \(T_{max}\) 的上界（即证明解爆破得很快），利用 \(1+y^4 \ge 1+y^2\)（当 \(|y| \ge 1\)）或者简单地利用 \(\frac{dy}{dt} = 1+y^4\) 的解析积分形式。
- 为了找 \(T_{max}\) 的下界（即证明解至少能存活一段时间），利用 \(1+y^4 \le (1+y^2)^2\) 或者利用 \(\frac{dy}{dt} \le C(1+y^4)\) 的放缩？
- 推荐思路：
- 因为 \(y' > y^2\) (当 \(y\) 较大且 \(y(0)\) 调整后) -> 导致快速爆破。
- 因为 \(y' = 1+y^4 < (1+y)^4\) -> 比较方程 \(u' = (1+u)^4\)。解出 \(u(t)\) 的爆破时间 \(T_1\)，那么原方程的 \(T_{max}\) 肯定比 \(T_1\) 长（因为原方程长得慢）。
教师寄语：
做这些题目时，请记住核心心法：
- 寻找控制函数：对于 \(y' = f(t,y)\)，试着找 \(g(y)\) 和 \(h(y)\) 使得 \(g(y) \le f(t,y) \le h(y)\)。
- 解的寿命：
- 如果你被一个超线性（如 \(y^2, y^3\)）的方程从下方顶住，你往往会死掉（爆破）。
- 如果你被一个线性（如 \(ky\)）的方程从上方盖住，你往往能长生不老（整体存在）。
希望这些题目能帮你巩固知识！如果你对其中某道题卡住了，随时来问我具体的解答步骤。
下面给出一组“用比较定理/微分不等式/构造上下解/Grönwall不等式/能量估计”等方式来**估计解的存在区间（或保证不爆破的时间下界）**的题目。每题都刻意让方法侧重点不同；难度中等，避免一眼秒杀，但不需要很深的技巧。默认讨论经典解，且解在其最大存在区间内满足方程。
题目1（用上解导出爆破时间下界：与可解方程比较）
初值问题
\[y'=y^2+\sin t,\qquad y(0)=1.\]- 证明解在某个正时间区间上存在且满足一个显式上界。
- 利用比较的方法给出解在 \(t>0\) 方向的存在时间下界（即保证至少存在到某个 \(T>0\)）。
提示：用 \(\sin t\le 1\)，与 \(u'=u^2+1\) 或 \(u'=u^2\) 比较，得到可积分的不等式。
题目2（用下解保证单调性 + 再比较控制增长）
\[y' = y^2 - t,\qquad y(0)=0.\]- 证明在某个区间内 \(y(t)\le 0\)。
- 进一步用比较定理估计最大存在区间至少包含 \([0,T]\) 的某个显式 \(T>0\)。
提示：先从 \(y'=y^2-t\le y^2\) 入手不够，需要结合 \(t\ge 0\) 产生的“负反馈”；可尝试证明 \(y(t)\ge -\sqrt t\) 之类的粗界，或构造合适的比较函数。
题目3（用Grönwall：把非线性压成线性不等式）
\[y' = (1+t)\,y + y^2,\qquad y(0)=\tfrac12.\]证明存在 \(T>0\) 使得解在 \([0,T]\) 上存在，且给出一个显式的 \(T\) 下界。
提示：在短时间内先假设 \(|y|\le 1\) 并自洽验证；或写成
\[y' \le (1+t)y + y^2 \le (1+t)y + y\]（在 \(y\le 1\) 时），然后用Grönwall得到先验界，再回代闭合估计。
题目4（用积分不等式 + 反证法控制不爆破）
\[y' = e^{t}\,(1+y^2),\qquad y(0)=0.\]- 用比较定理证明 \(y(t)\) 在正向会爆破。
- 估计爆破时间 \(T^*\) 的上下界（给出两个显式数 \(0<T_-\le T^*\le T_+\)）。
提示：与 \(u' = e^{t}u^2\) 或令 \(\arctan y\) 变换：\((\arctan y)' = e^t\)。
（这一题会用到“变换 + 比较”，和前面方法不同。）
题目5（用“夹逼比较”：上下两侧均比较，估计存在区间）
\[y' = y^3 + \cos t,\qquad y(0)=0.\]证明解在 \([0,T]\) 上存在并给出一个显式 \(T>0\)，并且在该区间给出 \(y(t)\) 的显式夹逼
\[\ell(t)\le y(t)\le u(t).\]提示：\(\cos t\in[-1,1]\)，分别与
\[u' = u^3 + 1,\quad \ell' = \ell^3 - 1\]比较。注意三次方程可分离变量积分，至少能得到“存在到某个时间”的显式下界。
题目6（用能量型估计：二阶方程转为一阶不等式）
考虑二阶初值问题
\[x'' + x = (x')^2,\qquad x(0)=0,\ x'(0)=1.\]- 证明解在 \(t\ge 0\) 的某个区间存在且 \(x'(t)>0\)。
- 通过比较或构造能量函数估计其最大存在时间的一个正下界。
提示：设 \(v=x'\)，得到系统；或考虑 \(v' = v^2 - x\)，再用 \(x\ge 0\)（需证明）得到 \(v'\le v^2\) 之类的控制。
题目7（用“延拓判据 + 有界性”：证明全局存在的粗方法）
\[y' = \frac{y}{1+t^2} - y^3,\qquad y(0)=1.\]证明解在 \([0,\infty)\) 上全局存在，并给出一个简单的上界（如 \(y(t)\le 1\) 或更好）。
提示：证明 \(y\) 不会变得无界：当 \(|y|\) 大时 \(-y^3\) 主导；可用比较 \(y'\le \frac{y}{1+t^2}\) 得上界，再用延拓定理（解在有限时间爆破必须无界）。
题目8（用局部Lipschitz给局部存在，再用比较估计“能延多远”）
\[y' = (t+1)\sqrt{1+y^2},\qquad y(0)=0.\]- 证明该问题在全实轴上解唯一。
- 用比较定理或变量代换估计 \(y(t)\) 的增长并证明解在 \([0,\infty)\) 全局存在。
提示：与 \(\sqrt{1+y^2}\le 1+|y|\) 比较可导出线性不等式；或用 \(\operatorname{arsinh}(y)\) 变换：\((\operatorname{arsinh} y)' = t+1\)。
题目9（用“分段比较”：先控在小区间，再扩大区间）
\[y' = y^2 + t^2,\qquad y(0)=0.\]给出一个分段的存在时间下界：先证明在 \([0,1]\) 上 \(y(t)\le C\)，再由此推出可至少延拓到 \([0,2]\)（或给出一个更一般的迭代延拓策略）。
提示：先用 \(t^2\le 1\) 在 \([0,1]\) 比较；得到 \(y\) 在 \([0,1]\) 上有界后，再利用 \(t^2\le 4\) 在 \([0,2]\) 上比较，形成“逐段”存在估计。
题目10（用反向比较：从下界推出有限时间爆破）
\[y' = y^2 + e^{-t},\qquad y(0)=1.\]证明解在正向有限时间爆破，并给出爆破时间的一个显式上界。
提示：因为 \(e^{-t}\ge 0\)，有 \(y'\ge y^2\)，与 \(u'=u^2\) 比较即可给出爆破上界；但要注意严格比较条件与初值一致。
如果你愿意，我可以把这些题按课堂作业/小测/期末复习三个梯度重新编排；也可以从中任选 3–5 题给出“标准解答”（详细写出比较函数选择、适用的比较定理版本、以及存在区间如何从先验界推出）。你更希望偏向“证明全局存在”还是“估计爆破时间/存在时间下界”？
你好！很高兴能扮演你的常微分方程教师。针对数学专业的同学，掌握各类非线性方程的变换技巧以及高阶变系数线性方程的解法是核心基本功。
这里我为你精心挑选了 6道 题目。这些题目涵盖了你要求的Bernoulli方程、Ricatti方程、Euler方程，同时也加入了一些需要技巧的全微分方程和高阶降阶问题。难度适中，旨在考察对标准解法和变换技巧的熟练度，计算量我都控制在合理范围内。
请准备好纸笔，我们开始吧！
第一组：一阶非线性方程的转化
题目 1：Bernoulli 方程 求下列方程的通解： \(y' + \frac{y}{x} = x^2 y^4, \quad (x > 0)\)
考察点： 这是标准的Bernoulli方程 (\(y' + P(x)y = Q(x)y^n\))。 提示： 注意 \(n=4\)。通常令 \(z = y^{1-n}\) 将其转化为线性方程。
题目 2：Riccati 方程 已知方程 \(y' = y^2 - \frac{2}{x^2}\) 有一个特解形式为 \(y = \frac{A}{x}\)（其中 \(A\) 为常数），试求该方程的通解。
考察点： Riccati方程一般不可积，但已知一个特解后可以通过变换 \(y = y^* + \frac{1}{u}\) 化为线性方程，或者 \(y = y^* + z\) 化为Bernoulli方程。 提示： 先代入形式解求出 \(A\)（会有两个值，取简单的即可），然后做变换。
第二组：全微分方程与积分因子
题目 3：分组凑全微分 求解方程： \((x^2 + y^2 + x) dx + y dy = 0\)
考察点： 并不是所有的全微分方程都要死板地求积分因子。数学专业的学生应具备敏锐的“凑微分”直觉。 提示： 观察 \(x^2+y^2\) 和 \(xdx + ydy\) 之间的关系。尝试将 \(xdx+ydy\) 写成 \(d(\dots)\) 的形式。
第三组：高阶变系数线性方程
题目 4：Euler 方程 求下列方程在 \(x>0\) 时的通解： \(x^2 y'' - 3xy' + 4y = x^3\)
考察点： Euler方程是变系数转化为常系数的经典例子。 提示： 做变换 \(x = e^t\)。注意二阶导数 \(y''\) 变换后的形式不仅仅是 \(\frac{d^2y}{dt^2}\)，还有一项一阶导数。
题目 5：不显含因变量的降阶法 求解初值问题： \((1+x^2)y'' + (y')^2 + 1 = 0, \quad y(0)=0, \ y'(0)=1\)
考察点： 方程中不显含 \(y\)，属于 \(F(x, y', y'')=0\) 类型。 提示： 令 \(p = y'\)，降阶为关于 \(p\) 的一阶方程。这是一个分离变量方程。
第四组：常系数线性方程系统的解结构
题目 6：待定系数法求特解（共振情形） 求如下二阶常系数非齐次线性方程的通解： \(y'' - 2y' + y = \frac{e^x}{x}\)
考察点： 虽然是常系数，但右端项 \(f(x)\) 不是常见的多项式、指数或三角函数形式，或者是“共振”情况下的特殊形式。 提示： 注意特征根的情况。对于 \(e^x/x\) 这种形式，常规待定系数法失效，请使用 常数变易法 (Variation of Parameters)。
教师寄语与简要参考思路
做完之后，可以对照一下简单的思路，看看你的方向是否正确：
- Bernoulli方程：令 \(z = y^{-3}\)，化简为 \(z' - \frac{3}{x}z = -3x^2\)。
- Riccati方程：代入 \(A/x\) 可得 \(A=2\) 或 \(A=-1\)。取 \(y_1 = \frac{2}{x}\)，令 \(y = \frac{2}{x} + \frac{1}{u}\)，化为关于 \(u\) 的一阶线性方程。
- 全微分：方程可重写为 \((x^2+y^2)dx + (xdx+ydy) = 0 \Rightarrow (x^2+y^2)dx + \frac{1}{2}d(x^2+y^2) = 0\)。两边同除以 \((x^2+y^2)\) 即可分离变量。
- Euler方程：特征方程为 \(r(r-1) - 3r + 4 = 0 \Rightarrow (r-2)^2=0\)。齐次解 \(C_1 x^2 + C_2 x^2 \ln x\)。特解设为 \(Ax^3\) 代入求系数。
- 降阶法：令 \(y'=p\)，则 \((1+x^2)p' = -(1+p^2)\)。分离变量解出 \(p(x)\)，再积分一次求 \(y\)。
- 常数变易法：齐次解为 \((C_1 + C_2 x)e^x\)。设特解 \(y^* = u_1(x)e^x + u_2(x)xe^x\)，利用Wronskian行列式求解 \(u_1', u_2'\)。
祝你解题愉快！如果有哪一道题卡住了，随时来问我具体的步骤。