Persamaan Diferensial: paket kilat

Lembar rumus, langkah-langkah metode utama untuk setiap topik, dan 8 soal latihan, dengan pembahasan lengkap di bagian belakang.

Soal baruKembali ke mata kuliah

Di jendela cetak, pilih “Simpan sebagai PDF” untuk menyimpan salinan di perangkat Anda.

1Rumus dan metode

Persamaan diferensial

Pertumbuhan/peluruhan eksponensial
dy/dt = ky → y = y0ekt
Waktu paruh
k = ln 2 / t½
Hukum pendinginan Newton
T(t) = Ta + (T0 − Ta)e−kt
Faktor integrasi
y′ + p(x)y = q(x): μ = e∫p dx, (μy)′ = μq
Tangki pencampuran
dm/dt = r·c − (r/V)·m

Integral

Aturan pangkat
∫xⁿ dx = xⁿ⁺¹/(n + 1) + C
Integral tentu
∫ab f dx = F(b) − F(a)
Integral umum
∫eˣ dx = eˣ, ∫1/x dx = ln|x|, ∫sin x dx = −cos x
Luas daerah antara kurva
∫ab (f − g) dx

Metode

  1. Tulis ulang integran sebagai jumlah suku-suku sederhana bila memungkinkan.
  2. Terapkan aturan pangkat: ∫xⁿ dx = xⁿ⁺¹/(n + 1), dan integral standar untuk ex, 1/x, sin, dan cos.
  3. Gunakan substitusi saat Anda melihat suatu fungsi bersama turunannya.
  4. Untuk integral tentu, hitung F(b) − F(a).

Turunan

Aturan pangkat
(xⁿ)′ = n·xⁿ⁻¹
Aturan perkalian
(uv)′ = u′v + uv′
Aturan hasil bagi
(u/v)′ = (u′v − uv′)/v²
Aturan rantai
f(g(x))′ = f′(g(x))·g′(x)
Turunan umum
(eˣ)′ = eˣ, (ln x)′ = 1/x, (sin x)′ = cos x, (cos x)′ = −sin x

Metode

  1. Kenali strukturnya: jumlah, hasil kali, hasil bagi, atau komposisi.
  2. Terapkan aturan pangkat pada setiap pangkat x.
  3. Gunakan aturan hasil kali, hasil bagi, atau rantai bila diperlukan.
  4. Sederhanakan, dan hitung nilainya di titik tersebut jika diminta.

Metode numerik

Bisection
c = (a + b)/2, width after n steps (b − a)/2ⁿ
Newton–Raphson
xn+1 = xₙ − f(xₙ)/f′(xₙ)
Interpolasi linear
y = y₁ + (y₂ − y₁)(x − x₁)/(x₂ − x₁)
Kaidah trapesium
T = (h/2)[f₀ + 2f₁ + … + 2fn−1 + fₙ], h = (b − a)/n
Kaidah Simpson (n genap)
S = (h/3)[f₀ + 4f₁ + 2f₂ + … + 4fn−1 + fₙ]
Batas galat trapesium
|E| ≤ (b − a)h²M/12, M = max|f″|
Euler
yn+1 = yₙ + h·f(tₙ, yₙ)
Heun
yn+1 = yₙ + (h/2)(k₁ + k₂), k₂ = f(tₙ + h, yₙ + hk₁)
RK4
yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄)

Metode

  1. Untuk akar: periksa pergantian tanda pada bisection, atau hitung f dan f′ untuk Newton, dan pertahankan enam angka atau lebih pada tabel iterasi.
  2. Untuk nilai tabel: interpolasikan secara linear di antara dua titik yang bertetangga.
  3. Untuk integral: cari h = (b − a)/n, daftarkan f di setiap simpul, terapkan bobot (1, 2, …, 2, 1 untuk trapesium; 1, 4, 2, …, 4, 1 untuk Simpson dengan n genap) dan kalikan dengan h/2 atau h/3.
  4. Untuk PD: melangkahlah dari t₀ dengan yₙ₊₁ = yₙ + h·f(tₙ, yₙ) (Euler) atau kemiringan rata-rata (Heun), satu baris per langkah, sampai mencapai t yang diminta.

2Soal latihan

1.Persamaan diferensialSedang
t (min)T (°C)14Ta​204 °C123 °C

Sebuah komponen baja tempa bersuhu 340 °C dibiarkan mendingin di bengkel bersuhu 23 °C. Setelah 14 menit suhunya 204 °C. Benda tersebut mengikuti hukum pendinginan Newton, dT/dt = −k·(T − Tₐ). Berapa lama setelah awal, t, yang diperlukan untuk mencapai 123 °C?

