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

 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...

★ 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!