Differensial tenglamalarni yechish bo`yicha
MatLab dasturining funksiyalari
Differensial tenglamalar va ularning sistemalarini yechish uchun MATLAB paketida quyidagi funksiyalar tashkil qilingan:
Bu funksiyalarning kirish parametrlari:
- vektor funksiya bo`lib, tenglamani hisoblash uchun qo`llanilgan;
- boshlang’ich shart vektori;
- ikkita sondan iborat massiv bo`lib, differensial tenglama yoki sistemaning integrallash intervalini aniqlaydi;
- differensial tenglama yoki ularning sistemalarini yechishning borishini boshqarish parametri.
Barcha funksiyalar quyidagi natijalar chiqaradi:
T massiv – yechim izlanayotgan to`rning koordinatalari.
X matritsa – i – ustuni yechim vektorining Ti bo`lakdagi qiymati.
funksiyada to`rtinchi-beshinchi tartibli Runge-Kutta usuli, da ikkinchi – uchinchi tartibli Runge-Kutta usuli, funksiyasida esa Adams usuli kiritilgan.
Qattiq sistemalarni yechishga mo`ljallangan funsiyalar , ya’ni bu funksiyada Gir usuli kiritilgan. Rozenbrok usuli funksiyasida, qattiq sistemaning yanada yuqori aniqlikdagi yechimini olish uchun funksiyasini qo`llash mumkin.
5. MatLab dasturida xususiy hosilali differensial tenglamalarni yechish
Differensial tenglamalar sistemasining “qattiq sistema” bo`lish ta’rifini keltiramiz. - tartibli differensial tenglamalar sistemasi
(5.1)
qattiq sistema deyiladi [7], agar quyidagi shart o`rinli bo`lsa:
B matrisa barcha xos sonlarining haqiqiy qismi musbat bo`lsa:
Sistemaning qattiqlik soni deb ataluvchi son, katta bo`lsa.
(5.2)
, ,
(5.2) chiziqsiz differensial tenglamalar sistemasini qattiqlikka tekshirishda matrisa rolida xususiy hosilalar matrisasi ishlatiladi.
Uncha katta bo`lmagan qattiqlik soni bilan berilgan sistemalarni yechish uchun ode23t, shunga o`xshash sistemalarni baholash uchun ode23tb, funksiyalari xizmat qiladi.
Bu funksiyalarning qo`llanilishini aniq misollarda ko`ramiz.
5.1.-masala. Quyidagi chegaraviy masalani [2,25; 2] intervalda yeching:
(5.3)
MATLAB funksiyalaridan foydalanish mumkin bo`lishi uchun tenglamani sistemaga keltiramiz. Buning uchun almashtirish bajaramiz va
(5.4)
tenglamalar sistemasiga ega bo`lamiz. Sistema uchun quyidagi
(5.5)
boshlang`ich shart o`rinli bo`lsin.
(5.4) sistemani hisoblash funksiyasini tuzamiz (2.8-listing). 2.9-listing da (5.4) tenglamani ode45 funksiyasi yordamida yechish tasvirlangan, yechim grafigi 32-rasmda keltirilgan.
5.1-listing.
function F=FF(t,x)
F=[-4*x(1)-13*x(2)+exp(t); x(1)];
end
5.2-listing.
% boshlang`ich shart vektorini hosil qilamiz
x0=[1,-1];
% Integrallash intervalini, ya’ni ikki sonli massivni
% hosil qilamiz
interval=[0.25 2];
% ode45 funksiyasiga murojaat qilamiz
[T,X]=ode45(@FF, interval, x0);
% grafik yechimni chiqarish
plot(T,X(:,1),’:’,T,X(:,2),’-’);
legend(“y”, ‘x - Yechim’);
grid on;
1-rasm. (5.4) sistemaning grafik yechimi.
Differensial tenglamalar va ularning sistemalarini yechish uchun mo`ljallangan boshqa funksiyalarga ham shu tarzda murojaat qilish mumkin. Differensial tenglamalarni yechishda qo`llaniladigan MATLAB funksiyalarini izchil o`rganish uchun paketning ma’lumotlar tizimiga [4] murojaat qilish zarur.
Xulosa
Differensial tenglamalarning yechimlari aniq (analitik) va taqribiy (sonli) bo`lishi mumkin. Ba’zi differensial tenglamalarni aniq yechish mumkin bo`lsa, amaliyotda shunday tenglamalar, ayniqsa, ularning shunday sistemalari mavjudki, ularning aniq yechimlarini topib bo`lmaydi. Hattoki, analitik yechimga ega bo`lgan tenglamalar uchun ham ba`zi hollarda oldindan berilgan qiymatlardagi sonli yechimlarni topishga to`g`ri keladi. Shuning uchun ham oddiy differensial tenglamalarni sonli yechish usullari rivoj topdi.
Ushbu referatda hisoblash matematikasining amaliyotda ko`p uchraydigan va kompyuterda hisoblashlari zaruriyati yuqori bo`lgan masalalarni MatLab matematik paketida yechish uslubiyotini tadqiq qilish maqsad qilib qo`yilgan edi. Boshqa matematik paketlarning ichida aynan MatLab paketining tanlab olinishi, unda dasturlash imkoniyatining mavjudligi va hisoblash jarayonini to`liq kuzatish va boshqarish mumkinligidadir.
MatLab (Matrix Laboratory) tizimi Amerikaning MathWorks firmasi mahsuloti bo’lib, bu tizim katta imkoniyatlarga ega bo’lgan dasturiy mahsulotdir [5]. Uning birinchi versiyasi 1970 yilda foydalanuvchilarga havola etilgan. U ilmiy va muhandislik masalalarini yechuvchi ko’plab maxsus dasturlardan tashkil topgan. Uning asosiy elementi - bu MatLab sistemasining yadrosi. Bunga qo’shimcha tarzda unda 60 ga yaqin buyruqlar kompleksi ("Toolboxes") biriktirilgan. U Curve Fitting Toolbox, Optimization Toolbox, Partial Differential Equation Toolbox, Statistics Toolbox, Symbolic Math Toolbox va boshqa amaliy dasturlar paketlaridir. MatLab tizimining boshqa kompyuter algebrasi tizimlariga nisbatan yana bir muhim tomoni shundaki, unda dasturlash imkoniyatining mavjudligi va hisoblash jarayonini boshqarish hamda kuzatish mumkinligidadir.
Ushbu referat “MatLab dasturida xususiy hosilalali differensial tenglamalarni yechish” mavzusiga bag`ishlangan bo`lib, bunda differensial tenglamalar haqida umumiy tushunchalar berilgan, xususiy hosilali differensial tenglamalar va ularni yechish usullari keltirilgan. Shuningdek, differensial tenglamalarni yechish bo`yicha MatLab dasturining funksiyalari tadqiq qilingan, hamda MatLab dasturida xususiy hosilali differensial tenglamalarni yechish misollar orqali batafsil yoritilgan.
Do'stlaringiz bilan baham: |