Diferencialinės lygtys: pasiruošimo paketas

Formulių lapai, pagrindiniai kiekvienos temos metodų žingsniai ir 8 uždavinių, o pilni sprendimai pateikti gale.

Nauji uždaviniaiAtgal į kursą

Spausdinimo lange pasirinkite „Išsaugoti kaip PDF“, kad nuosavą kopiją išsaugotumėte savo įrenginyje.

1Formulės ir metodai

Diferencialinės lygtys

Eksponentinis augimas / mažėjimas
dy/dt = ky → y = y0ekt
Pusėjimo trukmė
k = ln 2 / t½
Niutono aušimo dėsnis
T(t) = Ta + (T0 − Ta)e−kt
Integruojantysis daugiklis
y′ + p(x)y = q(x): μ = e∫p dx, (μy)′ = μq
Maišymo talpa
dm/dt = r·c − (r/V)·m

Integralai

Laipsnio išvestinė
∫xⁿ dx = xⁿ⁺¹/(n + 1) + C
Apibrėžtinis integralas
∫ab f dx = F(b) − F(a)
Pagrindiniai integralai
∫eˣ dx = eˣ, ∫1/x dx = ln|x|, ∫sin x dx = −cos x
Plotas tarp kreivių
∫ab (f − g) dx

Metodas

  1. Kai įmanoma, pointegralinę funkciją perrašykite kaip paprastų narių sumą.
  2. Taikykite laipsninės funkcijos taisyklę: ∫xⁿ dx = xⁿ⁺¹/(n + 1), o ex, 1/x, sin ir cos – lentelinius integralus.
  3. Keitinį naudokite, kai matote funkciją kartu su jos išvestine.
  4. Apibrėžtiniam integralui apskaičiuokite F(b) − F(a).

Išvestinės

Laipsnio išvestinė
(xⁿ)′ = n·xⁿ⁻¹
Sandaugos taisyklė
(uv)′ = u′v + uv′
Dalmens išvestinė
(u/v)′ = (u′v − uv′)/v²
Grandinės taisyklė
f(g(x))′ = f′(g(x))·g′(x)
Pagrindinės išvestinės
(eˣ)′ = eˣ, (ln x)′ = 1/x, (sin x)′ = cos x, (cos x)′ = −sin x

Metodas

  1. Nustatykite struktūrą: suma, sandauga, dalmuo ar sudėtinė funkcija.
  2. Kiekvienam x laipsniui taikykite laipsninės funkcijos išvestinės taisyklę.
  3. Kur reikia, taikykite sandaugos, dalmens arba sudėtinės funkcijos išvestinės taisyklę.
  4. Suprastinkite ir, jei prašoma, apskaičiuokite reikšmę taške.

Skaitiniai metodai

Pusiaukirta
c = (a + b)/2, width after n steps (b − a)/2ⁿ
Niutono–Rafsono metodas
xn+1 = xₙ − f(xₙ)/f′(xₙ)
Tiesinė interpoliacija
y = y₁ + (y₂ − y₁)(x − x₁)/(x₂ − x₁)
Trapecijų formulė
T = (h/2)[f₀ + 2f₁ + … + 2fn−1 + fₙ], h = (b − a)/n
Simpsono formulė (n lyginis)
S = (h/3)[f₀ + 4f₁ + 2f₂ + … + 4fn−1 + fₙ]
Trapecijų paklaidos įvertis
|E| ≤ (b − a)h²M/12, M = max|f″|
Oileris
yn+1 = yₙ + h·f(tₙ, yₙ)
Hoinas
yn+1 = yₙ + (h/2)(k₁ + k₂), k₂ = f(tₙ + h, yₙ + hk₁)
RK4
yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄)

Metodas

  1. Šakniai: pusiaukirtai patikrinkite ženklo pokytį, Niutono metodui apskaičiuokite f ir f′, o iteracijų lentelėje palikite šešis ar daugiau skaitmenų.
  2. Lentelės reikšmei: tiesiškai interpoliuokite tarp dviejų gretimų taškų.
  3. Integralui: raskite h = (b − a)/n, išvardykite f kiekviename mazge, pritaikykite svorius (1, 2, …, 2, 1 trapecijų; 1, 4, 2, …, 4, 1 Simpsono, kai n lyginis) ir padauginkite iš h/2 arba h/3.
  4. Diferencialinei lygčiai: žengkite nuo t₀ su yₙ₊₁ = yₙ + h·f(tₙ, yₙ) (Oileris) arba vidurkintu nuolydžiu (Hoinas), po eilutę kiekvienam žingsniui, kol pasieksite prašomą t.

2Uždaviniai

1.Diferencialinės lygtysVidutinis
tyy(0) = 2y(2) = ?

Išspręskite y″ + 6y′ + 9y = 0, kai y(0) = 2 ir y′(0) = 2. Raskite y(2).

y″ + 6y′ + 9y = 0y(0) = 2y′(0) = 2
y(2) =
2.Diferencialinės lygtysSunkus
t (min)m (kg)15 minm = ?

