Metode koefisien tak tentu

Persamaan diferensial
Persamaan diferensial Navier–Stokes digunakan untuk mensimulasikan aliran udara di sekitar suatu rintangan
Persamaan diferensial Navier–Stokes digunakan untuk mensimulasikan aliran udara di sekitar suatu rintangan
Ruang lingkup
Bidang
  • Ilmu alam
  • Perekayasaan
  • Astronomi
  • Fisika
  • Kimia

  • Biologi
  • Geologi
Matematika terapan
  • Mekanika kontinuum
  • Teori kekacauan
  • Sistem dinamis
Ilmu sosial
  • Ekonomi
  • Dinamika populasi

Daftar persamaan diferensial bernama
Klasifikasi
Jenis
  • Biasa
  • Linear
  • Non-linear
  • Parsial
  • Pecahan
  • Diferensial-aljabar
  • Integro-diferensial
Berdasarkan jenis variabel
  • Variabel bebas dan terikat
  • Autonom
  • Coupled / Decoupled
  • Eksak
  • Homogen / Non-homogen
Sifat-sifat
Hubungan dengan proses
  • Persamaan beda (versi diskret)
  • Stokastik
  • Stokastik parsial
  • Tunda
Solusi
Kewujudan dan ketunggalan
  • Teorema Picard–Lindelöf
  • Teorema kewujudan Peano
  • Teorema kewujudan Carathéodory
  • Teorema Cauchy–Kowalevski
Topik umum
Metode menemukan solusi
  • Inspeksi
  • Metode karakteristik

  • Euler
  • Rumus respons eksponen
  • Beda hingga (Crank–Nicolson)
  • Elemen hingga
    • Elemen takhingga
  • Volume hingga
  • Galerkin
    • Petrov–Galerkin
  • Fungsi Green
  • Faktor integrasi
  • Transformasi integral
  • Teori perturbasi
  • Runge–Kutta
Tokoh
Daftar
  • l
  • b
  • s

Dalam matematika, metode koefisien tak tentu adalah suatu pendekatan untuk mencari solusi khusus dari suatu persamaan diferensial biasa nonhomogen dan relasi perulangan nonhomogen.

Metode koefisien tak tentu tidak seumum metode variasi parameter, sebab metode ini hanya berlaku untuk persamaan diferensial yang memiliki bentuk tertentu. Untuk persamaan yang kompleks, pencarian solusi menggunakan metode variasi parameter akan memakan waktu yang lebih sedikit.[1]

Deskripsi Metode

Perhatikan persamaan diferensial biasa linear nonhomogen dengan bentuk umum sebagai berikut

