basic method

25 Jan 2022 • viridi | history

Berdasarkan kebiasaan metode-metode untuk menyelesaikan masalah dalam fisika matematika dapat dibagi menjadi dua kelas, metode analitik (analytical methods) dan metode hampiran (aproximate methods) [1]. Dalam matematika beberapa masalah dapat dipecahkan secara analitik dan numerik, dengan metoda analitik lebih disukai karena lebih cepat dan solusinya eksak [2]. Solusi eksak dalam fisika digunakan untuk solusi yang menangkap keseluruhan aspek fisika dan matematika suatu masalah [3], karena tidak menggunakan aproksimasi [4] dan merupakan representasi simbolik [5], yang bila cukup dengan jumlah data berhingga akan menjadi bentuk tertutup atau closed form [6]. Dalam komputasi matematika metode iteratif (iterative method) merupakan prosedur matematika yang menggunakan nilai awal untuk menghasilkan rangkaian hampiran solusi yang lebih baik dengan kriteria penghentian tertentu, di mana metode ini berlawanan dengan metode langsung (direct method) yang menyelesaikan masalah dengan sejumlah berhingga rangkaian operasi [7].

polynomial roots

Akar suatu persamaan polinomial dapat dicari dengan menggunakan metode langsung dan iteratif [8], dengan contoh metode langsungnya adalah metode kuadrat akar Graeffe [9] dan metode iteratifnya adalah metode bisection dan iterasi titik tetap [10]. Untuk persamaan kuadrat atau polinomial orde dua, formula kuadrat termasuk dalam metode langsung karena jumlah langkah operasi, yang meliputi operasi penjumlahan, pengurangan, pembagian dan perkalian dapat ditentukan, akan tetapi untuk polinomial dengan orde lebih tinggi metode ini sulit untuk digunakan, walapun ada [11].

linear system of equations

Sistem persamaan linier, yang merupakan terjemahan dari system of linear equations (SLE) dan linear system of equations [12], dapat dipecahkan menggunakan metode langsung seperti eliminasi Gauss dan faktorisasi LU [13] atau metode iteratif seperti metode Gauss-Seidel dan Successive Over-Relaxation (SOR) [14].

Tabel 1 Perbandingan metode langsung dan iteratif untuk pemecahan SLE.

Metode Langsung Iteratif
Langkah Tertentu bergantung
ukuran matriks
Fleksibel dengan
memilih antara akurasi
atau waktu komputasi
Waktu komputasi Dapat diperkirakan
berdasarkan jumlah
langkah
Bergantung akurasi
yang diinginkan
Tepat untuk Matriks kecil Matriks besar
dengan banyak nol
Penghentian Harus diselesaikan
atau solusi belum
diperoleh
Minimal satu siklus,
lalu sesuai akurasi
yang diinginkan

Metode langsung menggunakan sumber daya komputasi secara besar-besaran, akan tetapi amat bagus untuk menyelesaikan matriks berukuran kecil, handal dan cepat [15].

example

Sebagai contoh akan dipilih mencari satu akar dari suatu persamaan kuadrat dengan diskriminan $D > 0$.

import math

a = 1
b = -100
c = 196

D = b*b - 4*a*c
x = (-b - math.sqrt(D)) / (2*a)
print("x = ", x, sep='')

Program di atas merupakan penerapan metode langsung dengan menggunakan formula kuadrat. Kode dapat dicoba di OneCompiler 3xrdp6ggx. Bandingkan dengan kode berikut ini.

def f(x):
  a = 1
  b = -100
  c = 196
  y = a*x*x + b*x + c
  return y

def fx(x):
  a = 2
  b = -100
  y = a*x + b
  return y

err = 10
eps = 1E-20
x = 0
i = 0
print(i, ": x = ", x, sep='')

while err > eps:
  xnew = x - f(x)/fx(x)
  err = abs(xnew - x)
  x = xnew
  i = i + 1
  print(i, ": x = ", x, sep='')

Program kedua di atas merupakan penerapan dari metode iteratif yang dikenal sebagai metode Newton-Raphson. Kode dapat dicoba di OneCompiler 3xrdpfxwq.

Tabel 2 Mencari satu akar dari persamaan kuadrat dengan metode langsung dan iteratif.

Kelompok
Metode
Langsung Iteratif
Metode Formula Kuadrat Newton-Raphson
Fungsi $f(x) = ax^2 + bx + c$ $f(x) = ax^2 + bx + c$
Turunan
fungsi
Tidak perlu $f’(x) = 2ax + b$
Persamaan $\displaystyle x = \frac{-b - \sqrt{b^2 - 4ac}}{2a}$ $\displaystyle x_{i + 1} = x_i - \frac{f(x_i)}{f’(x_i)}$
Tebakan
awal
Tidak perlu $x = 0$
Keluaran x = 2.0 0: x = 0
1: x = 1.96
2: x = 1.9999833472106578
3: x = 1.9999999999971112
4: x = 1.9999999999999998
5: x = 2.0
6: x = 2.0
Batas
kesalahan
perhitungan
$\varepsilon$
Tidak perlu $10^{-20}$

