yj+1 = yj + h(k1,j + 2k2,j + 2k3,j + k4,j)
bu yerda
k1, j
f (t j , yi ),
k2, j
k3, j
f (t j
f (t j
h , y
2 j
h , y
2 j
k1, j ),
k2, j ),
k4, j
f (t j
h, y j
hk3, j ),
Runge-Kutta usulining qo‟llanilish dasturi matematik paketda ko‟zda tutilgan, shuning uchun undan dasturda ko‟rsatilgan tartibda foydalanilgan. Runge-Kutta usulining blok-sxemasi 1.1-rasmda tasvirlangan.
Maple da differensial tenglamalarning analitik yechimlarini topish uchun quyidagi komanda ishlitiladi:
dsolve(eq,var,options),
bu yerda
eq – differensial tenglama; var – noaniq funksiyalar; options – parametrlar.
Parametrlar masalaning yechilish metodini ko‟rsatishi mumkin, masalan, jimlik qoidasi bo‟yicha analitik yechim quyidagicha izlanadi: type=exact.
1.1-rasm. Runge-Kutta usulining blok-sxemasi.
Differensial tenglamani kiritishda hosilani bildirish uchun diff komanda ishlatiladi, masalan,
y''+y=x
differensial tenglama quyidagi ko‟rinishda yoziladi:
diff(y(x),x$2)+y(x)=x.
Differensial tenglamalarning umumiy yechimi soni differensial tenglamaning tartibiga bog‟liq bo‟lgan ixriyoriy o‟zgarmaslarga bog‟liqdir. Maple da bunday o‟zgarmaslar qoida bo‟yicha _C1, _C2, va h.k.lar bilan belgilanadi.
Bir jinsli bo‟lmagan chiziqli differensial tenglamaning umumiy yechimi hamma vaqt shunday chiqariladiki, ushbu yechimning strukturasi aniq ko‟rinadi. Shu bilan birga bir jinsli bo‟lmagan chiziqli differensial tenglamaning umumiy yechimi unga mos keluvchi bir jinsli differensial tenglamaning umumiy yechim- lari yig‟indisiga hamda berilgan bir jinsli bo‟lmagan differensial tenglamaning xususiy yechimiga teng. Shuning uchun ham bir jinsli bo‟lmagan chiziqli differensial tenglamaning yechimini chiqarish satri hamma vaqt ixtiyoriy o‟zgarmaslarni o‟z ichiga olgan qo‟shiluvchilardan iborat (bu mos keluvchi differensial tenglamaning umumiy yechimi) va ixtiyoriy o‟zgarmaslarsiz bo‟lgan yig‟indidan iborat (bu bir turli bo‟lmagan differensial tenglamaning xususiy yechimi) bo‟lishi mumkin.
dsolve komanda differensial tenglamaning yechimini hisoblanmaydigan shaklda beradi. Hosil bo‟lgan yechim ustidan keyinchalik ishlash uchun (masalan, yechim grafigini yasash) hosil bo‟lgan yechimning o‟ng tomonini rhs(%) komanda bilan ajratish kerak.
Muammoni oydinlashtirishni mashqlarda bajarib ko‟raylik va quyidagi tadbiqlarni bajaraylik :
misol. Ushbu
y'+ycosx=sinxcosx
birinchi tartibli differensial tenglamaning umumiy yechimini toping.
Yechish:
restart;
de:=diff(y(x),x)+y(x)*cos(x)=sin(x)*cos(x);
de:=
y(x)
y(x) cos(x)
sin( x) cos(x)
y( x)
sin( x) 1
e( sin( x)) _ C 1
Demak, izlanayotgan tenglamaning yechim funksiyasi
y(x)
sin( x) 1
e( sin( x)) _ C 1.
Eslatma: Maple da differensial tenglamaning yechimini chiqarish satrida ixtiyoriy konstanta _C1 kabi belgilanadi.
misol. Ushbu
y'' 2 y'+ y=sin x+ e x
ikkinchi tartibli differensial tenglamaning umumiy yechimini toping.
Yechish:
restart;
deq:=diff(y(x),x$2)-2*diff(y(x),x)+y(x)
=sin(x)+exp(-x);
deq:=
y( x)
2 y( x) x
y( x)
sin( x)
e( x)
y( x)
_ C1ex
_ C2ex x
1 cos(x)
2
1 e( x)
4
Eslatma: berilgan tenglama ikkinchi tartibli bo‟lganligi sababli olingan natijada ikkita ixtiyoriy konstantalar mavjud, ular Maple da _ C1 va _ C2 kabi belgilanadi. Yechimda birinchi ikkita qo‟shiluvchilar berilgan bir jinsli differensial tenglamaning umumiy yechimi, qolgan ikkitasi esa bir jinsli bo‟lmagan diffe- rensial tenglamaning xususiy yechimidir.
misol. Ushbu
y''+ k2y=sin( qx) tartibda berilgan differensial tenglamaning q umumiy yechimini toping.
Yechish:
restart; de:=diff(y(x),x$2)+k^2*y(x)=sin(q*x);
va q= k (rezonans) ikki holda
de:=
y( x)
k 2y( x)
sin(
qx)
1 cos(( k
q) x)
1 cos(( k
q) x)
sin( kx)
y( x)
2 k q
2 k q k
1 sin(( k q) x)
sin(( k q)x)
cos( kx)
2 k q
2 k q k
_ C1sin( kx)
_ C2 cos(kx)
Endi yechimni rezonans holatda izlaymiz. Buning uchun esa dsolve
komandani chaqirishdan oldin q= k deb olish kerak.
y( x)
_ C1sin( kx)
_ C2 cos( kx)
Eslatma: bu ikki holda ham bir jinsli bo‟lmagan differensial tenglamaning ixtiyoriy o‟zgarmaslarni o‟z ichiga olgan xususiy hamda umumiy yechimlar alohida qo‟shiluvchilar ko‟rinishida chiqariladi.
Dostları ilə paylaş: |