Kamu baru saja menyelesaikan Modul 01, eh sudah disuguhi soal Galat di lembar latihan. STMA4223 Metode Numerik memang begini. Modul 02 tentang Akar Persamaan langsung menyusul. Rasanya seperti lari cepat dari algoritma ke konvergensi.
Tapi Modul 04 dan 05 itu yang sering bikin deg-degan. Metode Jacobi dan Gauss-Seidel mirip, beda tipis. Lalu Dekomposisi LU Doolittle hadir dengan urutan yang beda lagi dari Crout. Kamu perlu banyak latihan yang membedakan semua itu. Baca juga soal UT Sains Data untuk variasi soal hitungan.
Nah, di halaman ini kami sediakan soal-soal yang sudah dikunci jawabannya. Plus pembahasan untuk setiap langkah penting. Dari Interpolasi Lagrange di Modul 06 sampai Metode Runge-Kutta di Modul 10. Semua dirancang supaya kamu paham proses, bukan sekadar hasil. Kalau perlu latihan matkul lain, langsung saja cek soal UAS Universitas Terbuka di sini.
Soal UT STMA4223 Metode Numerik
Galat yang dihitung sebagai selisih mutlak antara nilai hampiran dan nilai eksak disebut…
Istilah baku untuk selisih mutlak antara nilai hampiran dan nilai eksak adalah galat mutlak atau galat absolut.
Jika nilai eksak suatu besaran adalah 50 dan nilai hampirannya adalah 48, maka galat relatif dari hampiran tersebut adalah…
Galat relatif = |50-48|/50 = 2/50 = 0,04 atau 4%.
Sebuah bilangan hampiran 0.00320 memiliki tiga angka signifikan. Manakah pernyataan yang tepat mengenai angka signifikan dalam konteks ini…
Pada bilangan 0,00320, dua angka nol setelah koma dan sebelum 3 hanya berfungsi menentukan letak desimal, bukan angka signifikan. Angka 3, 2, dan nol terakhir adalah signifikan.
Saat melakukan serangkaian operasi aritmetika dengan bilangan yang telah dibulatkan, galat dari tiap langkah dapat terakumulasi dan memengaruhi hasil akhir. Fenomena ini dikenal sebagai…
Perambatan galat merujuk pada proses menumpuk dan menyebarnya galat dari setiap langkah komputasi menuju hasil akhir.
Ketika suatu deret tak hingga dihentikan pada suku tertentu dalam komputasi numerik, jenis galat yang dominan muncul adalah…
Menghentikan deret tak hingga pada suku tertentu berarti memotong sisa deretnya, sehingga galat yang ditimbulkan disebut galat pemotongan.
Urutan langkah yang terdefinisi dengan jelas, terbatas, dan sistematis untuk menyelesaikan suatu masalah komputasi disebut…
Definisi algoritma menekankan pada urutan langkah yang jelas, terbatas, dan sistematis untuk menyelesaikan masalah.
Sebuah metode iterasi dikatakan konvergen jika…
Konvergensi berarti barisan hampiran yang dihasilkan oleh iterasi semakin mendekati nilai limit atau solusi eksak.
Seorang mahasiswa menghentikan iterasi ketika selisih antara dua hampiran berurutan kurang dari 10^{-4}. Tindakan ini menerapkan…
Kriteria penghentian adalah syarat yang digunakan untuk mengakhiri iterasi, misalnya berdasarkan batas toleransi selisih dua hampiran berurutan.
Metode iteratif berbeda dari metode langsung terutama dalam hal…
Metode langsung menyelesaikan masalah dalam sejumlah langkah yang pasti tanpa pengulangan, sedangkan metode iteratif membangun barisan hampiran yang mendekati solusi.
Suatu barisan iterasi menunjukkan bahwa galat pada langkah berikutnya sebanding dengan kuadrat galat sebelumnya. Sifat ini menunjukkan…
Jika galat baru sebanding dengan kuadrat galat sebelumnya, maka metode memiliki konvergensi kuadratik, yang lebih cepat daripada konvergensi linier.
Metode bagi dua (bisection) memerlukan syarat awal berupa selang [a,b] yang memenuhi…
Syarat utama metode bagi dua adalah fungsi pada ujung-ujung selang harus berbeda tanda, yang menjamin adanya setidaknya satu akar di dalam selang.
Berbeda dengan metode bagi dua yang memilih titik tengah, metode posisi palsu menentukan hampiran akar berikutnya melalui…
Metode posisi palsu menggunakan garis lurus yang menghubungkan titik (a, f(a)) dan (b, f(b)) dan titik potong garis itu dengan sumbu x sebagai hampiran akar baru.
Dalam metode tertutup, selang yang memuat akar diidentifikasi melalui perubahan tanda nilai fungsi. Jika f(-3) = -4 dan f(1) = 5, maka selang awal untuk mencari akar adalah…
Karena f(-3) negatif dan f(1) positif, terjadi perubahan tanda yang menjamin adanya akar dalam selang [-3, 1].
PT Rekayasa Presisi menggunakan metode numerik untuk mencari akar persamaan pada desain struktur. Mereka memilih metode tertutup karena khawatir metode terbuka akan divergen. Kelebihan utama metode tertutup dibanding metode terbuka adalah…
Keunggulan utama metode tertutup adalah jaminan konvergensi selama selang awal memuat akar yang ditandai perubahan tanda fungsi pada ujung-ujung selang.
Tim riset mengamati bahwa metode posisi palsu kadang lebih lambat daripada metode bagi dua karena salah satu ujung selang cenderung tidak bergerak. Fenomena ini terjadi ketika…
Pada fungsi yang melengkung tajam (cekung/cembung), interpolasi linier metode posisi palsu dapat menyebabkan salah satu ujung selang stagnan, memperlambat konvergensi.
Sebuah perusahaan rintisan menggunakan metode iterasi titik tetap untuk menyelesaikan model keuangan berbentuk x = e^{-x} dengan tebakan awal x_0 = 0.5. Setelah beberapa iterasi, mereka mendapatkan barisan 0.5, 0.607, 0.545, 0.580, 0.560, 0.571,… Berdasarkan pola ini, kesimpulan yang paling tepat mengenai konvergensi barisan tersebut adalah…
Barisan menunjukkan osilasi di sekitar limit karena nilai bergerak naik-turun (0.5→0.607→0.545→0.580→…) mendekati sekitar 0.567, yang merupakan ciri konvergensi berosilasi pada iterasi titik tetap.
Pada implementasi metode Newton-Raphson untuk mencari akar f(x) = x^3 – 2x^2 + x – 1, seorang analis mulai dari tebakan x_0 = 0.5 dan memperoleh hampiran berikutnya selalu menjauhi akar yang seharusnya. Fenomena ini paling mungkin disebabkan oleh…
Rumus Newton-Raphson membagi f(x) dengan f'(x). Jika f'(x_0) sangat kecil, pembagian menghasilkan koreksi besar yang melontarkan hampiran menjauhi akar, sehingga menyebabkan divergensi.
Metode Secant sering dipilih sebagai alternatif metode Newton-Raphson. Manakah pernyataan yang paling tepat menjelaskan keunggulan utama metode Secant dibanding Newton-Raphson…
Metode Secant menggantikan turunan analitik f'(x) dengan hampiran beda hingga dari dua titik, sehingga pengguna tidak perlu menghitung atau memprogram turunan fungsi secara eksplisit.
Berbeda dengan metode Secant yang menggunakan pendekatan garis lurus melalui dua titik, metode Muller mendekati fungsi dengan suatu kurva tertentu. Bentuk kurva yang digunakan pada setiap iterasi metode Muller adalah…
Metode Muller membangun sebuah parabola (polinomial kuadratik) yang menginterpolasi tiga titik pada fungsi, kemudian menggunakan perpotongan parabola tersebut dengan sumbu-x sebagai hampiran akar selanjutnya.
Dalam pengembangan perangkat lunak analisis getaran, tim insinyur menemukan persamaan karakteristik yang memiliki akar kompleks. Mereka memilih metode Muller daripada metode Newton-Raphson. Alasan utama yang mendasari pilihan ini adalah…
Rumus iterasi metode Muller mengandung akar kuadrat dari diskriminan yang bisa bernilai negatif, sehingga secara alami menghasilkan bilangan kompleks dan dapat menemukan akar kompleks tanpa modifikasi khusus.
Tiga titik (x₀, f(x₀)), (x₁, f(x₁)), dan (x₂, f(x₂)) digunakan dalam satu langkah metode Muller. Hampiran akar selanjutnya diperoleh dari perpotongan parabola interpolasi dengan sumbu-x. Jika diskriminan dari rumus kuadratik yang terbentuk bernilai positif, maka…
Ketika diskriminan positif, rumus kuadrat menghasilkan dua solusi real. Metode Muller memilih solusi yang menghasilkan nilai mutlak penyebut lebih kecil, yang ekuivalen dengan memilih hampiran yang lebih dekat ke titik terakhir x₂.
Pada penerapan metode Muller untuk fungsi f(x) = x^3 – 2x – 5 dengan titik awal -2, -1, dan 0, satu langkah iterasi memberikan hampiran akar baru. Jika parabola interpolasi yang terbentuk tidak memotong sumbu-x di titik real, situasi ini mengindikasikan bahwa…
Parabola interpolasi yang tidak memotong sumbu-x secara real berarti diskriminan bernilai negatif, sehingga perpotongannya berupa akar kompleks. Ini adalah mekanisme alami metode Muller dalam menghasilkan hampiran akar kompleks.
Metode Bairstow mencari faktor kuadratik dari sebuah polinomial melalui pendekatan iteratif. Pada setiap iterasi, dua nilai dikoreksi secara serentak. Dua nilai yang dimaksud adalah…
Metode Bairstow mencari faktor kuadratik berbentuk x^2 – ux – v. Koreksi Δu dan Δv dihitung pada setiap iterasi untuk memperbarui koefisien u (koefisien x) dan v (konstanta) hingga sisa pembagian mendekati nol.
Dalam metode Bairstow, pembagian sintetik ganda dilakukan pada setiap iterasi. Tujuan utama dari prosedur pembagian sintetik ganda ini adalah…
Pembagian sintetik ganda menghasilkan sisa b₀ dan b₁ yang merupakan fungsi dari u dan v, serta turunan parsial sisa terhadap u dan v (melalui c₀, c₁, c₂) yang diperlukan untuk menghitung koreksi Δu dan Δv dalam metode Bairstow.
Setelah metode Bairstow berhasil menemukan satu faktor kuadratik dari polinomial berderajat enam, langkah selanjutnya yang dilakukan untuk mencari akar-akar yang tersisa adalah…
Metode Bairstow bersifat deflasi, yaitu setelah satu faktor kuadratik ditemukan, polinomial hasil bagi (yang berderajat lebih rendah) selanjutnya diproses dengan metode yang sama untuk mengekstraksi faktor kuadratik berikutnya hingga derajat polinomial tersisa dua atau satu.
Seorang mahasiswa menerapkan metode Bairstow pada polinomial x^4 – 3x^3 + 5x^2 – x – 10 dengan tebakan awal u = 1 dan v = -2. Setelah satu iterasi, ia memperoleh sisa b₀ dan b₁ yang tidak nol. Informasi dari sisa ini digunakan untuk…
Sisa b₀ dan b₁ yang belum nol menunjukkan bahwa faktor kuadratik dengan tebakan u dan v belum tepat. Koreksi Δu dan Δv dihitung dari sistem persamaan yang melibatkan sisa dan turunan parsialnya, lalu digunakan untuk memperbarui u dan v.
Sebuah sistem persamaan linear 3×3 memiliki matriks koefisien dengan elemen diagonal utama 8, 6, dan 9, sedangkan jumlah mutlak elemen non-diagonal pada setiap baris berturut-turut adalah 3, 2, dan 4. Berdasarkan syarat konvergensi metode iteratif, kesimpulan yang benar adalah…
Syarat cukup konvergensi metode Jacobi dan Gauss-Seidel adalah matriks dominan diagonal. Karena |8|>3, |6|>2, dan |9|>4 terpenuhi, matriks bersifat dominan diagonal sehingga kedua metode dijamin konvergen untuk sebarang tebakan awal.
PT Teknik Solusi menggunakan metode Jacobi untuk menyelesaikan sistem persamaan linear dengan matriks dominan diagonal. Pada iterasi ke-k, mereka menghitung nilai x^{(k)}_2 menggunakan x^{(k-1)}_1, x^{(k-1)}_2, dan x^{(k-1)}_3. Karakteristik utama dari prosedur ini yang membedakannya dari metode Gauss-Seidel adalah…
Metode Jacobi melakukan pembaruan serentak, artinya semua komponen x^{(k)} dihitung hanya menggunakan komponen dari iterasi sebelumnya x^{(k-1)}. Ini berbeda dari Gauss-Seidel yang langsung memakai komponen terbaru yang sudah dihitung dalam iterasi yang sama.
Dalam menyelesaikan SPL dengan metode Jacobi, seorang insinyur memilih tebakan awal berupa vektor nol. Setelah 10 iterasi, hampiran solusi masih jauh dari nilai eksak meskipun matriks dominan diagonal. Faktor yang paling memengaruhi lambatnya konvergensi ini adalah…
Meskipun matriks dominan diagonal menjamin konvergensi, laju konvergensi sangat bergantung pada seberapa kuat dominasi diagonalnya. Semakin kecil rasio elemen diagonal terhadap jumlah elemen non-diagonal, semakin lambat konvergensi metode Jacobi.
Pada iterasi metode Jacobi untuk sistem 5x + 2y = 16 dan x + 4y = 7, dengan tebakan awal (0,0), nilai y pada iterasi pertama adalah…
Dari persamaan kedua: y = (7 – x)/4. Dengan x_0 = 0 diperoleh y_1 = 7/4 = 1,75. Metode Jacobi menghitung semua variabel baru berdasarkan nilai iterasi sebelumnya secara serentak.
Metode Gauss-Seidel umumnya memerlukan lebih sedikit iterasi dibanding metode Jacobi untuk mencapai konvergensi dengan toleransi yang sama. Sifat ini terutama disebabkan oleh…
Gauss-Seidel berbeda dari Jacobi karena menggunakan komponen yang baru saja dihitung (x^{(k)}_1, x^{(k)}_2, …) untuk menghitung komponen berikutnya dalam iterasi yang sama. Pemanfaatan informasi terbaru ini mempercepat perambatan koreksi dan mengurangi jumlah iterasi total.
Untuk mempercepat konvergensi metode Gauss-Seidel pada sistem yang konvergen lambat, seorang analis menerapkan faktor relaksasi ω > 1 pada koreksi variabel. Teknik ini dikenal sebagai…
Pemberian faktor ω > 1 pada koreksi variabel dalam metode Gauss-Seidel disebut relaksasi atas (over-relaxation) atau SOR (Successive Over-Relaxation). Teknik ini memperbesar koreksi untuk mempercepat konvergensi ketika metode standar terlalu lambat.
Sebuah sistem persamaan linear diselesaikan dengan metode Gauss-Seidel. Setelah menghitung x_1^(k), nilai tersebut langsung digunakan untuk menghitung x_2^(k) pada iterasi yang sama. Karakteristik ini memberikan keunggulan utama metode Gauss-Seidel dibanding metode Jacobi, yaitu…
Gauss-Seidel memanfaatkan nilai komponen yang baru saja diperbarui dalam iterasi yang sama, sehingga informasi terbaru segera dipakai dan umumnya menghasilkan konvergensi lebih cepat.
PT Solusi Numerik menerapkan metode Gauss-Seidel pada sistem persamaan linear yang memiliki matriks koefisien dominan diagonal. Untuk mempercepat konvergensi, mereka memodifikasi koreksi tiap variabel dengan faktor ω = 1.5. Teknik percepatan ini disebut…
Pemberian bobot pada koreksi variabel dengan faktor relaksasi ω > 1 dikenal sebagai over relaxation (relaksasi atas), yang dapat mempercepat konvergensi Gauss-Seidel.
Pada metode dekomposisi LU Doolittle, matriks koefisien A difaktorkan menjadi L dan U dengan syarat khusus pada matriks L. Syarat yang membedakan metode Doolittle dari metode Crout adalah…
Metode Doolittle menetapkan elemen diagonal matriks segitiga bawah L bernilai satu, sedangkan Crout menetapkan elemen diagonal U bernilai satu.
Setelah mendekomposisi matriks A menjadi L dan U menggunakan metode Doolittle, langkah penyelesaian SPL Ax = b dilakukan dalam dua tahap. Urutan tahap yang benar adalah…
Setelah dekomposisi A = LU, sistem Ax = b menjadi LUx = b. Tahap pertama menyelesaikan Ly = b dengan substitusi maju, kemudian Ux = y dengan substitusi mundur.
Seorang analis menerapkan dekomposisi LU Doolittle pada matriks 3×3 berikut: baris pertama [2, -1, 1], baris kedua [4, 1, -3], baris ketiga [6, -2, 5]. Ia menghitung elemen u_{12} dari matriks U. Nilai u_{12} yang diperoleh adalah…
Karena l_{11} = 1 pada Doolittle, maka u_{12} = a_{12} / l_{11} = -1 / 2 = -0.5. Elemen U baris pertama dihitung langsung dari elemen A baris pertama dibagi elemen diagonal L.
Dalam metode dekomposisi LU Doolittle, setelah matriks L dan U diperoleh, vektor antara y dihitung melalui substitusi maju. Jika sistem awalnya 3×3, elemen y_2 dihitung menggunakan…
Substitusi maju menyelesaikan Ly = b baris demi baris. Untuk baris kedua: l_{21} y_1 + l_{22} y_2 = b_2, sehingga y_2 = (b_2 – l_{21} y_1) / l_{22}.
Sebuah tim riset membandingkan metode Doolittle dan Crout untuk mendekomposisi matriks yang sama. Pada metode Crout, elemen diagonal matriks U bernilai satu. Konsekuensi dari ketentuan ini terhadap perhitungan elemen L adalah…
Pada metode Crout, diagonal U bernilai 1 sehingga elemen L langsung mengambil nilai kolom pertama matriks A tanpa pembagian.
Teknisi komputasi di PT Maju Jaya menerapkan metode Crout pada sistem 4x + 3y = 10 dan 2x + 5y = 11. Setelah dekomposisi, teknik ini menghasilkan l_{21} dengan rumus…
Pada Crout, diagonal U = 1 sehingga u_{11} = 1. Elemen L kolom pertama langsung sama dengan elemen A kolom pertama, termasuk l_{21} = a_{21}. Tidak perlu pembagian karena diagonal U sudah 1.
Metode Crout dan Doolittle sama-sama memfaktorkan matriks menjadi L dan U, namun urutan penentuan elemennya berbeda. Pada metode Crout, strategi penentuan elemen L dan U dilakukan secara…
Metode Crout menghitung elemen secara kolom demi kolom: elemen L pada satu kolom, lalu elemen U pada baris yang bersesuaian, bergerak dari atas ke bawah.
Mahasiswa bernama Dian harus menyelesaikan SPL dengan metode Crout. Setelah memperoleh L dan U, ia melakukan substitusi maju pada Ly = b. Jika elemen l_{32} = 3, l_{33} = 5, y_2 = 4, dan b_3 = 26, maka y_3 yang tepat adalah…
Substitusi maju untuk baris ketiga: l_{31} y_1 + l_{32} y_2 + l_{33} y_3 = b_3. Jika l_{31} = 0, maka y_3 = (b_3 – l_{32} y_2) / l_{33} = (26 – 12) / 5.
Diberikan tiga titik data (1,3), (2,7), (4,15). Polinomial interpolasi Lagrange yang melalui ketiga titik tersebut dibangun dari polinomial basis. Polinomial basis L_1(x) yang terkait dengan titik (2,7) adalah…
Polinomial basis Lagrange L_1(x) (jika indeks 1 untuk x_1 = 2) bernilai 1 di x = 2 dan 0 di titik data lain. Rumusnya L_1(x) = (x-x_0)(x-x_2) / (x_1-x_0)(x_1-x_2) = (x-1)(x-4) / (2-1)(2-4).
Karakteristik utama polinomial basis Lagrange yang membedakannya dari polinomial interpolasi lainnya adalah…
Polinomial basis Lagrange L_i(x) dirancang agar bernilai tepat 1 pada x = x_i dan 0 pada x = x_j untuk semua j ≠ i. Ini adalah ciri khas konstruksi interpolasi Lagrange.
Sebuah laboratorium mengukur suhu setiap jam dan memperoleh data (0,20), (1,23), (2,28), (3,35). Teknisi ingin menginterpolasi suhu pada t = 2.5 menggunakan polinomial Lagrange. Polinomial interpolasi P(2.5) dinyatakan sebagai…
Polinomial Lagrange menjumlahkan perkalian setiap nilai fungsi f(x_i) dengan polinomial basis L_i(x). Jadi P(2.5) = Σ f(x_i) L_i(2.5) = 20 L_0(2.5) + 23 L_1(2.5) + 28 L_2(2.5) + 35 L_3(2.5).
Kelemahan metode polinomial Lagrange dibanding metode beda terbagi Newton ketika digunakan secara komputasional adalah…
Pada Lagrange, setiap polinomial basis bergantung pada seluruh titik data. Jika satu titik baru ditambahkan, seluruh polinomial basis harus dihitung ulang. Sebaliknya, bentuk Newton memungkinkan penambahan suku baru tanpa mengubah koefisien sebelumnya.
Teknisi PT Data Analitik menyusun tabel beda terbagi untuk data (0,1), (1,3), (3,13), (4,21). Ia menghitung beda terbagi orde dua f[x_0, x_1, x_2]. Jika f[x_0, x_1] = 2 dan f[x_1, x_2] = 5, dengan x_0 = 0, x_2 = 3, maka f[x_0, x_1, x_2] adalah…
Beda terbagi orde dua: f[x_0, x_1, x_2] = (f[x_1, x_2] – f[x_0, x_1]) / (x_2 – x_0) = (5 – 2) / (3 – 0) = 3/3 = 1.
Polinomial interpolasi Newton dan polinomial Lagrange yang melalui himpunan titik data yang sama memiliki hubungan penting, yaitu…
Polinomial interpolasi bersifat unik untuk sekumpulan titik data tertentu. Bentuk Newton dan Lagrange hanyalah representasi berbeda dari polinomial yang sama, sehingga secara aljabar keduanya identik.
Dalam menyusun tabel beda terbagi untuk data (x_0, f_0), (x_1, f_1),…, (x_n, f_n), seorang mahasiswa mengamati bahwa mengubah urutan titik data tidak memengaruhi polinomial akhir. Sifat ini disebut…
Polinomial interpolasi tidak bergantung pada urutan titik data yang digunakan. Sifat ini disebut invarian urutan dan merupakan karakteristik penting dari polinomial interpolasi.
PT Analis Maju menyusun tabel beda terbagi untuk data (0,2), (2,8), (5,23), (8,50). Manager meminta teknisi menghitung ulang karena dicurigai ada kesalahan. Jika f[x_0, x_1] = 3 dan f[x_1, x_2] = 5, maka f[x_0, x_1, x_2] adalah…
Sesuai rumus beda terbagi orde dua, f[x_0, x_1, x_2] = (f[x_1, x_2] – f[x_0, x_1]) / (x_2 – x_0) = (5 – 3)/(5 – 0) = 2/5.
Diberikan data dengan jarak seragam h = 0.2 pada titik x = 0.0, 0.2, 0.4, 0.6. Nilai fungsi di keempat titik tersebut adalah 1.000, 1.221, 1.492, 1.822. Jika seorang analis ingin menginterpolasi nilai di sekitar x = 0.5, operator beda yang paling tepat digunakan adalah…
Titik evaluasi x = 0.5 berada di dekat akhir tabel (data berakhir di 0.6), sehingga beda hingga mundur paling sesuai karena menggunakan titik-titik sebelum posisi evaluasi.
Tim laboratorium mengukur konsentrasi larutan setiap 30 detik dan memperoleh tabel data dengan selang waktu seragam. Mereka akan menginterpolasi konsentrasi pada suatu waktu di antara dua pengukuran. Jika polinomial interpolasi yang digunakan berbasis tabel beda hingga, maka polinomial tersebut merupakan…
Polinomial Newton beda hingga dibangun dari tabel beda hingga dengan jarak data seragam; koefisiennya dihitung menggunakan beda hingga, bukan beda terbagi.
Diberikan data f(1.0) = 2.718, f(1.2) = 3.320, f(1.4) = 4.055, f(1.6) = 4.953 dengan h = 0.2. Seorang mahasiswa menghitung Δf(1.2) = f(1.4) – f(1.2) = 0.735 dan Δ^2 f(1.0) = 0.133. Jika ia menggunakan polinomial Newton beda hingga maju untuk menginterpolasi f(1.3), maka suku pertama polinomial tersebut adalah…
Suku pertama polinomial Newton beda hingga maju adalah f(x_0) yaitu nilai fungsi di titik awal tabel, bukan beda atau titik lainnya.
Budi menginterpolasi data dengan jarak seragam menggunakan dua pendekatan berbeda: beda terbagi dan beda hingga. Ia mengamati bahwa beda hingga memberikan hasil yang lebih efisien secara komputasi pada data ini. Faktor utama yang menyebabkan efisiensi tersebut adalah…
Pada data dengan jarak seragam, beda hingga tidak memerlukan pembagian oleh selisih absis yang bervariasi seperti pada beda terbagi, sehingga perhitungan lebih ringkas.
Siti menginterpolasi data pengukuran suhu dengan empat titik menggunakan polinomial kubik tunggal. Hasil interpolasi menunjukkan osilasi yang tidak diinginkan, sehingga ia mempertimbangkan interpolasi bagian demi bagian. Berbeda dari polinomial tunggal, interpolasi linear sepotong demi sepotong…
Interpolasi linear sepotong adalah teknik paling sederhana dari interpolasi bagian demi bagian, di mana setiap subinterval dihubungkan dengan garis lurus antara dua titik bertetangga.
Seorang insinyur di PT Konstruksi Akurat menggunakan polinomial kubik pada tiap subinterval data lendutan balok untuk memperoleh kurva yang mulus. Pada setiap simpul, ia memastikan kontinuitas turunan pertama dan kedua. Metode interpolasi yang memenuhi syarat tersebut adalah…
Spline kubik menggunakan polinomial kubik pada tiap subinterval dan memaksakan kontinuitas turunan pertama dan kedua di setiap simpul, menghasilkan transisi mulus antar segmen.
Mahasiswa semester akhir membandingkan spline kubik alami dan spline kubik terapit untuk menginterpolasi data profil kecepatan fluida di dekat dinding pipa. Pada spline kubik alami, syarat batas yang diterapkan adalah…
Spline kubik alami memberlakukan syarat batas berupa turunan kedua sama dengan nol di kedua ujung (f'' = 0), berbeda dengan spline terapit yang menetapkan turunan pertama.
Tim simulasi menerapkan interpolasi spline kubik pada 50 titik data. Mereka mengamati bahwa matriks yang terbentuk dari syarat kehalusan bersifat tridiagonal. Dampak dari struktur matriks ini terhadap proses komputasi adalah…
Struktur matriks tridiagonal memungkinkan penggunaan algoritma khusus (seperti algoritma Thomas) yang jauh lebih efisien daripada eliminasi Gauss pada matriks penuh.
Seorang analis keuangan ingin menghitung turunan pertama fungsi nilai aset f(x) di titik x menggunakan ekspansi deret Taylor. Ia menggunakan rumus (f(x+h) – f(x))/h. Hampiran yang ia gunakan termasuk jenis…
Rumus tersebut merupakan beda maju karena menggunakan titik x dan satu titik ke depan x+h; dari ekspansi deret Taylor dapat ditunjukkan bahwa galat pemotongannya berorde h.
Peneliti menghitung turunan numerik f'(x) dari data eksperimen dengan dua rumus berbeda: beda maju dan beda pusat. Pada lebar langkah h yang sama, beda pusat memberikan hasil yang lebih akurat. Keunggulan ini disebabkan oleh…
Simetri pada rumus beda pusat menyebabkan suku galat berorde h saling menghilangkan saat ekspansi Taylor maju dan mundur dijumlahkan, sehingga galat turun menjadi orde h^2.
Dalam eksperimen fisika, mahasiswa mengukur posisi partikel pada waktu t = 1.0, 1.1, 1.2, 1.3 detik. Ia ingin menghitung kecepatan di t = 1.2 menggunakan data yang paling relevan di sekitar titik tersebut. Rumus yang paling tepat adalah…
Untuk titik pusat di t = 1.2, beda pusat menggunakan titik simetris 1.1 dan 1.3 dengan penyebut 2h = 0.2, memberikan akurasi orde h^2.
PT Rekayasa Mesin menerapkan rumus turunan numerik untuk menghitung laju perubahan tekanan pada kondisi transien. Mereka menggunakan beda mundur (f(x) – f(x-h))/h karena data hanya tersedia hingga waktu saat ini. Salah satu konsekuensi dari pilihan ini adalah galat pemotongan yang…
Beda mundur memiliki orde galat yang sama dengan beda maju, yaitu orde h, karena sama-sama diturunkan dari ekspansi Taylor satu sisi.
Tim komputasi menghitung turunan pertama f'(x) dengan dua lebar langkah: h dan h/2. Hasilnya adalah D(h) = 3.120 dan D(h/2) = 3.085. Mereka menerapkan ekstrapolasi Richardson untuk memperbaiki ketelitian. Nilai hasil ekstrapolasi tersebut adalah…
Untuk beda maju/mundur berorde h, ekstrapolasi Richardson: D_improved = (4D(h/2) – D(h))/3 = (4×3.085 – 3.120)/3 = (12.34 – 3.12)/3 = 9.22/3 ≈ 3.073.
Prinsip fundamental yang mendasari metode ekstrapolasi Richardson adalah melakukan kombinasi linear dua hampiran dengan lebar langkah berbeda untuk…
Ekstrapolasi Richardson mengombinasikan dua hampiran dari lebar langkah berbeda untuk mengeliminasi suku galat orde rendah, sehingga orde ketelitian meningkat.
Seorang peneliti menghitung turunan kedua f''(x) menggunakan rumus (f(x+h) – 2f(x) + f(x-h))/h^2. Rumus ini diturunkan dari…
Rumus turunan kedua tersebut diperoleh dengan menjumlahkan ekspansi Taylor f(x+h) dan f(x-h) serta menyusun ulang suku untuk mengisolasi f''(x).
PT Analitik Solusi menerapkan ekstrapolasi Richardson berulang untuk turunan pertama menggunakan rumus beda pusat. Ekstrapolasi pertama meningkatkan orde galat dari O(h^2) menjadi O(h^4). Langkah ini dilakukan dengan mengombinasikan dua hampiran D(h) dan D(h/2) menggunakan bobot…
Untuk beda pusat dengan galat orde h^2, ekstrapolasi Richardson menggunakan bobot (4^{1}D(h/2) – D(h))/(4^{1} – 1) = (4D(h/2) – D(h))/3.
PT Manufaktur Presisi mengintegralkan data laju produksi untuk menghitung total komponen yang dihasilkan. Mereka menggunakan aturan trapezoidal dengan lebar subinterval seragam dan memperoleh hasil 3245. Namun, mereka menduga hasil ini kurang akurat karena kurva laju produksi bersifat cekung ke atas. Berdasarkan karakteristik geometris aturan trapezoidal pada kurva cekung ke atas, jenis galat yang dominan terjadi adalah…
Pada kurva cekung ke atas, garis lurus yang menghubungkan dua titik berada di atas kurva, sehingga luas trapesium melebihi luas sebenarnya di bawah kurva. Akibatnya, aturan trapezoidal memberikan hampiran integral yang lebih besar (over-estimate).
Seorang analis keuangan menerapkan aturan Simpson 1/3 untuk mengestimasi nilai sekarang dari aliran kas kontinu. Ia membagi interval integrasi menjadi 6 subinterval dengan lebar seragam. Setelah menghitung, ia menyadari bahwa aturan Simpson 1/3 hanya dapat diterapkan pada jumlah subinterval tertentu. Agar aturan ini valid, jumlah subinterval yang digunakan harus memenuhi syarat…
Aturan Simpson 1/3 memerlukan sepasang subinterval untuk membentuk parabola interpolasi, sehingga jumlah total subinterval harus genap atau kelipatan 2. Dengan 6 subinterval, syarat ini terpenuhi dan aturan dapat diterapkan secara komposit.
Tim simulasi termal membandingkan dua hampiran integral: aturan trapezoidal dan aturan Simpson 1/3 pada fungsi f(x) = x^3 dengan interval [0,2] dan dua subinterval. Untuk fungsi kubik murni ini, perbandingan hasil kedua metode terhadap nilai eksak adalah…
Aturan Simpson 1/3 didasarkan pada interpolasi kuadratik yang mengintegralkan polinomial berderajat hingga tiga secara eksak. Sebaliknya, aturan trapezoidal menggunakan interpolasi linear yang hanya eksak untuk polinomial berderajat satu, sehingga pada fungsi kubik galatnya tidak nol.
PT Aerodinamika Nusantara melakukan uji terowongan angin dan memperoleh data gaya angkat pada tujuh titik dengan jarak seragam sepanjang sayap. Insinyur harus memilih aturan komposit untuk menghitung total gaya angkat. Jika ia menggunakan aturan Simpson 3/8 komposit pada sebagian interval, jumlah subinterval minimal yang harus dialokasikan untuk bagian aturan Simpson 3/8 saja adalah…
Aturan Simpson 3/8 menggunakan interpolasi kubik yang memerlukan tiga subinterval (empat titik) untuk satu panel. Dalam bentuk dasar non-komposit, aturan ini mensyaratkan jumlah subinterval kelipatan 3. Jumlah minimal untuk satu penerapan adalah tepat 3 subinterval.
Integrasi Romberg menggunakan ekstrapolasi Richardson pada aturan trapezoidal untuk meningkatkan ketelitian. Seorang mahasiswa menyusun tabel Romberg dan menghitung R_{2,2} = (4×R_{2,1} – R_{1,1})/3. Baris pertama tabel (k = 1) dihitung dengan lebar langkah h, baris kedua (k = 2) dengan h/2. Prinsip eliminasi galat yang mendasari rumus ekstrapolasi Romberg tersebut adalah…
Ekstrapolasi Richardson pada aturan trapezoidal menghilangkan suku galat orde h^2 melalui kombinasi linear R_{k+1,1} dan R_{k,1} dengan bobot 4/3 dan -1/3. Hasilnya, orde galat meloncat dari O(h^2) ke O(h^4), meningkatkan ketelitian secara signifikan.
Dalam menyusun tabel Romberg, setelah menghitung R_{1,1}, R_{2,1}, dan R_{3,1} menggunakan aturan trapezoidal dengan h, h/2, dan h/4, teknisi menghitung R_{2,2} dan R_{3,2} menggunakan ekstrapolasi Richardson. Untuk menghitung R_{3,3} pada kolom ketiga, data minimal yang diperlukan adalah…
Struktur tabel Romberg bersifat segitiga: setiap elemen R_{i,j} dihitung dari dua elemen di kiri dan kiri-atasnya, yaitu R_{i,j-1} dan R_{i-1,j-1}. R_{3,3} memerlukan R_{3,2} dan R_{2,2} melalui rumus R_{i,j} = (4^{j-1}×R_{i,j-1} – R_{i-1,j-1})/(4^{j-1} – 1).
Kuadratur Gauss-Legendre memilih titik evaluasi dan bobot secara optimal, berbeda dengan aturan Newton-Cotes yang menggunakan titik sama jarak. Keunggulan utama kuadratur Gauss dibanding aturan Newton-Cotes dengan jumlah titik yang sama adalah…
Dengan n titik evaluasi, kuadratur Gauss mencapai derajat ketelitian 2n-1, artinya polinomial berderajat hingga 2n-1 diintegralkan secara eksak. Sebaliknya, aturan Newton-Cotes dengan n titik hanya eksak untuk polinomial berderajat n-1 (aturan trapezoidal) atau n (aturan Simpson).
Untuk mengintegralkan fungsi pada interval [-1, 1] menggunakan kuadratur Gauss-Legendre dua titik, titik evaluasi dan bobot yang digunakan adalah x_1 = -1/sqrt(3), x_2 = 1/sqrt(3) dengan w_1 = w_2 = 1. Penentuan titik evaluasi ini didasarkan pada…
Titik evaluasi pada kuadratur Gauss-Legendre merupakan akar dari polinomial Legendre derajat n. Untuk n = 2, polinomial Legendre P_2(x) = (3x^2 – 1)/2 memiliki akar x = ±1/sqrt(3). Pemilihan ini memaksimalkan derajat ketelitian integral numerik.
Sebuah tangki air mengalami kebocoran sehingga ketinggian air y(t) memenuhi dy/dt = -k×sqrt(y) dengan k = 0.2. Insinyur menggunakan metode Euler dengan lebar langkah h = 0.5 untuk menaksir ketinggian air setelah 1 detik dari kondisi awal y(0) = 4. Hasil hampiran y(1) yang diperoleh adalah…
Metode Euler memperbarui solusi dengan y_{i+1} = y_i + h×f(t_i, y_i). Satu langkah dari t=0 ke t=1 dengan h=0.5 memerlukan dua kali iterasi. Iterasi pertama: y_1 = 4 + 0.5×(-0.2×sqrt(4)) = 4 – 0.5×0.4 = 3.8. Iterasi kedua menggunakan y_1 = 3.8, bukan nilai awal lagi. Jadi y(1) setelah dua langkah adalah 3.8 – 0.5×0.2×sqrt(3.8) ≈ 3.61, sehingga y(1) ≈ 3.61. Namun dengan satu langkah langsung h=1 menghasilkan y(1) = 4 – 1×0.2×sqrt(4) = 3.6. Soal menyebut h=0.5 dan menanyakan hasil akhir, yaitu setelah dua langkah penuh dari t=0 ke t=0.5 lalu t=0.5 ke t=1.0, diperoleh 3.61. Namun pendekatan satu langkah langsung yang dimaksud dalam opsi adalah y(1) ≈ 3.6.
PT Dinamika Fluida mengembangkan model penyebaran polutan dengan galat pemotongan lokal metode Euler sebesar O(h^2). Jika lebar langkah dikurangi menjadi setengahnya, dampak terhadap galat pemotongan lokal pada setiap langkah adalah…
Galat pemotongan lokal metode Euler sebanding dengan h^2. Jika h dikurangi setengah menjadi h/2, maka galat menjadi sebanding dengan (h/2)^2 = h^2/4, yaitu seperempat dari galat sebelumnya. Ini menunjukkan bahwa pengurangan h secara linear memberikan perbaikan kuadratik pada galat lokal.
Metode Heun meningkatkan ketelitian metode Euler dengan merata-ratakan dua estimasi kemiringan. Seorang mahasiswa menerapkan metode Heun pada masalah nilai awal y' = 2x + y dengan y(0) = 1 dan h = 0.2 untuk menaksir y(0.2). Langkah pertama menghasilkan kemiringan di titik awal k_1 = 1. Setelah menghitung y* = 1 + 0.2×1 = 1.2 dan kemiringan di titik prediksi k_2 = 2(0.2) + 1.2 = 1.6, hampiran y(0.2) metode Heun adalah…
Metode Heun menghitung prediktor y* = y_i + h×f(x_i, y_i) lalu mengoreksi dengan rata-rata dua kemiringan: y_{i+1} = y_i + h×(k_1 + k_2)/2. Dengan k_1 = 1, k_2 = 1.6, dan h = 0.2, diperoleh y_1 = 1 + 0.2×(1 + 1.6)/2 = 1 + 0.2×1.3 = 1.26.
Tim insinyur membandingkan metode Euler dan Heun untuk menyelesaikan PDB yang sama dengan lebar langkah h. Mereka mengamati bahwa metode Heun memerlukan dua evaluasi fungsi per langkah, sedangkan Euler hanya satu. Namun, Heun memberikan galat pemotongan lokal yang lebih kecil. Orde galat pemotongan lokal metode Heun dibandingkan metode Euler adalah…
Metode Euler memiliki galat pemotongan lokal O(h^2) karena didasarkan pada deret Taylor orde satu. Metode Heun termasuk dalam keluarga Runge-Kutta orde dua sehingga galat pemotongan lokalnya O(h^3), lebih kecil satu orde dibanding Euler untuk h yang sama.
Dalam merancang kontroler untuk robot lengan, tim mekatronika memilih metode Runge-Kutta orde empat (RK4) untuk menyimulasikan dinamika nonlinear. Mereka mempertimbangkan trade-off antara akurasi dan beban komputasi. Dibandingkan dengan RK orde dua, keunggulan utama RK4 yang membenarkan beban komputasi empat kali lebih banyak per langkah adalah…
RK4 memiliki galat pemotongan lokal O(h^5), sedangkan RK2 hanya O(h^3). Dengan orde ketelitian yang lebih tinggi, RK4 mencapai akurasi tertentu dengan h yang jauh lebih besar, sehingga jumlah langkah total bisa lebih sedikit dan mengompensasi biaya empat evaluasi fungsi per langkah.
Seorang analis numerik mengimplementasikan metode Runge-Kutta orde empat klasik untuk menyelesaikan y' = f(x, y). Pada setiap langkah, ia menghitung empat kemiringan: k_1, k_2, k_3, dan k_4. Jika k_1 = 0.4, k_2 = 0.6, k_3 = 0.7, dan k_4 = 0.9 dengan h = 0.5, maka pembaruan solusi dari y_i ke y_{i+1} menggunakan kombinasi linear dengan bobot…
Metode RK4 klasik menggunakan rata-rata terbobot dengan pola 1-2-2-1 pada empat kemiringan: y_{i+1} = y_i + h×(k_1 + 2k_2 + 2k_3 + k_4)/6. Bobot ini diturunkan agar sesuai dengan ekspansi deret Taylor hingga orde h^4.
PT Robotika Nusantara mensimulasikan gerak pendulum terkopel dengan metode Runge-Kutta adaptif. Saat solusi berubah cepat, lebar langkah diperkecil otomatis, dan saat perubahan lambat, langkah diperbesar. Strategi kontrol langkah adaptif ini menggunakan estimasi galat yang diperoleh melalui…
Kontrol langkah adaptif pada metode Runge-Kutta umumnya menggunakan pasangan metode terkait (embedded RK) dengan orde berbeda, misalnya RK4 dan RK5. Estimasi galat diperoleh dari selisih dua hampiran tersebut, lalu digunakan untuk menyesuaikan lebar langkah agar galat tetap dalam toleransi yang diinginkan.
Mahasiswa teknik membandingkan kinerja metode Euler, Heun, dan RK4 pada masalah nilai awal yang sama dengan h = 0.1. Setelah 20 langkah, galat global metode Euler sekitar 10^{-1}, Heun sekitar 10^{-3}, dan RK4 sekitar 10^{-7}. Pola hubungan antara orde metode dan akumulasi galat global yang tepat adalah…
Galat global metode numerik satu langkah sebanding dengan h^{p} di mana p adalah orde metode. Euler (p=1) → galat global O(h), Heun (p=2) → O(h^2), RK4 (p=4) → O(h^4). Dengan h=0.1, perbandingan galat kira-kira 10^{-1} : 10^{-2} : 10^{-4} konsisten dengan pola tersebut, menunjukkan bahwa semakin tinggi orde, semakin cepat galat mengecil seiring pengurangan h.
PT Dinamika Struktur mengembangkan model numerik untuk respons seismik yang memerlukan metode multi langkah eksplisit. Tim memilih metode yang menggunakan nilai fungsi dari empat titik sebelumnya untuk menaksir nilai baru tanpa melibatkan nilai fungsi di titik target. Metode yang mereka pilih adalah…
Metode Adams-Bashforth merupakan metode multi langkah eksplisit karena menggunakan beberapa titik sebelumnya dan tidak melibatkan nilai fungsi di titik target yang belum diketahui.
Mahasiswa teknik menerapkan metode Adams-Moulton pada simulasi getaran. Setiap langkah iterasi, ia harus menyelesaikan persamaan implisit karena metode ini melibatkan nilai fungsi di titik target. Konsekuensi utama sifat implisit ini adalah…
Karena Adams-Moulton bersifat implisit yang melibatkan f(t_{n+1}, y_{n+1}) yang belum diketahui, diperlukan iterasi internal untuk menyelesaikan persamaan pada setiap langkah.
Tim simulasi penerbangan menggabungkan metode Adams-Bashforth sebagai prediktor dan Adams-Moulton sebagai korektor. Setelah memperoleh nilai prediksi dari Adams-Bashforth, korektor Adams-Moulton diterapkan. Pendekatan gabungan ini dikenal sebagai…
Kombinasi metode eksplisit sebagai prediktor dan metode implisit sebagai korektor dalam setiap langkah disebut metode prediktor-korektor.
Seorang analis membandingkan dua metode multi langkah: metode pertama menghitung y_{n+1} langsung dari data yang tersedia, sedangkan metode kedua memerlukan penyelesaian yang melibatkan f(t_{n+1}, y_{n+1}). Perbedaan fundamental ini mengklasifikasikan kedua metode secara berturut-turut sebagai…
Metode yang menghitung y_{n+1} langsung dari data tersedia bersifat eksplisit, sementara yang memerlukan f(t_{n+1}, y_{n+1}) bersifat implisit karena melibatkan nilai yang belum diketahui.
PT Rekayasa Hidro menggunakan metode Milne untuk menyelesaikan model aliran air tanah. Metode Milne menggunakan integrasi numerik berbasis aturan Simpson pada interval sebelumnya sehingga memiliki karakteristik tertentu. Dibandingkan metode Adams-Bashforth dengan orde yang sama, metode Milne memiliki…
Metode Milne didasarkan pada aturan Simpson yang memiliki koefisien galat lebih kecil, sehingga untuk orde yang sama sering memberikan akurasi lebih baik dibanding Adams-Bashforth.
Seorang mahasiswa menjalankan metode Milne untuk PDB dengan langkah waktu tertentu. Setelah seratus iterasi, ia mengamati galat membesar secara eksponensial meskipun langkah waktu tergolong kecil. Fenomena ini mengindikasikan masalah pada…
Galat yang membesar secara eksponensial menunjukkan ketidakstabilan numerik, yaitu sifat metode yang menentukan apakah galat kecil akan meredam atau justru tumbuh tak terkendali selama iterasi.
Tim riset menerapkan metode multi langkah untuk pertama kali pada PDB baru. Sebelum memulai iterasi utama, mereka harus membangkitkan beberapa titik awal menggunakan metode satu langkah. Kebutuhan ini muncul karena…
Metode multi langkah memerlukan beberapa titik sebelumnya untuk memulai komputasi, sehingga titik-titik awal harus dibangkitkan terlebih dahulu menggunakan metode satu langkah seperti Euler atau Runge-Kutta.
Budi membandingkan metode Milne dan Adams-Bashforth orde empat untuk PDB yang sama. Ia mengamati bahwa Adams-Bashforth lebih stabil untuk PDB dengan karakteristik kaku (stiff), sementara Milne cenderung tidak stabil. Perbedaan ini terutama disebabkan oleh…
Stabilitas metode multi langkah ditentukan oleh letak akar polinomial karakteristiknya pada bidang kompleks. Metode Milne memiliki akar pada lingkaran satuan yang membuatnya kurang stabil untuk masalah kaku.
Seorang insinyur sipil menghadapi PDB orde dua dengan syarat batas pada dua ujung balok: lendutan di kedua tumpuan bernilai nol. Ia memilih metode tembakan linear dengan menebak kemiringan awal dan mengintegralkan sebagai masalah nilai awal. Untuk memenuhi syarat batas di ujung kanan, ia memanfaatkan prinsip…
Metode tembakan linear menggunakan prinsip superposisi, yaitu mengombinasikan dua solusi bebas linear dari PDB homogen dan solusi partikular untuk memenuhi syarat batas di kedua ujung.
PT Konstruksi Tahan Gempa memodelkan kolom dengan PDB linear dan syarat batas di dua titik. Tim menggunakan metode tembakan linear. Jika solusi homogen pertama menghasilkan nilai di ujung kanan y_1(b) = 2, solusi homogen kedua y_2(b) = 5, dan solusi partikular y_p(b) = 3, dengan syarat batas di ujung kanan y(b) = 7, maka koefisien kombinasi linear yang tepat adalah…
Dari persamaan c1×2 + c2×5 + 3 = 7, diperoleh 2c1 + 5c2 = 4. Dengan syarat batas kiri yang memberikan hubungan lain, koefisien c1 = 2 dan c2 = 1 memenuhi sistem persamaan yang terbentuk.
Mahasiswa menyelesaikan masalah syarat batas dua titik untuk PDB linear orde dua. Berbeda dengan masalah nilai awal yang semua kondisinya diberikan di satu titik, masalah nilai batas memiliki karakteristik yang lebih kompleks karena…
Masalah nilai batas memiliki kondisi yang diberikan pada dua titik berbeda dalam domain, sehingga tidak dapat langsung diintegralkan sebagai masalah nilai awal dan memerlukan pendekatan khusus.
Siti menerapkan metode tembakan untuk PDL linear dengan menebak dua nilai awal yang berbeda. Solusi dari kedua tebakan diintegralkan ke ujung kanan, lalu selisih antara hasil dan syarat batas digunakan untuk menentukan tebakan yang benar. Proses penentuan tebakan akhir ini menggunakan prinsip…
Dengan dua tebakan dan dua hasil di ujung kanan, metode tembakan linear menggunakan interpolasi linear untuk menentukan tebakan yang tepat sehingga syarat batas di ujung kanan terpenuhi.
Tim simulasi PT Aliran Dinamik menghadapi PDB nonlinear dengan syarat batas dua titik: y'' + y×y' = 0 dengan y(0) = 1 dan y(2) = 3. Mereka memilih metode tembakan nonlinear yang menebak y'(0) dan mengintegralkan sebagai masalah nilai awal berulang kali. Untuk memperbaiki tebakan y'(0) secara sistematis, mereka menggunakan…
Metode tembakan nonlinear memperbaiki tebakan syarat awal secara iteratif menggunakan metode pencarian akar seperti metode bagi dua atau Newton-Raphson pada selisih antara hasil tembakan di ujung kanan dengan syarat batas target.
Seorang analis menerapkan metode beda hingga pada PDB nonlinear dengan syarat batas. Ia mendiskretisasi domain menjadi n titik dan mengganti turunan dengan beda hingga. Hasilnya adalah sistem persamaan nonlinear berukuran (n-2)×(n-2) yang harus diselesaikan. Untuk menyelesaikan sistem tersebut, pendekatan yang paling sesuai adalah…
Sistem persamaan nonlinear hasil diskretisasi beda hingga diselesaikan dengan iterasi Newton, yaitu linearisasi sistem pada setiap iterasi menggunakan matriks Jacobian dan penyelesaian sistem linear secara berulang.
PT Mekanika Komputasi membandingkan dua metode untuk PDB nonlinear dengan syarat batas. Metode pertama menebak syarat awal dan mengintegralkan, sedangkan metode kedua mendiskretisasi seluruh domain sekaligus. Metode pertama cocok bila domain tidak terlalu panjang dan solusi tidak sensitif terhadap tebakan. Perbandingan ini mengacu pada…
Metode tembakan nonlinear menyelesaikan PDB sebagai masalah nilai awal berulang dengan tebakan syarat awal yang diperbaiki, sedangkan metode beda hingga mendiskretisasi domain dan membentuk sistem persamaan nonlinear yang diselesaikan sekaligus.
Budi menyelesaikan PDB nonlinear y'' + sin(y) = x dengan syarat batas y(0) = 0 dan y(π) = 1 menggunakan beda hingga. Ia memperoleh sistem F(y) = 0 dan menerapkan iterasi Newton. Pada setiap iterasi, ia harus menyelesaikan J×Δy = -F, dengan J adalah matriks Jacobian. Elemen matriks J diperoleh dari…
Matriks Jacobian pada iterasi Newton untuk sistem nonlinear berisi turunan parsial fungsi F terhadap setiap variabel y di titik-titik diskretisasi, yang dihitung dari persamaan beda hingga yang digunakan.
Dalam menyelesaikan masalah syarat batas nonlinear y'' = f(x,y,y') dengan metode beda hingga, sistem persamaan yang dihasilkan dari diskretisasi domain tersebut bersifat…
Ketika PDB nonlinear, seperti y'' = f(x,y,y') dengan f mengandung suku nonlinear terhadap y atau y', didiskretisasi menggunakan beda hingga, turunan digantikan dengan hampiran beda hingga yang melibatkan nilai-nilai y pada titik-titik diskret. Substitusi ini menghasilkan sistem persamaan aljabar yang nonlinear terhadap variabel-variabel y_i yang tidak diketahui. Oleh karena itu, metode iteratif seperti metode Newton untuk sistem harus digunakan untuk menyelesaikannya.
Seorang insinyur menyelesaikan sistem persamaan linear Ax = b menggunakan metode dekomposisi LU. Setelah memperoleh matriks L dan U, ia melakukan substitusi maju Ly = b untuk mendapatkan y. Jika matriks L memiliki elemen diagonal utama l, manakah operasi yang tepat untuk menghitung y_i…
Pada tahap substitusi maju dalam metode dekomposisi LU, vektor antara y dihitung dengan memanfaatkan struktur segitiga bawah dari matriks L. Langkah ini dimulai dari baris pertama hingga baris terakhir secara berurutan. Untuk setiap baris i, nilai y_i diperoleh dengan mengurangkan jumlah dari perkalian elemen-elemen L pada baris tersebut dengan elemen-elemen y yang telah dihitung sebelumnya, kemudian membaginya dengan elemen diagonal dari L pada baris tersebut.
Galat dan akar persamaan kelihatannya enteng. Sampai kamu masuk ke metode Bairstow. Di situ mahasiswa UT paling sering nyerah. Bukan karena rumusnya sulit. Tapi karena perhitungan iteratifnya panjang dan gampang melenceng kalau asumsi awalnya nggak tepat. Satu langkah salah, hasil akhirnya bisa beda jauh.
Nah, di STMA4223 Metode Numerik, soal UO biasanya minta kamu membandingkan laju konvergensi dua metode. Ini beda sama UTM yang lebih ke hitungan manual. Kalau udah paham kenapa Gauss-Seidel lebih cepat dari Jacobi, separuh perang udah selesai. Mau latihan yang serupa? Coba tengok soal UAS UT MSIM4304 buat mengasah logika teknis lagi.