Dapat terlihat secara tersirat pada Tabel 2 bahwa jumlah langkah dengan metode langsung dapat diperkirakan, akan tetapi pada metode iteratif masih bergantung pada batas kesalahan perhitungan $\varepsilon$ yang diberikan.

Coba modifikasi program yang diberikan untuk fungsi kuadrat $f(x) = x^2 - 100x + 1875$ bagaimanakah hasilnya? Jangan gunakan alat bantu seperti Google [16]. :smile:

exer

  1. Dengan mengacu pada Tabel 2 berapakah nilai akar bila digunakan $\varepsilon = 0.001$?
  2. Dengan metode iteratif pada Tabel 2 bagaimana untuk mendapatkan akar keduanya? Dan berapa akarnya? Perlu berapa iterasi? Gunakan $\varepsilon = 10^{-20}$.
  3. Dengan metode langsung pada Tabel 2 bagaimana untuk mendapatkan akar keduanya? Berapa akarnya?

note

  1. V. K. Andreev, “Basic Method for Solving Equations of Mathematical Physics”, Computational Methods and Algorithms, Volume 1, Vladimir V. Shaidurov, Olivier Pironneau (Eds.), Encyclopedia of Life Support Systems (EOLSS), url https://www.eolss.net/sample-chapters/c02/E6-04-01.pdf [20220125].
  2. Jason Brownlee, “”, Machinery Learning Mastery, 13 Apr 2018, url https://machinelearningmastery.com/analytical-vs-numerical-solutions-in-machine-learning/ [20220125].
  3. Eric W. Weisstein, “Exact Solution”, from MathWorld–A Wolfram Web Resource, Wolfram Research, 21 Jan 2022, url https://mathworld.wolfram.com/ExactSolution.html [20220125].
  4. Pengwuino, “Answer to ‘What does an exact solution to a pde mean?’”, Physics Forum, 30 Nov 2011, url https://www.physicsforums.com/threads/what-does-an-exact-solution-to-a-pde-mean.555623/#post-3643407 [20220125].
  5. Alex Jones, “Asnwer to ‘What is the main difference between the exact and numerical solution?’”, Quora, 29 Jul 2022, url https://www.quora.com/What-is-the-main-difference-between-the-exact-and-numerical-solution/answer/Alex-Jones-1227 [20220125].
  6. Mark van Hoeij, “Closed Form Solutions”, Florida State University, ISSAC 2017, 2 Aug 2017, url https://www.math.fsu.edu/~hoeij/issac2017.pdf [20220125].
  7. Wikipedia contributors, “Iterative method”, Wikipedia, The Free Encyclopedia, 3 January 2022, 05:59 UTC, url https://en.wikipedia.org/w/index.php?oldid=1063462370 [20220125].
  8. Wikipedia contributors, “Numerical analysis”, Wikipedia, The Free Encyclopedia, 27 December 2021, 12:57 UTC, url https://en.wikipedia.org/w/index.php?oldid=1062275306 [20220125].
  9. Sanyasiraju V. S. S. Yedida, “Graeffe’s Root SquaringMethod”, Transcendental Equations, Department of Mathematics, Indian Institute Of Technology Madras, Chennai, url https://math.iitm.ac.in/public_html/sryedida/caimna/transcendental/polynomial%20methods/direct.html [20220125].
  10. “Iterations Methods for Locating a Root”, Solution of Non-linear Equations in One Variable, IGNOU, Indira Gandhi National Open University, 2017, url https://www.egyankosh.ac.in/bitstream/123456789/18065/1/Unit-2.pdf [20220125].
  11. FaadooEngineers, “Methods of solving roots of polynomial equations”, FaaDoO Engineers, url http://www.faadooengineers.com/online-study/post/cse/numerical-analysis/1427/methods-of-solving-roots-of-polynomial-equations [20220125].
  12. “Systems of Linear Equation”, Science Direct, 2022, url https://www.sciencedirect.com/topics/mathematics/systems-of-linear-equation [20220125].
  13. Chen Greif, “Numerical Solution of Linear Systems”, Department of Computer Science, The University of British Columbia, Vancouver B.C. 17 Dec 2008, url https://www.cs.tau.ac.il/~dcor/Graphics/adv-slides/Solving.pdf [20220125].
  14. Yousef Saad, “Iterative methods for linear systems of equations: A brief historical journey”, arXiv:1908.01083v1 [math.HO], Fri, 2 Aug 2019 22:42:47 UTC, url https://arxiv.org/abs/1908.01083v1 [20220125].
  15. StudySession, “Direct Vs Iterative Numerical Methods”, YouTube, 28.07.2020, url https://www.youtube.com/watch?v=CHUlnCSkg5k [20220125].
  16. url https://www.google.com/search?q=solving+x%5E2+-+100x+%2B+1875&oq=solving+x%5E2+-+100x+%2B+1875&aqs=chrome..69i57.2854j0j9&sourceid=chrome&ie=UTF-8 [20220125].

comments

#bug0420

— Sparisoma Viridi (@6unpnp) January 25, 2022

 

1) $x = 1.9999999999971112$;   2) syarat awal $x = 100$ memberikan akar $x = 98.0$ dalam 5 iterasi;   3) gunakan $\displaystyle x = \frac{-b + \sqrt{b^2 - 4ac}}{2a}$, diperoleh $x = 98.0$;