解法选择先于代数运算
同是一阶方程,右端的结构不同,可靠的变形路线也不同。开始积分前先完成三项检查:确定自变量和未知函数,写出方程与初值共同有效的区域,再辨认它是否可分离、线性或恰当。一个方程可能同时属于多类,此时应选步骤更短且不易丢解的方法。
常见识别信号如下:
y ′ = g ( t ) h ( y ) y'=g(t)h(y) y ′ = g ( t ) h ( y ) 具有自变量与状态变量的乘积结构,是可分离方程;
y ′ + p ( t ) y = q ( t ) y'+p(t)y=q(t) y ′ + p ( t ) y = q ( t ) 对 y y y 与 y ′ y' y ′ 都是一次,是一阶线性方程;
M ( t , y ) d t + N ( t , y ) d y = 0 M(t,y)\,\mathrm dt+N(t,y)\,\mathrm dy=0 M ( t , y ) d t + N ( t , y ) d y = 0 若来自某个势函数的全微分,就是恰当方程;
非线性式 y ′ + p ( t ) y = q ( t ) y α y'+p(t)y=q(t)y^\alpha y ′ + p ( t ) y = q ( t ) y α 可用幂代换化为线性方程;
y ′ = F ( y / t ) y'=F(y/t) y ′ = F ( y / t ) 在 t ≠ 0 t\ne0 t = 0 的区间上可用比例 v = y / t v=y/t v = y / t 降为可分离方程。
所有方法最后都要做同一件事:把结果代回原方程和初值,并给出包含初始时刻的连通定义区间。积分常数正确不等于区间自动正确。
可分离方程:除法前先保存平衡解
可分离方程
若方程可写成
y ′ = g ( t ) h ( y ) , y'=g(t)h(y), y ′ = g ( t ) h ( y ) , 则称它可分离。在 h ( y ) ≠ 0 h(y)\ne0 h ( y ) = 0 的区域内可写为
d y h ( y ) = g ( t ) d t , \frac{\mathrm dy}{h(y)}=g(t)\,\mathrm dt, h ( y ) d y = g ( t ) d t , 从而得到隐式关系
∫ d y h ( y ) = ∫ g ( t ) d t + C . \int\frac{\mathrm dy}{h(y)}
=\int g(t)\,\mathrm dt+C. ∫ h ( y ) d y = ∫ g ( t ) d t + C . 每个满足 h ( y ∗ ) = 0 h(y_*)=0 h ( y ∗ ) = 0 的常值函数 y ≡ y ∗ y\equiv y_* y ≡ y ∗ 都应在除以 h ( y ) h(y) h ( y ) 前单独记录。
把 d y \mathrm dy d y 和 d t \mathrm dt d t 分列两侧是一种简洁记号,其合法依据可由链式法则说明。若
H ′ ( y ) = 1 / h ( y ) H'(y)=1/h(y) H ′ ( y ) = 1/ h ( y ) 、G ′ ( t ) = g ( t ) G'(t)=g(t) G ′ ( t ) = g ( t ) ,那么沿非平衡解
d d t H ( y ( t ) ) = y ′ ( t ) h ( y ( t ) ) = g ( t ) = G ′ ( t ) , \frac{\mathrm d}{\mathrm dt}H(y(t))
=\frac{y'(t)}{h(y(t))}=g(t)=G'(t), d t d H ( y ( t )) = h ( y ( t )) y ′ ( t ) = g ( t ) = G ′ ( t ) ,
故 H ( y ( t ) ) − G ( t ) H(y(t))-G(t) H ( y ( t )) − G ( t ) 为常数。这个论证也明确了 h ( y ) ≠ 0 h(y)\ne0 h ( y ) = 0 的适用范围。
例 1:Logistic 方程中的平衡解与分支
求解
y ′ = y ( 1 − y ) , y ( 0 ) = 1 2 . y'=y(1-y),\qquad y(0)=\frac12. y ′ = y ( 1 − y ) , y ( 0 ) = 2 1 . 先记录平衡解 y ≡ 0 y\equiv0 y ≡ 0 与 y ≡ 1 y\equiv1 y ≡ 1 。给定初值位于两者之间,唯一性保证对应解不会穿过平衡线,于是可在 0 < y < 1 0<y<1 0 < y < 1 中分离:
1 y ( 1 − y ) d y = d t . \frac{1}{y(1-y)}\,\mathrm dy=\mathrm dt. y ( 1 − y ) 1 d y = d t . 利用
1 y ( 1 − y ) = 1 y + 1 1 − y \frac1{y(1-y)}=\frac1y+\frac1{1-y} y ( 1 − y ) 1 = y 1 + 1 − y 1 积分得
ln y 1 − y = t + C . \ln\frac{y}{1-y}=t+C. ln 1 − y y = t + C . 初值给出 C = 0 C=0 C = 0 ,所以
y ( t ) = e t 1 + e t = 1 1 + e − t . y(t)=\frac{e^t}{1+e^t}=\frac1{1+e^{-t}}. y ( t ) = 1 + e t e t = 1 + e − t 1 . 分母始终为正,且 0 < y ( t ) < 1 0<y(t)<1 0 < y ( t ) < 1 ,故最大存在区间为全部实数。求导可得
y ′ = e − t / ( 1 + e − t ) 2 = y ( 1 − y ) y'=e^{-t}/(1+e^{-t})^2=y(1-y) y ′ = e − t / ( 1 + e − t ) 2 = y ( 1 − y ) ,完成回代。
如果初值恰是零或一,分离时的除法会排除该点,但原方程仍有对应常值解。若初值在其他区域,积分后的绝对值和代数分支也要结合初值选择,不能在最后任意添加正负号。
一阶线性方程:积分因子来自乘积求导
一阶线性方程与积分因子
在 p , q p,q p , q 连续的区间 I I I 上,标准形式
y ′ + p ( t ) y = q ( t ) y'+p(t)y=q(t) y ′ + p ( t ) y = q ( t ) 的一阶线性方程可取积分因子
μ ( t ) = exp ( ∫ p ( t ) d t ) . \mu(t)=\exp\left(\int p(t)\,\mathrm dt\right). μ ( t ) = exp ( ∫ p ( t ) d t ) . 它满足 μ ′ = p μ \mu'=p\mu μ ′ = p μ ,因而
( μ y ) ′ = μ y ′ + μ ′ y = μ ( y ′ + p y ) = μ q . (\mu y)'=\mu y'+\mu' y
=\mu(y'+py)=\mu q. ( μ y ) ′ = μ y ′ + μ ′ y = μ ( y ′ + p y ) = μ q . 积分后得到
y ( t ) = 1 μ ( t ) ( C + ∫ μ ( t ) q ( t ) d t ) . y(t)=\frac{1}{\mu(t)}
\left(C+\int \mu(t)q(t)\,\mathrm dt\right). y ( t ) = μ ( t ) 1 ( C + ∫ μ ( t ) q ( t ) d t ) .
积分因子乘上任意非零常数不会改变最终解族。真正决定区间的是标准化之前是否除以了可能为零的系数,以及 p , q p,q p , q 在哪里连续。若原式为
a ( t ) y ′ + b ( t ) y = r ( t ) a(t)y'+b(t)y=r(t) a ( t ) y ′ + b ( t ) y = r ( t ) ,必须先限制到 a ( t ) ≠ 0 a(t)\ne0 a ( t ) = 0 的连通区间,再除以 a ( t ) a(t) a ( t ) 。
例 2:积分因子与初值常数
求解
y ′ + 2 y = e − t , y ( 0 ) = 0. y'+2y=e^{-t},\qquad y(0)=0. y ′ + 2 y = e − t , y ( 0 ) = 0. 取 μ = e 2 t \mu=e^{2t} μ = e 2 t ,则
( e 2 t y ) ′ = e 2 t e − t = e t . (e^{2t}y)'=e^{2t}e^{-t}=e^t. ( e 2 t y ) ′ = e 2 t e − t = e t . 积分后
e 2 t y = e t + C , y = e − t + C e − 2 t . e^{2t}y=e^t+C,
\qquad
y=e^{-t}+Ce^{-2t}. e 2 t y = e t + C , y = e − t + C e − 2 t . 由 y ( 0 ) = 0 y(0)=0 y ( 0 ) = 0 得 C = − 1 C=-1 C = − 1 ,故
y ( t ) = e − t − e − 2 t . y(t)=e^{-t}-e^{-2t}. y ( t ) = e − t − e − 2 t . 左侧回代为
y ′ + 2 y = ( − e − t + 2 e − 2 t ) + 2 ( e − t − e − 2 t ) = e − t . y'+2y=(-e^{-t}+2e^{-2t})+2(e^{-t}-e^{-2t})=e^{-t}. y ′ + 2 y = ( − e − t + 2 e − 2 t ) + 2 ( e − t − e − 2 t ) = e − t . 所有系数在实轴连续,因此最大存在区间为 ( − ∞ , ∞ ) (-\infty,\infty) ( − ∞ , ∞ ) 。
恰当方程:寻找等势曲线
把一阶方程写成微分形式
M ( t , y ) d t + N ( t , y ) d y = 0. M(t,y)\,\mathrm dt+N(t,y)\,\mathrm dy=0. M ( t , y ) d t + N ( t , y ) d y = 0.
若存在势函数 Φ ( t , y ) \Phi(t,y) Φ ( t , y ) 使
Φ t = M \Phi_t=M Φ t = M 、Φ y = N \Phi_y=N Φ y = N ,那么沿解曲线
d d t Φ ( t , y ( t ) ) = Φ t + Φ y y ′ = M + N y ′ = 0 , \frac{\mathrm d}{\mathrm dt}\Phi(t,y(t))
=\Phi_t+\Phi_y y'=M+Ny'=0, d t d Φ ( t , y ( t )) = Φ t + Φ y y ′ = M + N y ′ = 0 ,
所以解由等势曲线 Φ ( t , y ) = C \Phi(t,y)=C Φ ( t , y ) = C 描述。
恰当性检验与势函数
若 M , N M,N M , N 在单连通区域上具有连续一阶偏导,则微分形式恰当,当且仅当
恢复势函数时,可先对 M M M 关于 t t t 积分:
Φ ( t , y ) = ∫ M ( t , y ) d t + g ( y ) , \Phi(t,y)=\int M(t,y)\,\mathrm dt+g(y), Φ ( t , y ) = ∫ M ( t , y ) d t + g ( y ) , 再用 Φ y = N \Phi_y=N Φ y = N 求只依赖 y y y 的函数 g g g 。也可先积分 N N N ,两条路线应给出相差常数的势函数。
例 3:从交叉偏导恢复势函数
考虑
( 2 t y + y 2 ) d t + ( t 2 + 2 t y ) d y = 0 , y ( 1 ) = 1. (2ty+y^2)\,\mathrm dt+(t^2+2ty)\,\mathrm dy=0,
\qquad y(1)=1. ( 2 t y + y 2 ) d t + ( t 2 + 2 t y ) d y = 0 , y ( 1 ) = 1. 这里
M y = 2 t + 2 y = N t , M_y=2t+2y=N_t, M y = 2 t + 2 y = N t , 所以方程在整个平面恰当。对 M M M 关于 t t t 积分:
Φ = t 2 y + t y 2 + g ( y ) . \Phi=t^2y+ty^2+g(y). Φ = t 2 y + t y 2 + g ( y ) . 再求偏导得 Φ y = t 2 + 2 t y + g ′ ( y ) \Phi_y=t^2+2ty+g'(y) Φ y = t 2 + 2 t y + g ′ ( y ) 。与 N N N 比较可知 g ′ ( y ) = 0 g'(y)=0 g ′ ( y ) = 0 ,于是隐式解为
t 2 y + t y 2 = C . t^2y+ty^2=C. t 2 y + t y 2 = C . 代入 ( 1 , 1 ) (1,1) ( 1 , 1 ) 得 C = 2 C=2 C = 2 。在初始点
Φ y = N ( 1 , 1 ) = 3 ≠ 0 \Phi_y=N(1,1)=3\ne0 Φ y = N ( 1 , 1 ) = 3 = 0 ,隐式函数定理保证这条等势曲线在初始点附近可唯一写成 y ( t ) y(t) y ( t ) 。势函数方程本身可能包含多个代数分支,初值和局部唯一性负责选择正确分支。
如果 M y ≠ N t M_y\ne N_t M y = N t ,有时可乘积分因子使其恰当。只依赖 t t t 的积分因子存在的一种可检验情形是
M y − N t N = r ( t ) , \frac{M_y-N_t}{N}=r(t), N M y − N t = r ( t ) ,
此时可取 μ ( t ) = exp ( ∫ r ( t ) d t ) \mu(t)=\exp(\int r(t)\,\mathrm dt) μ ( t ) = exp ( ∫ r ( t ) d t ) ;只依赖 y y y 的情形相应要求
N t − M y M = s ( y ) , \frac{N_t-M_y}{M}=s(y), M N t − M y = s ( y ) ,
并取 μ ( y ) = exp ( ∫ s ( y ) d y ) \mu(y)=\exp(\int s(y)\,\mathrm dy) μ ( y ) = exp ( ∫ s ( y ) d y ) 。这些商式只在分母非零的区域内用于推导,得到积分因子后仍应直接检查新的交叉偏导。
代换的目标是回到已知类型
Bernoulli 方程
y ′ + p ( t ) y = q ( t ) y α , α ≠ 0 , 1 , y'+p(t)y=q(t)y^\alpha,
\qquad \alpha\ne0,1, y ′ + p ( t ) y = q ( t ) y α , α = 0 , 1 ,
看似非线性。在 y y y 的所选实数分支上令 v = y 1 − α v=y^{1-\alpha} v = y 1 − α ,则
v ′ + ( 1 − α ) p ( t ) v = ( 1 − α ) q ( t ) , v'+(1-\alpha)p(t)v=(1-\alpha)q(t), v ′ + ( 1 − α ) p ( t ) v = ( 1 − α ) q ( t ) ,
转化为一阶线性方程。若 y = 0 y=0 y = 0 在原方程中有意义,应在除以 y α y^\alpha y α 前单独检查;非整数幂还要求明确正负与实数定义域。
例 4:Bernoulli 方程化为线性方程
求解
y ′ + y = t y 2 , y ( 0 ) = 1. y'+y=ty^2,\qquad y(0)=1. y ′ + y = t y 2 , y ( 0 ) = 1. 在初值附近 y ≠ 0 y\ne0 y = 0 ,令 v = y − 1 v=y^{-1} v = y − 1 。因为 v ′ = − y ′ / y 2 v'=-y'/y^2 v ′ = − y ′ / y 2 ,原式除以 y 2 y^2 y 2 后化为
线性方程的积分因子为 e − t e^{-t} e − t :
( e − t v ) ′ = − t e − t . (e^{-t}v)'=-te^{-t}. ( e − t v ) ′ = − t e − t . 积分得到 e − t v = ( t + 1 ) e − t + C e^{-t}v=(t+1)e^{-t}+C e − t v = ( t + 1 ) e − t + C ,即
v = t + 1 + C e t v=t+1+Ce^t v = t + 1 + C e t 。初值给 v ( 0 ) = 1 v(0)=1 v ( 0 ) = 1 ,故 C = 0 C=0 C = 0 ,于是
y ( t ) = 1 t + 1 . y(t)=\frac1{t+1}. y ( t ) = t + 1 1 . 包含零且不穿过分母零点的最大区间为 ( − 1 , ∞ ) (-1,\infty) ( − 1 , ∞ ) 。回代有
y ′ + y = − 1 / ( t + 1 ) 2 + 1 / ( t + 1 ) = t / ( t + 1 ) 2 = t y 2 y'+y=-1/(t+1)^2+1/(t+1)=t/(t+1)^2=ty^2 y ′ + y = − 1/ ( t + 1 ) 2 + 1/ ( t + 1 ) = t / ( t + 1 ) 2 = t y 2 。
另一常见结构是 y ′ = F ( y / t ) y'=F(y/t) y ′ = F ( y / t ) 。在不含 t = 0 t=0 t = 0 的区间令 y = t v y=tv y = t v ,便有
y ′ = v + t v ′ y'=v+tv' y ′ = v + t v ′ ,从而
t v ′ = F ( v ) − v . tv'=F(v)-v. t v ′ = F ( v ) − v .
这通常可分离。代换并没有消除定义域问题:t = 0 t=0 t = 0 在原比例中就是特殊点,F ( v ) − v = 0 F(v)-v=0 F ( v ) − v = 0 对应的常比例解也必须在除法前保存。
方法选择与独立验算
实际方程不一定以标准形式出现。可靠的选择顺序是先做不改变解集的整理,再检查结构。若右端已经能分解为
g ( t ) h ( y ) g(t)h(y) g ( t ) h ( y ) ,可分离法通常最短;若展开后对 y y y 和 y ′ y' y ′ 都是一次,应优先标准化为线性方程;若方程自然给成
M d t + N d y = 0 M\,\mathrm dt+N\,\mathrm dy=0 M d t + N d y = 0 ,先比较交叉偏导往往比强行解出 y ′ y' y ′ 更有效。只有这些直接结构都不成立时,才寻找幂、比例、平移等代换。代换的价值必须由“新方程落入哪一种已知类型”说明,仅仅换了字母并不构成进展。
求出结果后,可以用两层检查减少符号错误。第一层是局部代数检查:对显式解求导,逐项代入原方程;对隐式解作隐式微分,确认重新得到原微分形式。第二层是问题条件检查:代入初值,确认积分常数和分支;列出对数真数、根式、分母、原方程系数与代换的非零条件;最后取包含初始时刻的最大连通区间。若解法中曾除以某个含 y y y 的量,还要返回原方程补查其零点对应的常值解。
不同方法也能互相核验。例如一阶线性方程乘积分因子后,本身变成恰当形式;可分离方程积分得到的关系
H ( y ) − G ( t ) = C H(y)-G(t)=C H ( y ) − G ( t ) = C 也是一个势函数。若同一方程既可分离又线性,两条路线所得公式应在重新命名常数后完全一致。出现不一致时,应优先检查部分分式符号、积分因子指数、初值代入和绝对值分支,而不是把差异解释成两个不同解族。
方程的量纲还能提供快速筛查。若 t t t 表示时间、y y y 表示浓度,则 y ′ y' y ′ 的量纲是浓度除以时间;方程右侧每一项都必须相同。积分因子本身应无量纲,因此指数
∫ p ( t ) d t \int p(t)\,\mathrm dt ∫ p ( t ) d t 必须无量纲。量纲一致不能证明推导正确,却能在回代前发现遗漏系数或把增长率误当总增量等错误。
四类高频错误
见到乘积就把变量分开
只有能写成 g ( t ) h ( y ) g(t)h(y) g ( t ) h ( y ) 的乘积才可直接分离。t + y t+y t + y 不是一个只含 t t t 的函数与一个只含 y y y 的函数之积;把
d y / ( t + y ) = d t \mathrm dy/(t+y)=\mathrm dt d y / ( t + y ) = d t 写成两边可独立积分是非法的。
除以含 y 的因子不会改变解集
除以 h ( y ) h(y) h ( y ) 会排除 h ( y ) = 0 h(y)=0 h ( y ) = 0 的状态,而这些状态常对应平衡解。正确顺序是先求零点并代回原方程,再在每个非零区域分离。
交叉偏导相等在任何区域都自动给全局势函数
单连通性排除了绕孔积分产生的全局障碍。在非单连通区域,交叉偏导相等只提供局部恰当性,还需检查沿闭曲线的积分或直接构造全局势函数。
显式公式可代回,不代表它跨奇点仍是一条解
y = 1 / ( t + 1 ) y=1/(t+1) y = 1/ ( t + 1 ) 在 t = − 1 t=-1 t = − 1 两侧分别满足同一微分方程,但它没有在该点定义。初值位于零时,只能选包含零的
( − 1 , ∞ ) (-1,\infty) ( − 1 , ∞ ) ,不能把两个分支合写成一个跨越奇点的解区间。
练习:分类、求解、回代、定区间
练习 1:反正切分支 标记完成
所属知识 可分离方程
难度 3/5 求解
y ′ = t ( 1 + y 2 ) , y ( 0 ) = 0 , y'=t(1+y^2),\qquad y(0)=0, y ′ = t ( 1 + y 2 ) , y ( 0 ) = 0 , 并给出最大存在区间。
查看提示 把
1 + y 2 1+y^2 1 + y 2 移到 dy 一侧,并用初值选择正切分支。
查看解答 方程可分离:
d y 1 + y 2 = t d t . \frac{\mathrm dy}{1+y^2}=t\,\mathrm dt. 1 + y 2 d y = t d t . 积分并代入初值得
arctan y = t 2 2 , y = tan t 2 2 . \arctan y=\frac{t^2}{2},
\qquad
y=\tan\frac{t^2}{2}. arctan y = 2 t 2 , y = tan 2 t 2 . 包含零的连续分支要求 t 2 / 2 < π / 2 t^2/2<\pi/2 t 2 /2 < π /2 ,即
∣ t ∣ < π |t|<\sqrt\pi ∣ t ∣ < π 。在端点处正切函数爆破,因此最大存在区间是
( − π , π ) (-\sqrt\pi,\sqrt\pi) ( − π , π ) 。求导得到
y ′ = t sec 2 ( t 2 / 2 ) = t ( 1 + y 2 ) y'=t\sec^2(t^2/2)=t(1+y^2) y ′ = t sec 2 ( t 2 /2 ) = t ( 1 + y 2 ) 。
练习 2:先选择系数连续区间 标记完成
所属知识 一阶线性方程
难度 3/5 在 t > 0 t>0 t > 0 上求初值问题
t y ′ + 2 y = t 2 , y ( 1 ) = 0. ty'+2y=t^2,\qquad y(1)=0. t y ′ + 2 y = t 2 , y ( 1 ) = 0. 查看提示 在 t>0 上除以 t,积分因子是一项幂函数。
查看解答 除以 t t t 得
y ′ + 2 t y = t . y'+\frac2t y=t. y ′ + t 2 y = t . 在 t > 0 t>0 t > 0 上取积分因子 μ = t 2 \mu=t^2 μ = t 2 ,于是
( t 2 y ) ′ = t 3 . (t^2y)'=t^3. ( t 2 y ) ′ = t 3 . 积分得
y = t 2 4 + C t 2 . y=\frac{t^2}{4}+\frac{C}{t^2}. y = 4 t 2 + t 2 C . 初值给 C = − 1 / 4 C=-1/4 C = − 1/4 ,所以
y ( t ) = t 2 4 − 1 4 t 2 . y(t)=\frac{t^2}{4}-\frac{1}{4t^2}. y ( t ) = 4 t 2 − 4 t 2 1 . 原方程在 t = 0 t=0 t = 0 退化,且解含 t − 2 t^{-2} t − 2 。包含初始点一的最大系数连续区间为
( 0 , ∞ ) (0,\infty) ( 0 , ∞ ) ;直接代入可验证等式。
练习 3:恢复势函数 标记完成
所属知识 恰当方程
难度 3/5 求经过 ( 1 , 0 ) (1,0) ( 1 , 0 ) 的隐式解:
( 2 t y + 3 t 2 ) d t + ( t 2 + 4 y ) d y = 0. (2ty+3t^2)\,\mathrm dt+(t^2+4y)\,\mathrm dy=0. ( 2 t y + 3 t 2 ) d t + ( t 2 + 4 y ) d y = 0. 查看提示 先比较
M y M_y M y 与
N t N_t N t ,再对 M 关于 t 积分。
查看解答 令 M = 2 t y + 3 t 2 M=2ty+3t^2 M = 2 t y + 3 t 2 、N = t 2 + 4 y N=t^2+4y N = t 2 + 4 y ,则
M y = 2 t = N t M_y=2t=N_t M y = 2 t = N t ,方程恰当。对 M M M 关于 t t t 积分:
Φ = t 2 y + t 3 + g ( y ) . \Phi=t^2y+t^3+g(y). Φ = t 2 y + t 3 + g ( y ) . 由 Φ y = t 2 + g ′ ( y ) = N = t 2 + 4 y \Phi_y=t^2+g'(y)=N=t^2+4y Φ y = t 2 + g ′ ( y ) = N = t 2 + 4 y 得
g ( y ) = 2 y 2 g(y)=2y^2 g ( y ) = 2 y 2 。因此解族是
t 2 y + t 3 + 2 y 2 = C . t^2y+t^3+2y^2=C. t 2 y + t 3 + 2 y 2 = C . 代入 ( 1 , 0 ) (1,0) ( 1 , 0 ) 得 C = 1 C=1 C = 1 。在该点
Φ y = 1 ≠ 0 \Phi_y=1\ne0 Φ y = 1 = 0 ,所以它在初始点附近确定唯一的函数分支 y ( t ) y(t) y ( t ) 。对隐式式求导可恢复
M + N y ′ = 0 M+Ny'=0 M + N y ′ = 0 。
练习 4:幂代换后的爆破时刻 标记完成
所属知识 Bernoulli 方程
难度 4/5 求解
y ′ − y = e t y 2 , y ( 0 ) = 1 , y'-y=e^t y^2,\qquad y(0)=1, y ′ − y = e t y 2 , y ( 0 ) = 1 , 并确定包含零的最大存在区间。
查看提示 令
v = 1 / y v=1/y v = 1/ y ,先求线性方程,再检查分母何时为零。
查看解答 初值非零,令 v = y − 1 v=y^{-1} v = y − 1 。原方程除以 y 2 y^2 y 2 后为
y ′ / y 2 − 1 / y = e t y'/y^2-1/y=e^t y ′ / y 2 − 1/ y = e t ,所以
v ′ + v = − e t . v'+v=-e^t. v ′ + v = − e t . 乘积分因子 e t e^t e t 得
( e t v ) ′ = − e 2 t . (e^t v)'=-e^{2t}. ( e t v ) ′ = − e 2 t . 因而
v = − 1 2 e t + C e − t . v=-\frac12e^t+Ce^{-t}. v = − 2 1 e t + C e − t . v ( 0 ) = 1 v(0)=1 v ( 0 ) = 1 给出 C = 3 / 2 C=3/2 C = 3/2 ,于是
y ( t ) = 2 e t 3 − e 2 t . y(t)=\frac{2e^t}{3-e^{2t}}. y ( t ) = 3 − e 2 t 2 e t . 分母在 t = 1 2 ln 3 t=\tfrac12\ln3 t = 2 1 ln 3 首次为零,向负无穷方向没有有限零点。故最大存在区间为
( − ∞ , 1 2 ln 3 ) . \left(-\infty,\frac12\ln3\right). ( − ∞ , 2 1 ln 3 ) . 将公式求导并代回原式可得等号成立,且 y ( 0 ) = 1 y(0)=1 y ( 0 ) = 1 。
概念连接与方法路线
常微分方程
提供阶数、线性与自治等基本分类语言。
初值问题
用初始状态确定积分常数、代数分支和最大区间。
积分与累积
支撑分离变量、积分因子和势函数恢复。
积分技巧
处理部分分式、换元和指数积分等实际计算。
导数与微分
为乘积求导、链式法则与全微分形式提供依据。
课程 · 2011 MIT 18.03SC Differential Equations Arthur Mattuck, Haynes Miller
为傅里叶方法、模态分解及其微分方程应用提供连续课程背景。
打开官方来源
MIT 18.03SC Differential Equations 的一阶方程单元提供方向场、可分离模型、线性方程和积分因子的连续课程脉络,适合用同一初值问题比较几何判断与解析解法。
书籍 · 2016 Calculus Volume 2 Gilbert Strang, Edwin Herman
用于核对数列与级数的收敛条件、积分技巧、反常积分和幂级数例题。
打开官方来源
OpenStax《Calculus Volume 2》的微分方程章节包含可分离方程、一阶线性方程与典型增长衰减模型,可用于补充计算练习并核对初值代入步骤。
一阶方程的稳定流程可以压缩为五步:识别结构,声明合法区域,执行变形,使用初值选分支,最后回代并确定最大区间。后续高阶线性方程与线性系统会扩大未知状态的维数,但“结构决定方法、初值选择解、区间约束结论”的原则保持不变。