T0 = 340 °CTa = 23 °CT(14) = 204 °CT = 123 °C
t = min
2.Persamaan diferensialSulit
xy13(1, 1)y(3) = ?

Selesaikan masalah nilai awal y′ + 2y/x = 3 untuk x > 0, dengan y(1) = 1. Cari y(3).

y′ + 2y/x = 3y(1) = 1
y(3) =
3.IntegralSedang
xy1y = x·cos x

Hitung I = ∫₀1 x·cos x dx (x dalam radian).

∫₀1 x·cos x dx
I =
4.IntegralSulit
xy1y = x²eˣ

Hitung I = ∫₀1 x²·eˣ dx.

∫₀1 x²eˣ dx
I =
5.TurunanSedang
xyx = −2

Misalkan f(x) = ln(2x2 + 4). Tentukan f′(−2).

f(x) = ln(2x2 + 4)x = −2
f′(−2) =
6.TurunanSulit
xyx = 3

Garis singgung f(x) = 3x2 + 3x − 6 di x = 3 memotong sumbu y di y = b. Tentukan b.

f(x) = 3x2 + 3x − 6x = 3
b =
7.Metode numerikSedang
xy00,33330,666711,3331,6672

Hampiri ∫02 1/(1 + x³) dx dengan kaidah Simpson memakai n = 6 pita. Berikan S.

f(x) = 1/(1 + x³)[a, b] = [0, 2]n = 6
S =
8.Metode numerikSulit
012012ty(0, 0,2)

Selesaikan y′ = y(2 − y), y(0) = 0,2 dengan metode Euler dan langkah h = 0,5 sampai t = 2,5. Solusi eksaknya y = 2y₀/(y₀ + (2 − y₀)e−2t). Carilah galat global |E| = |y(2,5) − y5|.

y′ = y(2 − y)y(0) = 0,2h = 0,5
|E| =

3Larutan

1.Persamaan diferensial

Jawaban: t = 28,82 min

Persamaan diferensial orde satu: dy/dt = k·y memiliki penyelesaian y = y₀·ekt (pertumbuhan untuk k > 0, peluruhan untuk k < 0); waktu paruh T memberikan k = ln 2/T. Hukum pendinginan Newton, dT/dt = −k·(T − Tₐ), memberikan T(t) = Tₐ + (T₀ − Tₐ)·e−kt. Persamaan linear y′ + p(x)·y = q(x) diselesaikan dengan faktor integrasi μ = e∫p dx: (μy)′ = μq, sehingga μy = ∫μq dx + C.

Konstanta pendinginan

k = ln(317/181)/14 = 0,04003 min⁻¹

Waktu untuk mencapai target

t = ln(317/100)/0,04003 = 28,82 min

Jawaban

t = 28,82 min
2.Persamaan diferensial

Jawaban: y(3) = 3

Persamaan diferensial orde satu: dy/dt = k·y memiliki penyelesaian y = y₀·ekt (pertumbuhan untuk k > 0, peluruhan untuk k < 0); waktu paruh T memberikan k = ln 2/T. Hukum pendinginan Newton, dT/dt = −k·(T − Tₐ), memberikan T(t) = Tₐ + (T₀ − Tₐ)·e−kt. Persamaan linear y′ + p(x)·y = q(x) diselesaikan dengan faktor integrasi μ = e∫p dx: (μy)′ = μq, sehingga μy = ∫μq dx + C.

Faktor integrasi

μ = e∫2/x dx = x2

Integralkan

(x2y)′ = 3x2 → x2y = 3x3/3 + C

Kondisi awal

1 = 3/3 + C → C = 0

Hitung nilainya

y(3) = (3·33/3 + 0)/32 = 3

Jawaban

y(3) = 3
3.Integral

Jawaban: I = 0,3818

Integrasi parsial: ∫u·dv = u·v − ∫v·du, dari aturan perkalian. Pilih u sebagai faktor yang menjadi lebih sederhana saat diturunkan (pangkat x, atau ln x), dan dv sebagai faktor yang dapat Anda integralkan (eˣ, sin x, cos x, atau x bila u = ln x). Untuk integral tentu, hitung u·v pada kedua batas: ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. Faktor x² memerlukan dua putaran.

Pilih u dan dv

u = x, du = dx; dv = cos x dx, v = sin x

Antiturunan

∫x cos x dx = x sin x − ∫sin x dx = x sin x + cos x

Limit

(1·sin 1 + cos 1) − (0 + 1) = 0,3818

Jawaban

I = 0,3818
4.Integral

Jawaban: I = 0,7183

