logo

Mapleda differensial tenglamalarni yechish funksiyalari

Загружено в:

08.08.2023

Скачано:

0

Размер:

101.95703125 KB
Mavzu: Mapleda differensial tenglamalarni yechish funksiyalari.
1. Nazariy qism
2. Berilgan topshiriqning bajarilish 
qismi
3.  2.1. Laboratoriya topshirig’i 
mazmuni
2.2.   Mustaqil ishning bajarilishi
4. Xulosa
5. Adabiyotlar
                                                                                                                 1 . Nazariy qism
Chiziqli differensial tenglamalar.
.Maple   tizimida,   DEtools   paketining   subsetini   ishlatgan   holatda   differensial
operatorlar bilan ishlash mumkin. Differensial operatorlar bu holda ko'phad obyektlar
bo'lib quyidagi shaklga ega bo'ladi:
L  :=  a  (  x  )  DF (
  n
  )
    +
.  .  . +
a( x )	
DF		a (  x  )
N
1 0
a (  x  ) koeffisiyentlari  soha  ustidagi ratsional funksiya  hisoblanadi.  Bu  yerda   D	
i
DF(  x  )  1    va	
DF( 	u v	 )	u DF( 	v )	DF( 	u ) 	v
. kabi xossalarni  qonotlantiradigan
obyekt hisoblanadi . Bu operatorlar ustida ko'paytirish, simmetrik ko'paytmalar hosil
qilish, bir tomonli eng katta umumiy bo'luvchilarni topish, ko'paytuvchilarga ajratish,
va boshqa funksiyalarni bajarish mumkin.
Mark van Xoyeyjga (Nijmegen Universiteti) xos bu funksiya chiziqli differensial
tenglamalarning yopiq shakldagi natijalarini topishda kelajakdagi rivojlanishlarga imkon
yaratadi. Quyida DEtools ning subpaketidagi amaliyotlar keltirilgan.
DEtools paketida biz with komandasidan komandalarning qisqa shaklini ishlatishga
imkon yaratishi uchun foydalanamiz.
> restart; > with(DEtools):
Differensial operatorlar bilan algebraik amallar bajarish . Differensial operatorlarni qo'shish, ko'paytirish, o'rniga qo'yish va boshqalarga 
o'xshash o'rin almashtirib bo'lmaydi (noncommutative domain) . C(x)[DF]dagi L
differensial operatori buL := a 	( x 	) DF	( n ) + . . . +
a	( x )	
DF		a( x ), bo'lib , bu yerda
n
1 0
a (  x  ) C(x)ning elementlari.  C(x)[Dx] dagi L elementi L( y(x) )=0 chiziqli bir xil	
i
differensial tenglamaga to'g'ri keladi.
C(x)[Dx] doirasida ko'paytirish differensial operatorlarni shakllantirishga 
to'g'ri keladi. Shunday qilib agar L = mult(f,g) bo'lsa, u holda L( y(x) ) = 
f(g( y(x) )). Xususan, mult(DF,x) = x*DF + 1.
Misol tariqasida qanday algebraik amallarni bajarish mumkinligini ko'rish uchun
quyidagi differensial operatorlarga e'tibor bering:
> L1 := x^2*DF^2 - x*DF + (a-x^2);
L1  :=  x 2
 DF 2
 x DF  a  x 2
> L2 := x*DF - (x^2-b);
L2  :=  x DF  x 2
 b
> L3 := DF^2 - x;
L3  :=  DF 2
 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  :=  x 3
 DF 3
 (    x 4
 x 2
 x 2
 b  )  DF 2
 (  a x  x b  x  4  x 3
  )  DF  a x 2
 x 4
 x 2
 b  a
