L2 := x DF - x2 + b
> L3 := DFA2 - x;
L3 := DF2 - x
Biz bu operatorlari ko'paytirib shuni eslatib o'tishimiz mumkinki, ularni o'rnini almashtirib bo'lmaydi:
L4 := collect( mult( L1, L2, [DF, x] ), DF );
L4 := x3 DF3 + ( -x4 + x2 + x2 b) DF2 + ( ax - x b - x - 4 x3) DF - ax2 + x4 - x2 b + a b
L5 := collect( mult( L2, L1, [DF, x] ), DF );
L5 : =
x3 DF + ( -x4 + x2 + x2 b ) DF2 + ( -x b - x + ax) DF + a b + x4 - 2 x2 - ax2 - x2 b
Argument [DF, x] mult komandasiga ko'paytirish DF va x belgilab bergan differensial sohasi ustida ekanligini ko'rsatadi, shuning uchun a va b o'zgaruvchilari o'zgarmas songa teng. Bu yana _Envdiffopdomain := [DF, x], muhit o'zrgaruvchisi orqali belgilanishi mumkin.
Bir tomonli eng kichik umumiy ko'paytuvchi va eng katta umumiy bo'luvchi tushunchasi shunday sohalarda mavjud bo'ladi. Misol tariqasida,
L6 := LCLM( L3, L2, [DF, x] );
(-b3 - x3 - 2 b - 3 x4 b + 3 x2 b2 + 3 x4 - 3 b2 + x3 b + x6 - x5) DF2
x ( b + b2 + x2 - 2 x2 b - x3 + x4 )
L6 := DF3
-b3 +
x DF
3 b - 3 x4 b + 3 x2 b2 + 2 x4 - 4 b2 + x3 b + x6 - x5 - x2 + 2 x2 b b + b2 + x2 - 2 x2 b - x3 + x4
va
> L7 := GCRD( L4, L6, [DF, x] );
L7 := DF -
x2 - b
x
Buni o'ng yoki chap tomonni ko'paytirishi orqali tekshirish mumkin. Bizning
masalamizda biz quyidagiga egamiz:
> rightdivision(L6,L7, [DF,x] );
" . (-x3 - 2 b + 2 x4 - 2 b2) DF
DF —
_ x ( b + b2 + x2 - 2 x2 b - x3 + x4 )
0
2 x3 + 4 xb - 6 b - x4 - 2 x2 + xb2 - 2 x3 b + x5
b + b2 + x2 - 2 x2 b - x3 + x4
> rightdivision(L4,L7, [DF,x] );
[ x3 DF2 + x2 DF - x3 + ax - x, 0]
Ikki holatda ham yuqoridagi amal 0 qoldiqli bo'linmani beradi. Hisobni to'g'riligi quyidagi ko'paytirish orqali tekshirish mumkin:
> collect( L4 - mult( %[1], L7, [DF,x] ), DF, normal );
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.
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)
iSx )
> 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+e~x 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:
Ml ^ ( Pj Л
^ y(x) -1|— y(x)|+ y(x) = sin(x) + e(-x)
ydx ) V«x )
> dsolve(deq,y(x));
y( x) = _ C1ex + _ C2exx + 1cos(x) +1 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 q±k va q=k (rezonans) ikki holda umumiy yechimini topish.
de:=
d ( л
2y( x)
dx 2
+ к 2y( x) = sin( qx)
> restart; de:=diff(y(x),x$2)+kA2*y(x)=sin(q*x);
> dsolve(deq,y(x));
cos((k + q)x) 1 cos((k - q)x) .
У( x) =
к + q 2 к - q
sin((k - q)x)
к - q
sin^ + q) x)
к + q к
cos (к») +
C1 sin^x) +
C2 cos(кx)
к
Endi yechimni rezonans holatda izlaymiz. Buning uchun esa dsolve komandani chaqirishdan oldin q=k deb olish kerak.
2
1 cos^x) sin(kx)
x =-2 к2 C1 sin(kx) + _ C2 cos(kx)
- ^cos^^m^x) +1 kx \cos(kx)
к
2
+
> q:=k: dsolve(de,y(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.
Yechimning fundamental (bazisli) sistema ko’rimshi.
dsolve komanda differensial tenglama yechimining fundamen-tal sistemali (bazisli funksiyalar) ko’rinishini topish imkonini beradi. Buning uchun dsolve komandaning parametrlarida output=basis ni ko’rsatish kerak.
Quyidagi berilgan differensial tenglamaning fundamental sistema yechimini toping: y(4)+2y"+y=0.
> de:=diff(y(x),x$4)+2*diff(y(x),x$2)+y(x)=0;
de: =
Л ( д2 + 2
-^y( x)
dx 2
+ y( x) = 0
> dsolve(de, y(x), output=basis);
[cos(x), sin( x), x cos(x), x sin( x)]
Koshi masalasi yoki chegaraviy masalaning yechilishi. dsolve komanda Koshi masalasi yoki chegaraviy masalaning yechimini topishi mumkin, agarda berilgan differensial tenglama uchun noaniq funksiyaning boshlang’ich hamda chegaraviy shartlari berilsa. Boshlang’ich yoki chegaraviy shartlarda hosilalarni belgi-lash uchun differensial operator D ishlatiladi, masalan, y''(0)=2 shartni (D@@2)(y)(0) = 2 kabi berishga to’g’ri keladi yoki y'(1)=0 shart-ni: D(y)(i) = 0. Eslatib o’tamiz, n-chi tartibli hosila (D@@n)(y) kabi yoziladi.
Koshi masalasining yechimini topish: y(4)+y"=2cosv, y(0)=-2, y'(0)=1, y"(0)=0, ym(0)=0.
de :=
( д4 / \
дй *x)
+
д f ч -^y( x)
дт2
= 2cos(x)
cond:=y(0)=-2, D(y)(0)=1, (D@@2)(y)(0)=0,
(D@@3)(y)(0)=0;
cond:=y(0)=-2, D(y)(0)=1, (D(2))(y)(0)=0, (D(3))(y)(0)=0
y(v)=-2cos(.x)-.xsin(.x)+.x
Quyidagi chegaraviy masalaning yechimini topamiz: y"+y = 2x -n , y(0) = о , y\n 1 = о. Yechim grafigini yasang.
> restart; de:=diff(y(x),x$2)+y(x)=2*x-Pi;
de:
f д2 л
-^y( x)
дх2 \дл
+ y( x) = 2x -n
> cond:=y(0)=0,y(Pi/2)=0;
( n
cond := y(0) = 0, y\ — | = 0
> dsolve({de,cond},y(x));
y(.x)=2.x-n+ncos(.x)
> de:=diff(y(x),x$4)+diff(y(x),x$2)=2*cos(x);
Berilgan topshiriqning bajarilish qismi
misol: Quyidagi differensial tenglamaning umumiy yechimini toping:
y(4)+2ym-2y"=2cosx
Yechish:
>with(DETools):
DifTeng:=diff(y(x),x$4)+3*diff(y(x),x$3)-
2 cos(x)
>dsolve(DifTeng,y(x));
y( x) =
2
- + -
2
3 V17 2 2
3 у 17 2 2
+ - cos( x)
1
3
sin( x ) + C3 x + C4
2*diff(y(x),x$2)=2*cos(x);
d4 , '
|
|
d3 '
|
|
f d2
|
^y( x ) j
|
+ 3
|
d?y( x),
|
- 2
|
v dx2
|
DifTeng :=
;y( x )
(-3/2 x + 1/2 x-(rf) ^ (-3/2 x - 1/2 x/l7)
C2 e C1 e 1
misol: Quyidagi Koshi masalasini yechimini aniqlang:
y(4)+2yy'"-2y"=2cosx, y(0)=-2, y’(0)=1, y"(0)=0, ym(0)=0.
Yechish:
> Order:=5;
Order := 5
>
dsolve({diff(y(x),x)=y(x)+x*exp(y(x)),y(0)=0},y(x),type=s
eries);
y( x ) = - x2 + - x3 + - x4 + O( x5)
misol: Quyidagi Koshi masalasini sonli yechimining grafigini [-2; 3] oraliqda DEplot komandasidan foydalanib chizing:
>with(DETools):
>DifTeng:=diff(y(x),x$3)+3*x*diff(y(x),x$2)- diff(y(x),x)+5*x=0;
DifTeng:
(dd t Л . f d2 Л
з y( x ) + 3 x 2 y( x )
l dx ) l dx )
dxy( x)
+ 5 x = 0
/W- У+5*=0, у(0)=1, у'(0)=-4, у" (0)=2.
>Bosh_shart:=[[y(0)=1,D(y)(0)=-4,D[1$2](y)(0)=2]];
Bosh_shart:= [[y(0) = 1, D(y)(0) = -4, (D(2))(y)(0) = 2]]
>DEplot( DifTeng, y(x), x=-2..3, Bosh shart );
2-misol.
Do'stlaringiz bilan baham: |