c n y ( n ) + c n 1 y ( n 1 ) + + c 2 y + c 1 y + c 0 y = F ( x ) {\displaystyle c_{n}\,y^{(n)}+c_{n-1}\,y^{(n-1)}+\,\ldots \,+c_{2}\,y''+c_{1}\,y'+c_{0}\,y=F(x)}
dimana

  • y ( i ) {\displaystyle y^{(i)}} menyatakan turunan ke- i {\displaystyle i} dari y {\displaystyle y}
  • c i {\displaystyle c_{i}} adalah suatu konstanta
  • c n 0 {\displaystyle c_{n}\neq 0}

Metode koefisien tak tentu menyediakan cara untuk memperoleh solusi dari PDB ini apabila fungsi F ( x ) {\displaystyle F(x)} merupakan:[2]

  1. fungsi polinomial
  2. fungsi eksponensial (dengan bentuk umum e m x {\displaystyle e^{mx}} )
  3. fungsi sinus dan kosinus (dengan bentuk umum sin ( ω x ) {\displaystyle \sin(\omega x)} atau cos ( ω x ) {\displaystyle \cos(\omega x)} )
  4. hasil perkalian dan penjumlahan berhingga dari (1), (2), dan (3) (misalnya e x + x 3 cos ( 2 x ) {\displaystyle e^{x}+x^{3}\,\cos(2x)} )

untuk suatu konstanta m {\displaystyle m} dan ω {\displaystyle \omega }

Misalkan y = g ( x ) {\displaystyle y=g(x)} adalah solusi persamaan diferensial linear yang ruas kanannya adalah G ( x ) {\displaystyle G(x)} , dan misalkan y = h ( x ) {\displaystyle y=h(x)} adalah solusi persamaan diferensial linear yang ruas kanannya adalah H ( x ) {\displaystyle H(x)} . Secara matematis, maka

  • c n d n g d x n + c n 1 d n 1 g d x n 1 + + c 2 d 2 g d x 2 + c 1 d g d x + c 0 g = G ( x ) {\displaystyle c_{n}\,{\dfrac {{\text{d}}^{n}g}{{\text{d}}x^{n}}}+c_{n-1}\,{\dfrac {{\text{d}}^{n-1}g}{{\text{d}}x^{n-1}}}+\,\ldots \,+c_{2}\,{\dfrac {{\text{d}}^{2}g}{{\text{d}}x^{2}}}+c_{1}\,{\dfrac {{\text{d}}g}{{\text{d}}x}}+c_{0}\,g=G(x)}
  • c n d n h d x n + c n 1 d n 1 h d x n 1 + + c 2 d 2 h d x 2 + c 1 d h d x + c 0 h = H ( x ) {\displaystyle c_{n}\,{\dfrac {{\text{d}}^{n}h}{{\text{d}}x^{n}}}+c_{n-1}\,{\dfrac {{\text{d}}^{n-1}h}{{\text{d}}x^{n-1}}}+\,\ldots \,+c_{2}\,{\dfrac {{\text{d}}^{2}h}{{\text{d}}x^{2}}}+c_{1}\,{\dfrac {{\text{d}}h}{{\text{d}}x}}+c_{0}\,h=H(x)}

Oleh karena operator diferensial bersifat linier, maka dengan menggunakan asas superposisi, didapatkan[3]

c n d n g d x n + c n 1 d n 1 g d x n 1 + + c 2 d 2 g d x 2 + c 1 d g d x + c 0 g = G ( x ) c n d n h d x n + c n 1 d n 1 h d x n 1 + + c 2 d 2 h d x 2 + c 1 d h d x + c 0 h = H ( x ) + c n d n ( g + h ) d x n + c n 1 d n 1 ( g + h ) d x n 1 + + c 2 d 2 ( g + h ) d x 2 + c 1 d ( g + h ) d x + c 0 ( g + h ) = G ( x ) + H ( x ) {\displaystyle {\begin{array}{cccccccccccclr}c_{n}\,{\dfrac {{\text{d}}^{n}g}{{\text{d}}x^{n}}}&+&c_{n-1}\,{\dfrac {{\text{d}}^{n-1}g}{{\text{d}}x^{n-1}}}&+&\ldots &+&c_{2}\,{\dfrac {{\text{d}}^{2}g}{{\text{d}}x^{2}}}&+&c_{1}\,{\dfrac {{\text{d}}g}{{\text{d}}x}}&+&c_{0}\,g&=&G(x)&\\c_{n}\,{\dfrac {{\text{d}}^{n}h}{{\text{d}}x^{n}}}&+&c_{n-1}\,{\dfrac {{\text{d}}^{n-1}h}{{\text{d}}x^{n-1}}}&+&\ldots &+&c_{2}\,{\dfrac {{\text{d}}^{2}h}{{\text{d}}x^{2}}}&+&c_{1}\,{\dfrac {{\text{d}}h}{{\text{d}}x}}&+&c_{0}\,h&=&H(x)&+\\\hline c_{n}\,{\dfrac {{\text{d}}^{n}(g+h)}{{\text{d}}x^{n}}}&+&c_{n-1}\,{\dfrac {{\text{d}}^{n-1}(g+h)}{{\text{d}}x^{n-1}}}&+&\ldots &+&c_{2}\,{\dfrac {{\text{d}}^{2}(g+h)}{{\text{d}}x^{2}}}&+&c_{1}\,{\dfrac {{\text{d}}(g+h)}{{\text{d}}x}}&+&c_{0}\,(g+h)&=&G(x)+H(x)&\\\end{array}}}

Dengan kata lain, jika fungsi F ( x ) {\displaystyle F(x)} di ruas kanan dapat dinyatakan sebagai G ( x ) + H ( x ) {\displaystyle G(x)+H(x)} , maka solusi akhirnya adalah jumlahan dari masing-masing solusi, yaitu y = g ( x ) + h ( x ) {\displaystyle y=g(x)+h(x)} . Jika G ( x ) = 0 {\displaystyle G(x)=0} (yang mengakibatkan H ( x ) = F ( x ) {\displaystyle H(x)=F(x)} ), maka fungsi g ( x ) {\displaystyle g(x)} disebut sebagai solusi umum, dan h ( x ) {\displaystyle h(x)} disebut sebagai solusi khusus.[butuh rujukan]

Metode ini secara umum terdiri dari dua bagian, yaitu:

  1. Pencarian solusi umum, yaitu suatu fungsi y u {\displaystyle y_{u}} sedemikian sehingga fungsi tersebut memenuhi
    c n y ( n ) + c n 1 y ( n 1 ) + + c 2 y + c 1 y + c 0 y = 0 {\displaystyle c_{n}\,y^{(n)}+c_{n-1}\,y^{(n-1)}+\,\ldots \,+c_{2}\,y''+c_{1}\,y'+c_{0}\,y=0}
    Dengan kata lain, persamaan diferensialnya dipandang sebagai persamaan diferensial homogen terlebih dahulu.
  2. Pencarian solusi khusus, yaitu suatu fungsi y k {\displaystyle y_{k}} sedemikian sehingga fungsi tersebut memenuhi
    c n y ( n ) + c n 1 y ( n 1 ) + + c 2 y + c 1 y + c 0 y = g ( x ) {\displaystyle c_{n}\,y^{(n)}+c_{n-1}\,y^{(n-1)}+\,\ldots \,+c_{2}\,y''+c_{1}\,y'+c_{0}\,y=g(x)}

Setelah diperoleh y u {\displaystyle y_{u}} dan y k {\displaystyle y_{k}} , maka solusi akhir dari persamaan diferensialnya ialah

y = y u + y k {\displaystyle y=y_{u}+y_{k}} [3]

Bentuk umum beserta solusinya

Untuk mendapatkan solusi persamaan diferensialnya, maka terlebih dahulu harus 'ditebak' bentuk umumnya, yang nantinya beberapa koefisien yang ada akan menjadi variabel, yang kemudian akan dicari nilainya. Berikut adalah beberapa jenis fungsi beserta bentuk umum solusinya.

Bentuk umum F ( x ) {\displaystyle F(x)} Bentuk umum dari y {\displaystyle y}
k e m x {\displaystyle ke^{mx}} c e m x {\displaystyle ce^{mx}}
k x n n = 0 , 1 , 2 , {\displaystyle kx^{n}\qquad n=0,\,1,\,2,\,\ldots } a 0 + a 1 x + a 2 x 2 + + a n x n {\displaystyle a_{0}+a_{1}x+a_{2}x^{2}+\,\ldots \,+a_{n}x^{n}}
k cos ( ω x ) {\displaystyle k\cos(\omega x)} λ 1 cos ( ω x ) + λ 2 sin ( ω x ) {\displaystyle \lambda _{1}\cos(\omega x)+\lambda _{2}\sin(\omega x)}
k sin ( ω x ) {\displaystyle k\sin(\omega x)}
k e m x cos ( ω x ) {\displaystyle ke^{mx}\cos(\omega x)} λ 1 e m x cos ( ω x ) + λ 2 e m x sin ( ω x ) {\displaystyle \lambda _{1}e^{mx}\cos(\omega x)+\lambda _{2}e^{mx}\sin(\omega x)}
k e m x sin ( ω x ) {\displaystyle ke^{mx}\sin(\omega x)}
( k 0 + k 1 x + k 2 x 2 + + k n x n ) sin ( ω x ) {\displaystyle (k_{0}+k_{1}x+k_{2}x^{2}+\,\ldots \,+k_{n}x^{n})\,\sin(\omega x)} ( a 0 + a 1 x + a 2 x 2 + + a n x n ) cos ( ω x ) {\displaystyle (a_{0}+a_{1}x+a_{2}x^{2}+\,\ldots \,+a_{n}x^{n})\,\cos(\omega x)}
+ ( b 0 + b 1 x + b 2 x 2 + + b n x n ) sin ( ω x ) {\displaystyle +\,(b_{0}+b_{1}x+b_{2}x^{2}+\,\ldots \,+b_{n}x^{n})\,\sin(\omega x)}
( k 0 + k 1 x + k 2 x 2 + + k n x n ) cos ( ω x ) {\displaystyle (k_{0}+k_{1}x+k_{2}x^{2}+\,\ldots \,+k_{n}x^{n})\,\cos(\omega x)}
( k 0 + k 1 x + k 2 x 2 + + k n x n ) e m x {\displaystyle (k_{0}+k_{1}x+k_{2}x^{2}+\,\ldots \,+k_{n}x^{n})\,e^{mx}} ( a 0 + a 1 x + a 2 x 2 + + a n x n ) e m x {\displaystyle (a_{0}+a_{1}x+a_{2}x^{2}+\,\ldots \,+a_{n}x^{n})\,e^{mx}}
( k 0 + k 1 x + k 2 x 2 + + k n x n ) e m x cos ( ω x ) {\displaystyle (k_{0}+k_{1}x+k_{2}x^{2}+\,\ldots \,+k_{n}x^{n})e^{mx}\cos(\omega x)} ( a 0 + a 1 x + a 2 x 2 + + a n x n ) e m x cos ( ω x ) {\displaystyle (a_{0}+a_{1}x+a_{2}x^{2}+\,\ldots \,+a_{n}x^{n})e^{mx}\cos(\omega x)}
+ ( b 0 + b 1 x + b 2 x 2 + + b n x n ) e m x sin ( ω x ) {\displaystyle +\,(b_{0}+b_{1}x+b_{2}x^{2}+\,\ldots \,+b_{n}x^{n})e^{mx}\sin(\omega x)}
( k 0 + k 1 x + k 2 x 2 + + k n x n ) e m x sin ( ω x ) {\displaystyle (k_{0}+k_{1}x+k_{2}x^{2}+\,\ldots \,+k_{n}x^{n})e^{mx}\sin(\omega x)}

Jika salah satu suku pada bentuk umum dari y {\displaystyle y} muncul pada solusi homogen, maka bentuk umumnya harus dikalikan dengan perpangkatan x {\displaystyle x} yang cukup besar agar solusinya menjadi bebas linier.[1]

Contoh pengerjaan

Contoh 1

Untuk mencari solusi dari persamaan diferensial

d y d t + y = t 3 {\displaystyle {\dfrac {{\text{d}}y}{{\text{d}}t}}+y=t^{3}}

maka perhatikan bahwa t 3 {\displaystyle t^{3}} adalah fungsi polinomial berderajat 3, sehingga solusi khususnya juga merupakan fungsi polinomial berderajat 3, dengan bentuk umum

y k = a t 3 + b t 2 + c t + d {\displaystyle y_{k}=at^{3}+bt^{2}+ct+d}

yang mengakibatkan

d y d t = 3 a t 2 + 2 b t + c {\displaystyle {\dfrac {{\text{d}}y}{{\text{d}}t}}=3at^{2}+2bt+c}

Substitusikan hasil di atas pada persamaan diferensial di awal, maka didapatkan

d y d t + y = t 3 ( 3 a t 2 + 2 b t + c ) + ( a t 3 + b t 2 + c t + d ) = ( 1 ) t 3 + ( 0 ) t 2 + ( 0 ) t + 0 ( a ) t 3 + ( 3 a + b ) t 2 + ( 2 b + c ) t + ( c + d ) = ( 1 ) t 3 + ( 0 ) t 2 + ( 0 ) t + 0 {\displaystyle {\begin{aligned}{\dfrac {{\text{d}}y}{{\text{d}}t}}+y&=t^{3}\\\left(3at^{2}+2bt+c\right)+\left(at^{3}+bt^{2}+ct+d\right)&=(1)t^{3}+(0)t^{2}+(0)t+0\\(a)t^{3}+(3a+b)t^{2}+(2b+c)t+(c+d)&=(1)t^{3}+(0)t^{2}+(0)t+0\\\end{aligned}}}

sehingga diperoleh sistem persamaan

a = 1 3 a + b = 0 2 b + c = 0 c + d = 0 {\displaystyle {\begin{array}{rcrcrcrcl}a&&&&&&&=&1\\3a&+&b&&&&&=&0\\&&2b&+&c&&&=&0\\&&&&c&+&d&=&0\end{array}}}

yang solusinya ialah [ a b c d ] = [ 1 3 6 6 ] {\displaystyle {\begin{bmatrix}a\\b\\c\\d\end{bmatrix}}={\begin{bmatrix}1\\-3\\6\\-6\end{bmatrix}}} . Sehingga, didapatkan

y k = t 3 3 t 2 + 6 t 6 {\displaystyle y_{k}=t^{3}-3t^{2}+6t-6}

Oleh karena solusi umum dari persamaan diferensial

d y d t + y = 0 {\displaystyle {\dfrac {{\text{d}}y}{{\text{d}}t}}+y=0}

adalah y u = k e t {\displaystyle y_{u}=ke^{-t}} , untuk sembarang konstanta k {\displaystyle k} , maka solusi akhir dari persamaan diferensial

d y d t + y = t 3 {\displaystyle {\dfrac {{\text{d}}y}{{\text{d}}t}}+y=t^{3}}

adalah y = k e t + t 3 3 t 2 + 6 t 6 {\displaystyle y=ke^{-t}+t^{3}-3t^{2}+6t-6}

Contoh 2

Perhatikan persamaan diferensial linier nonhomogen berikut:

y + 4 y = cos ( t ) {\displaystyle y''+4y=\cos(t)}

Oleh karena bagian nonhomogen dari persamaan di atas adalah cos ( t ) {\displaystyle \cos(t)} , maka solusi khususnya akan memiliki bentuk umum

y k = a cos ( t ) + b sin ( t ) {\displaystyle y_{k}=a\cos(t)+b\sin(t)}

Substitusikan bentuk di atas ke persamaan diferensial di awal, maka didapatkan

y + 4 y = cos ( t ) ( a cos ( t ) b sin ( t ) ) + 4 ( a cos ( t ) + b sin ( t ) ) = cos ( t ) 3 a cos ( t ) + 3 b sin ( t ) = ( 1 ) cos ( t ) + ( 0 ) sin ( t ) {\displaystyle {\begin{aligned}y''+4y&=\cos(t)\\\left(-a\cos(t)-b\sin(t)\right)+4\left(a\cos(t)+b\sin(t)\right)&=\cos(t)\\3a\cos(t)+3b\sin(t)&=(1)\cos(t)+(0)\sin(t)\end{aligned}}}

Dengan membandingkan koefisien pada kedua ruas, maka diperoleh

3 a = 1 a = 1 3 3 b = 0 b = 0 {\displaystyle {\begin{aligned}3a&=1\\a&={\tfrac {1}{3}}\\\\3b&=0\\b&=0\end{aligned}}}

Sehingga, solusi khusus dari persamaan diferensial tersebut adalah

y k = 1 3 cos ( t ) {\displaystyle y_{k}={\frac {1}{3}}\cos(t)}

Dengan menggunakan informasi bahwa y u = p cos ( 2 t ) + q sin ( 2 t ) {\displaystyle y_{u}=p\cos(2t)+q\sin(2t)} adalah solusi dari persamaan diferensial linier homogen

y + 4 y = 0 {\displaystyle y''+4y=0}

untuk sembarang konstanta p {\displaystyle p} dan q {\displaystyle q} , maka solusi akhir dari persamaan diferensialnya ialah

y = p cos ( 2 t ) + q sin ( 2 t ) + 1 3 cos ( t ) {\displaystyle y=p\cos(2t)+q\sin(2t)+{\frac {1}{3}}\cos(t)}

Contoh 3

Untuk mencari solusi dari persamaan diferensial linier nonhomogen

d y d x 5 y = 6 e 5 x {\displaystyle {\dfrac {{\text{d}}y}{{\text{d}}x}}-5y=6e^{5x}}

maka perhatikan bahwa y = c e 5 x {\displaystyle y=ce^{5x}} adalah solusi umum dari persamaan diferensial linier homogen

d y d x 5 y = 0 {\displaystyle {\dfrac {{\text{d}}y}{{\text{d}}x}}-5y=0}

untuk sembarang konstanta c {\displaystyle c} . Akan tetapi, fungsi e 5 x {\displaystyle e^{5x}} juga muncul pada bagian nonhomogen dari persamaan diferensial yang diberikan (bagian ruas kanan), yang membuat solusi umumnya tidak bebas linier dengan bentuk umum solusi khususnya (yaitu y = λ e 5 x {\displaystyle y=\lambda e^{5x}} ). Alhasil, bentuk umum dari solusi khususnya harus dikalikan dengan perpangkatan x {\displaystyle x} yang cukup besar agar solusinya menjadi bebas linier. Dalam kasus ini, bentuk umum solusi khususnya menjadi

y k = λ x e 5 x {\displaystyle y_{k}=\lambda xe^{5x}}

Apabila fungsi tersebut (beserta turunannya) disubstitusikan ke persamaan diferensial yang diberikan, maka nilai λ {\displaystyle \lambda } dapat diperoleh sebagai berikut

d y d x 5 y = 6 e 5 x ( λ e 5 x + 5 λ x e 5 x ) 5 ( λ x e 5 x ) = 6 e 5 x λ e 5 x = 6 e 5 x λ = 6 {\displaystyle {\begin{aligned}{\dfrac {{\text{d}}y}{{\text{d}}x}}-5y&=6e^{5x}\\\left(\lambda e^{5x}+5\lambda xe^{5x}\right)-5\left(\lambda xe^{5x}\right)&=6e^{5x}\\\lambda e^{5x}&=6e^{5x}\\\lambda &=6\end{aligned}}}

Maka dari itu, solusi akhir dari persamaan diferensialnya ialah

y = c e 5 x + 6 x e 5 x {\displaystyle y=ce^{5x}+6xe^{5x}}

Referensi

  1. ^ a b Ralph P. Grimaldi (2000). "Nonhomogeneous Recurrence Relations". Section 3.3.3 of Handbook of Discrete and Combinatorial Mathematics. Kenneth H. Rosen, ed. CRC Press. ISBN 0-8493-0149-1.
  2. ^ Zill, Dennis G.; Wright, Warren S. (2014). Advanced Engineering Mathematics [Matematika Teknik Lanjut] (dalam bahasa Inggris). Jones and Bartlett. hlm. 125. ISBN 978-1-4496-7977-4. 
  3. ^ a b Dennis G. Zill (14 May 2008). A First Course in Differential Equations [Kursus Pertama pada Persamaan Diferensial] (dalam bahasa Inggris). Cengage Learning. ISBN 978-0-495-10824-5. 

Bacaan lanjutan

  • Boyce, W. E.; DiPrima, R. C. (1986). Elementary Differential Equations and Boundary Value ProblemsPerlu mendaftar (gratis) [Persamaan Diferensial Elementer dan Masalah Nilai Batas] (dalam bahasa Inggris) (edisi ke-4th). John Wiley & Sons. ISBN 0-471-83824-1. 
  • Riley, K. F.; Bence, S. J. (2010). Mathematical Methods for Physics and EngineeringPerlu mendaftar (gratis) [Metode Matematis untuk Fisika dan Teknik] (dalam bahasa Inggris). Cambridge University Press. ISBN 978-0-521-86153-3. 
  • Tenenbaum, Morris; Pollard, Harry (1985). Ordinary Differential EquationsPerlu mendaftar (gratis) [Persamaan Diferensial Biasa] (dalam bahasa Inggris). Dover. ISBN 978-0-486-64940-5. 
  • de Oliveira, O. R. B. (2013). "A formula substituting the undetermined coefficients and the annihilator methods". Int. J. Math. Educ. Sci. Technol (dalam bahasa Inggris). 44 (3): 462–468. arXiv:1110.4425 alt=Dapat diakses gratis. Bibcode:2013IJMES..44..462R. doi:10.1080/0020739X.2012.714496.  Parameter |s2cid= yang tidak diketahui akan diabaikan (bantuan)

Templat:Topik persamaan diferensial