学径XUEJING · 个人学习平台通知设置

课程库 Calculus 2 S5 · 一阶常微分方程 › 第 8

混合问题 —— 水箱里的一阶线性方程

普通课45 分钟

本课导读

单元正课的收官战。混合问题是一阶 ODE 最经典的应用题型:一只水箱、盐水进进出出,问某时刻箱里有多少盐。它的珍贵之处不在水箱本身,而在同构——药物在血液里的代谢、房间空气的净化、湖泊污染的治理,全是这道题换了戏服。学会一次建模,四处免费复用。技术上它是积分因子的主场,还会客串一次"水量会变"的变系数好戏。预计 45 分钟。

说明

前置工具箱:积分因子四步流程(第 3 课,本课主力);平衡解与稳定性(第 5-6 课,用来"读"答案);浓度常识:浓度=溶质量体积\text{浓度} = \dfrac{\text{溶质量}}{\text{体积}}

🔧 30 秒回顾 · 本课要用的旧手艺(熟练可跳过)

① 积分因子四步(第 3 课):标准化 → μ=epdt\mu = e^{\int p\,dt} → 左边折叠成 (μA)(\mu A)' → 积分。微例:p=2100+tμ=e2ln(100+t)=(100+t)2p = \dfrac{2}{100+t} \Rightarrow \mu = e^{2\ln(100+t)} = (100+t)^2——例题 2 原式。⚠️ 翻车点:eklnu=uke^{k\ln u} = u^k 的化简。

② 指数衰减模板A=AτA=A0et/τA' = -\dfrac{A}{\tau} \Rightarrow A = A_0 e^{-t/\tau},与终点的剩余差距每 τln2\tau\ln 2 减半。微例:练习 2 的通风问题整题就是它。⚠️ 翻车点:τ\tau 在分母——时间常数越大衰减越慢。

③ 定性读答案(第 5-6 课):f(A)=0f(A^*) = 0 找平衡,f(A)<0f'(A^*) < 0 则稳定。例题 1 第 2 步用它不解方程先"预习"出 50 kg。

建模总纲:一条守恒式打天下

A(t)A(t)tt 时刻箱中溶质(盐)的总量。核心只有一条收支守恒

dAdt=cinrin流入速率    AVrout流出速率\frac{dA}{dt} = \underbrace{c_{\text{in}} \cdot r_{\text{in}}}_{\text{流入速率}} \;-\; \underbrace{\frac{A}{V} \cdot r_{\text{out}}}_{\text{流出速率}}
  • 流入:进水浓度 cinc_{\text{in}}(kg/L)× 进水流量 rinr_{\text{in}}(L/min)——两个已知数相乘,简单
  • 流出:出水浓度 × 出水流量。关键假设:搅拌均匀(well-mixed),出水浓度 = 箱内当前浓度 AV\dfrac{A}{V}——未知数 AA 由此进入方程,题目的全部难度都在这一项
  • 体积V(t)=V0+(rinrout)tV(t) = V_0 + (r_{\text{in}} - r_{\text{out}})\,t。进出流量相等时 VV 是常数(方程自治,好算);不等时 VV 随时间变(变系数线性方程,例题 2 的戏码)

提示

类比记忆:混合问题是记账——存量的变化率 = 进账速率 − 出账速率。进账是"外面给多少"(已知),出账是"按现有家底的比例往外流"(含未知数 AA)。所有同构问题(药物、污染、通风)都是这本账换个户名。

例题 1 · 标准水箱(水量不变)

100 L 水箱盛满清水。浓度 0.5 kg/L 的盐水以 2 L/min 流入,充分搅拌后混合液以 2 L/min 流出。求 A(t)A(t) 并回答:长期来看箱里会有多少盐?

第 1 步 · 建模

进出流量相等,V100V \equiv 100 L。记账:

