工学数学 2：直前対策パック

1 階微分方程式：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 階線形方程式：y″ + b·y′ + c·y = 0 では y = e^rt とおきます。これから特性方程式 r² + b·r + c = 0 が得られます。異なる 2 つの実根 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 に自由応答を加えたものです。

行列：2×2 行列の行列式は ad − bc です。3×3 では 1 つの行に沿って展開します：各成分に、その行と列を除いた 2×2 行列の行列式を掛け、符号は + − + と交互にします。固有値 λ は det(A − λI) = 0 を満たし、2×2 では λ² −（トレース）·λ + det = 0 です。

行列：2×2 行列の行列式は ad − bc です。3×3 では 1 つの行に沿って展開します：各成分に、その行と列を除いた 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.

定積分：∫ₐᵇ f(x) dx は、a から b までのグラフと x 軸の間の符号付き面積です。原始関数 F（F′ = f）を求め、F(b) − F(a) を計算します。べき乗の法則：∫xⁿ dx = xⁿ⁺¹/(n + 1)。2 曲線の間の面積は ∫(上 − 下) dx です。

回転体の体積：x = a から x = b までの y = f(x) の下の領域を x 軸のまわりに回転させると、x での断面が半径 f(x) の円板である立体ができます。円板を足し合わせると V = π∫ₐᵇ f(x)² dx になります。領域が外側の曲線 R(x) と内側の曲線 r(x) の間にあるときは、各スライスはワッシャー（円環）で、V = π∫ₐᵇ (R(x)² − r(x)²) dx です。引く前に半径を 2 乗しましょう。

勾配：∇f = (f_x, f_y) は最も急に増加する方向を向き、|∇f| は変化率の最大値です。方向 v に沿った方向微分は、単位ベクトル u = v/|v| を用いて D_uf = ∇f·u です。点 (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² を用いた 2 階微分判定法：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 に近づかなければ、級数は発散します（第 n 項判定法）。Σ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 は a と x の間での |f⁽ⁿ⁺¹⁾| の最大値です。

ラプラス変換：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) は LTI システムを完全に表します：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)| をもちます。