Talpoje yra 500 L gryno vandens. Sūrymas, kurio litre yra 0,05 kg druskos, įteka 2 L/min greičiu, o gerai išmaišytas tirpalas tokiu pat greičiu išteka. Kiek druskos m bus talpoje po 15 min?

V = 500 Lm0 = 0 kgc = 0,05 kg/Lr = 2 L/min
m = kg
3.IntegralaiVidutinis
xy11,5y = x ln x

Apskaičiuokite I = ∫₁1,5 x·ln x dx.

∫₁1,5 x ln x dx
I =
4.IntegralaiSunkus
xy0,5y = x²eˣ

Apskaičiuokite I = ∫₀0,5 x²·eˣ dx.

∫₀0,5 x²eˣ dx
I =
5.IšvestinėsVidutinis
xyx = 2

Tegu f(x) = (4x + 1)4. Raskite f′(2).

f(x) = (4x + 1)4x = 2
f′(2) =
6.IšvestinėsSunkus
xy

Tegu f(x) = x3 + 9x2 + 24x + 2. Raskite lokaliojo minimumo x koordinatę x.

f(x) = x3 + 9x2 + 24x + 2
x =
7.Skaitiniai metodaiVidutinis
xy01,53

Simpsono formule, kai n = 2, apytiksliai apskaičiuokite ∫03 sin x dx ir palyginkite su tikslia reikšme. Raskite absoliučiąją paklaidą |E| = |S − I|.

f(x) = sin x[a, b] = [0, 3]n = 2
|E| =
8.Skaitiniai metodaiSunkus
00,050,100,050,1ty(0, 0)

Atlikite vieną klasikinio Rungės–Kutos metodo (RK4) žingsnį su h = 0,1, kai y′ = e−ty, y(0) = 0. Raskite y₁ ≈ y(0,1).

y′ = e−tyy(0) = 0h = 0,1
y1 =

3Tirpalai

1.Diferencialinės lygtys

Atsakymas: y(2) = 0,04462

Antros eilės tiesinės lygtys: lygčiai y″ + b·y′ + c·y = 0 bandykite y = ert: gaunama charakteristinė lygtis r² + b·r + c = 0. Dvi skirtingos realiosios šaknys r₁ ≠ r₂ duoda y = A·er₁t + B·er₂t (viršslopinamas atvejis); dviguba šaknis r duoda y = (A + Bt)·ert (kritinis slopinimas); kompleksinės šaknys −α ± jω duoda y = e−αt(A cos ωt + B sin ωt) (nepakankamas slopinimas, virpesiai). Pradinės vertės y(0) ir y′(0) nustato A ir B. Užrašius y″ + 2ζωny′ + ωn²y = 0, ωn yra savasis dažnis, o ζ – slopinimo santykis: kai ζ < 1, atsiranda virpesių. Esant pastoviam įėjimui, y″ + by′ + cy = c·K, sprendinys yra nusistovėjusi vertė K plius laisvoji dedamoji.

Dviguba šaknis

(r + 3)² = 0 → r = −3

Konstantos

A = y(0) = 2
y′(0) = B − 3A → B = 2 + 6 = 8

Apskaičiuokite

y(2) = (2 + 8·2)·e−6 = 0,04462

Atsakymas

y(2) = 0,04462
2.Diferencialinės lygtys

Atsakymas: m = 1,456 kg

Pirmos eilės diferencialinės lygtys: lygties dy/dt = k·y sprendinys yra y = y₀·ekt (augimas, kai k > 0, mažėjimas, kai k < 0); pusėjimo trukmė T duoda k = ln 2/T. Niutono aušimo dėsnis dT/dt = −k·(T − Tₐ) duoda T(t) = Tₐ + (T₀ − Tₐ)·e−kt. Tiesinė lygtis y′ + p(x)·y = q(x) sprendžiama su integruojančiuoju daugikliu μ = e∫p dx: (μy)′ = μq, todėl μy = ∫μq dx + C.

Sudarykite lygtį

dm/dt = 2·0,05 − (2/500)·m = 0,1 − 0,004·m

Išspręskite

Integruojantysis daugiklis yra ert/V; sprendinys artėja prie pusiausvyrinės reikšmės c·V.

m(t) = 25 + (−25)·e−t/250

Įstatykite t

m(15) = 25 + (−25)·0,9418 = 1,456 kg

Atsakymas

m = 1,456 kg
3.Integralai

Atsakymas: I = 0,1436

Integravimas dalimis: ∫u·dv = u·v − ∫v·du, išplaukia iš sandaugos diferencijavimo taisyklės. Kaip u rinkitės dauginamąjį, kuris diferencijuojant tampa paprastesnis (x laipsnį ar ln x), o kaip dv – dauginamąjį, kurį galite suintegruoti (eˣ, sin x, cos x arba x, kai u = ln x). Apibrėžtiniam integralui u·v apskaičiuokite abiejose ribose: ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. Dauginamajam x² reikia dviejų ciklų.

Pasirinkite u ir dv

