Metode Galerkin: Perbedaan antara revisi

Konten dihapus Konten ditambahkan
Boulevard (bicara | kontrib)
Tidak ada ringkasan suntingan
k melakukan perapian bagian; mengubah "simetrik" menjadi simetris"
 
(34 revisi perantara oleh 6 pengguna tidak ditampilkan)
Baris 1:
{{rapikan}}
 
Dalam [[matematika]], khususnya bidang [[analisis numerik]], '''metode galerkin''' merupakan metode yang digunakan untuk mengubah masalah operasioperator kontinu (seperti [[persamaan differensial]]) ke masalah diskret. Dalam prinsipnya, metode ini penerapannya mirip penerapannya dengan ''[[metode variasi]]'' ke ruang fungsi dengan merubahmengubah parsamaannya ke [[formulasi lemah]]. TipeYang yangsecara pertamakhusus menggunakanmenerapkan beberapa batasan pada ruang fungsi untuk menggolongkan ruang denganpada suatu himpunan terbatas dari basis fungsi. Seringkali ketikapada penggunaannya, metode Galerkin yang pertama menyajikan juga metode approksimasi yang biasa digunakan pada umumnya, seperti metode [[Petrov-Galerkin]] atau [[metode Ritz-Galerkin]].
Metode Galerkin
 
Pendekatan yang berharga untukoleh matematikawan Rusia [[Boris Galerkin]].
Dalam [[matematika]], khususnya bidang [[analisis numerik]], '''metode galerkin''' merupakan metode untuk mengubah masalah operasi kontinu (seperti [[persamaan differensial]]) ke masalah diskret. Dalam prinsipnya, metode ini penerapannya mirip dengan ''[[metode variasi]]'' ke ruang fungsi dengan merubah parsamaannya ke [[formulasi lemah]]. Tipe yang pertama menggunakan beberapa batasan pada ruang fungsi untuk menggolongkan ruang dengan suatu himpunan terbatas dari basis fungsi. Seringkali ketika penggunaannya, metode Galerkin yang pertama menyajikan juga metode approksimasi yang biasa digunakan, seperti metode Petrov-Galerkin atau metode Ritz-Galerkin.
 
Sejak keindahan metode Galerikin terungkap dalam cara yang sangat abstrak dari studi mereka, maka pertama kali kita akan memberikan abstrak turunannya. Pada akhirnya, kita akan memberikan contoh untuk penggunaannya.
Pendekatan yang berharga untuk matematikawan Rusia Boris Galerkin.
 
Sejak keindahan metode Galerikin terungkap dalam cara yang sangat abstrak dari studi mereka, pertama kali kita akan memberikan abstrak turunannya. Pada akhirnya, kita akan memberikan contoh untuk penggunaannya.
Contoh-contoh metode Galerkin adalah:
1. # [[Metode elemen berhingga]]
2. # [[Metode elemen pembatas]] untuk menyelesaikammenyelesaikan persamaan integral
3. # [[Metode subruang Kyrlov]]
 
== Pengenalan Masalahmasalah abstrak Abstrak==
 
a. === Masalah padadalam formulasi lemah ===
Misalkan kita memperkenalkan metode Galerkin dengan sebuah masalah abstrak yang merupakan suatu formulasi lemah pada ruang Hilbert yaitu ''V'', jika diketahui ''u''\in''V''dan untuk setiap''v''\in''V'' maka
a(u,v)=f(v)
adalah benar. a(. . . , . . . )Sekarang adalah bentuk bilinear (penjelasan yang eksak atas a(. . . , . . . ) akan ditentukan selanjutnya) dan ''f'' adalah operator linear terbatasi pada ''V''.
 
Misalkan kita memperkenalkanmemasukkan metode Galerkin denganpada sebuah masalah abstrak yang merupakan suatu [[formulasi lemah]] pada [[ruang Hilbert]] yaitu ''V'', jika diketahui ''<math> u''\in'' V''dan </math> sehingga untuk setiap'' <math> v''\in'' V'' </math> maka
b. Diskretisasi Galerkin
<center><math> a(u,v) = f(v) </math>.</center>
Dipilih subruangdengan dimensi yang lebih kecil (sebenarnya, kita akan mengasumsikan bahwa indeks n menujukkan dimensinya) dan memecahkan masalah yang perhitungkan.
adalah benar. Sekarang <math>a(. . . \cdots, . . .\cdots )Sekarang</math> adalah bentuk [[bilinear]] (penjelasan yang eksak atas <math>a(. . . \cdots, . . .\cdots )</math> akan ditentukan selanjutnya) dan ''f'' adalah operator linear terbatasipembatas pada ''V''.
Jika diketahui dan untuk setiap maka
Kita akan menyebut persamaan ini sebagai persamaan Galerkin. Dengan catatan bahwa persamaan ini tidak dapat dirubah dan hanya ruangnya yang dapat dirubah.
 
