Differensial tenglamalarning umumiy yechimi.
Maple da differensial tenglamalarning analitik yechimlarini topish uchun
quyidagi komanda ishlitiladi:
dsolve(eq,var,options),
bu yerda eq – differensial tenglama, var – noaniq funkslar, options – parametrlar.
Parametrlar masalaning yechilish metodini ko’rsatishi mumkin, masalan, jimlik qoidasi
bo’yicha analitik yechim quyidagicha izlanadi: type=exact. 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’zgar-maslarga bog’liqdir. Maple da bunday
o’zgarmaslar qoida bo’yicha _S1, _S2, va h.k.lar bilan belgilanadi.
L6 DF3 (b3x32 b3 x4 b3 x2 b23 x43 b2x3 bx6x5 ) DF2
:= x (bb2x22 x2 bx3x4 ) x DF
b33 b3 x4 b3 x2 b22 x44 b2x3 bx6x5x22 x2 b
bb2x22 x2 bx3x4
L7 := DFx2b
x
DF2 (x32 b2 x42 b2 ) DF
x (bb2x22 x2 bx3x4 )
2 x34 x b6 bx42 x2x b22 x3 bx5
bb2x22 x2 bx3x4 , 0
[x3 DF2x2 DFx3a xx, 0]
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’lma-gan chiziqli differensial tenglamaning umumiy yechimi unga
mos keluvchi bir jinsli differensial tenglamaning umumiy yechim-lari yig’indisiga
hamda berilgan bir jinsli bo’lmagan diffe-rensial tenglamaning xususiy yechimiga teng.
Shuning uchun ham bir jinsli bo’lmagan chiziqli differensial tenglamaning yechi-mini
chiqarish satri hamma vaqt ixtiyoriy o’zgarmaslarni o’z ichi-ga 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 teng-lamaning xususiy yechimi) bo’lishi mumkin.
dsolve komanda differensial tenglamaning yechimini hi-soblanmaydigan shaklda
beradi. Hosil bo’lgan yechim ustidan ke-yinchalik ishlash uchun (masalan, yechim
grafigini yasash) hosil bo’lgan yechimning o’ng tomonini rhs(%)komanda bilan
ajratish kerak.
y'+ycosx=sinxcosx differensial tenglamaning umumiy yechimini topish.
> 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)
x
> dsolve(de,y(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 _S1 kabi belgilanadi.
y''2y'+y=sinx+ex ikkinchi tartibli differensial tenglamaning umumiy yechimini
toping.
> restart;
> deq:=diff(y(x),x$2)-2*diff(y(x),x)+y(x)
=sin(x)+exp(-x);
deq:= 2 ( )
2
y( ) 2 y(x) y(x) sin(x) e x
x
x
x
> dsolve(deq,y(x));
( )
1 4
cos( )
1 2
y(x) _ C1e x _ C2e x x x e x
Eslatma: berilgan tenglama ikkinchi tartibli bo’lganligi sa-babli olingan natijada
ikkita ixtiyoriy konstantalar mavjud, ular Maple da _S1 i _S2 kabi balgilanadi.
Yechimda birinchi ikkita qo’shiluvchilar berilgan bir jinsli differensial tenglamaning
umumiy yechimi, qolgan ikkitasi esa bir jinsli bo’lmagan differensial tenglamaning
xususiy yechimidir.
y''+k2y=sin(qx) tartibda berilgan differensial tenglamaning qk va q=k (rezonans)
ikki holda umumiy yechimini topish.
> restart; de:=diff(y(x),x$2)+k^2*y(x)=sin(q*x);
de:= y( ) 2y( ) sin( )
2
2
x k x qx
x
> dsolve(deq,y(x));
_ sin( ) _ cos( )
sin(( ) ) cos( )
1 2
sin(( ) )
1 2
cos(( ) ) sin( )
1 2
cos(( ) )
1 2
y( )
C1 kx C2 kx
k
kx
k q
k q x
k q
k q x
k
kx
k q
k q x
k q
k q x
x
Endi yechimni rezonans holatda izlaymiz. Buning uchun esa dsolve komandani
chaqirishdan oldin q=k deb olish kerak.
> q:=k: dsolve(de,y(x));
_ sin( ) _ cos( )
cos( )
1 2
cos( )sin( )
1 2
cos( ) sin( )
1 2
y( )
2 2
2
C1 kx C2 kx
k
kx kx kx kx
k
kx kx
x
Eslatma: bu ikki holda ham bir jinsli bo’lmagan differen-sial tenglamaning
ixtiyoriy o’zgarmaslarni o’z ichiga olgan xususiy hamda umumiy yechimlar alohida
qo’shiluvchilar ko’rinishida chiqa-riladi.
Do'stlaringiz bilan baham: |