Laboratoriniai darbai

Triįstrižainių lygčių sprendimas ir aproksimacija splainais

9.0   (3 atsiliepimai)
Triįstrižainių lygčių sprendimas ir aproksimacija splainais 1 puslapis
Triįstrižainių lygčių sprendimas ir aproksimacija splainais 2 puslapis
Triįstrižainių lygčių sprendimas ir aproksimacija splainais 3 puslapis
Triįstrižainių lygčių sprendimas ir aproksimacija splainais 4 puslapis
Triįstrižainių lygčių sprendimas ir aproksimacija splainais 5 puslapis
Triįstrižainių lygčių sprendimas ir aproksimacija splainais 6 puslapis
Triįstrižainių lygčių sprendimas ir aproksimacija splainais 7 puslapis
www.nemoku.lt
www.nemoku.lt
Aukščiau pateiktos peržiūros nuotraukos yra sumažintos kokybės. Norėdami matyti visą darbą, spustelkite peržiūrėti darbą.
Ištrauka

 Vilniaus Gedimino Technikos Universitetas Matematinio modeliavimo katedra Skaitiniai metodai ir algoritmai Laboratorinis darbas Nr. 1 Triįstrižainių lygčių sprendimas ir aproksimacija splainais Tikrino: O. Štikonienė Programos kodas: Triistrizaines sistemos sprendimas perkelties metodu Program perktr real*8,allocatable,dimension(:,:)::matr,Cd(:),D(:),r(:),x(:),a(:),b(:),c(:) integer n,i,j open(unit=1,file='duom.txt') open(unit=2,file='rezult.txt') read(1,*) n allocate(matr(n,n),Cd(n),D(n),r(n),x(n),a(n),b(n),c(n),stat=klaida) if(klaida/=0) then print*, 'Atmintis nepaskirta' stop end if ! Veiksmai su matrica read(1,*) ((matr(i,j),j=1,n),i=1,n) ! Nustatomi koeficientai a,b,c a(1)=0 do i=2,n a(i)=matr(i,i-1) end do do i=1,n b(i)=matr(i,i) end do c(n)=0 do i=2,n c(i-1)=matr(i-1,i) end do !Tikrinu ar pagr. istrizaine yra vyramoji do i=1,n do j=1,n if (DABS(a(i))+DABS(c(i))>DABS(b(i))) then print*, 'Pagr. istrizaines vyravimo salyga neispildyta' stop end if if(DABS(a(i))+DABS(c(i))f10.6)) 11 format (/25x, 'Reiksmiu stulpelis:'//(10x,f10.6)) 12 format (/30x, 'Sprendiniai:'//(25x, 'x(',I4,3x,')=', f10.6)) 13 format (/30x/(25x, 'x(', I4,3x, ')=', f10.6)) end program Kubinis splainas Program kubspl real,allocatable,dimension(:,:)::matr,x(:),y(:),g(:),Hd(:),Gd(:),e(:),a(:),b(:),c(:),Cd(:),D(:),r(:) real h,n,t integer k,i,j open(unit=1,file='rezult.txt') allocate (matr(9,9),x(11),y(11),g(11),Gd(10),Hd(10),e(10),a(9),b(9),c(9),Cd(9),D(9),r(9),stat=klaida) if (klaida/=0) then write(*,*) 'Atmintis nepaskirta!' stop end if ! Intervale [0,4] zingsnis h=4/10=0.4: h=0.4 k=11 do i=2,k x(1)=0 x(i)=x(i-1)+h end do write(1,*)'X:',(x(i), ';' ,i=1,k) ! Surandam funkcijos y(i) reiksmes intervale do i=1,k y(i)=sqrt(x(i))*cos(2*x(i)) end do write(1,*)'Y:',(y(i), ';' ,i=1,k) ! Sprendziant pekelties metodu do i=1,9 a(i)=h b(i)=4*h c(i)=h a(1)=0 c(9)=0 end do do i=1,9 do j=1,9 if (ABS(a(i))+ABS(c(i))>ABS(b(i))) then print*, 'Pagr. istrizaines vyravimo salyga neispildyta' stop end if end do end do ! d=6((y(i+2)-y(i+1)-y(i+1)-y(i))/h) do i=1,9 r(i)=6*((y(i+2)-y(i+1)-y(i+1)-y(i))/h) end do Cd(1)=-c(1)/b(1) D(1)=r(1)/b(1) Cd(9)=0 do i=2,8 Cd(i)=-c(i)/(Cd(i-1)*a(i)+b(i)) end do do i=2,9 D(i)=(r(i)-D(i-1)*a(i))/(Cd(i-1)*a(i)+b(i)) end do ! Gaunam g(i) reiksmes g=0 g(10)=D(9) do i=9,2,-1 g(i)=Cd(i-1)*g(i+1)+D(i-1) end do g(11)=0 g(1)=0 ! Koeficientai e, G, H do i=1,10 e(i)=(y(i+1)-y(i))/h -g(i+1)*h/6 - g(i)*h/3 Gd(i)=g(i)/2 Hd(i)=(g(i+1)-g(i))/6*h end do ! Sudarom kubini splaina write(1,*)' S(x):' write(*,*)' S(x):' do i=1,10 write(*,*) i,'=',Hd(i),'*(x-',x(i),')**3+',Gd(i),'*(x-',x(i),')**2+',e(i),'*(x-',x(i),')+', y(i) write(1,*) i,'=',Hd(i),'*(x-',x(i),')**3+',Gd(i),'*(x-',x(i),')**2+',e(i),'*(x-',x(i),')+', y(i) end do print*, 'Iveskite taska intervale[0;4]' print*, 'kuriame norite apskaiciuoti splaino funkcijos reiksme' read(*,*) t do i=1,11 if (t>x(i)) then if (t

Daugiau informacijos...

Šį darbą sudaro 699 žodžiai, tikrai rasi tai, ko ieškai!

★ Klientai rekomenduoja


Šį rašto darbą rekomenduoja mūsų klientai. Ką tai reiškia?

Mūsų svetainėje pateikiama dešimtys tūkstančių skirtingų rašto darbų, kuriuos įkėlė daugybė moksleivių ir studentų su skirtingais gabumais. Būtent šis rašto darbas yra patikrintas specialistų ir rekomenduojamas kitų klientų, kurie po atsisiuntimo įvertino šį mokslo darbą teigiamai. Todėl galite būti tikri, kad šis pasirinkimas geriausias!

Detali informacija
Darbo tipas
Lygis
Universitetinis
Failo tipas
Word failas (.doc)
Apimtis
7 psl., (699 ž.)
Darbo duomenys
  • Matematikos laboratorinis darbas
  • 7 psl., (699 ž.)
  • Word failas 86 KB
  • Lygis: Universitetinis
www.nemoku.lt Atsisiųsti šį laboratorinį darbą
Privalumai
Pakeitimo garantija Darbo pakeitimo garantija

Atsisiuntei rašto darbą ir neradai jame reikalingos informacijos? Pakeisime jį kitu nemokamai.

Sutaupyk 25% pirkdamas daugiau Gauk 25% nuolaidą

Pirkdamas daugiau nei vieną darbą, nuo sekančių darbų gausi 25% nuolaidą.

Greitas aptarnavimas Greitas aptarnavimas

Išsirink norimus rašto darbus ir gauk juos akimirksniu po sėkmingo apmokėjimo!

Atsiliepimai
www.nemoku.lt
Dainius Studentas
Naudojuosi nuo pirmo kurso ir visad randu tai, ko reikia. O ypač smagu, kad įdėjęs darbą gaunu bet kurį nemokamai. Geras puslapis.
www.nemoku.lt
Aurimas Studentas
Puiki svetainė, refleksija pilnai pateisino visus lūkesčius.
www.nemoku.lt
Greta Moksleivė
Pirkau rašto darbą, viskas gerai.
www.nemoku.lt
Skaistė Studentė
Užmačiau šią svetainę kursiokės kompiuteryje. :D Ką galiu pasakyti, iš kitur ir nebesisiunčiu, kai čia yra viskas ko reikia.
Palaukite! Šį darbą galite atsisiųsti visiškai NEMOKAMAI! Įkelkite bet kokį savo turimą mokslo darbą ir už kiekvieną įkeltą darbą būsite apdovanoti - gausite dovanų kodus, skirtus nemokamai parsisiųsti jums reikalingus rašto darbus.
Vilkti dokumentus čia:

.doc, .docx, .pdf, .ppt, .pptx, .odt