u = ln x, du = dx/x; dv = x dx, v = x²/2

Pirmykštė funkcija

∫x ln x dx = (x²/2)ln x − ∫(x/2) dx = (x²/2)ln x − x²/4

Ribos

(1,125·ln 1,5 − 0,5625) − (0 − 0,25) = 0,1436

Atsakymas

I = 0,1436
4.Integralai

Atsakymas: I = 0,0609

Integravimas dalimis: ∫u·dv = u·v − ∫v·du, išplaukia iš sandaugos diferencijavimo taisyklės. Kaip u rinkitės dauginamąjį, kuris diferencijuojant tampa paprastesnis (x laipsnį ar ln x), o kaip dv – dauginamąjį, kurį galite suintegruoti (eˣ, sin x, cos x arba x, kai u = ln x). Apibrėžtiniam integralui u·v apskaičiuokite abiejose ribose: ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. Dauginamajam x² reikia dviejų ciklų.

Pirmasis ciklas

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

Antrasis ciklas

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

Ribos

1,649·(1,25) − 1·2 = 0,0609

Atsakymas

I = 0,0609
5.Išvestinės

Atsakymas: f′(2) = 11664

Diferencijavimas: išvestinė f′(x) – tai grafiko krypties koeficientas taške x. Laipsnio taisyklė: (xⁿ)′ = n·xⁿ⁻¹. Sandaugos taisyklė: (uv)′ = u′v + uv′. Grandinės taisyklė: (f(g(x)))′ = f′(g(x))·g′(x). Kur f′(x) = 0, grafikas turi horizontalią liestinę: maksimumą, minimumą arba balno tašką.

Grandinės taisyklė

f′(x) = 4·(4x + 1)3·4

Įstatykite

f′(2) = 4·93·4 = 11664

Atsakymas

f′(2) = 11664
6.Išvestinės

Atsakymas: x = −2

Diferencijavimas: išvestinė f′(x) – tai grafiko krypties koeficientas taške x. Laipsnio taisyklė: (xⁿ)′ = n·xⁿ⁻¹. Sandaugos taisyklė: (uv)′ = u′v + uv′. Grandinės taisyklė: (f(g(x)))′ = f′(g(x))·g′(x). Kur f′(x) = 0, grafikas turi horizontalią liestinę: maksimumą, minimumą arba balno tašką.

Išvestinė

f′(x) = 3x2 + 18x + 24 = 3·(x2 + 6x + 8)

Išspręskite f′(x) = 0

x = −4 arba x = −2

Kuris yra kuris

f″(x) = 6x + 18: neigiama maksimume, teigiama minimume.

f″(−4) = −6, f″(−2) = 6
x = −2

Atsakymas

x = −2
7.Skaitiniai metodai

Atsakymas: |E| = 0,07556

Skaitinis integravimas: padalykite [a, b] į n juostų, kurių plotis h = (b − a)/n, su mazgais x₀ … xₙ. Trapecijų formulė: T = (h/2)[f₀ + 2f₁ + 2f₂ + … + 2fn−1 + fₙ]. Simpsono formulė (n lyginis): S = (h/3)[f₀ + 4f₁ + 2f₂ + 4f₃ + … + 4fn−1 + fₙ]. Trapecijų paklaida mažėja kaip h² (įvertis (b − a)h²·M/12, kur M = max|f″|), Simpsono – kaip h⁴, todėl su tais pačiais mazgais Simpsono formulė paprastai yra gerokai tikslesnė.

Tiksli reikšmė

I = 1,989992

Simpsono formulė

S = (1,5/3)·(0 + 4·0,997495 + 0,14112) = 2,06555

Paklaida

|E| = |2,06555 − 1,989992| = 0,07556

Atsakymas

|E| = 0,07556

Santykinė paklaida: 3,797 %.

8.Skaitiniai metodai

Atsakymas: y1 = 0,09967

Skaitinis y′ = f(t, y) sprendimas: Oilerio metodas žengia išilgai nuolydžio kiekvieno žingsnio pradžioje: yn+1 = yₙ + h·f(tₙ, yₙ). Hoino metodas (patobulintas Oilerio) vidurkina nuolydį pradžioje ir Oilerio prognozės taške: k₁ = f(tₙ, yₙ), k₂ = f(tₙ + h, yₙ + h·k₁), yn+1 = yₙ + (h/2)(k₁ + k₂). Klasikinis Rungės–Kutos metodas (RK4) sujungia keturis nuolydžius: yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄). Oilerio globalioji paklaida proporcinga h, Hoino – h², RK4 – h⁴.

Keturi nuolydžiai

k₁ = f(0, 0) = 1
k₂ = f(0,05, 0,05) = 0,997503
k₃ = f(0,05, 0,0498752) = 0,997509
k₄ = f(0,1, 0,0997509) = 0,990074

Sujunkite

y₁ = 0 + (0,1/6)(1 + 2·0,997503 + 2·0,997509 + 0,990074) = 0,09967

Atsakymas

y1 = 0,09967

ExamProblems · examproblems.com/courses/differential-equations