b
> L5 := collect( mult( L2, L1, [DF, x] ), DF );
L5  :=
x 3
 DF 3
 (    x 4
 x 2
 x 2
 b  )  DF 2
 (    x b  x  a x  )  DF  a b  x 4
 2  x 2
 a x 2
 x 2
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] ); va
> L7 := GCRD( L4, L6, [DF, x] );
L7  :=  DF    x
2 
x b
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 '+ y cos x =sin x cos x  differensial tenglamaning umumiy yechimini topish.
> restart;
> de:=diff(y(x),x)+y(x)*cos(x)=sin(x)*cos(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 ''  2 y '+ y =sin x + 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);
>  dsolve(deq,y(x));
y(  x )    _  C1e   x
    _  C2e   x
  x      1
2  cos( x )     1
4   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.  
>  restart; de:=diff(y(x),x$2)+k^2*y(x)=sin(q*x);
 := de 


 



d d 2
x 2 ( )y x k 2
( )y x ( )sin q x
Endi  yechimni  rezonans  holatda  izlaymiz.  Buning  uchun  esa  dsolve  komandani 
chaqirishdan oldin q=k deb olish kerak.
                  >  #Quyidagi Koshi masalasini yechimini toping  Xursanov Akrom 316-
guruh#
>  de:=diff(y(x),x$3)-diff(y(x),x)=tan(x);
 := de 


 



d d3
x
3 ( )y x 

 


d d
x ( )y x ( )tan x
>  cond:=y(0)=0,D(y)(0)=-1,(D@@2)(y)(0)=1;	
 := 	cond	,	,		(	)	y	0	0		(	)	(	)	D	y	0	-1		(	)	(	)	(	)	D()2	y	0	1
>  dsolve({de,cond},y(x));
( )y x d





	
0
x   e	
(	)	_ z1 1
2 d

	
0
_ z1 e	(	)	_ z1 ( )tan _z1 _z1 e	_ z1 1
2 e	(	)	_ z1 d

	
0
_ z1 e	_ z1 ( )tan _z1 _z1 _z1
                                            Foydalanilgan adabiyotlar
1. Матросов А. Решение задач математики и механики в среде Maple 6. СПб.:
Питер, 2000.
2. В.З. АЛАДЬЕВ. Основы программирования в Maple. Таллинн, 2006.
3. Основы использования математического пакета Maple в 
моделировании:
Учебное пособие / Международный институт компьютерных 
технологий.
Липецк, 2006. 119с.
4.  Дьяконов В. Maple 6. Учебный ,рс СПб.: Питер, 2001.
5. Аладьев В.З., Лиопо В.А., Никитин А.В. Математический пакет Maple в
физическом   моделировании.-   Гродно:   Гродненский   госу-дарственный
университет им. Янки Купалы, 2002, 416 с.
6. O’runbayev E., Murodov F. Kompyuter algebrasi tizimlarining 
amaliy tadbiqlari. –SamDU nashri – Samarqand, 2003, 96 s.
7. Аладьев В.З., Богдявичюс М.А .   Решение физико-технических и  
математических задач с пакетом  Maple V .- Вильнюс: Изд-во 
Техника , 1999, 686 c., ISBN 9986-05-398-6.
8. Аладьев В.З., Богдявичюс М.А .  Maple 6:  Решение математичес-ких, статистических и инженерно-физических задач.- Москва: 
Лаборатория   Базовых Знаний , 2001, 850   с. + CD-ROM, ISBN 5-
93308-085-X.
9. Математика на компьютере: Maple 8. — М.: СОЛОН-Пресс, 2003.176 с:

Mavzu: Mapleda differensial tenglamalarni yechish funksiyalari. 1. Nazariy qism 2. Berilgan topshiriqning bajarilish qismi 3. 2.1. Laboratoriya topshirig’i mazmuni 2.2. Mustaqil ishning bajarilishi 4. Xulosa 5. Adabiyotlar

1 . Nazariy qism Chiziqli differensial tenglamalar. .Maple tizimida, DEtools paketining subsetini ishlatgan holatda differensial operatorlar bilan ishlash mumkin. Differensial operatorlar bu holda ko'phad obyektlar bo'lib quyidagi shaklga ega bo'ladi: L := a ( x ) DF ( n ) + . . . + a( x ) DF  a ( x ) N 1 0 a ( x ) koeffisiyentlari soha ustidagi ratsional funksiya hisoblanadi. Bu yerda D i DF( x )  1 va DF( u v ) u DF( v ) DF( u ) v . kabi xossalarni qonotlantiradigan obyekt hisoblanadi . Bu operatorlar ustida ko'paytirish, simmetrik ko'paytmalar hosil qilish, bir tomonli eng katta umumiy bo'luvchilarni topish, ko'paytuvchilarga ajratish, va boshqa funksiyalarni bajarish mumkin. Mark van Xoyeyjga (Nijmegen Universiteti) xos bu funksiya chiziqli differensial tenglamalarning yopiq shakldagi natijalarini topishda kelajakdagi rivojlanishlarga imkon yaratadi. Quyida DEtools ning subpaketidagi amaliyotlar keltirilgan. DEtools paketida biz with komandasidan komandalarning qisqa shaklini ishlatishga imkon yaratishi uchun foydalanamiz. > restart;

> with(DEtools): Differensial operatorlar bilan algebraik amallar bajarish .

Differensial operatorlarni qo'shish, ko'paytirish, o'rniga qo'yish va boshqalarga o'xshash o'rin almashtirib bo'lmaydi (noncommutative domain) . C(x)[DF]dagi L differensial operatori buL := a ( x ) DF ( n ) + . . . + a ( x ) DF  a( x ), bo'lib , bu yerda n 1 0 a ( x ) C(x)ning elementlari. C(x)[Dx] dagi L elementi L( y(x) )=0 chiziqli bir xil i differensial tenglamaga to'g'ri keladi. C(x)[Dx] doirasida ko'paytirish differensial operatorlarni shakllantirishga to'g'ri keladi. Shunday qilib agar L = mult(f,g) bo'lsa, u holda L( y(x) ) = f(g( y(x) )). Xususan, mult(DF,x) = x*DF + 1. Misol tariqasida qanday algebraik amallarni bajarish mumkinligini ko'rish uchun quyidagi differensial operatorlarga e'tibor bering: > L1 := x^2*DF^2 - x*DF + (a-x^2); L1 := x 2 DF 2  x DF  a  x 2 > L2 := x*DF - (x^2-b); L2 := x DF  x 2  b > L3 := DF^2 - x; L3 := DF 2  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 := x 3 DF 3  (  x 4  x 2  x 2 b ) DF 2  ( a x  x b  x  4 x 3 ) DF  a x 2  x 4  x 2 b  a b > L5 := collect( mult( L2, L1, [DF, x] ), DF ); L5 := x 3 DF 3  (  x 4  x 2  x 2 b ) DF 2  (  x b  x  a x ) DF  a b  x 4  2 x 2  a x 2  x 2 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] );