Untuk punca: semak perubahan tanda untuk pembahagi dua, atau hitung f dan f′ untuk Newton, dan simpan enam digit atau lebih dalam jadual lelaran.
Untuk nilai jadual: interpolasi secara linear antara dua titik jiran.
Untuk kamiran: cari h = (b − a)/n, senaraikan f pada setiap nod, gunakan pemberat (1, 2, …, 2, 1 untuk trapezium; 1, 4, 2, …, 4, 1 untuk Simpson dengan n genap) dan darab dengan h/2 atau h/3.
Untuk PPB: melangkah dari t₀ dengan yₙ₊₁ = yₙ + h·f(tₙ, yₙ) (Euler) atau cerun purata (Heun), satu baris bagi setiap langkah, sehingga anda sampai ke t yang diminta.
2Soalan latihan
1.Persamaan pembezaanSederhana
Sekaleng minuman bergas pada 3 °C dikeluarkan daripada peti sejuk ke dalam bilik pada 24 °C. Selepas 6 min suhunya ialah 12 °C. Ia mengikut hukum penyejukan Newton, dT/dt = −k·(T − Tₐ). Berapa lama selepas permulaan, t, untuk mencapai 7 °C?
T0 = 3 °CTa = 24 °CT(6) = 12 °CT = 7 °C
t = min
2.Persamaan pembezaanSukar
Satu sistem mematuhi y″ + y′ + 9.25y = 9.25 untuk t ≥ 0, bermula dari pegun: y(0) = 0, y′(0) = 0. Cari y(0.25); sudut dalam radian.
y″ + y′ + 9.25y = 9.25y(0) = 0y′(0) = 0
y(0.25) =
3.KamiranSederhana
Nilaikan I = ∫26 4/x dx.
f(x) = 4/x[2, 6]
I =
4.KamiranSukar
Cari luas A kawasan antara parabola y = 4 − x² dan paksi-x.
y = 4 − x²
A =
5.TerbitanSederhana
Biar f(x) = ln(2x2 + 5). Cari f′(−2).
f(x) = ln(2x2 + 5)x = −2
f′(−2) =
6.TerbitanSukar
Biar f(x) = x3 − 4.5x2 + 6x + 6. Cari koordinat-x x bagi minimum setempat.
f(x) = x3 − 4.5x2 + 6x + 6
x =
7.Kaedah berangkaSederhana
Anggarkan ∫01 eˣ dx dengan petua Simpson, n = 2, dan bandingkan dengan nilai tepat. Cari ralat mutlak |E| = |S − I|.
f(x) = eˣ[a, b] = [0, 1]n = 2
|E| =
8.Kaedah berangkaSukar
Selesaikan y′ = −y, y(0) = 2 dengan kaedah Euler dan langkah h = 0.1 sehingga t = 0.4. Penyelesaian tepat ialah y = y₀e−t. Cari ralat global |E| = |y(0.4) − y4|.
y′ = −yy(0) = 2h = 0.1
|E| =
3Larutan
1.Persamaan pembezaan
Jawapan: t = 2.266 min
Persamaan pembezaan peringkat pertama: dy/dt = k·y mempunyai penyelesaian y = y₀·ekt (pertumbuhan bagi k > 0, pereputan bagi k < 0); separuh hayat T memberikan k = ln 2/T. Hukum penyejukan 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 pengamiran μ = e∫p dx: (μy)′ = μq, maka μy = ∫μq dx + C.
Pemalar penyejukan
k = ln(-21/-12)/6 = 0.09327 min⁻¹
Masa untuk mencapai sasaran
t = ln(-21/(−17))/0.09327 = 2.266 min
Jawapan
t = 2.266 min
2.Persamaan pembezaan
Jawapan: y(0.25) = 0.254
Persamaan linear tertib kedua: bagi y″ + b·y′ + c·y = 0 cuba y = ert: ini memberikan persamaan ciri r² + b·r + c = 0. Dua punca nyata r₁ ≠ r₂ memberikan y = A·er₁t + B·er₂t (redaman lampau); punca berganda r memberikan y = (A + Bt)·ert (redaman genting); punca kompleks −α ± jω memberikan y = e−αt(A cos ωt + B sin ωt) (redaman kurang, berayun). Nilai awal y(0) dan y′(0) menetapkan A dan B. Ditulis sebagai y″ + 2ζωny′ + ωn²y = 0, ωn ialah frekuensi tabii dan ζ nisbah redaman: ζ < 1 berayun. Dengan input malar, y″ + by′ + cy = c·K, penyelesaiannya ialah nilai mantap K tambah sambutan bebas.
Kamiran tentu: ∫ₐᵇ f(x) dx ialah luas bertanda antara graf dan paksi-x dari a hingga b. Cari antiterbitan F (dengan F′ = f) dan nilaikan F(b) − F(a). Petua kuasa: ∫xⁿ dx = xⁿ⁺¹/(n + 1). Luas antara dua lengkung ialah ∫(atas − bawah) dx.
Antiterbitan
F(x) = 4·ln x
Nilaikan pada had
I = 4·(ln 6 − ln 2) = 4·ln 3 = 4.394
Jawapan
I = 4.394
4.Kamiran
Jawapan: A = 10.67
Kamiran tentu: ∫ₐᵇ f(x) dx ialah luas bertanda antara graf dan paksi-x dari a hingga b. Cari antiterbitan F (dengan F′ = f) dan nilaikan F(b) − F(a). Petua kuasa: ∫xⁿ dx = xⁿ⁺¹/(n + 1). Luas antara dua lengkung ialah ∫(atas − bawah) dx.
Sifar
4 − x² = 0 → x = ±2
Kamirkan
A = [4x − x³/3]−22 = 5.333 − (−5.333) = 10.67
Jawapan
A = 10.67
5.Terbitan
Jawapan: f′(−2) = −0.6154
Pembezaan: terbitan f′(x) ialah kecerunan graf pada x. Petua kuasa: (xⁿ)′ = n·xⁿ⁻¹. Petua hasil darab: (uv)′ = u′v + uv′. Petua rantai: (f(g(x)))′ = f′(g(x))·g′(x). Di mana f′(x) = 0, graf mempunyai tangen mengufuk: maksimum, minimum atau titik pelana.
Petua rantai untuk logaritma
f′(x) = 4x/(2x2 + 5)
Gantikan
f′(−2) = −8/13 = −0.6154
Jawapan
f′(−2) = −0.6154
6.Terbitan
Jawapan: x = 2
Pembezaan: terbitan f′(x) ialah kecerunan graf pada x. Petua kuasa: (xⁿ)′ = n·xⁿ⁻¹. Petua hasil darab: (uv)′ = u′v + uv′. Petua rantai: (f(g(x)))′ = f′(g(x))·g′(x). Di mana f′(x) = 0, graf mempunyai tangen mengufuk: maksimum, minimum atau titik pelana.
Terbitan
f′(x) = 3x2 − 9x + 6 = 3·(x2 − 3x + 2)
Selesaikan f′(x) = 0
x = 1 atau x = 2
Yang mana satu
f″(x) = 6x − 9: negatif pada maksimum, positif pada minimum.
f″(1) = −3, f″(2) = 3
x = 2
Jawapan
x = 2
7.Kaedah berangka
Jawapan: |E| = 0.0005793
Pengamiran berangka: bahagikan [a, b] kepada n jalur selebar h = (b − a)/n dengan nod x₀ … xₙ. Petua trapezium: T = (h/2)[f₀ + 2f₁ + 2f₂ + … + 2fn−1 + fₙ]. Petua Simpson (n genap): S = (h/3)[f₀ + 4f₁ + 2f₂ + 4f₃ + … + 4fn−1 + fₙ]. Ralat trapezium mengecil seperti h² (had (b − a)h²·M/12 dengan M = maks|f″|), Simpson seperti h⁴, maka Simpson biasanya jauh lebih tepat untuk nod yang sama.
Nilai tepat
I = 1.718282
Petua Simpson
S = (0.5/3)·(1 + 4·1.648721 + 2.718282) = 1.718861
Ralat
|E| = |1.718861 − 1.718282| = 0.0005793
Jawapan
|E| = 0.0005793
Ralat relatif: 0.03372%.
8.Kaedah berangka
Jawapan: |E| = 0.02844
Penyelesaian berangka bagi y′ = f(t, y): kaedah Euler melangkah di sepanjang cerun pada permulaan setiap langkah: yn+1 = yₙ + h·f(tₙ, yₙ). Kaedah Heun (Euler yang dipertingkat) mengambil purata cerun pada permulaan dan pada ramalan Euler: k₁ = f(tₙ, yₙ), k₂ = f(tₙ + h, yₙ + h·k₁), yn+1 = yₙ + (h/2)(k₁ + k₂). Kaedah Runge–Kutta klasik (RK4) menggabungkan empat cerun: yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄). Ralat global Euler berkadar dengan h, Heun dengan h², RK4 dengan h⁴.
Langkah Euler
y4 = 1.3122
Nilai tepat
y(0.4) = 1.34064
Ralat global
|E| = |1.34064 − 1.3122| = 0.02844
Graf
Jawapan
|E| = 0.02844
Membahagi dua h kira-kira membahagi dua ralat kaedah Euler.