微分方程式
- 指数関数的な増加・減少
- 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 問の演習問題。完全な解答は巻末にあります。
y(0) = 3、y′(0) = 0 として y″ + 5y′ + 4y = 0 を解きなさい。y(1.5) を求めなさい。
系が t ≥ 0 で y″ + 3y′ + 18.25y = 182.5 に従い、静止状態 y(0) = 0、y′(0) = 0 から始まります。ステップ応答のピーク値 ymax を求めなさい。
示された行列について C = AB とします。成分 c12(第 1 行、第 2 列)を求めなさい。
示された行列 A の大きい方の固有値 λ を求めなさい。
ベクトル a = (−6, 4, 2) と b = (1, 2, 1) のなす角 θ を度で求めよ。
力 F = (−40, −60, 50) N が点 A(−3, −3, 3) に作用している。座標の単位はメートルである。 点 O(−2, 1, −1) まわりの力のモーメントの z 成分 Mz を求めよ。
I = ∫₀1 x·cos x dx(x はラジアン)を計算しなさい。
放物線 y = 4 − x² と x 軸で囲まれた部分の面積 A を求めなさい。
(x, y) = (3, −2) の点での z = −4xy³ + 4x² − 2y² − 3xy の接平面は、z = px + qy + c と書けます。c を求めてください。
f(x, y) = x³ − 3x + y² の臨界点を求めて分類し、極小での f の値 fmin を答えてください。
和 S = Σk=1∞ 5/(k(k + 4)) を求めてください。
f(x) = ln(1 + x) を、x = 0.2 で 2 次のマクローリン多項式で近似します。M = 2! = 2 としてラグランジュの誤差評価 M·|x|ⁿ⁺¹/(n + 1)! を使い、誤差の最大値 |R₂| を求めてください。
F(s) = (4s + 28)/(s² + 2s + 1) の f(t) を求め、f(0.6) を計算せよ。
ラプラス変換を用いて、y(0) = 1、y′(0) = −15 の条件で y″ + 4y′ + 3y = 0 を解き、y(0.25) を求めよ。
入力 x(t) = 1·e−2tu(t) が、h(t) = e−1tu(t) のシステムに加えられます。t = 2 s での出力 y(2) を求めなさい。
パルス列の高さは 5、パルス幅は 0.3 ms、周期は 1 ms です。1000 Hz にある高調波次数 1 の振幅 c1 を求めなさい。
答え: y(1.5) = 0.89
2 階線形方程式:y″ + b·y′ + c·y = 0 では y = ert とおきます。これから特性方程式 r² + b·r + c = 0 が得られます。異なる 2 つの実根 r₁ ≠ r₂ のときは y = A·er₁t + B·er₂t(過減衰)、重根 r のときは y = (A + Bt)·ert(臨界減衰)、複素根 −α ± 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 に自由応答を加えたものです。
答え: ymax = 13.08
2 階線形方程式:y″ + b·y′ + c·y = 0 では y = ert とおきます。これから特性方程式 r² + b·r + c = 0 が得られます。異なる 2 つの実根 r₁ ≠ r₂ のときは y = A·er₁t + B·er₂t(過減衰)、重根 r のときは y = (A + Bt)·ert(臨界減衰)、複素根 −α ± 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 に自由応答を加えたものです。
答え: c12 = 15
行列:2×2 行列の行列式は ad − bc です。3×3 では 1 つの行に沿って展開します:各成分に、その行と列を除いた 2×2 行列の行列式を掛け、符号は + − + と交互にします。固有値 λ は det(A − λI) = 0 を満たし、2×2 では λ² −(トレース)·λ + det = 0 です。
答え: λ = 2
行列:2×2 行列の行列式は ad − bc です。3×3 では 1 つの行に沿って展開します:各成分に、その行と列を除いた 2×2 行列の行列式を掛け、符号は + − + と交互にします。固有値 λ は det(A − λI) = 0 を満たし、2×2 では λ² −(トレース)·λ + det = 0 です。
答え: θ = 77.4 °
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.
答え: Mz = −100 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.
答え: I = 0.3818
部分積分:∫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² があるときは 2 回行います。
答え: A = 10.67
定積分:∫ₐᵇ f(x) dx は、a から b までのグラフと x 軸の間の符号付き面積です。原始関数 F(F′ = f)を求め、F(b) − F(a) を計算します。べき乗の法則:∫xⁿ dx = xⁿ⁺¹/(n + 1)。2 曲線の間の面積は ∫(上 − 下) dx です。
答え: c = −334
勾配:∇f = (fx, fy) は最も急に増加する方向を向き、|∇f| は変化率の最大値です。方向 v に沿った方向微分は、単位ベクトル u = v/|v| を用いて Duf = ∇f·u です。点 (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。
答え: fmin = −2
臨界点と二重積分:臨界点では fx = fy = 0 です。D = fxxfyy − fxy² を用いた 2 階微分判定法:D > 0 かつ fxx > 0 なら極小、D > 0 かつ fxx < 0 なら極大、D < 0 なら鞍点、D = 0 なら何も決まりません。一般の領域上の二重積分では、内側の積分の限界を領域を囲む曲線にします。円板や扇形では極座標 x = r cos θ、y = r sin θ、dA = r dr dθ を使います(係数 r を忘れないでください)。
fx = 0 と fy = 0 を同時に解きましょう。
答え: S = 2.604
収束:項が 0 に近づかなければ、級数は発散します(第 n 項判定法)。Σ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 に減少していく交代級数は収束します。
答え: |R₂| = 0.002667
テイラー多項式: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(0.2) − P₂(0.2)| = 0.002322。
答え: f(0.6) = 10.1
ラプラス変換: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) = −1.366
ラプラス変換: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 = 1。
答え: y(2) = 0.117
連続たたみ込み: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) です。
答え: c1 = 2.575
インパルス応答とフーリエ級数:インパルス応答 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)| をもちます。
ExamProblems · examproblems.com/courses/engineering-mathematics-2