Ritzova metoda

Iz varijacijskog principa slijedi da se umjesto rješavanja jednadžbe, rješenje rubnog problema može dobiti rješavanjem problema minimizacije odgovarajućeg funkcionala. Ima raznih metoda koje koriste varijacijsku formulaciju. Jedna od njih je Ritzova. Pogledajmo kako ona funkcionira kad se radi o rubnom problemu #math2783#

#tex2html_wrap_indisplay38592#

gdje je #tex2html_wrap_inline38594# za svaki #tex2html_wrap_inline38596# Izaberemo #tex2html_wrap_inline38598# linearno nezavisnih funkcija #math2784##tex2html_wrap_inline38600# koje zadovoljavaju Dirichletov homogeni rubni uvjet. Rješenje se pretpostavi u obliku #math2785#

#tex2html_wrap_indisplay38602#

i neodređeni koeficijenti se odrede iz uvjeta da #tex2html_wrap_inline38604# minimizira funkcional #math2786#

#tex2html_wrap_indisplay38606#

Tako dobiveni #tex2html_wrap_inline38608# leži u vektorskom prostoru razapetom s funkcijama #math2787##tex2html_wrap_inline38610# tj. u vektorskom prostoru svih linearnih kombinacija funkcija #math2788##tex2html_wrap_inline38612# Rješenje ne mora biti linearna kombinacija tih funkcija, pa u tom slučaju #tex2html_wrap_inline38614# nije točno već samo približno rješenje. No, grešku koju tom pretpostavkom činimo, možemo umanjiti uzimanjem većeg #tex2html_wrap_inline38616# Prvi problem s kojim se susrećemo kod Ritzove metode je određivanje funkcija #tex2html_wrap_inline38618# koje se često zovu koordinatne funkcije. One se obično biraju u skladu s problemom koji se rješava. Na pr. ako iz fizikalnih razloga očekujemo periodičko rješenje, onda ćemo takvima pretpostaviti i funkcije #tex2html_wrap_inline38620# Inače se za koordinatne funkcije mogu koristiti polinomi. Nakon što smo izabrali funkcije #tex2html_wrap_inline38622# pretpostavljeno rješenje uvrstimo u funkcional #math2789#

#tex2html_wrap_indisplay38624#

#math2790#

#tex2html_wrap_indisplay38626#

#math2791#

#tex2html_wrap_indisplay38628#

#tex2html_wrap_inline38630# je derivabilna funkcija od #tex2html_wrap_inline38632# varijabli #math2792##tex2html_wrap_inline38634# pa jednadžbe #math2793#

#tex2html_wrap_indisplay38636#

za #math2794##tex2html_wrap_inline38638# predstavljaju nužan uvjet za ekstrem funkcije #tex2html_wrap_inline38640# u točki #math2795##tex2html_wrap_inline38642# Ovo je sustav od #tex2html_wrap_inline38644# linearnih algebarskih jednadžbi s #tex2html_wrap_inline38646# nepoznanica. Ako stavimo #math2796#

#tex2html_wrap_indisplay38648#

dobivamo sustav jednadžbi #math2797#

#tex2html_wrap_indisplay38650#

On se može matrično zapisati

#math2798#
#tex2html_wrap_indisplay38652# (3.26)

gdje je #math2799##tex2html_wrap_inline38654# Matrica #tex2html_wrap_inline38656# koja se inače zove matrica krutosti, je očito simetrična, jer je #math2800#

#tex2html_wrap_indisplay38658#

Ona je i pozitivno definitna. Doista, <#12120#> #math2801#

#tex2html_wrap_indisplay38660#

#math2802#

#tex2html_wrap_indisplay38662#

<#12120#> Budući da je #tex2html_wrap_inline38664# za svaki #tex2html_wrap_inline38666# ovaj integral se može poništavati samo tako da bude #math2803#

#tex2html_wrap_indisplay38668#

odnosno #math2804#

#tex2html_wrap_indisplay38670#

a to se može dogoditi jedino u slučaju da je #tex2html_wrap_inline38672# za svaki #tex2html_wrap_inline38674# tj. #math2805##tex2html_wrap_inline38676# Dakle #math2806##tex2html_wrap_inline38678# za #math2807##tex2html_wrap_inline38680# pa je matrica #tex2html_wrap_inline38682# pozitivno definitna. Kao što smo vidjeli, takva matrica se može dijagonalizirati i ima pozitivne vlastite vrijednosti. To znači da je regularna. Zaista, neka su njezine vlastite vrijednosti #math2808##tex2html_wrap_inline38684# Tada je #math2809#

#tex2html_wrap_indisplay38686#

Odatle slijedi da jednadžba (#eq:ritzmatr#10184>) ima jedno i samo jedno rješenje. Nedostak ove metode je u tome što je matrica #tex2html_wrap_inline38688# puna matrica, tj. općenito je svaki njezin element različit od nule.

Primjer 3.20   Riješimo Ritzovom metodom sljedeći rubni problem

  #tex2html_wrap_indisplay38690# ;SPMnbsp;;SPMnbsp;;SPMnbsp;
  #tex2html_wrap_indisplay38691# ;SPMnbsp;;SPMnbsp;;SPMnbsp;

Rješenje. Za koordinatne funkcije uzmimo polinome. Budući da na lijevom rubu imamo homogen Dirichletov uvjet, polinomi se moraju poništavati u nuli. Neka su, dakle, koordinatne funkcije #math2810#

#tex2html_wrap_indisplay38693#

Rješenje pretpostavljamo u obliku #math2811#

#tex2html_wrap_indisplay38695#

Kad s #tex2html_wrap_inline38697# uđemo u funkcional, i njegove derivacije po koeficijentima #math2812##tex2html_wrap_inline38699# izjednačimo s nulom, dobivamo sustav jednadžbi
#math2813#
#tex2html_wrap_indisplay38702# #tex2html_wrap_indisplay38704# #tex2html_wrap_indisplay38706#  
#tex2html_wrap_indisplay38708# #tex2html_wrap_indisplay38710# #tex2html_wrap_indisplay38712#  
#tex2html_wrap_indisplay38714# #tex2html_wrap_indisplay38716# #tex2html_wrap_indisplay38718#  
#tex2html_wrap_indisplay38720# #tex2html_wrap_indisplay38722# #tex2html_wrap_indisplay38724#  

Odavde dobijemo za koeficijente #math2814#

#tex2html_wrap_indisplay38726#

pa je rješenje #math2815#

#tex2html_wrap_indisplay38728#

Točno rješenje je #math2816#

#tex2html_wrap_indisplay38730#

Greška računata u točkama #math2817##tex2html_wrap_inline38732# iznosi #math2818#

#tex2html_wrap_indisplay38734#

#math2819#

#tex2html_wrap_indisplay38736#