580 likes | 790 Views
Lietojumi. Jebkurš uzdevums, kurā līdz ar kāda lieluma x maiņu ir jāņem vērā arī šīs maiņas ātrums x ’, var tikt aprakstīts ar diferenciālvienādojumu. Daudzos uzdevumos vienādojuma sastādīšanas pamatprincips ir:
E N D
Lietojumi Jebkurš uzdevums, kurā līdz ar kāda lieluma x maiņu ir jāņem vērā arī šīs maiņas ātrums x’, var tikt aprakstīts ar diferenciālvienādojumu. Daudzos uzdevumos vienādojuma sastādīšanas pamatprincips ir: “lieluma maiņas ātrums ir vienāds ar tā palielināšanās ātrumu mīnus samazināšanās ātrums” Ja uzdevumā jāievēro arī paātrinājums, kā tas notiek daudzos mehānikas uzdevumos, procesu modelē ar otrās kārtas vienādojumu. Šādos uzdevumos bieži izmanto Ņūtona otro likumu F=ma.
Uzdevumi par ceļu un ātrumu x – pārvietojums, t – laiks, x’=v - ātrums Ja ātrums ir laika funkcija v=f (t), rodas diferenciālvienādojums: Ja zināms kustīgā punkta stāvoklis fiksētā sākuma momentā atrisinot Košī problēmu, viennozīmīgi atrodam stāvokli patvaļīgā momentā t:
Ja f (t)=const, vienādojums apraksta vienmērīgu kustību, f (t)=at, kustība ir vienmērīgi paātrināta. Analoģiska Košī problēma apraksta kustību arī gadījumā, ja ātrums ir ne tikai laika, bet arī noietā ceļa funkcija v=f (t,x): Piemērs: ātrums proporcionāli noietajam ceļam samazinās, “gājējs nogurst”:
Sastādot modeli, būtiska loma ir vienādojumā ieejošo parametru noteikšanai. • Eksperimentālo datu izmantošana. • Piemērs. Sprinta modelis: (1973.g. J.B.Keller) • http://www.math.duke.edu/education/ccp/materials/diffeq/sprints/sprints1.html , dv/dt = A - v/b A = 12.2 m/sec2 , b = 0.892 sec.
Radioaktīvā sabrukšana Radioaktīvās vielas sabrukšanas ātrums ir proporcionāls vielas esošajam daudzumam m(0)=m0 Pussabrukšanas periods
Piemērs. Metodes lietojums. Radioaktīvā oglekļa C14 relatīvais daudzums katrā dzīvā organismā (arī augos) ir tāds pats kā apkārtējā gaisā. Kad organisms mirst, ogļskābās gāzes uzņemšana beidzas, turpinās tikai radioaktīvā sabrukšana. C14 pussabrukšanas periods ir Dzīvā organismā Geigera skaitītājs uzrāda 13.5 sabrukšanas minūtē uz gramu vielas. Zinot Geigera skaitītāja rādījumu pētāmās vielas paraugā, var aprēķināt tās vecumu. Uzdevums. Noteikt Francijas aizvēsturisko alu gleznojumu vecumu, ja atrastajam oganiskā materiāla paraugam Geigera skaitītājs uzrāda 1.69 sabrukšanas minūtē uz gramu vielas.
Pieņemsim t=0 šobrīd, T<0 momentā, kad organiskais paraugs gāja bojā, q(t) oglekļa saturs paraugā momentā t.
Tā kā sabrukšanas ātrumu raksturo Geigera skaitītāja rādījumi, var atrast
Bioloģijas piemēri. Baktēriju vairošanās. y’=y, y(0)=y0 Ja dzīves apstākļi baktērijām (vai citām būtnēm) ir ļoti labi, var gadīties: Ja y(0)=y0, šāds likums ātri noved pie katastrofas, jo ,
Populāciju augšanas uzdevums… Kenijas iedzīvotāju skaits (miljonos) 40 gadu laika posmā no 1950. līdz 1990. gadam statistiski dots tabulā
Ar lineāro interpolāciju rēķinot, iegūst, ka dubultošanās periods ir apmēram 22,30049 gadi. Modelējot ar lineāro diferenciālvienādojumu un atrisinājumu ar eksponentfunkciju, dabū
Salīdzinājumam grafikā statistiskie dati un eksponentfunkcijas vērtības Faktiski iedzīvotāju skaits aug ātrāk nekā eksponentfunkcijas vērtība
Logistiskais likums Situācijā, kad apkārtējās vides resursi ir ierobežoti, notiek konkurence par tiem, populācijas augšanas ātrumu ierobežo savstarpējo kontaktu biežums, jo kontaktējoties sugas īpatņi viens otru var iznīcināt. populācija pieaug populācija samazinās
Ir cerība, ka Zemeslodes iedzīvotāju skaits pieaug pēc logistiskā likuma
Ja populācijas īpatņi migrē prom no dotā areāla vai arī tiek rūpnieciski izmantoti ar konstantu ātrumu, pieaugšanas ātrums samazinās: B pieaugot…
yn+1=pyn, Pietiekoši lielam qN:
Divu populāciju mijiedarbība Plēsoņu (y) - upuru (x) izturēšanās modelis • x’=ax-bxy • y’=-cy+kxy
Piemērs. Volterra – Lotkas sistēma. Eksperimentālie dati. Novērojumi veikti ilgstošā periodā. http://www.math.montana.edu/frankw//ccp/modeling/continuous/twovars/body.htm http://www.biology.ualberta.ca/courses.hp/bio331.hp/lectures/lect22/ PredatorPreyDynamics.htm
Kolmogorova pieeja. Par konkrētām parametru vērtībām neinteresējamies, uzsverot tikai proporcionalitāti, lielāks – mazāks, pozitīvs – negatīvs utml. īpašības. Jebkurā gadījumā ir svarīgi maksimāli samazināt ieejošo parametru skaitu, atstājot tikai būtiskos. Izvēlēties svarīgākos parametrus ir viena no modelēšanas mākslām.
Piemērs.Parametru skaita samazināšana Volterrra – Lotkas sistēmā. Reparametrizācija. Substitūcija: Izvēlamies koeficientus A, B, T, lai sistēma būtu pēc iespējas vienkārša:
Aprēķins Formulas: Sistēmā paliek viens būtisks parametrs Pāreja uz bezdimensionāliem mainīgajiem jāveic katrā uzdevumā!
A modelī raksturo upuru vairošanos B raksturo upuru apēšanas ātrumu, B(x,y)=B(x,.) B(.,y) B(x,.) nosaka apēšanas ātruma atkarību no upuru populācijas blīvuma fiksētam y. B(x,.)=kx.šādam modelimplēsoņa nekad nav paēdis. B(x,.)=kx x <x*,tad iestājas piesātinājums un x>x* B(x,.)=kx*
plēsoņa nevar apēst vairāk kā upurus, plēsoņam tāpat iestājas piesātinājums, bet mazam upuru populācijas blīvumam apēšanas ātrums ir mazākskā iepriekšējā gadījumā
Ja neņem vērā plēsoņu konkurenci par upuriem, tad B(.,y)=y Ar konkurenci C(y)=-cy-hy2 Konkrēts uzdevums no mikrobioloģijas u'(T) = (1 - u(T)) - Mu(T) v(T) / (A + u(T)), v'(T) = Mu(T) v(T) / (A + u(T)) - v(T), u ir proporcionāls substrāta koncentrācijai hemostatā, v mikroorganismu kultūras apjomam
Epidemioloģija 1. I(t) inficēto skaits momentā t, r – inficēšanās ātrums (cik cilvēkus dienā inficētā persona inficē no jauna), a - ātrums ar kādu atbrīvojas no infekcijas (atveseļojas vai nomirst) -1/dienu skaitu, kurās cilvēks ir inficēts 2. SIS modelis. S- uzņēmīgi, I - inficētie
4. SEIR modelis latentais periods infekcijas periods
http://www.muk.uni-hannover.de/~jansen/master_hjansen.pdf http://www.me.ucsb.edu/~moehlis/APC514/tutorials/tutorial_seasonal/node5.html
Statistikas dati SARS 2003. Mirušie Saslimušie Izdzīvojušie
Ķīmijas uzdevumi i0, i=1,...,k, i0,i=k+1,...,n. Stohiometriskie koeficienti
Piemēri. nA+nX=C
nA+nC=a, nB-nA=b Iespējamas divas līdzsvara vērtības!
Mehānika Piemērs Masas punkts ar masu m=1 kustas ārējā spēka F un berzes spēka iespaidā. Berze ļoti liela Pārejam uz citu laika skalu: Praktiski pirmās kārtas vienādojums
Pieņemsim Vienādojums Elastības spēks! Viens pats līdzsvara stāvoklis x=0
Pieliekam vēl konstantu ārējo spēku Saglabājas viens pats (stabils) līdzsvara stāvoklis
II kārtas vienādojumu lietojumiMehānika, elektrība u.c. Lineāri vienādojumi. Mehāniskās svārstības. Ideālās svārstības - bez berzes.
Trajektorijas fāzu plaknē
Svārstības ar lineāru berzi Maza berze