dAdt=0.5×2A100×2=1A50,A(0)=0\frac{dA}{dt} = 0.5 \times 2 - \frac{A}{100} \times 2 = 1 - \frac{A}{50}, \qquad A(0) = 0
第 2 步 · 先定性“预习”答案

自治方程!平衡点:1A50=0A=501 - \dfrac{A^*}{50} = 0 \Rightarrow A^* = 50 kg;f(A)=150<0f'(A) = -\dfrac{1}{50} < 0稳定。不解方程已知结局:盐量将单调爬向 50 kg。合理性核对:50 kg / 100 L = 0.5 kg/L——箱内浓度终将追平进水浓度,直觉完全认账(一直灌 0.5 的水,最后整箱都是 0.5)。

第 3 步 · 积分因子求精确解

标准化:A+150A=1A' + \dfrac{1}{50}A = 1μ=et/50\mu = e^{t/50}(et/50A)=et/50(e^{t/50}A)' = e^{t/50},积分得 et/50A=50et/50+Ce^{t/50}A = 50e^{t/50} + C,通解 A=50+Cet/50A = 50 + Ce^{-t/50}。初值 A(0)=0A(0)=0C=50C = -50

A(t)=50(1et/50)A(t) = 50\left(1 - e^{-t/50}\right)

验证:A=et/50A' = e^{-t/50},而 1A50=et/501 - \dfrac{A}{50} = e^{-t/50} ✓;A(0)=0A(0) = 0 ✓;tt \to \inftyA50A \to 50,与定性预习互相盖章。指数逼近的节奏和咖啡如出一辙——它本来就是咖啡方程的亲戚。

例题 2 · 水量会变的水箱(变系数登场)

100 L 水箱含 10 kg 盐。浓度 1 kg/L 的盐水以 3 L/min 流入,混合液以 2 L/min 流出。求 A(t)A(t)

第 1 步 · 建模(体积在涨)

净流入 32=13 - 2 = 1 L/min,故 V(t)=100+tV(t) = 100 + t。记账:

dAdt=1×3A100+t×2A+2100+tA=3,A(0)=10\frac{dA}{dt} = 1 \times 3 - \frac{A}{100 + t} \times 2 \quad\Rightarrow\quad A' + \frac{2}{100+t}A = 3, \qquad A(0) = 10

变系数一阶线性——不再自治(规则看钟表:同样的 AA,水多时流失浓度低),定性捷径失效,正面强攻。

第 2 步 · 积分因子(变系数照样四步)μ=e2100+tdt=e2ln(100+t)=(100+t)2\mu = e^{\int \frac{2}{100+t}dt} = e^{2\ln(100+t)} = (100+t)^2

折叠积分:((100+t)2A)=3(100+t)2\big((100+t)^2 A\big)' = 3(100+t)^2,两边积分得 (100+t)2A=(100+t)3+C(100+t)^2 A = (100+t)^3 + C

A=(100+t)+C(100+t)2A = (100+t) + \frac{C}{(100+t)^2}
第 3 步 · 定 C,读答案

A(0)=10A(0) = 10100+C104=10C=9×105100 + \dfrac{C}{10^4} = 10 \Rightarrow C = -9 \times 10^5

A(t)=(100+t)900000(100+t)2A(t) = (100 + t) - \frac{900000}{(100+t)^2}

核对 A(0)=10090=10A(0) = 100 - 90 = 10 ✓。读答案:tt 大时修正项消失,A100+t=V(t)A \approx 100 + t = V(t),即浓度 AV1\dfrac{A}{V} \to 1 kg/L——又追平进水浓度。两道例题殊途同归的这条规律,值得放进直觉库:只进一种水,箱子终成那种水。(现实提醒:若水箱有容量上限,解只在灌满前有效——第 2 课"解的有效期"再现。)

✍️ 练习

练习 1 · 例题 1 的追问

例题 1 中,盐量达到极限值的一半(25 kg)需要多久?

✅ 参考答案

