Interpolasi lagrange

Interpolasi lagrange

Apa interpolasi Lagrange?

Interpolasi LaGrange adalah metode numerik perkiraan fungsi, yang memanfaatkan polinomial yang melewati titik -titik fungsi tertentu yang diketahui yang dimaksudkan untuk mendekati perkiraan.

Jika fungsi perkiraan lunak, bahkan di luar nilai yang diberikan atau diketahui, polinomial mengambil nilai dekat dengan fungsi yang menarik, terutama jika nilai -nilai ini berada di antara titik yang diberikan. Itulah sebabnya polinomial dianggap sebagai pendekatan yang baik untuk fungsi tersebut.

Gambar 1.- Formula untuk membangun polinomial lagrange. Sumber: f. Zapata.

Sekarang, misalkan Anda ingin memperkirakan suatu fungsi f (x) yang hanya diketahui nilainya X-yo-, dengan yo dari 0 sampai N-1. Yaitu, mereka saling kenal N poin (X-yo, Danyo) dengan Danyo = f (xyo), Dimana indeks yo Pergi dari 0 sampai N-1.

Dalam metode interpolasi lagrange, polinomial yang mendekati fungsi f (x) Itu adalah polinomial P (x) gelar N-1, dibangun oleh kombinasi linier N Polinomial Lyo(X) gelar N-1. Ini adalah Polinomial lagrange, yang dinyatakan sebagai berikut:

Nilai -nilai Danyo Mereka mewakili ordinat yang sesuai dengan absis Xyo Dimana fungsinya f (x) Diketahui, yaitu: Danyo = f (xyo).

Polinomial lagrange

Melalui kombinasi linier di antara mereka, polinomial lagrange bertindak sebagai dasar untuk pembangunan polinomial kelas N -1 yang akan berfungsi untuk interpolasi N poin yang diketahui.

Notasi untuk polinomial adalah lyo(x), dengan indeks I dalam kisaran dari 0 hingga n-1. Formula untuk membangun polinomial lagrange adalah sebagai berikut:

Simbol yang ditunjukkan menunjukkan bahwa produktori monomial n -1 harus dilakukan, mulai dari polinomial j = 0.

Karakteristik polinomial lagrange

1.- Polinomial lagrange persis sama dengan unit ketika dievaluasi di absis yang sesuai dengan indeks mereka, yaitu:

Lyo(Xyo) = 1

2.- Mereka dibatalkan di absis titik interpolasi dengan indeks berbeda dari polinomial yang sama:

Dapat melayani Anda: statistik deskriptif: sejarah, karakteristik, contoh, konsep

Lyo(XJ) = 0, dengan i ≠ j.

3.- Mengambil nilai absis lain yang berbeda dari titik interpolasi, polinomial lagrange memperoleh nilai antara -1 dan +1.

4.- Untuk mendapatkan polinomial lagrange, hanya perlu mengetahui absis titik ke interpoch.

Polinomial Lagrange Kedua

Polinomial Lagrange Kedua adalah yang paling sering digunakan saat Anda ingin melakukan interpolasi tiga poin.

Misalkan fungsi interpolar dikenal dalam tiga titik, yaitu:

(X0,Dan0); (X1, Dan1); (X2, Dan2)

Kemudian polinomial lagrange yang sesuai L0, L1 Dan L2 Mereka menjadi seperti ini:

L0(x) = [(x - x1) / (X0 - X1)] [(x - x2) / (X0 - X2)]

L1(x) = [(x - x0) / (X1 - X0)] [(x - x2) / (X1 - X2)]

L2(x) = [(x - x0) / (X2 - X0)] [(x - x1) / (X2 - X1)]

Perlu dicatat itu L0(X0) = L1(X1) = L2(X2) = 1, ketika Lyo(XJ) = 0 selama yo≠ j.

Polinomial interpolasi dari tingkat kedua