Integrasi parsial: ∫u·dv = u·v − ∫v·du, dari aturan perkalian. Pilih u sebagai faktor yang menjadi lebih sederhana saat diturunkan (pangkat x, atau ln x), dan dv sebagai faktor yang dapat Anda integralkan (eˣ, sin x, cos x, atau x bila u = ln x). Untuk integral tentu, hitung u·v pada kedua batas: ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. Faktor x² memerlukan dua putaran.

Putaran pertama

∫x²eˣ dx = x²eˣ − ∫2x eˣ dx

Putaran kedua

∫2x eˣ dx = 2x eˣ − 2eˣ
∫x²eˣ dx = eˣ(x² − 2x + 2)

Limit

2,718·(1) − 1·2 = 0,7183

Jawaban

I = 0,7183
5.Turunan

Jawaban: f′(−2) = −0,6667

Diferensial: turunan f′(x) adalah kemiringan grafik di x. Aturan pangkat: (xⁿ)′ = n·xⁿ⁻¹. Aturan perkalian: (uv)′ = u′v + uv′. Aturan rantai: (f(g(x)))′ = f′(g(x))·g′(x). Di mana f′(x) = 0, grafik memiliki garis singgung mendatar: maksimum, minimum, atau titik pelana.

Aturan rantai untuk logaritma

f′(x) = 4x/(2x2 + 4)

Substitusikan

f′(−2) = −8/12 = −0,6667

Jawaban

f′(−2) = −0,6667
6.Turunan

Jawaban: b = −33

Diferensial: turunan f′(x) adalah kemiringan grafik di x. Aturan pangkat: (xⁿ)′ = n·xⁿ⁻¹. Aturan perkalian: (uv)′ = u′v + uv′. Aturan rantai: (f(g(x)))′ = f′(g(x))·g′(x). Di mana f′(x) = 0, grafik memiliki garis singgung mendatar: maksimum, minimum, atau titik pelana.

Titik dan kemiringan

f(3) = 30
f′(x) = 6x + 3 → f′(3) = 21

Garis singgung di x = 0

b = f(x₀) − f′(x₀)·x₀ = 30 − 21·3 = −33

Jawaban

b = −33
7.Metode numerik

Jawaban: S = 1,091

Integrasi numerik: bagi [a, b] menjadi n pita selebar h = (b − a)/n dengan simpul x₀ … xₙ. Kaidah trapesium: T = (h/2)[f₀ + 2f₁ + 2f₂ + … + 2fn−1 + fₙ]. Kaidah Simpson (n genap): S = (h/3)[f₀ + 4f₁ + 2f₂ + 4f₃ + … + 4fn−1 + fₙ]. Galat trapesium mengecil seperti h² (batas (b − a)h²·M/12 dengan M = maks|f″|), Simpson seperti h⁴, jadi Simpson biasanya jauh lebih akurat untuk simpul yang sama.

ixᵢf(xᵢ)berat badan001110,3333330,964286420,6666670,7714292310,5441,3333330,296703251,6666670,1776324620,1111111

Simpul dan nilai

h = 0,3333
Σ wᵢfᵢ = 1 + 4·0,964286 + 2·0,771429 + 4·0,5 + 2·0,296703 + 4·0,177632 + 0,111111 = 9,815044

Kaidah Simpson

S = (0,3333/3)·9,815044 = 1,091

Jawaban

S = 1,091

Integral ini tidak memiliki antiturunan elementer, dan justru pada saat seperti itulah kaidah numerik diperlukan.

8.Metode numerik

Jawaban: |E| = 0,04568

Penyelesaian numerik y′ = f(t, y): metode Euler melangkah menurut kemiringan di awal setiap langkah: yn+1 = yₙ + h·f(tₙ, yₙ). Metode Heun (Euler yang diperbaiki) merata-ratakan kemiringan di awal dan pada prediksi Euler: k₁ = f(tₙ, yₙ), k₂ = f(tₙ + h, yₙ + h·k₁), yn+1 = yₙ + (h/2)(k₁ + k₂). Metode Runge–Kutta klasik (RK4) menggabungkan empat kemiringan: yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄). Galat global Euler sebanding dengan h, Heun dengan h², RK4 dengan h⁴.

ntₙyₙf(tₙ, yₙ)yₙ₊₁000,20,360,3810,50,380,61560,6878210,68780,9025311,13906631,51,1390660,9806611,629396421,6293960,6038611,931326

Langkah Euler

y5 = 1,931326

Nilai eksak

y(2,5) = 1,885651

Galat global

|E| = |1,885651 − 1,931326| = 0,04568
012012ty

Grafik

Jawaban

|E| = 0,04568

Membagi dua h kira-kira membagi dua galat metode Euler.

ExamProblems · examproblems.com/courses/differential-equations