c. Ortogonalitas=== Diskretisasi Galerkin ===
 
Hal ini adalah sifat mendasar yang membuat analisis matematika dari metode Galerkin sangat jelas. Karena , kita dapat menggunakan sebagai vector dalam persamaan awal. Substitusi persamaan yang kedua, kita dapati ortogonalitas Galerkin untuk galat
DipilihPilih subruangdengansubruang <math> v_n \subset V</math> dengan dimensi yang lebih kecil (sebenarnya, kita akan mengasumsikan bahwa indeks ''n'' menujukkan dimensinya) dan memecahkan masalah yang perhitungkan.
Jika diketahui <math> u_n \in V_n </math> dan untuk setiap <math> v_n \in V_n </math> maka
Sekarang adalah galat antara solusi masalah awal dan persamaan Galerkin secara berturut-turut.
<center><math> a(u_n,v_n) = f(v_n) </math>.</center>
Kita akan menyebut persamaan ini sebagai persamaan Galerkin. Dengan catatan bahwa persamaan ini tidak dapat dirubahdiubah dan hanya ruangnya yang dapat dirubahdiubah.
 
=== Ortogonalitas Galerkin ===
 
Hal ini adalahmerupakan sifat mendasar yang membuat analisis matematika dari metode Galerkin sangat jelas. Karena <math>v_n \subset V</math>, kita dapat menggunakan <math> v_n </math> sebagai vector dalam persamaan awal. Substitusi persamaan yang kedua, kita dapati ortogonalitas Galerkin untuk galat
<center><math> a(e_n, v_n) = a(u,v_n) - a(u_n, v_n) = f(v_n) - f(v_n) = 0 </math>.</center>
Sekarang, <math>e_n</math> = ''u'' – <math>u_n</math> adalah galat antara solusi masalah awal ''u'' dan persamaan Galerkin <math>u_n</math> secara berturut-turut.
 
d. === Bentuk Matriks ===
 
d. Bentuk Matriks
Karena tujuan dari metode Galerkin adalah membentuk sistem persamaan linear, maka kita membangun bentuk matriksnya, sehingga dapat digunakan untuk menghitung solusi dengan program computer.
Misal <math>e_1, e_2, \cdots,e_n </math> basis untuk <math> v_n </math>. Maka hal ini cukup untuk menguji coba persamaan Galerkin, sebagai contoh:
Diketahui <math> u_n \in V_n </math> sehingga
<center><math> a(u_n, e_i) = f(e_i) </math>.</center>
Kita akan mengembangkan <math> u_n </math> menjadi basis seperti ini, <math> u_n = \sum_{j=1}^n u_j e_j</math> dan masukkanmemasukkannya kedalam persamaan di atas, sehingga diperoleh
<center><math> a(\sum_{j=1}^n u_j e_j, e_i) = \sum_{j=1}^n u_j a(e_j,e_i) = f(e_i) </math> untuk <math>i = 1, \cdots, n</math>.</center>
Dalam persamaan sebelumnya, sebenarnya merupakan sistem persamaan linear <math> A_u = f </math>, dimana
<math> a_ij = a(e_j, e_i) </math> dengan <math> f_i = f(e_i) </math>
 
e. === Matriks SimetrikSimetris ===
 
Dalam kaitannya dengan definisi dari matriks entry, matriks dari persamaan Galerkin adalah simetriksimetris jika dan hanya jika bentuk [[bilinear]] <math>a( \cdots, \cdots )</math> adalah simetriksimetris.
 
== Analisis dari Metode Galerkin ==
 
Sekarang, kita akan membatasi diri kita pada bentuk bilinear simetriksimetris, yaitu:
<center><math> a(u,v) = a(u,v) </math>.</center>
BilamanaKarena ini tidakbukan benar-benar sebuah pembatasbatas dari metode Galerkin, aplikasi dari teori standar ini menjadi sangat mudah. Selanjutnya, metode [[Petrov-Galerkin]] dibutuhkan dalam kasus non-simetriksimetris.
Analisis dari hasil metode ini dihasilkan dalam dua langkah. Yang pertama, kita akan menunjukkan bahwa persamaan Galerkin adalah [[well-posed problem]] dalam pengertianmenurut [[Hadamard]] dan oleh karena itu kita mengakui persamaan ini sebagai solusi yang tunggal. DalamPada langkah kedua, kita mempelajari pendekatan kualitassifat dari solusi Galerkin <math>u_n</math> .
 
