考研数学二·零基础讲义 · 第 14 讲(二)|常微分方程:一阶线性方程与应用
学完本篇,你能推导并熟练使用一阶线性方程的通解公式,会"把 x x x 当未知函数"的技巧,了解伯努利方程的换元思路,并能把积分方程、几何与冷却问题化成一阶方程求解。预计用时:3 小时。
本讲地图#
第 14 讲拆成五篇:(一)基本概念、可分离变量与齐次方程;(二)一阶线性方程与应用;(三)可降阶方程与解的结构;(四)常系数齐次方程与非齐次特解(类型 Ⅰ);(五)非齐次特解(类型 Ⅱ)、欧拉方程与应用。
一阶线性方程 :y ′ + P ( x ) y = Q ( x ) y'+P(x)y=Q(x) y ′ + P ( x ) y = Q ( x ) ,通解公式及推导——大纲要求"掌握",是一阶方程中最重要的内容。
伯努利方程(选学) :说明:根据近年的数学二考试大纲原文,伯努利方程只写在数学一的要求里,数学二大纲没有单列 ;它的解法只是"换元 + 一阶线性",这里作为选学简要介绍。
应用 :积分方程、几何问题、冷却问题——大纲要求"会用微分方程解决一些简单的应用问题"。
一阶方程类型识别表 :五种类型怎么认。
前置知识#
第 14 讲(一):基本概念、可分离变量方程、齐次方程。
分部积分(见第 09 讲)、变限积分求导 ( ∫ 0 x f ( t ) d t ) ′ = f ( x ) \left(\int_0^xf(t)\,dt\right)'=f(x) ( ∫ 0 x f ( t ) d t ) ′ = f ( x ) (见第 10 讲)。
一句话心态:一阶线性方程只有一个公式,但化标准形 和符号 是两大雷区。每次套公式前先写出 P ( x ) P(x) P ( x ) 、Q ( x ) Q(x) Q ( x ) 各是什么。
一、一阶线性微分方程(本篇重点)#
1.1 大白话#
"线性"的意思是:方程里 y y y 和 y ′ y' y ′ 都只以一次方 出现,并且没有 y y ′ yy' y y ′ 、sin y \sin y sin y 这种"y y y 被加工过"的项。一阶线性方程长这样:
y ′ + P ( x ) y = Q ( x ) . y'+P(x)\,y=Q(x). y ′ + P ( x ) y = Q ( x ) .
Q ( x ) ≡ 0 Q(x)\equiv0 Q ( x ) ≡ 0 时叫齐次 线性方程(它其实可以分离变量);Q ( x ) ≢ 0 Q(x)\not\equiv0 Q ( x ) ≡ 0 时叫非齐次 线性方程。
它有一个现成的通解公式 。考研里这个公式必须背熟,而且要知道它是怎么来的(推导过程本身就是一种常考的思路)。
1.2 通解公式#
y = e − ∫ P ( x ) d x [ ∫ Q ( x ) e ∫ P ( x ) d x d x + C ] . y=e^{-\int P(x)\,dx}\left[\int Q(x)\,e^{\int P(x)\,dx}\,dx+C\right]. y = e − ∫ P ( x ) d x [ ∫ Q ( x ) e ∫ P ( x ) d x d x + C ] .
1.3 逐行拆解:公式是怎么来的#
方法一:积分因子法(推荐理解)
方程两边同乘 e ∫ P d x e^{\int P\,dx} e ∫ P d x (叫积分因子 ):
e ∫ P d x y ′ + P e ∫ P d x y = Q e ∫ P d x . e^{\int P\,dx}\,y'+P\,e^{\int P\,dx}\,y=Q\,e^{\int P\,dx}. e ∫ P d x y ′ + P e ∫ P d x y = Q e ∫ P d x .
左边恰好是乘积 ( y ⋅ e ∫ P d x ) \left(y\cdot e^{\int P\,dx}\right) ( y ⋅ e ∫ P d x ) 的导数!(验证:( y e ∫ P ) ′ = y ′ e ∫ P + y ⋅ P e ∫ P \left(ye^{\int P}\right)'=y'e^{\int P}+y\cdot Pe^{\int P} ( y e ∫ P ) ′ = y ′ e ∫ P + y ⋅ P e ∫ P 。)所以
( y e ∫ P d x ) ′ = Q e ∫ P d x . \left(y\,e^{\int P\,dx}\right)'=Q\,e^{\int P\,dx}. ( y e ∫ P d x ) ′ = Q e ∫ P d x .
两边积分:y e ∫ P d x = ∫ Q e ∫ P d x d x + C y\,e^{\int P\,dx}=\displaystyle\int Qe^{\int P\,dx}dx+C y e ∫ P d x = ∫ Q e ∫ P d x d x + C ,再除过去,就得到通解公式。
方法二:常数变易法(了解思路)
先解齐次方程 y ′ + P y = 0 y'+Py=0 y ′ + P y = 0 :分离变量得 y = C e − ∫ P d x y=Ce^{-\int P\,dx} y = C e − ∫ P d x 。然后把常数 C C C "变成"函数 C ( x ) C(x) C ( x ) ,设 y = C ( x ) e − ∫ P d x y=C(x)e^{-\int P\,dx} y = C ( x ) e − ∫ P d x ,代入原方程,可以算出 C ′ ( x ) = Q e ∫ P d x C'(x)=Qe^{\int P\,dx} C ′ ( x ) = Q e ∫ P d x ,积分得 C ( x ) C(x) C ( x ) ,结果和公式一样。
使用公式的注意事项 :
先把方程化成标准形 y ′ + P ( x ) y = Q ( x ) y'+P(x)y=Q(x) y ′ + P ( x ) y = Q ( x ) :y ′ y' y ′ 的系数必须是 1,P P P 包括前面的符号。
公式里的 ∫ P d x \int P\,dx ∫ P d x 只取一个原函数 ,不用加 C C C ;而且出现 ln ∣ x ∣ \ln|x| ln ∣ x ∣ 时通常可以直接写 ln x \ln x ln x (考虑 x > 0 x>0 x > 0 的情形,结果形式相同)。
整个公式里只有最后那个 C C C 是任意常数。
结构理解:C e − ∫ P Ce^{-\int P} C e − ∫ P 是对应齐次方程的通解,剩下那部分是非齐次方程的一个特解——"非齐次通解 = 齐次通解 + 非齐次特解",第 14 讲(三)讲解的结构时还会遇到。
1.4 例题#
例 1 求 y ′ + y x = sin x x y'+\dfrac yx=\dfrac{\sin x}{x} y ′ + x y = x sin x 的通解。
解 :
① 标准形已经是了:P = 1 x P=\dfrac1x P = x 1 ,Q = sin x x Q=\dfrac{\sin x}{x} Q = x sin x 。
② ∫ P d x = ln x \int P\,dx=\ln x ∫ P d x = ln x ,e ∫ P = x e^{\int P}=x e ∫ P = x ,e − ∫ P = 1 x e^{-\int P}=\dfrac1x e − ∫ P = x 1 。
③ ∫ Q e ∫ P d x = ∫ sin x x ⋅ x d x = ∫ sin x d x = − cos x \displaystyle\int Qe^{\int P}dx=\int\frac{\sin x}{x}\cdot x\,dx=\int\sin x\,dx=-\cos x ∫ Q e ∫ P d x = ∫ x sin x ⋅ x d x = ∫ sin x d x = − cos x 。
④ 通解:
y = 1 x ( − cos x + C ) = C − cos x x . y=\frac1x\left(-\cos x+C\right)=\frac{C-\cos x}{x}. y = x 1 ( − cos x + C ) = x C − cos x .
例 2 求 y ′ − 2 y x + 1 = ( x + 1 ) 5 2 y'-\dfrac{2y}{x+1}=(x+1)^{\frac52} y ′ − x + 1 2 y = ( x + 1 ) 2 5 的通解。
解 :
① P = − 2 x + 1 P=-\dfrac{2}{x+1} P = − x + 1 2 (注意负号 ),Q = ( x + 1 ) 5 / 2 Q=(x+1)^{5/2} Q = ( x + 1 ) 5/2 。
② ∫ P d x = − 2 ln ( x + 1 ) \int P\,dx=-2\ln(x+1) ∫ P d x = − 2 ln ( x + 1 ) ,所以 e ∫ P = ( x + 1 ) − 2 e^{\int P}=(x+1)^{-2} e ∫ P = ( x + 1 ) − 2 ,e − ∫ P = ( x + 1 ) 2 e^{-\int P}=(x+1)^2 e − ∫ P = ( x + 1 ) 2 。
③ ∫ ( x + 1 ) 5 / 2 ( x + 1 ) − 2 d x = ∫ ( x + 1 ) 1 / 2 d x = 2 3 ( x + 1 ) 3 / 2 \displaystyle\int(x+1)^{5/2}(x+1)^{-2}dx=\int(x+1)^{1/2}dx=\frac23(x+1)^{3/2} ∫ ( x + 1 ) 5/2 ( x + 1 ) − 2 d x = ∫ ( x + 1 ) 1/2 d x = 3 2 ( x + 1 ) 3/2 。
④ 通解:
y = ( x + 1 ) 2 [ 2 3 ( x + 1 ) 3 2 + C ] . y=(x+1)^2\left[\frac23(x+1)^{\frac32}+C\right]. y = ( x + 1 ) 2 [ 3 2 ( x + 1 ) 2 3 + C ] .
例 3 求 y ′ + y = e − x y'+y=e^{-x} y ′ + y = e − x 满足 y ( 0 ) = 1 y(0)=1 y ( 0 ) = 1 的特解。
解 :P = 1 P=1 P = 1 ,Q = e − x Q=e^{-x} Q = e − x ;e ∫ P = e x e^{\int P}=e^x e ∫ P = e x 。
y = e − x [ ∫ e − x e x d x + C ] = e − x ( x + C ) . y=e^{-x}\left[\int e^{-x}e^{x}dx+C\right]=e^{-x}(x+C). y = e − x [ ∫ e − x e x d x + C ] = e − x ( x + C ) .
代入 y ( 0 ) = 1 y(0)=1 y ( 0 ) = 1 :C = 1 C=1 C = 1 。特解 y = ( x + 1 ) e − x y=(x+1)e^{-x} y = ( x + 1 ) e − x 。
例 4 求 y ′ + y tan x = cos x y'+y\tan x=\cos x y ′ + y tan x = cos x 的通解(− π 2 < x < π 2 -\frac\pi2<x<\frac\pi2 − 2 π < x < 2 π )。
解 :P = tan x P=\tan x P = tan x ,∫ tan x d x = − ln cos x \int\tan x\,dx=-\ln\cos x ∫ tan x d x = − ln cos x ,所以 e ∫ P = 1 cos x e^{\int P}=\dfrac{1}{\cos x} e ∫ P = cos x 1 ,e − ∫ P = cos x e^{-\int P}=\cos x e − ∫ P = cos x 。
∫ cos x ⋅ 1 cos x d x = x , y = cos x ( x + C ) . \int\cos x\cdot\frac{1}{\cos x}dx=x ,\qquad y=\cos x\,(x+C). ∫ cos x ⋅ cos x 1 d x = x , y = cos x ( x + C ) .
例 5 求 y ′ − y = e 2 x y'-y=e^{2x} y ′ − y = e 2 x 满足 y ( 0 ) = 3 y(0)=3 y ( 0 ) = 3 的特解。
解 :P = − 1 P=-1 P = − 1 ,e ∫ P = e − x e^{\int P}=e^{-x} e ∫ P = e − x ,e − ∫ P = e x e^{-\int P}=e^x e − ∫ P = e x 。
y = e x [ ∫ e 2 x e − x d x + C ] = e x ( e x + C ) = e 2 x + C e x . y=e^x\left[\int e^{2x}e^{-x}dx+C\right]=e^x(e^x+C)=e^{2x}+Ce^x . y = e x [ ∫ e 2 x e − x d x + C ] = e x ( e x + C ) = e 2 x + C e x .
y ( 0 ) = 1 + C = 3 y(0)=1+C=3 y ( 0 ) = 1 + C = 3 ,C = 2 C=2 C = 2 。特解 y = e 2 x + 2 e x y=e^{2x}+2e^x y = e 2 x + 2 e x 。
例 6(把 x x x 看成 y y y 的函数) 求 d y d x = y x + y 3 \dfrac{dy}{dx}=\dfrac{y}{x+y^3} d x d y = x + y 3 y 的通解(y > 0 y>0 y > 0 )。
解 :关于 y y y 看,这个方程既不能分离,也不是线性的(y y y 出现在分母的 y 3 y^3 y 3 里)。换个角度 :把它倒过来,
d x d y = x + y 3 y = x y + y 2 , 即 d x d y − 1 y x = y 2 . \frac{dx}{dy}=\frac{x+y^3}{y}=\frac xy+y^2 ,\quad\text{即}\quad \frac{dx}{dy}-\frac1y\,x=y^2 . d y d x = y x + y 3 = y x + y 2 , 即 d y d x − y 1 x = y 2 .
这是以 y y y 为自变量、x x x 为未知函数 的一阶线性方程!P = − 1 y P=-\dfrac1y P = − y 1 ,Q = y 2 Q=y^2 Q = y 2 ,e ∫ P d y = 1 y e^{\int P\,dy}=\dfrac1y e ∫ P d y = y 1 ,e − ∫ P d y = y e^{-\int P\,dy}=y e − ∫ P d y = y 。
x = y [ ∫ y 2 ⋅ 1 y d y + C ] = y ( y 2 2 + C ) . x=y\left[\int y^2\cdot\frac1y\,dy+C\right]=y\left(\frac{y^2}{2}+C\right). x = y [ ∫ y 2 ⋅ y 1 d y + C ] = y ( 2 y 2 + C ) .
当方程关于 y y y 很别扭、关于 x x x 却是一次的时候,试试"把 x x x 当未知函数 "。
1.5 易错点 ⚠️#
不化标准形就套公式 是头号错误。比如 x y ′ + y = sin x xy'+y=\sin x x y ′ + y = sin x 要先除以 x x x 。
P ( x ) P(x) P ( x ) 带符号:y ′ − 2 y = … y'-2y=\ldots y ′ − 2 y = … 中 P = − 2 P=-2 P = − 2 。
e ∫ P e^{\int P} e ∫ P 和 e − ∫ P e^{-\int P} e − ∫ P 别放反:括号里面 乘的是 e + ∫ P e^{+\int P} e + ∫ P ,括号外面 是 e − ∫ P e^{-\int P} e − ∫ P 。
C C C 在括号里面,乘开后是 C e − ∫ P Ce^{-\int P} C e − ∫ P ,不是简单加在最后。
二、伯努利方程(选学)#
2.1 大白话#
伯努利方程比一阶线性方程多了一个 y n y^n y n :y ′ + P ( x ) y = Q ( x ) y n y'+P(x)y=Q(x)y^n y ′ + P ( x ) y = Q ( x ) y n 。它不是线性的,但用一个巧妙的换元就能变成线性方程。再次说明:数学二大纲没有单列伯努利方程 (数学一大纲有),这里只是作为"变量代换"思路的练习,时间紧可以跳过。
2.2 方法#
y ′ + P ( x ) y = Q ( x ) y n ( n ≠ 0 , 1 ) . y'+P(x)\,y=Q(x)\,y^n\qquad(n\ne0,1). y ′ + P ( x ) y = Q ( x ) y n ( n = 0 , 1 ) .
两边除以 y n y^n y n :y − n y ′ + P y 1 − n = Q y^{-n}y'+P\,y^{1-n}=Q y − n y ′ + P y 1 − n = Q ;
令 z = y 1 − n z=y^{1-n} z = y 1 − n ,则 z ′ = ( 1 − n ) y − n y ′ z'=(1-n)y^{-n}y' z ′ = ( 1 − n ) y − n y ′ ;
代入得一阶线性方程:
z ′ + ( 1 − n ) P ( x ) z = ( 1 − n ) Q ( x ) ; z'+(1-n)P(x)\,z=(1-n)Q(x); z ′ + ( 1 − n ) P ( x ) z = ( 1 − n ) Q ( x ) ;
用通解公式求 z z z ,再换回 y = z 1 1 − n y=z^{\frac{1}{1-n}} y = z 1 − n 1 。
2.3 逐行拆解#
n = 0 n=0 n = 0 时就是线性方程,n = 1 n=1 n = 1 时可以分离变量,所以排除这两种。
最常见的是 n = 2 n=2 n = 2 :令 z = 1 y z=\dfrac1y z = y 1 。
y = 0 y=0 y = 0 往往也是一个解(除以 y n y^n y n 时被丢掉了)。
2.4 例题#
例 7 求 y ′ + y = x y 2 y'+y=xy^2 y ′ + y = x y 2 的通解。
解 :n = 2 n=2 n = 2 ,令 z = y − 1 z=y^{-1} z = y − 1 ,则 z ′ = − y − 2 y ′ z'=-y^{-2}y' z ′ = − y − 2 y ′ 。原方程两边除以 y 2 y^2 y 2 :y − 2 y ′ + y − 1 = x y^{-2}y'+y^{-1}=x y − 2 y ′ + y − 1 = x ,即 − z ′ + z = x -z'+z=x − z ′ + z = x ,整理成标准形:
z ′ − z = − x . z'-z=-x . z ′ − z = − x .
P = − 1 P=-1 P = − 1 ,e ∫ P = e − x e^{\int P}=e^{-x} e ∫ P = e − x :
z = e x [ ∫ ( − x ) e − x d x + C ] . z=e^{x}\left[\int(-x)e^{-x}dx+C\right]. z = e x [ ∫ ( − x ) e − x d x + C ] .
分部积分:∫ ( − x ) e − x d x = x e − x + e − x \displaystyle\int(-x)e^{-x}dx=xe^{-x}+e^{-x} ∫ ( − x ) e − x d x = x e − x + e − x (验证:对它求导得 e − x − x e − x − e − x = − x e − x e^{-x}-xe^{-x}-e^{-x}=-xe^{-x} e − x − x e − x − e − x = − x e − x ✓)。所以
z = x + 1 + C e x , y = 1 x + 1 + C e x . z=x+1+Ce^x,\qquad y=\frac{1}{x+1+Ce^x}. z = x + 1 + C e x , y = x + 1 + C e x 1 .
例 8 求 y ′ − y x = − y 2 y'-\dfrac yx=-y^2 y ′ − x y = − y 2 的通解。
解 :n = 2 n=2 n = 2 ,令 z = 1 y z=\dfrac1y z = y 1 。除以 y 2 y^2 y 2 得 y − 2 y ′ − 1 x y − 1 = − 1 y^{-2}y'-\dfrac1x y^{-1}=-1 y − 2 y ′ − x 1 y − 1 = − 1 ,即 − z ′ − z x = − 1 -z'-\dfrac zx=-1 − z ′ − x z = − 1 ,标准形 z ′ + z x = 1 z'+\dfrac zx=1 z ′ + x z = 1 。
e ∫ 1 x d x = x e^{\int\frac1x dx}=x e ∫ x 1 d x = x :z = 1 x [ ∫ x d x + C ] = x 2 + C x z=\dfrac1x\left[\displaystyle\int x\,dx+C\right]=\dfrac x2+\dfrac Cx z = x 1 [ ∫ x d x + C ] = 2 x + x C 。所以
y = 1 x 2 + C x = 2 x x 2 + 2 C . y=\frac{1}{\frac x2+\frac Cx}=\frac{2x}{x^2+2C}. y = 2 x + x C 1 = x 2 + 2 C 2 x .
2.5 易错点 ⚠️#
换元后得到的是关于 z z z 的方程,系数要乘 ( 1 − n ) (1-n) ( 1 − n ) ,符号别丢。
最后要换回 y y y 。
三、一阶方程的简单应用#
3.1 大白话#
应用题的核心是把文字翻译成方程 :
"变化率""速度""斜率" → 导数;
"与……成正比" → 乘一个比例系数 k k k ;
题目里含 ∫ 0 x f ( t ) d t \int_0^xf(t)\,dt ∫ 0 x f ( t ) d t 的等式(积分方程)→ 两边求导 变成微分方程,同时代入特殊点 得到初始条件。
3.2 例题#
例 9(积分方程) 设 f ( x ) f(x) f ( x ) 连续,且 f ( x ) = ∫ 0 x f ( t ) d t + x f(x)=\displaystyle\int_0^xf(t)\,dt+x f ( x ) = ∫ 0 x f ( t ) d t + x ,求 f ( x ) f(x) f ( x ) 。
解 :
① 求导:两边对 x x x 求导(变限积分求导):f ′ ( x ) = f ( x ) + 1 f'(x)=f(x)+1 f ′ ( x ) = f ( x ) + 1 。
② 初始条件:令 x = 0 x=0 x = 0 ,积分为 0,得 f ( 0 ) = 0 f(0)=0 f ( 0 ) = 0 。(这一步最容易漏 。)
③ 解方程:f ′ − f = 1 f'-f=1 f ′ − f = 1 ,一阶线性,P = − 1 P=-1 P = − 1 :f = e x [ ∫ e − x d x + C ] = e x ( − e − x + C ) = C e x − 1 f=e^{x}\left[\displaystyle\int e^{-x}dx+C\right]=e^x(-e^{-x}+C)=Ce^x-1 f = e x [ ∫ e − x d x + C ] = e x ( − e − x + C ) = C e x − 1 。
④ 代 f ( 0 ) = 0 f(0)=0 f ( 0 ) = 0 :C = 1 C=1 C = 1 。所以 f ( x ) = e x − 1 f(x)=e^x-1 f ( x ) = e x − 1 。
验算:∫ 0 x ( e t − 1 ) d t + x = ( e x − 1 − x ) + x = e x − 1 \int_0^x(e^t-1)dt+x=(e^x-1-x)+x=e^x-1 ∫ 0 x ( e t − 1 ) d t + x = ( e x − 1 − x ) + x = e x − 1 ✓。
例 10(积分方程) 设 f ( x ) f(x) f ( x ) 连续,且 f ( x ) = x 2 + ∫ 0 x f ( t ) d t f(x)=x^2+\displaystyle\int_0^xf(t)\,dt f ( x ) = x 2 + ∫ 0 x f ( t ) d t ,求 f ( x ) f(x) f ( x ) 。
解 :求导得 f ′ = 2 x + f f'=2x+f f ′ = 2 x + f ;令 x = 0 x=0 x = 0 得 f ( 0 ) = 0 f(0)=0 f ( 0 ) = 0 。
f ′ − f = 2 x f'-f=2x f ′ − f = 2 x :f = e x [ ∫ 2 x e − x d x + C ] f=e^x\left[\displaystyle\int2xe^{-x}dx+C\right] f = e x [ ∫ 2 x e − x d x + C ] 。分部积分 ∫ 2 x e − x d x = − 2 x e − x − 2 e − x \int2xe^{-x}dx=-2xe^{-x}-2e^{-x} ∫ 2 x e − x d x = − 2 x e − x − 2 e − x ,所以
f = e x ( − 2 x e − x − 2 e − x + C ) = C e x − 2 x − 2. f=e^x\left(-2xe^{-x}-2e^{-x}+C\right)=Ce^x-2x-2 . f = e x ( − 2 x e − x − 2 e − x + C ) = C e x − 2 x − 2.
f ( 0 ) = C − 2 = 0 f(0)=C-2=0 f ( 0 ) = C − 2 = 0 ,C = 2 C=2 C = 2 。f ( x ) = 2 e x − 2 x − 2 f(x)=2e^x-2x-2 f ( x ) = 2 e x − 2 x − 2 。
例 11(几何问题) 一条曲线过点 ( 1 , 2 ) (1,2) ( 1 , 2 ) ,且曲线上任一点 ( x , y ) (x,y) ( x , y ) 处切线的斜率等于 2 y x \dfrac{2y}{x} x 2 y ,求这条曲线。
解 :切线斜率就是 y ′ y' y ′ ,所以 y ′ = 2 y x y'=\dfrac{2y}{x} y ′ = x 2 y ,y ( 1 ) = 2 y(1)=2 y ( 1 ) = 2 。
分离:d y y = 2 d x x \dfrac{dy}{y}=\dfrac{2dx}{x} y d y = x 2 d x ,ln ∣ y ∣ = 2 ln ∣ x ∣ + C 0 \ln|y|=2\ln|x|+C_0 ln ∣ y ∣ = 2 ln ∣ x ∣ + C 0 ,y = C x 2 y=Cx^2 y = C x 2 。由 y ( 1 ) = 2 y(1)=2 y ( 1 ) = 2 得 C = 2 C=2 C = 2 。曲线 y = 2 x 2 y=2x^2 y = 2 x 2 。
例 12(冷却问题) 牛顿冷却定律:物体温度的变化率与"物体温度和环境温度之差"成正比。一杯 100 ∘ C 100^\circ\mathrm{C} 10 0 ∘ C 的水放在 20 ∘ C 20^\circ\mathrm{C} 2 0 ∘ C 的房间里,10 分钟后降到 60 ∘ C 60^\circ\mathrm{C} 6 0 ∘ C 。问再过多久降到 30 ∘ C 30^\circ\mathrm{C} 3 0 ∘ C ?
解 :
① 建方程 :设 t t t 分钟时温度为 T ( t ) T(t) T ( t ) 。"变化率"是 d T d t \dfrac{dT}{dt} d t d T ;温度在下降,所以
d T d t = − k ( T − 20 ) , k > 0 , T ( 0 ) = 100. \frac{dT}{dt}=-k(T-20),\quad k>0,\quad T(0)=100 . d t d T = − k ( T − 20 ) , k > 0 , T ( 0 ) = 100.
② 解 :分离 d T T − 20 = − k d t \dfrac{dT}{T-20}=-k\,dt T − 20 d T = − k d t ,积分 ln ( T − 20 ) = − k t + C 0 \ln(T-20)=-kt+C_0 ln ( T − 20 ) = − k t + C 0 ,T = 20 + C e − k t T=20+Ce^{-kt} T = 20 + C e − k t 。由 T ( 0 ) = 100 T(0)=100 T ( 0 ) = 100 得 C = 80 C=80 C = 80 :T = 20 + 80 e − k t T=20+80e^{-kt} T = 20 + 80 e − k t 。
③ 定 k k k :T ( 10 ) = 60 T(10)=60 T ( 10 ) = 60 :80 e − 10 k = 40 80e^{-10k}=40 80 e − 10 k = 40 ,e − 10 k = 1 2 e^{-10k}=\frac12 e − 10 k = 2 1 ,k = ln 2 10 k=\dfrac{\ln2}{10} k = 10 ln 2 。于是 T = 20 + 80 ⋅ 2 − t / 10 T=20+80\cdot2^{-t/10} T = 20 + 80 ⋅ 2 − t /10 。
④ 求时间 :T = 30 T=30 T = 30 :80 ⋅ 2 − t / 10 = 10 80\cdot2^{-t/10}=10 80 ⋅ 2 − t /10 = 10 ,2 − t / 10 = 1 8 = 2 − 3 2^{-t/10}=\frac18=2^{-3} 2 − t /10 = 8 1 = 2 − 3 ,t = 30 t=30 t = 30 。
从开始算是 30 分钟,所以再过 20 分钟 降到 30 ∘ C 30^\circ\mathrm{C} 3 0 ∘ C 。
冷却问题:温差按指数衰减,30 分钟时降到 30 ℃
t
T
10
30
60
20
60
100
30
O
(0, 100)
(10, 60)
(30, 30)
室温 20 ℃(T 无限接近但达不到)
T = 20 + 80·2^(−t/10)
每 10 分钟温差减半:80 → 40 → 20 → 10
图 1:例 12 冷却曲线 T = 20 + 80·2^(−t/10):每 10 分钟温差减半,t = 30 时 T = 30 ℃
3.3 易错点 ⚠️#
积分方程题:求导 与代特殊点找初始条件 两步缺一不可。
"再过多久"要减去已经过去的时间。
方程里的符号要和实际一致:温度下降,导数为负。
四、一阶方程类型识别表#
拿到一个一阶方程,按这个顺序检查:
顺序
检查
是的话
解法
1
能写成 y ′ = f ( x ) g ( y ) y'=f(x)g(y) y ′ = f ( x ) g ( y ) ?
可分离变量
分离后两边积分
2
y y y 、y ′ y' y ′ 都是一次,能写成 y ′ + P ( x ) y = Q ( x ) y'+P(x)y=Q(x) y ′ + P ( x ) y = Q ( x ) ?
一阶线性
通解公式
3
右边能写成 φ ( y x ) \varphi\left(\frac yx\right) φ ( x y ) ?
齐次方程
令 u = y x u=\frac yx u = x y
4
形如 y ′ + P y = Q y n y'+Py=Qy^n y ′ + P y = Q y n ?
伯努利(选学)
令 z = y 1 − n z=y^{1-n} z = y 1 − n
5
关于 y y y 都不行,但关于 x x x 是一次的?
以 x x x 为未知函数的线性方程
倒过来写 d x d y \frac{dx}{dy} d y d x
一个方程可能同时属于几种类型(比如 y ′ = 2 x y y'=2xy y ′ = 2 x y 既可分离又是线性),选你觉得最快的那种。
公式/结论速记卡#
类型
标准形
解法要点
一阶线性
y ′ + P ( x ) y = Q ( x ) y'+P(x)y=Q(x) y ′ + P ( x ) y = Q ( x )
y = e − ∫ P d x [ ∫ Q e ∫ P d x d x + C ] y=e^{-\int Pdx}\left[\int Qe^{\int Pdx}dx+C\right] y = e − ∫ P d x [ ∫ Q e ∫ P d x d x + C ]
积分因子
—
两边乘 e ∫ P e^{\int P} e ∫ P ,左边变成 ( y e ∫ P ) ′ \left(ye^{\int P}\right)' ( y e ∫ P ) ′
换未知函数
关于 x x x 一次
写成 d x d y + P ( y ) x = Q ( y ) \frac{dx}{dy}+P(y)x=Q(y) d y d x + P ( y ) x = Q ( y )
伯努利(选学)
y ′ + P y = Q y n y'+Py=Qy^n y ′ + P y = Q y n
z = y 1 − n z=y^{1-n} z = y 1 − n ,z ′ + ( 1 − n ) P z = ( 1 − n ) Q z'+(1-n)Pz=(1-n)Q z ′ + ( 1 − n ) P z = ( 1 − n ) Q
积分方程
含 ∫ 0 x f ( t ) d t \int_0^xf(t)dt ∫ 0 x f ( t ) d t
求导得方程 + 代 x = 0 x=0 x = 0 得初始条件
自测练习#
1. 求 x y ′ = y + x xy'=y+x x y ′ = y + x 的通解。
答案与提示
化标准形:y ′ − 1 x y = 1 y'-\dfrac1xy=1 y ′ − x 1 y = 1 。P = − 1 x P=-\frac1x P = − x 1 ,e ∫ P = 1 x e^{\int P}=\frac1x e ∫ P = x 1 ,e − ∫ P = x e^{-\int P}=x e − ∫ P = x 。
y = x [ ∫ 1 ⋅ 1 x d x + C ] = x ( ln ∣ x ∣ + C ) y=x\left[\displaystyle\int1\cdot\frac1xdx+C\right]=x(\ln|x|+C) y = x [ ∫ 1 ⋅ x 1 d x + C ] = x ( ln ∣ x ∣ + C ) 。
(也可看成齐次方程 y ′ = y x + 1 y'=\frac yx+1 y ′ = x y + 1 来解,结果相同。)
2. 求 y ′ + 2 x y = 2 x e − x 2 y'+2xy=2xe^{-x^2} y ′ + 2 x y = 2 x e − x 2 的通解。
答案与提示
P = 2 x P=2x P = 2 x ,e ∫ P = e x 2 e^{\int P}=e^{x^2} e ∫ P = e x 2 。y = e − x 2 [ ∫ 2 x e − x 2 e x 2 d x + C ] = e − x 2 ( x 2 + C ) y=e^{-x^2}\left[\displaystyle\int2xe^{-x^2}e^{x^2}dx+C\right]=e^{-x^2}(x^2+C) y = e − x 2 [ ∫ 2 x e − x 2 e x 2 d x + C ] = e − x 2 ( x 2 + C ) 。
3. 求 y ′ − 2 y = e x y'-2y=e^x y ′ − 2 y = e x 满足 y ( 0 ) = 0 y(0)=0 y ( 0 ) = 0 的特解。
答案与提示
P = − 2 P=-2 P = − 2 ,e ∫ P = e − 2 x e^{\int P}=e^{-2x} e ∫ P = e − 2 x 。y = e 2 x [ ∫ e x e − 2 x d x + C ] = e 2 x ( − e − x + C ) = C e 2 x − e x y=e^{2x}\left[\displaystyle\int e^xe^{-2x}dx+C\right]=e^{2x}(-e^{-x}+C)=Ce^{2x}-e^x y = e 2 x [ ∫ e x e − 2 x d x + C ] = e 2 x ( − e − x + C ) = C e 2 x − e x 。
y ( 0 ) = C − 1 = 0 y(0)=C-1=0 y ( 0 ) = C − 1 = 0 ,C = 1 C=1 C = 1 。特解 y = e 2 x − e x y=e^{2x}-e^x y = e 2 x − e x 。
4. 求 y ′ − y x = x 2 y'-\dfrac yx=x^2 y ′ − x y = x 2 的通解(x > 0 x>0 x > 0 )。
答案与提示
P = − 1 x P=-\frac1x P = − x 1 ,e ∫ P = 1 x e^{\int P}=\frac1x e ∫ P = x 1 ,e − ∫ P = x e^{-\int P}=x e − ∫ P = x 。y = x [ ∫ x 2 ⋅ 1 x d x + C ] = x ( x 2 2 + C ) y=x\left[\displaystyle\int x^2\cdot\frac1x\,dx+C\right]=x\left(\frac{x^2}{2}+C\right) y = x [ ∫ x 2 ⋅ x 1 d x + C ] = x ( 2 x 2 + C ) 。
5. (选学)求伯努利方程 y ′ + y x = y 2 y'+\dfrac yx=y^2 y ′ + x y = y 2 的通解(x > 0 x>0 x > 0 )。
答案与提示
令 z = 1 y z=\frac1y z = y 1 :除以 y 2 y^2 y 2 得 y − 2 y ′ + 1 x y − 1 = 1 y^{-2}y'+\frac1xy^{-1}=1 y − 2 y ′ + x 1 y − 1 = 1 ,即 − z ′ + z x = 1 -z'+\frac zx=1 − z ′ + x z = 1 ,标准形 z ′ − 1 x z = − 1 z'-\frac1xz=-1 z ′ − x 1 z = − 1 。
e ∫ P = 1 x e^{\int P}=\frac1x e ∫ P = x 1 :z = x [ ∫ ( − 1 x ) d x + C ] = x ( C − ln x ) z=x\left[\displaystyle\int\left(-\frac1x\right)dx+C\right]=x(C-\ln x) z = x [ ∫ ( − x 1 ) d x + C ] = x ( C − ln x ) 。
所以 y = 1 x ( C − ln x ) y=\dfrac{1}{x(C-\ln x)} y = x ( C − ln x ) 1 。
6. 设 f ( x ) f(x) f ( x ) 连续,且 f ( x ) = e x + ∫ 0 x f ( t ) d t f(x)=e^x+\displaystyle\int_0^xf(t)\,dt f ( x ) = e x + ∫ 0 x f ( t ) d t ,求 f ( x ) f(x) f ( x ) 。
答案与提示
求导:f ′ = e x + f f'=e^x+f f ′ = e x + f ;令 x = 0 x=0 x = 0 :f ( 0 ) = 1 f(0)=1 f ( 0 ) = 1 。
f ′ − f = e x f'-f=e^x f ′ − f = e x :f = e x [ ∫ e x e − x d x + C ] = e x ( x + C ) f=e^x\left[\displaystyle\int e^xe^{-x}dx+C\right]=e^x(x+C) f = e x [ ∫ e x e − x d x + C ] = e x ( x + C ) 。f ( 0 ) = C = 1 f(0)=C=1 f ( 0 ) = C = 1 。所以 f ( x ) = ( x + 1 ) e x f(x)=(x+1)e^x f ( x ) = ( x + 1 ) e x 。
下一讲预告#
第 14 讲(三)进入高阶方程:三类可降阶方程,以及线性微分方程解的结构——为后面的二阶常系数方程打基础。
评论
还没有评论。 来写第一条。