Penting untuk dicatat bahwa dalam polinomial interpolasi Lagrange, ordinat titik interpolasi adalah faktor polinomial Lagrange.

Dengan cara ini, begitu polinomial diperoleh untuk nilai -nilai absis tertentu, mereka berfungsi untuk menghitung polinomial interpolasi dari berbagai fungsi, asalkan yang dipesan dalam absis yang sebelumnya diketahui sebelumnya diketahui diketahui.

Dalam kasus interpolasi kelas dua:

P (x) = f (x0) L0(x) + f (x1) L1(x) + f (x2) L2(X)

Dan p (x) mendekati fungsi f (x) dalam interval (X0, X2).

Gambar 2.- Gambar ini menunjukkan cara mendapatkan polinomial lagrange untuk tiga titik interpolasi dan dari mereka, polinomial interpoling. Sumber: f. Zapata.

Contoh

Contoh 1

Temukan polinomial lagrange yang sesuai dengan tiga titik absis X0= 0, X1= 1 Dan X2= 2.

Seperti yang terlihat pada bagian sebelumnya, polinomial ini adalah:

Dapat melayani Anda: Fungsi Overjektif: Definisi, Properti, Contoh

L0(x) = [(x - 1) / (0 - 1)] [(x - 2) / (0 - 2)] = - (x -1) ⋅ (-½) (x - 2) = ½ X2 - 3x + 2)

L1(x) = [(x - 0) / (1 - 0)] [(x - 2) / (1 - 2)] = x ⋅ (-1) (x - 2) = - x2 + 2x

L2(x) = [(x - 0) / (2 - 0)] [(x - 1) / (2 - 1)] = (½) x ⋅ (x - 1) = (½) (x2 - X)

Gambar 3. Polinomial Lagrange untuk Nilai Absis 0, 1 dan 2. Sumber: f. Zapata.

Contoh 2

Anda ingin memperkirakan fungsinya f (x) = arcan (x) Dalam interval [0, 2]. Fungsi ini hanya nilainya yang diketahui X0= 0, X1= 1 Dan X2= 2, yang masing -masing Dan0= 0, Dan1= π/4 = 0,785 Dan Dan2= 1.107.

Oleh karena itu Anda harus menemukan polinomial interpoling P (x) mendekat f (x) Dalam interval yang ditunjukkan.

Dalam Contoh 1, polinomial lagrange telah ditentukan untuk nilai absis yang ditunjukkan dalam pernyataan ini, sehingga tidak perlu mengulangi perhitungan. Polinomial interpoling sekarang akan menjadi:

P (x) = f (x0) L0(x) + f (x1) L1(x) + f (x2) L2(X)

Yang setara dengan:

P (x) = y0 L0(x) + dan1 L1(x) + dan2 L2(X)

Dalam kasus khusus ini adalah:

P (x) = 0 ∙ (½) (x2 - 3x + 2) + 0.785 ∙ (- x2 + 2x) + 1.107 ∙ (½) (x2 - X)

Di atas disederhanakan untuk:

P (x) = 0,785 ∙ (- x2 + 2x) + 1.107 ∙ (½) (x2 - X)

Dan akhirnya tetap:

P (x) = -0.2315 ∙ x2 + 1.0165 ∙ x

Gambar 4. Interpoling polinomial yang diperoleh melalui polinomial lagrange yang mendekati fungsi busur-tangent dalam interval (0, 2). Titik interpolasi juga ditampilkan. Sumber: f. Zapata.

Latihan

Latihan 1

Dapatkan polinomial lagrange yang memadai untuk memiliki pendekatan terhadap fungsi:

f (x) = sin (x)

Dalam interval [0, π] dan dengan lima titik interpolasi.

Larutan

Pertama -tama, absis titik interpolasi ditentukan, yang dipilih sama dan termasuk ujung interval perkiraan. Dengan ini Anda punya:

X0= 0; X1= π/4; X2= π/2; X3= 3 π/4; X4= π.