f. Analisis dari Metode Galerkin
Sekarang, kita akan membatasi diri kita pada bentuk bilinear simetrik, yaitu:
Bilamana ini tidak benar-benar sebuah pembatas dari metode Galerkin, aplikasi dari teori standar menjadi sangat mudah. Selanjutnya, metode Petrov-Galerkin dibutuhkan dalam kasus non-simetrik.
Analisis dari hasil metode ini dalam dua langkah. Yang pertama, kita akan menunjukkan bahwa persamaan Galerkin adalah well-posed problem dalam pengertian Hadamard dan oleh karena itu kita mengakui persamaan ini sebagai solusi tunggal. Dalam langkah kedua, kita mempelajari pendekatan kualitas dari solusi Galerkin .
Analisi ini kebanyakan akan mengacu pada dua sifat dari bentuk bilinear, yakni:
* Pembatasan: untuk setiap <math> u, v \in V</math> adalah benar bahwa
<center><math> a(u,v) \le C \lVert u \rVert \lVert v \rVert </math> untuk konstanta C > 0.</center>
* Eliptisitas: untuk setiap setiap <math> u \in V</math> adalah benar bahwa
<center><math> a(u,v) \ge c\lVert u\rVert ^2 </math> untuk konstanta c > 0 .</center>
 
Menurut teorema Lax-Milgram, ada dua kondisi implikatif well-posedness dari masalah awal dalam weak formulation ([[formulasi lemah)]]. Semua kaidah dalam bagian berikut ini akan dinormalisasikan untuk pertidaksamaan benar di atas (kaidah ini sering disebut juga kaidah energy).
 
g. === Well-posedness dari metode Galerkin ===
 
Karena <math> ,V_n \subset V </math> pembatasan dan eliptisitas dari bentuk bilinear berlaku bagi <math> V_n </math>. Oleh karena itu, Well-posedness dari metode Galerkin sebenarnya diturunkan dari Well-posedness dari masalah awal.
 
h. === Pendekatan Quasi-Best (Lemma Cèa) ===
 
Galat <math>e_n</math> = ''u'' – <math>u_n</math> antara solusi awal dan solusi Galerkin mengenal estimasi sbb:
 
<center><math>\lVert e_n \rVert</math> <math>\le</math> <math>\frac{C}{c}</math> <math>\overset {inf} {v_n \in V_n}</math> <math>\lVert u - v_n \rVert </math>.</center>
 
Ini artinya, bahwa sesuai dengan konstanta <math>\frac{C}{c}</math>, solusi Galerkin <math> u_n </math> adalah mendekati solusi awal ''u'' sebagai vector lainnya dalam <math> V_n </math> . Faktanya, hal ini cukup untuk mempelajari pendekatan dengan ruang <math> V_n </math>, dengan sepenuhnya melupakan tentang persamaan yang ssedang diselesaikan.
 
=== Bukti ===
Menurut teorema Lax-Milgram, ada dua kondisi implikatif well-posedness dari masalah awal dalam weak formulation (formulasi lemah). Semua kaidah dalam bagian berikut ini akan dinormalisasikan untuk pertidaksamaan benar di atas (kaidah ini sering disebut juga kaidah energy).
 
Karena buktinya sangat sederhana dan prinsip dasar dibalik semua metode Galerkin yaitu eliptisitas dan pembatasan pada bentuk [[bilinear]](pertidaksamaan) dan ortogonalitas Galerkin, kita punya <math> v_n \in V_n</math> sehingga:
g. Well-posedness dari metode Galerkin
Karena , pembatasan dan eliptisitas dari bentuk bilinear berlaku bagi . Oleh karena itu, Well-posedness dari metode Galerkin sebenarnya diturunkan dari Well-posedness dari masalah awal.
 
<center><math> c\lVert u\rVert ^2 \le a(e_n, e_n) = a(e_n, u-v_n) \le C\lVert e_n \rVert \lVert u-v_n \rVert </math>.</center>
h. Pendekatan Quasi-Best (Lemma Cèa)
Galat antara solusi awal dan solusi Galerkin mengenal estimasi sbb:
Ini artinya, bahwa sesuai dengan konstanta , solusi Galerkin adalah mendekati solusi awal u sebagai vector lainnya dalam . Faktanya, hal ini cukup untuk mempelajari pendekatan dengan ruang , dengan sepenuhnya melupakan tentang persamaan yang ssedang diselesaikan.
 
Bagi dengan <math> c\lVert e_n \rVert</math> dan ambil semua kemungkinan hasil akhir infimum lemma <math>v_h</math>.
i. Bukti
Karena buktinya sangat sederhana dan prinsip dasar dibalik semua metode Galerkin yaitu eliptisitas dan pembatasan pada bentuk bilinear(pertidaksamaan) dan ortogonalitas Galerkin, kita punya sehingga:
Bagi dengan dan ambil semua kemungkinan hasil akhir infimum lemma .