微分方程
- 指数增长/衰减
- dy/dt = ky → y = y0ekt
- 半衰期
- k = ln 2 / t½
- 牛顿冷却定律
- T(t) = Ta + (T0 − Ta)e−kt
- 积分因子
- y′ + p(x)y = q(x): μ = e∫p dx, (μy)′ = μq
- 混合槽
- dm/dt = r·c − (r/V)·m
公式表、每个专题的关键方法步骤和 16 道练习题,完整解答附在最后。
一杯 80 °C 的咖啡放在 23 °C 的房间里。 它遵循牛顿冷却定律 dT/dt = −k·(T − Tₐ)。 已知 k = 0.04 min⁻¹,求 15 min 后的温度 T。
求解初值问题 y′ + 2y/x = 5x2,x > 0,y(1) = −2。求 y(2)。
求所示 3×3 矩阵 A 的行列式 det A。
解方程组 −3x + 5y = 8 和 −2x − y = −7。求 x。
求 a = (2, 5, 4) 与 b = (−7, 1, −1) 的点积 a·b。
一个力 F = (−50, 0, 40) N 作用在点 A(−2, −1, 4)。坐标单位为米。 求该力对 原点 O(0, 0, 0) 的力矩的 y 分量 My。
将 y = √(4x) 之下、从 x = 0 到 x = 2 的区域绕 x 轴旋转。求该旋转体(抛物面体)的体积 V。
计算 I = ∫₀1.5 x²·eˣ dx。
z = 3 ln(x² + 4y²),其中 x = −t − 2,y = −t²。用链式法则求 t = 0 时的 dz/dt。
f(x, y) = x³ − 12x + y² − 4y 在 (2, 2) 处有一个临界点。计算该点的判别式 D = fxxfyy − fxy²。
求和 S = Σk=1∞ 2/(k(k + 3))。
用 f(x) = sin x 的 3 次麦克劳林多项式近似 f(−1.5)。给出 P₃(−1.5)。
用部分分式法求 F(s) = 20/(s(s + 3)(s + 6)) 对应的 f(t),并计算 f(0.25)。
某系统满足 y″ + 3y′ + 2y = 6·u(t),y(0) = 5,y′(0) = 0。在 t = 0 时施加阶跃输入 u(t) = 10。用拉普拉斯变换求 y(0.25)。
输入 x(t) = 4·e−6tu(t) 加到 h(t) = e−1tu(t) 的系统上。求 t = 1.5 s 时的输出 y(1.5)。
输入 x(t) = 2·cos(3.75t) 加到 h(t) = 5·e−5tu(t) 的系统上。求稳态输出相对于输入的相位 φ,单位为度。
答案: T = 54.28 °C
一阶微分方程:dy/dt = k·y 的解为 y = y₀·ekt(k > 0 时增长,k < 0 时衰减);半衰期 T 对应 k = ln 2/T。牛顿冷却定律 dT/dt = −k·(T − Tₐ) 给出 T(t) = Tₐ + (T₀ − Tₐ)·e−kt。线性方程 y′ + p(x)·y = q(x) 用积分因子 μ = e∫p dx 求解:(μy)′ = μq,所以 μy = ∫μq dx + C。
答案: y(2) = 7.25
一阶微分方程:dy/dt = k·y 的解为 y = y₀·ekt(k > 0 时增长,k < 0 时衰减);半衰期 T 对应 k = ln 2/T。牛顿冷却定律 dT/dt = −k·(T − Tₐ) 给出 T(t) = Tₐ + (T₀ − Tₐ)·e−kt。线性方程 y′ + p(x)·y = q(x) 用积分因子 μ = e∫p dx 求解:(μy)′ = μq,所以 μy = ∫μq dx + C。
答案: det A = −47
矩阵:2×2 矩阵的行列式为 ad − bc。对于 3×3,按某一行展开:每个元素乘以划去其所在行和列后剩下的 2×2 矩阵的行列式,符号按 + − + 交替。特征值 λ 满足 det(A − λI) = 0;对于 2×2,即 λ² −(迹)·λ + det = 0。
答案: x = 2.077
矩阵:2×2 矩阵的行列式为 ad − bc。对于 3×3,按某一行展开:每个元素乘以划去其所在行和列后剩下的 2×2 矩阵的行列式,符号按 + − + 交替。特征值 λ 满足 det(A − λI) = 0;对于 2×2,即 λ² −(迹)·λ + det = 0。
答案: a·b = −13
Dot product: a·b = axbx + ayby + azbz = |a|·|b|·cos θ, a scalar. It gives the angle between two vectors, cos θ = a·b/(|a|·|b|), and it is zero when they are perpendicular. The scalar projection of a onto b is a·b/|b|, and a constant force F moving through a displacement d does work W = F·d.
为负:两向量的夹角大于 90°。
答案: My = −120 N·m
Cross product: a × b = (aybz − azby, azbx − axbz, axby − aybx) is perpendicular to both vectors, and |a × b| = |a|·|b|·sin θ is the area of the parallelogram they span (the triangle has half of it). Order matters: b × a = −(a × b). The moment of a force F applied at point A about point O is MO = r × F, with r = A − O.
答案: V = 25.13
旋转体体积:把 y = f(x) 之下、从 x = a 到 x = b 的区域绕 x 轴旋转,得到一个旋转体,它在 x 处的横截面是半径为 f(x) 的圆盘。把圆盘相加得 V = π∫ₐᵇ f(x)² dx。当区域位于外曲线 R(x) 与内曲线 r(x) 之间时,每个薄片是一个垫圈:V = π∫ₐᵇ (R(x)² − r(x)²) dx。相减之前要先把半径平方。
答案: I = 3.602
分部积分:∫u·dv = u·v − ∫v·du,由乘积法则得出。把 u 选为求导后会变简单的因式(x 的幂或 ln x),把 dv 选为你能积分的因式(eˣ、sin x、cos x,或 u = ln x 时的 x)。对定积分,要在两个上下限处计算 u·v:∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du。含 x² 的因式需要做两轮。
答案: dz/dt = 3
梯度:∇f = (fx, fy) 指向最陡上升的方向,|∇f| 是最大变化率。沿方向 v 的方向导数为 Duf = ∇f·u,其中单位向量 u = v/|v|。在 (a, b) 处的切平面为 z = f(a, b) + fx(a, b)(x − a) + fy(a, b)(y − b);它给出 f 在该点附近的线性近似。链式法则:dz/dt = fx·dx/dt + fy·dy/dt。
答案: D = 24
临界点与二重积分:临界点满足 fx = fy = 0。用 D = fxxfyy − fxy² 做二阶导数检验:D > 0 且 fxx > 0 为局部极小值,D > 0 且 fxx < 0 为局部极大值,D < 0 为鞍点,D = 0 无法判断。对于一般区域上的二重积分,内层积分限取为围成该区域的曲线;对于圆盘和扇形,使用极坐标 x = r cos θ,y = r sin θ,dA = r dr dθ(不要忘了因子 r)。
答案: S = 1.222
收敛性:如果各项不趋于 0,级数发散(通项判别法)。Σ1/kᵖ 仅当 p > 1 时收敛,Σrᵏ 仅当 |r| < 1 时收敛。当一个级数的项在 k 很大时与这些级数的项相近时,可与之比较。比值判别法:L = lim |ak+1/ak|;根值判别法:L = lim |ak|1/k。L < 1 收敛,L > 1 发散,L = 1 无法判断。幂级数 Σck(x − x₀)ᵏ 在 |x − x₀| < R 时收敛,其中 R = lim |ck/ck+1|;项递减趋于 0 的交错级数收敛。
答案: P₃(−1.5) = −0.9375
泰勒多项式:在 x = a 附近,光滑函数可用 Pₙ(x) = f(a) + f′(a)(x − a) + f″(a)(x − a)²/2! + … + f⁽ⁿ⁾(a)(x − a)ⁿ/n! 近似;当 a = 0 时这就是麦克劳林多项式。标准级数:eˣ = Σxᵏ/k!,sin x = x − x³/3! + x⁵/5! − …,cos x = 1 − x²/2! + x⁴/4! − …,ln(1 + x) = x − x²/2 + x³/3 − …,1/(1 − x) = Σxᵏ。拉格朗日误差界:|f(x) − Pₙ(x)| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!,其中 M 是 |f⁽ⁿ⁺¹⁾| 在 a 与 x 之间的最大值。
精确值为 f(−1.5) = −0.997495;取更多项或取更靠近 0 的点时,误差 −0.059995 会变小。
答案: f(0.25) = 0.3093
拉普拉斯变换:F(s) = ∫₀∞ f(t)·e−st dt 把导数变成乘以 s:ℒ{f′} = sF(s) − f(0),ℒ{f″} = s²F(s) − s·f(0) − f′(0)。常用变换对:1 ↔ 1/s,e−at ↔ 1/(s + a),t·e−at ↔ 1/(s + a)²,sin ωt ↔ ω/(s² + ω²),cos ωt ↔ s/(s² + ω²);平移 s → s + σ 相当于把时间函数乘以 e−σt。要对有理函数 F(s) 求逆变换,先分解成部分分式:单极点 p 的留数为 A = [(s − p)·F(s)] 在 s = p 处的值。终值定理:f(∞) = lims→0 s·F(s),仅当 s·F(s) 的所有极点都在左半平面时成立。初值定理:f(0⁺) = lims→∞ s·F(s)。
答案: y(0.25) = 6.223
拉普拉斯变换:F(s) = ∫₀∞ f(t)·e−st dt 把导数变成乘以 s:ℒ{f′} = sF(s) − f(0),ℒ{f″} = s²F(s) − s·f(0) − f′(0)。常用变换对:1 ↔ 1/s,e−at ↔ 1/(s + a),t·e−at ↔ 1/(s + a)²,sin ωt ↔ ω/(s² + ω²),cos ωt ↔ s/(s² + ω²);平移 s → s + σ 相当于把时间函数乘以 e−σt。要对有理函数 F(s) 求逆变换,先分解成部分分式:单极点 p 的留数为 A = [(s − p)·F(s)] 在 s = p 处的值。终值定理:f(∞) = lims→0 s·F(s),仅当 s·F(s) 的所有极点都在左半平面时成立。初值定理:f(0⁺) = lims→∞ s·F(s)。
检验:y(0) = A + B + C = 5。
答案: y(1.5) = 0.1784
连续卷积:y(t) = ∫ x(τ)·h(t − τ) dτ。把 h 翻转,平移到位置 t,对重叠部分的乘积积分。对因果信号,积分从 0 到 t。有用的结论:y 的面积等于 x 的面积乘以 h 的面积;若 x 持续时间从 t₁ 到 t₂,h 从 t₃ 到 t₄,则 y 持续时间从 t₁ + t₃ 到 t₂ + t₄。对于 h(t) = e−atu(t) 和高度为 A 的阶跃,y(t) = (A/a)(1 − e−at)。
答案: φ = −36.87 °
冲激响应与傅里叶级数:冲激响应 h(t) 完全描述一个线性时不变系统:对于 y′ + a·y = b·x,h(t) = b·e−atu(t);对于 y[n] = a·y[n − 1] + b·x[n],h[n] = b·aⁿu[n]。正弦信号 cos ωt 通过系统时,增益为 |H(jω)|,相位为 ∠H(jω),其中 H 是 h 的傅里叶变换。周期信号是频率为 k/T 的各次谐波之和:幅值 ±A 的方波只含奇次谐波,幅值为 4A/(kπ);峰值为 A 的三角波的谐波幅值为 8A/(π²k²);高度为 A、占空比为 D 的脉冲串,其平均值为 A·D,谐波为 (2A/(kπ))·|sin(kπD)|。
ExamProblems · examproblems.com/courses/engineering-mathematics-2