Dapat melayani Anda: Ketimpangan segitiga: demonstrasi, contoh, latihan terpecahkan

Karena f (x) dibatalkan pada titik -titik ekstrem, tidak perlu mendapatkan polinomial lagrange l0 dan saya4.

Polinomial l1, L2 dan saya3 adalah:

L1 = [(x - x0) / (X1 - X0)] [(x - x2) / (X1 - X2)] [(x - x3) / (X1 - X3)] [(x - x4) / (X1 - X4)]

L2 = [(x - x0) / (X2 - X0)] [(x - x1) / (X2 - X1)] [(x - x3) / (X2 - X3)] [(x - x4) / (X2 - X4)]

L3 = [(x - x0) / (X3 - X0)] [(x - x1) / (X3 - X1)] [(x - x2) / (X3 - X2)] [(x - x4) / (X3 - X4)]

Sekarang kami mengganti nilai absis:

L1 = [(x - 0)/(π/4 - 0)] [(x - π/2)/(π/4 - π/2)] [(x - 3 π/4)/(π/4 - 3 π/4)] [(x - π)/(π/4 - π)]

L2 = [(x - 0)/(π/2 - 0)] [(x - π/4)/(π/2 - π/4)] [(x - 3 π/4)/(π/2 - 3 π/4)] [(x - π)/(π/2 - π)]

L3 = [(x - 0)/(3 π/4 - 0)] [(x - π/4)/(3 π/4 - π/4)] [(x - π/2)/(3 π/ 4 - π/2)] [(x - π)/(3 π/4 - π)]

Denominator diselesaikan:

L1 = [x/π/4] [(x - π/2)/( - π/4)] [(x - 3 π/4)/( - π/2)] [(x - π)/( - 3π/4)]

L2 = [x/π/2] [(x - π/4)/(π/4)] [(x - 3 π/4)/( - π/4)] [(x - π)/( - π /2)]

L3 = [x/(3 π/4)] [(x - π/4)/(π/2)] [(x - π/2)/(π/4)] [(x - π)/( - π/4)]

Disederhanakan dan dikelompokkan kembali untuk mendapatkan:

L1 = x (x - π/2) (x - 3 π/4) (x - π)/( - 3 π 4/128)

L2 = x (x - π/4) (x - 3 π/4) (x - π)/(π 4/64)

L3 = x (x - π/4) (x - π/2) (x - π)/( - 3 π 4/128)

Latihan 2

Dapatkan polinomial interpolasi yang mendekati fungsi sen (x) dalam interval [0, π] dengan lima titik interpolasi yang dipilih dalam Latihan 1 dan polinomial lagrange masing -masing.

Larutan

Polinomial interpolasi adalah:

P (x) = sin (0) * l0 + Sen (π/4) * l1 + Sen (π/2) * l2 + Sen (3π/4) * l3 + Sen (π) * l4

Mengevaluasi fungsi sinus dan pengalikan adalah:

P (x) = (√2/2) l1 + 1 * l2 + (-√2/2) l3

Setelah karya aljabar yang sulit, polinomial interpolasi adalah:

P (x) = 2. 7481 x4 -limabelas. 138 x3 +23. 467 x2 - 9. 5236 x

Referensi

  1. Goodman, a. L. H. seribu sembilan ratus sembilan puluh enam. Aljabar dan trigonometri dengan geometri analitik. Pendidikan Pearson.
  2. Harpe, hlm. D. (2000). Topik dalam teori kelompok geometris. University of Chicago Press.
  3. Hazewinkel, m. (2001). Interpolasi linier ", ensiklopedia matematika.
  4. Hoffmann, e. (2002). Hingga kronologi interpolasi: dari astronomi kuno hingga sinyal dan pemrosesan gambar modern. Prosiding IEEE.
  5. Wikipedia. Interpolasi polinomial lagrange. Pulih dari: wikipedia.com