25=50(1et/50)et/50=12t=50ln234.725 = 50(1 - e^{-t/50}) \Rightarrow e^{-t/50} = \dfrac{1}{2} \Rightarrow t = 50\ln 2 \approx 34.7 分钟。指数逼近的通用节奏:每过 50ln250\ln 2 分钟,与终点的剩余差距减半——和放射性半衰期是同一个数学。

练习 2 · 换戏服:会议室通风

60 m³ 的会议室散会后残留 300 g CO₂ 超标部分。新风系统以 3 m³/min 送入洁净空气(不含超标气体),混匀空气以 3 m³/min 排出。求超标量 A(t)A(t) 及降到 30 g(10%)所需时间。

✅ 参考答案

记账:流入为 0(洁净),流出 A60×3=A20\dfrac{A}{60}\times 3 = \dfrac{A}{20}

dAdt=A20A=300et/20\frac{dA}{dt} = -\frac{A}{20} \quad\Rightarrow\quad A = 300\,e^{-t/20}

(可分离送分题——混合问题不总需要积分因子。)令 A=30A = 30et/20=0.1e^{-t/20} = 0.1t=20ln1046t = 20\ln 10 \approx 46 分钟。同构清单打卡:这就是"进水浓度为零"的水箱,也是静脉停药后的血药浓度衰减。

练习 3 · 只建模,不求解

500 L 水箱含 40 kg 盐。浓度 0.2 kg/L 盐水以 5 L/min 流入,混合液以 5 L/min 流出。只写出初值问题,并用定性分析预言长期盐量(不许解方程)。

✅ 参考答案dAdt=1A100,A(0)=40\frac{dA}{dt} = 1 - \frac{A}{100}, \qquad A(0) = 40

平衡点 A=100A^* = 100 kg,f=1100<0f' = -\dfrac{1}{100} < 0 稳定:盐量将从 40 kg 单调爬向 100 kg(= 0.2×5000.2 \times 500,浓度追平进水,直觉盖章)。建模 + 定性 = 一滴积分不用,报告已经能交——这正是本单元想给你的完整工作流。

📝 随堂测验

随堂测验0 / 5 题正确

1. 混合问题中,流出速率写成 $\frac{A}{V} \cdot r_{out}$ 依赖的关键假设是?

2. 进出流量相等的标准水箱(例题 1),长期行为是?

3. 进出流量不等时(例题 2),方程发生的变化是?

4. 练习 2 的会议室通风问题与水箱问题的关系是?

5. 解混合问题的推荐工作流是?

本课小结

  • 建模一条式:dAdt=cinrinAVrout\dfrac{dA}{dt} = c_{\text{in}}r_{\text{in}} - \dfrac{A}{V}r_{\text{out}}——进账已知、出账按家底比例;地基假设是充分搅拌
  • 流量相等:VV 恒定、方程自治,定性能"预习"答案;流量不等:V(t)V(t) 线性变化、变系数线性方程,积分因子照旧
  • 通用直觉:只进一种水,箱子终成那种水(箱内浓度追平进水浓度,稳定平衡)
  • 指数逼近节奏:剩余差距按固定周期减半(t1/2=τln2t_{1/2} = \tau\ln 2)——与咖啡、半衰期同宗
  • 同构清单:药物代谢、房间通风、湖泊治污——一次建模,四处复用

延伸阅读(可选,跳过不影响主线)

  • Paul's Online Notes · Modeling with First Order DEs——混合、人口、落体三类建模一网打尽,混合问题部分有"进出流量不等 + 水箱溢出截断"的完整例题,正好接住本课例题 2 留的现实提醒(英文,跳过不影响主线)

八课正课全部收官。下一课是单元关卡:15 题综合测验,把解法三连、定性双雄和两大应用串成一张网——考的不是记忆,是"看到方程知道抄哪把兵器"的判断力。

划选正文任意文字可高亮、批注或加入复习卡

讨论

载入中…