工程数学 2：速成包

二阶线性方程：对 y″ + b·y′ + c·y = 0，试 y = e^rt：得到特征方程 r² + b·r + c = 0。两个不等实根 r₁ ≠ r₂ 给出 y = A·e^r₁t + B·e^r₂t（过阻尼）；二重根 r 给出 y = (A + Bt)·e^rt（临界阻尼）；复根 −α ± jω 给出 y = e^−αt(A cos ωt + B sin ωt)（欠阻尼，振荡）。初值 y(0) 和 y′(0) 决定 A 和 B。写成 y″ + 2ζω_ny′ + ω_n²y = 0 时，ω_n 是固有频率，ζ 是阻尼比：ζ < 1 时振荡。在常值输入下，y″ + by′ + cy = c·K 的解是稳态值 K 加上自由响应。

一阶微分方程：dy/dt = k·y 的解为 y = y₀·e^kt（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。

矩阵：2×2 矩阵的行列式为 ad − bc。对于 3×3，按某一行展开：每个元素乘以划去其所在行和列后剩下的 2×2 矩阵的行列式，符号按 + − + 交替。特征值 λ 满足 det(A − λI) = 0；对于 2×2，即 λ² −（迹）·λ + det = 0。

矩阵：2×2 矩阵的行列式为 ad − bc。对于 3×3，按某一行展开：每个元素乘以划去其所在行和列后剩下的 2×2 矩阵的行列式，符号按 + − + 交替。特征值 λ 满足 det(A − λI) = 0；对于 2×2，即 λ² −（迹）·λ + det = 0。

Dot product: a·b = a_xb_x + a_yb_y + a_zb_z = |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.

Cross product: a × b = (a_yb_z − a_zb_y, a_zb_x − a_xb_z, a_xb_y − a_yb_x) 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 M_O = r × F, with r = A − O.

旋转体体积：把 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。相减之前要先把半径平方。

分部积分：∫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² 的因式需要做两轮。

梯度：∇f = (f_x, f_y) 指向最陡上升的方向，|∇f| 是最大变化率。沿方向 v 的方向导数为 D_uf = ∇f·u，其中单位向量 u = v/|v|。在 (a, b) 处的切平面为 z = f(a, b) + f_x(a, b)(x − a) + f_y(a, b)(y − b)；它给出 f 在该点附近的线性近似。链式法则：dz/dt = f_x·dx/dt + f_y·dy/dt。

临界点与二重积分：临界点满足 f_x = f_y = 0。用 D = f_xxf_yy − f_xy² 做二阶导数检验：D > 0 且 f_xx > 0 为局部极小值，D > 0 且 f_xx < 0 为局部极大值，D < 0 为鞍点，D = 0 无法判断。对于一般区域上的二重积分，内层积分限取为围成该区域的曲线；对于圆盘和扇形，使用极坐标 x = r cos θ，y = r sin θ，dA = r dr dθ（不要忘了因子 r）。

收敛性：如果各项不趋于 0，级数发散（通项判别法）。Σ1/kᵖ 仅当 p > 1 时收敛，Σrᵏ 仅当 |r| < 1 时收敛。当一个级数的项在 k 很大时与这些级数的项相近时，可与之比较。比值判别法：L = lim |a_k+1/a_k|；根值判别法：L = lim |a_k|^1/k。L < 1 收敛，L > 1 发散，L = 1 无法判断。幂级数 Σc_k(x − x₀)ᵏ 在 |x − x₀| < R 时收敛，其中 R = lim |c_k/c_k+1|；项递减趋于 0 的交错级数收敛。

泰勒多项式：在 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(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(∞) = lim_s→0 s·F(s)，仅当 s·F(s) 的所有极点都在左半平面时成立。初值定理：f(0⁺) = lim_s→∞ s·F(s)。

拉普拉斯变换：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(∞) = lim_s→0 s·F(s)，仅当 s·F(s) 的所有极点都在左半平面时成立。初值定理：f(0⁺) = lim_s→∞ s·F(s)。

连续卷积：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)。

冲激响应与傅里叶级数：冲激响应 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)|。