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)35Ta​141 °C321 °C

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

T0 = 500 °CTa = 23 °CT(35) = 141 °CT = 321 °C
t = min
2.Persamaan diferensialSulit
tyK = 5ymax​ = ?

Sebuah sistem memenuhi y″ + 4y′ + 8y = 40 untuk t ≥ 0, dimulai dari keadaan diam: y(0) = 0, y′(0) = 0. Cari nilai puncak ymax dari respons langkahnya.

y″ + 4y′ + 8y = 40y(0) = 0y′(0) = 0
ymax =
3.IntegralSedang
xy02

Hitung I = ∫02 2x·(x² + 1)1 dx.

f(x) = 2x·(x² + 1)1[0, 2]
I =
4.IntegralSulit
xyy = 0,5x2​y = x

Tentukan luas A daerah yang dibatasi oleh y = 0,5x2 dan y = x.

y = 0,5x2y = x
A =
5.TurunanSedang
xyx = 0

Misalkan f(x) = (−2x − 1)5. Tentukan f′(0).

f(x) = (−2x − 1)5x = 0
f′(0) =
6.TurunanSulit
xy

Misalkan f(x) = x3 − 4,5x2 + 6x + 6. Tentukan koordinat x x dari minimum lokal.

f(x) = x3 − 4,5x2 + 6x + 6
x =
7.Metode numerikSedang
xy00,20,40,60,81

Hampiri ∫01 √x·e−x dx dengan kaidah trapesium memakai n = 5 pita. Berikan T.

f(x) = √x·e−x[a, b] = [0, 1]n = 5
T =
8.Metode numerikSulit
00,5112ty(0, 0,5)

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

y′ = t + yy(0) = 0,5h = 0,25
|E| =

3Larutan

1.Persamaan diferensial

Jawaban: t = 11,79 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(477/118)/35 = 0,03991 min⁻¹

Waktu untuk mencapai target

t = ln(477/298)/0,03991 = 11,79 min

Jawaban

t = 11,79 min
2.Persamaan diferensial

Jawaban: ymax = 5,216

Persamaan linear orde dua: untuk y″ + b·y′ + c·y = 0 coba y = ert: ini memberi persamaan karakteristik r² + b·r + c = 0. Dua akar real r₁ ≠ r₂ memberi y = A·er₁t + B·er₂t (redaman lebih); akar ganda r memberi y = (A + Bt)·ert (redaman kritis); akar kompleks −α ± jω memberi y = e−αt(A cos ωt + B sin ωt) (redaman kurang, berosilasi). Nilai awal y(0) dan y′(0) menentukan A dan B. Ditulis sebagai y″ + 2ζωny′ + ωn²y = 0, ωn adalah frekuensi alami dan ζ rasio redaman: ζ < 1 berosilasi. Dengan masukan konstan, y″ + by′ + cy = c·K, penyelesaiannya adalah nilai tunak K ditambah respons bebas.

Nilai tunak dan akar-akar

yss = 5, r = −2 ± 2j

Konstanta dari nilai awal

y(0) = 5 + A = 0 → A = −5
y′(0) = −2A + 2B = 0 → B = −5

Puncak

t = π/2 = 1,571
ymax = 5·(1 + e−2·π/2) = 5,216

Jawaban

ymax = 5,216
3.Integral

Jawaban: I = 12

Integral tentu: ∫ₐᵇ f(x) dx adalah luas bertanda antara grafik dan sumbu x dari a ke b. Cari antiturunan F (dengan F′ = f) dan hitung F(b) − F(a). Aturan pangkat: ∫xⁿ dx = xⁿ⁺¹/(n + 1). Luas antara dua kurva adalah ∫(atas − bawah) dx.

Substitusi

u = x² + 1, du = 2x dx
I = ∫15 u1 du

Integralkan dan substitusikan batas

I = [u2/2]15 = (25 − 1)/2 = 12

Jawaban

I = 12
4.Integral

Jawaban: A = 0,6667

Integral tentu: ∫ₐᵇ f(x) dx adalah luas bertanda antara grafik dan sumbu x dari a ke b. Cari antiturunan F (dengan F′ = f) dan hitung F(b) − F(a). Aturan pangkat: ∫xⁿ dx = xⁿ⁺¹/(n + 1). Luas antara dua kurva adalah ∫(atas − bawah) dx.

Titik potong

0,5x2 = x → x = 0, x = 2

Integralkan atas dikurangi bawah

A = [0,5x2 − 0,1667x3]02 = 0,6667

Jawaban

A = 0,6667
5.Turunan

Jawaban: f′(0) = −10

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

f′(x) = 5·(−2x − 1)4·(−2)

Substitusikan

f′(0) = 5·(−1)4·(−2) = −10

Jawaban

f′(0) = −10
6.Turunan

Jawaban: x = 2

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.

Turunan

f′(x) = 3x2 − 9x + 6 = 3·(x2 − 3x + 2)

Selesaikan f′(x) = 0

x = 1 atau x = 2

Mana yang mana

f″(x) = 6x − 9: negatif di maksimum, positif di minimum.

f″(1) = −3, f″(2) = 3
x = 2

Jawaban

x = 2
7.Metode numerik

Jawaban: T = 0,3602

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 badan000110,20,366148220,40,423948230,60,425108240,80,4018922510,3678791

Simpul dan nilai

h = 0,2
Σ wᵢfᵢ = 0 + 2·0,366148 + 2·0,423948 + 2·0,425108 + 2·0,401892 + 0,367879 = 3,602069

Kaidah trapesium

T = (0,2/2)·3,602069 = 0,3602

Jawaban

T = 0,3602

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

8.Metode numerik

Jawaban: |E| = 0,4153

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,50,50,62510,250,6250,8750,8437520,50,843751,343751,17968830,751,1796881,9296881,662109

Langkah Euler

y4 = 1,662109

Nilai eksak

y(1) = 2,077423

Galat global

|E| = |2,077423 − 1,662109| = 0,4153
00,5112ty

Grafik

Jawaban

|E| = 0,4153

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

ExamProblems · examproblems.com/courses/differential-equations