Metode Karakteristik

Tentang unit ini. Unit kumulatif ini menambahkan bagian Persamaan Diferensial Parsial Orde Pertama Kuasilinear dari Bab 1 pada terjemahan Prakata yang telah diverifikasi. Urutan, rumus, bukti, contoh, dan gambar sumber dipertahankan pada batas sumber yang dipin.

Persamaan Diferensial Parsial Orde Pertama Kuasilinear

Tinjau persamaan diferensial parsial

a(x,y,u)ux+b(x,y,u)uy=c(x,y,u),a(x,y,u)\,\frac{\partial{u}}{\partial x} + b(x,y,u)\,\frac{\partial{u}}{\partial y} = c(x,y,u) \ , dengan a,b,c:3\displaystyle a,b,c:\mathbb{R}^3\rightarrow \mathbb{R} merupakan tiga fungsi kelas C1C^1 yang diberikan, sedangkan u:2\displaystyle u:\mathbb{R}^2\rightarrow \mathbb{R} merupakan fungsi terdiferensialkan yang belum diketahui.

Persamaan yang berbentuk Persamaan (1.1) disebut persamaan diferensial parsial kuasilinear . Jika aa dan bb tidak bergantung pada uu, dan c(x,y,u)=p(x,y)u+q(x,y)c(x,y,u) = p(x,y) u + q(x,y) untuk fungsi-fungsi terdiferensialkan p:2\displaystyle p:\mathbb{R}^2\rightarrow \mathbb{R} dan q:2\displaystyle q:\mathbb{R}^2\rightarrow \mathbb{R}, maka persamaan diferensial parsial tersebut disebut linear.

Tujuannya adalah mencari solusi terdiferensialkan, yang juga disebut permukaan, u=u(x,y)u=u(x,y) yang memenuhi Persamaan (1.1) dan memuat suatu kurva Γ\Gamma yang diberikan oleh Γ={(x,y,u)=(x0(s),y0(s),u0(s)):sI}\Gamma = \left\{ (x,y,u) = \left(x_0(s), y_0(s),u_0(s)\right) : s\in I \right\} untuk fungsi-fungsi terdiferensialkan x0x_0, y0y_0, dan u0u_0 yang didefinisikan pada suatu interval terbuka II. Masalah yang solusinya (jika ada) bagi persamaan diferensial parsial tersebut harus memuat suatu kurva yang diberikan disebut masalah Cauchy .

Definisi. Suatu permukaan u=u(x,y)u=u(x,y) yang memenuhi Persamaan (1.1) dan memuat kurva Γ\Gamma yang diberikan disebut permukaan integral bagi Persamaan (1.1).

Kita dapat menulis ulang Persamaan (1.1) sebagai (a,b,c)(ux,uy,1)=0\displaystyle(a,b,c) \cdot \left(\frac{\partial{u}}{\partial x}, \frac{\partial{u}}{\partial y}, -1 \right) = 0. Vektor (ux,uy,1)\displaystyle\left(\frac{\partial{u}}{\partial x}, \frac{\partial{u}}{\partial y}, -1 \right) yang dievaluasi di (x,y)(x,y) tegak lurus terhadap permukaan u=u(x,y)u=u(x,y) di (x,y,u(x,y))\displaystyle\left(x,y,u(x,y)\right). Oleh karena itu, vektor (a,b,c)(a,b,c) yang dievaluasi di (x,y,u(x,y))(x,y,u(x,y)) terletak pada bidang singgung permukaan integral u=u(x,y)u=u(x,y) di titik (x,y,u(x,y))\displaystyle\left(x,y,u(x,y)\right).

Jika γ={(x(t),y(t),u(t)):tJ}\displaystyle\gamma = \left\{ (x(t),y(t),u(t)) : t \in J \right\} merupakan kurva pada permukaan integral u=u(x,y)u=u(x,y), maka (x(t),y(t),u(t))\left(x'(t),y'(t),u'(t)\right) terletak pada bidang singgung (lebih tepatnya, sejajar dengan bidang singgung) terhadap permukaan integral di titik (x(t),y(t),u(t))\left(x(t),y(t),u(t)\right), dengan JJ suatu interval tempat kurva tersebut terdefinisi.

Dari pembahasan sebelumnya, kita perlu meninjau keluarga kurva berikut. Untuk setiap nilai tetap sIs \in I, misalkan γs\gamma_s adalah kurva yang didefinisikan oleh γs={(x,y,u)=(x(t,s),y(t,s),u(t,s)):tJs},\gamma_s = \left\{ (x,y,u) = \left(x(t,s), y(t,s), u(t,s) \right) : t \in J_s \right\} \ , dengan JsJ_s interval maksimal yang memuat 00 tempat solusi berikut terdefinisi,

xt(t,s)=a(x(t,s),y(t,s),u(t,s))yt(t,s)=b(x(t,s),y(t,s),u(t,s))ut(t,s)=c(x(t,s),y(t,s),u(t,s))\begin{split} \frac{\partial{x}}{\partial t}(t,s) &= a(x(t,s),y(t,s),u(t,s)) \\ \frac{\partial{y}}{\partial t}(t,s) &= b(x(t,s),y(t,s),u(t,s)) \\ \frac{\partial{u}}{\partial t}(t,s) &= c(x(t,s),y(t,s),u(t,s)) \end{split} dan

(x(0,s),y(0,s),u(0,s))=(x0(s),y0(s),u0(s)).\left(x(0,s), y(0,s), u(0,s)\right) = \left(x_0(s), y_0(s), u_0(s)\right) \ . Untuk ss tetap, Persamaan (1.2) merupakan sistem persamaan diferensial biasa yang bergantung pada parameter ss dan memenuhi syarat awal Persamaan (1.3). Oleh karena itu, dalam literatur, turunan parsial sering digantikan dengan turunan dxdt\displaystyle\frac{\mathrm{d}{x}}{\mathrm{d}{t}}, dydt\displaystyle\frac{\mathrm{d}{y}}{\mathrm{d}{t}}, dan dudt\displaystyle\frac{\mathrm{d}{u}}{\mathrm{d}{t}}.

Definisi. Persamaan-persamaan dalam Persamaan (1.2) disebut persamaan karakteristik dan, untuk setiap ss, kurva γs\gamma_s disebut kurva karakteristik .

Dengan demikian, kita memperoleh suatu keluarga kurva {γs}sI\displaystyle\left\{ \gamma_s \right\}_{s\in I}. Jika kita mengasumsikan bahwa

(x,y)(t,s)|(t,s)=det(xtxsytys)|(t,s)0,\frac{\partial(x,y)}{\partial(t,s)}\bigg|_{(t,s)} = \det \begin{pmatrix} \displaystyle\frac{\partial{x}}{\partial t} & \displaystyle\frac{\partial{x}}{\partial s} \\[0.7em] \displaystyle\frac{\partial{y}}{\partial t} & \displaystyle\frac{\partial{y}}{\partial s} \end{pmatrix}\bigg|_{(t,s)} \neq 0 \ , maka, menurut Teorema Pemetaan Invers, ss dan tt dapat dinyatakan (secara lokal) sebagai fungsi dari xx dan yy. Dengan mensubstitusikan ungkapan ss dan tt dalam xx dan yy tersebut ke u=u(t,s)u=u(t,s), diperoleh persamaan suatu permukaan u=u(x,y)u=u(x,y) (Gambar 1).

Untuk menyederhanakan notasi, kita menggunakan uu untuk menyatakan fungsi dari tt dan ss maupun fungsi dari xx dan yy. Konteks akan menentukan fungsi yang dimaksud. Kita mempunyai u(x,y)=u(t,s)u(x,y) = u(t,s) dengan x=x(t,s)x=x(t,s) dan y=y(t,s)y=y(t,s).

Representasi skematis kurva-kurva karakteristik dan permukaan integral.

Teorema. Misalkan Γ\Gamma adalah kurva yang didefinisikan oleh Γ={(x,y,u)=(x0(s),y0(s),u0(s)):sI},\Gamma = \left\{ (x,y,u) = \left(x_0(s), y_0(s),u_0(s)\right) : s \in I \right\} \ , dengan II suatu interval terbuka dan x0,y0,u0Cb1(I)\displaystyle x_0, y_0, u_0 \in C_b^1(I). Untuk sIs \in I, misalkan γs={(x(t,s),y(t,s),u(t,s)):tJs}\gamma_s = \left\{ \left(x(t,s), y(t,s), u(t,s) \right) : t \in J_s \right\} merupakan orbit dari solusi t(x(t,s),y(t,s),u(t,s))t \mapsto \left(x(t,s), y(t,s), u(t,s) \right) bagi Persamaan (1.2) dengan syarat awal Persamaan (1.3), dengan JsJ_s \subset \mathbb{R} merupakan domain solusi tersebut.1Jika, untuk suatu s*Is_\ast\in I, det(a(x(0,s*),y(0,s*),u(0,s*))xs(0,s*)b(x(0,s*),y(0,s*),u(0,s*))ys(0,s*))0\det \begin{pmatrix} a\left(x(0,s_\ast),y(0,s_\ast),u(0,s_\ast)\right) & \displaystyle\frac{\partial{x}}{\partial s}(0,s_\ast) \\[0.7em] b\left(x(0,s_\ast),y(0,s_\ast),u(0,s_\ast)\right) & \displaystyle\frac{\partial{y}}{\partial s}(0,s_\ast) \end{pmatrix} \neq 0 maka terdapat lingkungan terbuka I*II_\ast\subset I dari s*s_\ast dan J*J_\ast\subset\mathbb{R} dari 00 sehingga keluarga kurva {γs}sI*\displaystyle\left\{ \gamma_s \right\}_{s\in I_\ast}, untuk tJ*t\in J_\ast, mendefinisikan suatu permukaan integral lokal bagi Persamaan (1.1) yang memuat bagian kurva Γ\Gamma dengan sI*s\in I_\ast. Permukaan integral lokal ini tunggal di antara permukaan C1C^1 yang memuat bagian kurva awal yang sama.

Bukti. Tahap i. Dari teorema eksistensi, ketunggalan, dan ketergantungan terdiferensialkan solusi persamaan diferensial biasa terhadap data awal, terdapat interval terbuka IbII_b\subset I yang memuat s*s_\ast dan JbJ_b\subset\mathbb{R} yang memuat 00 sedemikian sehingga solusi xx, yy, dan uu terdefinisi dan terdiferensialkan secara kontinu pada Jb×IbJ_b\times I_b. Menurut Teorema Pemetaan Invers pada (0,s*)(0,s_\ast), kita dapat memperkecil kedua interval sehingga pemetaan (t,s)(x(t,s),y(t,s))(t,s)\mapsto(x(t,s),y(t,s)) merupakan difeomorfisme dari Jb×IbJ_b\times I_b ke citranya. Khususnya, determinannya tidak nol di seluruh persegi panjang lokal ini. Definisikan Sb={(x(t,s),y(t,s),u(t,s)):(t,s)Jb×Ib}S_b = \left\{ \big( x(t,s),y(t,s),u(t,s)\big) : (t,s) \in J_b \times I_b \right\}

Tahap ii. Keberadaan permukaan terdiferensialkan u=u(x,y)u=u(x,y) diperoleh secara lokal dari Teorema Pemetaan Invers karena (x,y)(t,s)|(t,s)=det(xtxsytys)|(t,s)=det(a(x(t,s),y(t,s),u(t,s))xs(t,s)b(x(t,s),y(t,s),u(t,s))ys(t,s))0\frac{\partial(x,y)}{\partial(t,s)}\bigg|_{(t,s)} = \det \begin{pmatrix} \displaystyle\frac{\partial{x}}{\partial t} & \displaystyle\frac{\partial{x}}{\partial s} \\[0.7em] \displaystyle\frac{\partial{y}}{\partial t} & \displaystyle\frac{\partial{y}}{\partial s} \end{pmatrix}\bigg|_{(t,s)} = \det \begin{pmatrix} a\left(x(t,s),y(t,s),u(t,s)\right) & \displaystyle\frac{\partial{x}}{\partial s}(t,s) \\[0.7em] b\left(x(t,s),y(t,s),u(t,s)\right) & \displaystyle\frac{\partial{y}}{\partial s}(t,s) \end{pmatrix} \neq 0 di Jb×IbJ_b\times I_b. Dengan mensubstitusikan fungsi balik tunggal pemetaan tersebut ke u=u(t,s)u=u(t,s), diperoleh grafik lokal u=u(x,y)u=u(x,y) yang mendeskripsikan SbS_b.

Tahap iii. Permukaan u=u(x,y)u=u(x,y) ini memuat bagian lokal Γ\Gamma berdasarkan konstruksinya. Sesungguhnya,
Γ|Ib={(x(0,s),y(0,s),u(0,s)):sIb}\Gamma\big|_{I_b} = \{ (x(0,s),y(0,s),u(0,s)) : s \in I_b \}. Selain itu, titik-titik pada permukaan SbS_b berbentuk (x,y,u)=(x(t,s),y(t,s),u(t,s))=(x(t,s),y(t,s),u(x(t,s),y(t,s)))(x,y,u) = \left(x(t,s),y(t,s), u(t,s)\right) = \left(x(t,s),y(t,s), u\big(x(t,s),y(t,s)\big)\right) untuk suatu sIbs\in I_b dan tJbt\in J_b. Oleh karena itu, c(x,y,u)=c(x(t,s),y(t,s),u(t,s))=tu(t,s)=tu(x(t,s),y(t,s))=ux(x(t,s),y(t,s))xt(t,s)+uy(x(t,s),y(t,s))yt(t,s)=a(x(t,s),y(t,s),u(t,s))ux(x(t,s),y(t,s))+b(x(t,s),y(t,s),u(t,s))uy(x(t,s),y(t,s))=a(x,y,u)ux(x,y)+b(x,y,u)uy(x,y).\begin{align*} c\left(x,y, u\right) &= c\left(x(t,s),y(t,s),u(t,s)\right) = \frac{\partial}{\partial t}{u(t,s)} = \frac{\partial}{\partial t}{u(x(t,s),y(t,s))} \\ &= \frac{\partial{u}}{\partial x}\left(x(t,s),y(t,s)\right) \frac{\partial{x}}{\partial t}\left(t,s\right) + \frac{\partial{u}}{\partial y}\left(x(t,s),y(t,s)\right) \frac{\partial{y}}{\partial t}\left(t,s\right) \\ &= a(x(t,s),y(t,s),u(t,s)) \,\frac{\partial{u}}{\partial x}\left(x(t,s),y(t,s)\right) \\ & \qquad + b(x(t,s),y(t,s),u(t,s)) \,\frac{\partial{u}}{\partial y}\left(x(t,s),y(t,s)\right) \\ &= a(x,y,u) \,\frac{\partial{u}}{\partial x}\left(x,y\right) + b(x,y,u) \,\frac{\partial{u}}{\partial y}\left(x,y\right) \ . \end{align*} Jadi, Persamaan (1.1) dipenuhi.

Tahap iv. Untuk membuktikan ketunggalan lokal, kita menunjukkan bahwa kurva karakteristik yang berawal pada permukaan integral u=v(x,y)u=v(x,y) bagi Persamaan (1.1) dan memuat bagian kurva awal yang sama tetap berada pada permukaan integral tersebut. Berdasarkan ketunggalan kurva karakteristik, grafik u=v(x,y)u=v(x,y) dan grafik u=u(x,y)u=u(x,y) yang telah kita temukan berimpit pada lingkungan yang cukup kecil.

Misalkan g(t,s)=u(t,s)v(x(t,s),y(t,s)),g(t,s) = u(t,s) - v(x(t,s),y(t,s)) \ , dengan γs={(x,y,u)=(x(t,s),y(t,s),u(t,s)):tJb}\gamma_s = \left\{ (x,y,u) = (x(t,s), y(t,s) ,u(t,s)) : t \in J_b \right\} merupakan bagian kurva karakteristik, dengan sIbs\in I_b dianggap tetap. Fungsi xx, yy, dan uu memenuhi Persamaan (1.2) dan Persamaan (1.3). Oleh karena itu,

gt(t,s)=ut(t,s)vx(x(t,s),y(t,s))xt(t,s)vy(x(t,s),y(t,s))yt(t,s)=c(x(t,s),y(t,s),u(t,s))a(x(t,s),y(t,s),u(t,s))vx(x(t,s),y(t,s))b(x(t,s),y(t,s),u(t,s))vy(x(t,s),y(t,s))=c(x(t,s),y(t,s),g(t,s)+v(x(t,s),y(t,s)))a(x(t,s),y(t,s),g(t,s)+v(x(t,s),y(t,s)))vx(x(t,s),y(t,s))b(x(t,s),y(t,s),g(t,s)+v(x(t,s),y(t,s)))vy(x(t,s),y(t,s)).\begin{align} \frac{\partial{g}}{\partial t}(t,s) &= \frac{\partial{u}}{\partial t}(t,s) - \frac{\partial{v}}{\partial x}\left(x(t,s),y(t,s)\right) \frac{\partial{x}}{\partial t}(t,s) - \frac{\partial{v}}{\partial y}\left(x(t,s),y(t,s)\right) \frac{\partial{y}}{\partial t}(t,s) \nonumber \\ &= c(x(t,s),y(t,s),u(t,s)) - a(x(t,s),y(t,s),u(t,s))\, \frac{\partial{v}}{\partial x}\left(x(t,s),y(t,s)\right) \nonumber \\ &\qquad - b(x(t,s),y(t,s),u(t,s))\,\frac{\partial{v}}{\partial y}\left(x(t,s),y(t,s)\right) \nonumber \\ &= c(x(t,s),y(t,s),g(t,s)+v(x(t,s),y(t,s))) \nonumber \\ &\qquad - a(x(t,s),y(t,s),g(t,s)+v(x(t,s),y(t,s)))\, \frac{\partial{v}}{\partial x}(x(t,s),y(t,s)) \nonumber \\ &\qquad - b(x(t,s),y(t,s),g(t,s)+v(x(t,s),y(t,s)))\,\frac{\partial{v}}{\partial y}(x(t,s),y(t,s)) \ . \end{align} Ini merupakan persamaan diferensial biasa orde pertama bagi gg dengan syarat awal g(0,s)=u(0,s)v(x(0,s),y(0,s))=u0(s)v(x0(s),y0(s))=0g(0,s) = u(0,s) - v(x(0,s),y(0,s)) = u_0(s) - v(x_0(s),y_0(s)) = 0 karena u=v(x,y)u=v(x,y) merupakan permukaan yang memuat bagian kurva awal Γ|Ib\Gamma\big|_{I_b}. Berdasarkan ketunggalan solusi persamaan diferensial biasa dan fakta bahwa vv memenuhi Persamaan (1.1), diperoleh g(t,s)=0g(t,s)=0 untuk tJbt\in J_b dan sIbs\in I_b karena fungsi tersebut merupakan solusi Persamaan (1.5). Dengan demikian, v(x(t,s),y(t,s))=u(t,s)=u(x(t,s),y(t,s))v(x(t,s),y(t,s))=u(t,s)=u(x(t,s),y(t,s)) pada persegi panjang lokal yang sedang ditinjau; yakni, v=uv=u di citranya. \square

Catatan. Penting untuk diingat bahwa, untuk setiap kurva karakteristik
γs={(x(t,s),y(t,s),u(t,s)):tJs}\displaystyle \gamma_s = \left\{ \left(x(t,s),y(t,s),u(t,s)\right) : t \in J_s\right\}, vektor (xt(t,s),yt(t,s),ut(t,s))\displaystyle \left(\frac{\partial{x}}{\partial t}(t,s), \frac{\partial{y}}{\partial t}(t,s), \frac{\partial{u}}{\partial t}(t,s)\right) berada pada bidang singgung permukaan integral u=u(x,y)u=u(x,y) (lebih tepatnya, sejajar dengan bidang singgung) di titik (x(t,s),y(t,s),u(t,s))\left(x(t,s),y(t,s),u(t,s)\right).

Catatan. Jika Persamaan (1.4) tidak dipenuhi pada (t,s)=(0,s̃)(t,s) = (0,\tilde{s}), maka

0=(x,y)(t,s)|(0,s̃)=det(xtxsytys)|(0,s̃)=det(a(x(0,s̃),y(0,s̃),u(0,s̃))xs(0,s̃)b(x(0,s̃),y(0,s̃),u(0,s̃))ys(0,s̃))=det(a(x0(s̃),y0(s̃),u0(s̃))x0(s̃)b(x0(s̃),y0(s̃),u0(s̃))y0(s̃))=a(x0(s̃),y0(s̃),u0(s̃))y0(s̃)b(x0(s̃),y0(s̃),u0(s̃))x0(s̃).\begin{align} 0 &= \frac{\partial(x,y)}{\partial(t,s)}\bigg|_{(0,\tilde{s})} = \det \begin{pmatrix} \displaystyle\frac{\partial{x}}{\partial t} & \displaystyle\frac{\partial{x}}{\partial s} \\[0.7em] \displaystyle\frac{\partial{y}}{\partial t} & \displaystyle\frac{\partial{y}}{\partial s} \end{pmatrix}\bigg|_{(0,\tilde{s})} = \det \begin{pmatrix} a(x(0,\tilde{s}), y(0,\tilde{s}), u(0,\tilde{s})) & \displaystyle\frac{\partial{x}}{\partial s}(0,\tilde{s}) \\[0.7em] b(x(0,\tilde{s}), y(0,\tilde{s}), u(0,\tilde{s})) & \displaystyle\frac{\partial{y}}{\partial s}(0,\tilde{s}) \end{pmatrix} \nonumber \\ &= \det \begin{pmatrix} a(x_0(\tilde{s}), y_0(\tilde{s}), u_0(\tilde{s})) & \displaystyle x_0'(\tilde{s}) \\ b(x_0(\tilde{s}), y_0(\tilde{s}), u_0(\tilde{s})) & \displaystyle y_0'(\tilde{s}) \end{pmatrix} \nonumber \\ &= a(x_0(\tilde{s}), y_0(\tilde{s}), u_0(\tilde{s}))y_0'(\tilde{s}) - b(x_0(\tilde{s}), y_0(\tilde{s}),u_0(\tilde{s})) x_0'(\tilde{s}) \ .\end{align}

Jika Γ\Gamma berada pada permukaan integral u=u(x,y)u=u(x,y), maka u0(s)=u(x0(s),y0(s))u_0(s) = u(x_0(s), y_0(s)) sehingga

u0(s)=ux(x0(s),y0(s))x0(s)+uy(x0(s),y0(s))y0(s).u_0'(s) = \frac{\partial{u}}{\partial x}(x_0(s), y_0(s))x_0'(s) +\frac{\partial{u}}{\partial y}(x_0(s), y_0(s))y_0'(s) \ . Selain itu, dari Persamaan (1.1) kita memperoleh

a(x,y,u)ux(x,y)+b(x,y,u)uy(x,y)=c(x,y,u)a(x, y, u)\frac{\partial{u}}{\partial x}(x, y) + b(x, y, u)\frac{\partial{u}}{\partial y}(x, y) = c(x, y, u) pada (x,y,u)=(x0(s),y0(s),u0(s))(x,y,u)= (x_0(s), y_0(s), u_0(s)) untuk sIs \in I.

Oleh karena itu, dari Persamaan (1.7) dan Persamaan (1.8), kita memperoleh

b(x0(s̃),y0(s̃),u0(s̃))u0(s̃)c(x0(s̃),y0(s̃),u0(s̃))y0(s̃)=b(x0(s̃),y0(s̃),u0(s̃))(ux(x0(s̃),y0(s̃))x0(s̃)+uy(x0(s̃),y0(s̃))y0(s̃))(a(x,y,u)ux(x,y)+b(x,y,u)uy(x,y))|(x,y,u)=(x0(s̃),y0(s̃),u0(s̃))y0(s̃)=(b(x0(s̃),y0(s̃),u0(s̃))x0(s̃)a(x0(s̃),y0(s̃),u0(s̃))y0(s̃))ux(x0(s̃),y0(s̃))=0,\begin{align} & b(x_0(\tilde{s}), y_0(\tilde{s}),u_0(\tilde{s})) u_0'(\tilde{s}) - c(x_0(\tilde{s}), y_0(\tilde{s}),u_0(\tilde{s})) y_0'(\tilde{s}) \nonumber \\ &\qquad = b(x_0(\tilde{s}), y_0(\tilde{s}),u_0(\tilde{s})) \left(\frac{\partial{u}}{\partial x}(x_0(\tilde{s}), y_0(\tilde{s}))x_0'(\tilde{s}) +\frac{\partial{u}}{\partial y}(x_0(\tilde{s}), y_0(\tilde{s}))y_0'(\tilde{s})\right) \nonumber \\ &\qquad \quad- \left(a(x, y, u)\frac{\partial{u}}{\partial x}(x,y) + b(x, y, u) \frac{\partial{u}}{\partial y}(x,y) \right)\bigg|_{(x,y,u)= (x_0(\tilde{s}), y_0(\tilde{s}), u_0(\tilde{s}))} y_0'(\tilde{s}) \nonumber \\ &\qquad = \bigg( b(x_0(\tilde{s}), y_0(\tilde{s}),u_0(\tilde{s}))x_0'(\tilde{s}) - a(x_0(\tilde{s}), y_0(\tilde{s}), u_0(\tilde{s}))y_0'(\tilde{s}) \bigg) \frac{\partial{u}}{\partial x}(x_0(\tilde{s}), y_0(\tilde{s})) \nonumber \\ &\qquad = 0 \ ,\end{align} dengan identitas terakhir berasal dari Persamaan (1.6). Serupa dengan itu,

a(x0(s̃),y0(s̃),u0(s̃))u0(s̃)c(x0(s̃),y0(s̃),u0(s̃))x0(s̃)=a(x0(s̃),y0(s̃),u0(s̃))(ux(x0(s̃),y0(s̃))x0(s̃)+uy(x0(s̃),y0(s̃))y0(s̃))(a(x,y,u)ux(x,y)+b(x,y,u)uy(x,y))|(x,y,u)=(x0(s̃),y0(s̃),u0(s̃))x0(s̃)=(a(x0(s̃),y0(s̃),u0(s̃))y0(s̃)b(x0(s̃),y0(s̃),u0(s̃))x0(s̃))uy(x0(s̃),y0(s̃))=0.\begin{align} & a(x_0(\tilde{s}), y_0(\tilde{s}),u_0(\tilde{s})) u_0'(\tilde{s}) - c(x_0(\tilde{s}), y_0(\tilde{s}),u_0(\tilde{s})) x_0'(\tilde{s}) \nonumber \\ &\qquad = a(x_0(\tilde{s}), y_0(\tilde{s}),u_0(\tilde{s})) \left(\frac{\partial{u}}{\partial x}(x_0(\tilde{s}), y_0(\tilde{s}))x_0'(\tilde{s}) +\frac{\partial{u}}{\partial y}(x_0(\tilde{s}), y_0(\tilde{s}))y_0'(\tilde{s})\right) \nonumber \\ &\qquad \quad- \left(a(x, y, u)\frac{\partial{u}}{\partial x}(x,y) + b(x, y, u) \frac{\partial{u}}{\partial y}(x,y) \right)\bigg|_{(x,y,u)= (x_0(\tilde{s}), y_0(\tilde{s}),u_0(\tilde{s}))} x_0'(\tilde{s}) \nonumber \\ &\qquad = \bigg( a(x_0(\tilde{s}), y_0(\tilde{s}),u_0(\tilde{s}))y_0'(\tilde{s}) - b(x_0(\tilde{s}), y_0(\tilde{s}), u_0(\tilde{s}))x_0'(\tilde{s}) \bigg) \frac{\partial{u}}{\partial y}(x_0(\tilde{s}), y_0(\tilde{s})) \nonumber \\ &\qquad = 0 \ .\end{align} Dengan demikian, dari Persamaan (1.6), Persamaan (1.9), dan Persamaan (1.10), kita menyimpulkan bahwa (a(x,y,u),b(x,y,u),c(x,y,u))|(x,y,u)=(x0(s̃),y0(s̃),u0(s̃))×(x0(s̃),y0(s̃),u0(s̃))=0\left(a(x,y,u), b(x,y,u), c(x,y,u)\right) \big|_{(x,y,u)=(x_0(\tilde{s}), y_0(\tilde{s}), u_0(\tilde{s}))} \times \left(x_0'(\tilde{s}), y_0'(\tilde{s}), u_0'(\tilde{s}) \right) = 0 sehingga kedua vektor tersebut sejajar.

Jika Persamaan (1.4) tidak dipenuhi pada (t,s)=(0,s̃)(t,s)=(0,\tilde{s}), kesimpulannya bersifat lokal pada titik tersebut. Jika Persamaan (1.4) tidak dipenuhi dan data awal memenuhi syarat kompatibilitas di atas dan arah kurva Γ\Gamma berimpit dengan arah karakteristik, ketunggalan lokal dapat hilang: data yang kompatibel dapat termuat dalam lebih dari satu permukaan integral lokal. Jika syarat kompatibilitas tidak dipenuhi, tidak ada permukaan integral C1C^1 yang memuat Γ\Gamma di sekitar titik itu. Kegagalan pada satu nilai s̃\tilde{s} tidak dengan sendirinya menjadikan seluruh kurva Γ\Gamma karakteristik.

Contoh 1. Selesaikan persamaan diferensial parsial ux+uy=2\frac{\partial{u}}{\partial x} + \frac{\partial{u}}{\partial y} = 2 dengan syarat u(x,0)=x2\displaystyle u(x,0)=x^2.

Persamaan karakteristiknya adalah x(t)=1x'(t) = 1, y(t)=1y'(t) = 1, dan u(t)=2u'(t) = 2. Kurva awal diberikan oleh x(0)=sx(0) = s, y(0)=0y(0) = 0, dan u(0)=s2\displaystyle u(0) = s^2.

Persamaan x(t)=1x'(t) = 1 dengan x(0)=sx(0)=s memberikan x=x(t,s)=t+sx = x(t,s) = t + s, persamaan y(t)=1y'(t) = 1 dengan y(0)=0y(0)=0 memberikan y=y(t,s)=ty = y(t,s) = t, dan persamaan u(t)=2u'(t) = 2 dengan u(0)=s2\displaystyle u(0)=s^2 memberikan u=u(t,s)=2t+s2\displaystyle u = u(t,s) = 2t+s^2. Kita dapat menyelesaikan x=t+sx = t+s dan y=ty=t terhadap ss dan tt sehingga t=yt=y dan s=xys= x-y. Jika nilai-nilai ss dan tt ini kita substitusikan ke ekspresi untuk uu, kita memperoleh permukaan integral u=u(x,y)=2y+(xy)2\displaystyle u= u(x,y) = 2y + (x-y)^2.

Contoh 2. Selesaikan persamaan diferensial parsial xux+uy=ux \frac{\partial{u}}{\partial x} + \frac{\partial{u}}{\partial y} = u dengan syarat u(x,0)=ln(x)u(x,0)=\ln(x) untuk x>0x>0.

Persamaan karakteristiknya adalah x(t)=x(t)x'(t) = x(t), y(t)=1y'(t) = 1, dan u(t)=u(t)u'(t) = u(t). Kurva awal diberikan oleh x(0)=sx(0) = s, y(0)=0y(0) = 0, dan u(0)=ln(s)u(0) = \ln(s).

Persamaan x(t)=x(t)x'(t) = x(t) dengan x(0)=sx(0)=s memberikan x=x(t,s)=set\displaystyle x = x(t,s) = s e^t, persamaan y(t)=1y'(t) = 1 dengan y(0)=0y(0)=0 memberikan y=y(t,s)=ty = y(t,s) = t, dan persamaan u(t)=u(t)u'(t) = u(t) dengan u(0)=ln(s)u(0)=\ln(s) memberikan u=u(t,s)=ln(s)et\displaystyle u = u(t,s) = \ln(s)e^t. Kita dapat menyelesaikan x=set\displaystyle x = se^t dan y=ty=t terhadap ss dan tt sehingga t=yt=y dan s=xey\displaystyle s= xe^{-y}. Jika nilai-nilai ss dan tt ini kita substitusikan ke ekspresi untuk uu, kita memperoleh permukaan integral u=u(x,y)=ln(xey)ey=ln(x)eyyey.u= u(x,y) = \ln(xe^{-y}) e^y = \ln(x) e^y -y e^y \ .

Contoh 3. Selesaikan persamaan diferensial parsial ux+uy=u\frac{\partial{u}}{\partial x} + \frac{\partial{u}}{\partial y} = u dengan syarat u(x,0)=cos(x)u(x,0)=\cos(x).

Persamaan karakteristiknya adalah x(t)=1x'(t) = 1, y(t)=1y'(t) = 1, dan u(t)=u(t)u'(t) = u(t). Kurva awal diberikan oleh x(0)=sx(0) = s, y(0)=0y(0) = 0, dan u(0)=cos(s)u(0) = \cos(s).

Persamaan x(t)=1x'(t) = 1 dengan x(0)=sx(0)=s memberikan x=x(t,s)=t+sx = x(t,s) = t + s, persamaan y(t)=1y'(t) = 1 dengan y(0)=0y(0)=0 memberikan y=y(t,s)=ty = y(t,s) = t, dan persamaan u(t)=u(t)u'(t) = u(t) dengan u(0)=cos(s)u(0)=\cos(s) memberikan u=u(t,s)=cos(s)et\displaystyle u = u(t,s) = \cos(s)e^t. Kita dapat menyelesaikan x=t+sx = t+s dan y=ty=t terhadap ss dan tt sehingga t=yt=y dan s=xys= x-y. Jika nilai-nilai ss dan tt ini kita substitusikan ke ekspresi untuk uu, kita memperoleh permukaan integral u=u(x,y)=cos(xy)ey\displaystyle u= u(x,y) = \cos(x-y) e^y.

Contoh 4. Selesaikan persamaan diferensial parsial x2ux+y2uy=u2x^2 \frac{\partial{u}}{\partial x} + y^2 \frac{\partial{u}}{\partial y} = u^2 dengan syarat u(x,2x)=1u(x,2x)=1 untuk x0x\neq0.

Persamaan karakteristiknya adalah x(t)=x2(t)\displaystyle x'(t) = x^2(t), y(t)=y2(t)\displaystyle y'(t) = y^2(t), dan u(t)=u2(t)\displaystyle u'(t) = u^2(t). Kurva awal diberikan oleh x(0)=sx(0) = s, y(0)=2sy(0) = 2s, dan u(0)=1u(0) = 1 untuk s0s\neq0.

Persamaan x(t)=x2(t)\displaystyle x'(t) = x^2(t) dengan x(0)=sx(0)=s memberikan

x=x(t,s)=s1st.x = x(t,s) = \frac{s}{1-st} \ . Serupa dengan itu, persamaan y(t)=y2(t)\displaystyle y'(t) = y^2(t) dengan y(0)=2sy(0)=2s memberikan

y=y(t,s)=2s12st.y = y(t,s) = \frac{2s}{1-2st} \ . Terakhir, persamaan u(t)=u2(t)\displaystyle u'(t) = u^2(t) dengan u(0)=1u(0)=1 memberikan u=u(t,s)=11t\displaystyle u = u(t,s) = \frac{1}{1-t}. Kita dapat menyelesaikan Persamaan (1.11) dan Persamaan (1.12) terhadap ss dan tt dalam variabel xx dan yy sehingga diperoleh t=1x2y\displaystyle t = \frac{1}{x} - \frac{2}{y} dan s=xy2y2x\displaystyle s = \frac{xy}{2y-2x}. Jika nilai-nilai ss dan tt ini kita substitusikan ke ekspresi untuk uu, kita memperoleh permukaan integral u=u(x,y)=xyxyy+2x.u= u(x,y) = \frac{xy}{xy-y+2x} \ . Rumus ini berlaku secara lokal di sekitar masing-masing dari kedua bagian kurva awal y=2xy=2x, x0x\neq0, pada daerah tempat xyy+2x0xy-y+2x\neq0 dan solusi karakteristik tetap terdefinisi.

Metode Lagrange

Metode Lagrange menghasilkan solusi umum (yakni, suatu keluarga permukaan integral yang tidak diharuskan memuat kurva awal Γ\Gamma) dari suatu persamaan diferensial parsial kuasilinear orde pertama seperti Persamaan (1.1). Kita mengasumsikan bahwa kurva-kurva karakteristik diberikan oleh perpotongan dua keluarga permukaan yang berbentuk ϕ(x,y,u)=αdanψ(x,y,u)=β.\phi(x,y,u) = \alpha \quad \text{dan} \quad \psi(x,y,u) = \beta \ . Kurva karakteristik yang berkaitan dengan syarat awal (x0(s),y0(s),u0(s))Γ(x_0(s),y_0(s),u_0(s)) \in \Gamma merupakan perpotongan permukaan ϕ(x,y,u)=αs\phi(x,y,u) = \alpha_s dan ψ(x,y,u)=βs\psi(x,y,u) = \beta_s untuk suatu pasangan (αs,βs)(\alpha_s, \beta_s), sebagaimana diperlihatkan pada Gambar 2.

Perpotongan dua permukaan tingkat menghasilkan kurva karakteristik yang melalui data awal.

Himpunan {(αs,βs):sI}\displaystyle\left\{ (\alpha_s,\beta_s) : s \in I \right\} merupakan suatu kurva pada bidang α\alpha,β\beta. Jika kurva ini diberikan oleh F(α,β)=0F(\alpha,\beta)=0, permukaan integral persamaan kuasilinear yang memuat Γ\Gamma diberikan oleh F(ϕ(x,y,u),ψ(x,y,u))=0.F(\phi(x,y,u), \psi(x,y,u)) = 0 \ . Masalahnya adalah menemukan ϕ\phi dan ψ\psi. Dalam beberapa kasus yang jarang terjadi, kita dapat menemukan dua fungsi bernilai vektor v1:33\displaystyle v_1:\mathbb{R}^3\rightarrow \mathbb{R}^3 dan v2:33\displaystyle v_2:\mathbb{R}^3\rightarrow \mathbb{R}^3 sedemikian sehingga, pada setiap titik di daerah yang sedang ditinjau, v1(x,y,u)v_1(x,y,u) dan v2(x,y,u)v_2(x,y,u) bebas secara linear dan ortogonal terhadap (a(x,y,u),b(x,y,u),c(x,y,u))\big(a(x,y,u),b(x,y,u),c(x,y,u)\big) untuk semua (x,y,u)(x,y,u), serta ×(v1)=×(v2)=0\nabla\mathbin{\times}(v_1) = \nabla\mathbin{\times}(v_2)=0. Karena v1v_1 dan v2v_2 irrotasional, terdapat fungsi-fungsi terdiferensialkan ϕ:3\displaystyle\phi:\mathbb{R}^3\rightarrow \mathbb{R} dan ψ:3\displaystyle\psi:\mathbb{R}^3\rightarrow \mathbb{R} sedemikian sehingga ϕ=v1\nabla\phi = v_1 dan ψ=v2\nabla\psi = v_2. Fungsi ϕ\phi dan ψ\psi konstan sepanjang kurva-kurva karakteristik karena tϕ(x(t,s),y(t,s),u(t,s))=ϕ(x(t,s),y(t,s),u(t,s))(xt(t,s),yt(t,s),ut(t,s))=v1(x(t,s),y(t,s),u(t,s))(a(x,y,u),b(x,y,u),c(x,y,u))|(x,y,u)=(x(t,s),y(t,s),u(t,s))=0\begin{multline*} \frac{\partial}{\partial t}{ \phi(x(t,s),y(t,s),u(t,s)) } = \nabla\phi(x(t,s),y(t,s),u(t,s)) \cdot \left(\frac{\partial{x}}{\partial t}(t,s), \frac{\partial{y}}{\partial t}(t,s), \frac{\partial{u}}{\partial t}(t,s) \right) \\ = v_1(x(t,s),y(t,s),u(t,s)) \cdot (a(x,y,u),b(x,y,u),c(x,y,u))\bigg|_{(x,y,u)=(x(t,s),y(t,s),u(t,s))} = 0 \end{multline*} dan, serupa dengan itu, tψ(x(t,s),y(t,s),u(t,s))=0\frac{\partial}{\partial t}{ \psi(x(t,s),y(t,s),u(t,s)) } = 0 untuk semua tt dan ss.

Contoh 5. Temukan permukaan integral persamaan diferensial parsial yux+xuy=0-y \frac{\partial{u}}{\partial x} + x \frac{\partial{u}}{\partial y} = 0 dengan syarat u(x,x)=ex\displaystyle u(x,x)=e^x untuk x0x\neq0.

Kita mempunyai (a(x,y,u),b(x,y,u),c(x,y,u))=(y,x,0).\left(a(x,y,u),b(x,y,u),c(x,y,u)\right) = (-y,x,0) \ . Pada daerah x2+y2>0x^2+y^2>0, vektor v1(x,y,u)=(x,y,0)v_1(x,y,u)=(x,y,0) dan v2(x,y,u)=(0,0,1)v_2(x,y,u)=(0,0,1) bebas secara linear dan memenuhi ×(v1)=×(v2)=0\nabla\mathbin{\times}(v_1) = \nabla\mathbin{\times}(v_2) = 0 serta ortogonal terhadap
(a(x,y,u),b(x,y,u),c(x,y,u))(a(x,y,u),b(x,y,u),c(x,y,u)) untuk semua (x,y,u)(x,y,u). Kita mempunyai bahwa ϕ(x,y,u)=(x2+y2)/2\displaystyle\phi(x,y,u)=(x^2+y^2)/2 memenuhi ϕ=v1\nabla\phi = v_1 dan ψ(x,y,u)=u\psi(x,y,u)=u memenuhi ψ=v2\nabla\psi = v_2. Permukaan integral persamaan diferensial parsial berbentuk

F(x2+y2,u)=0F(x^2+y^2, u) = 0 untuk fungsi-fungsi F:2\displaystyle F:\mathbb{R}^2 \rightarrow \mathbb{R}.

Misalkan Fβ(α,β)0\displaystyle\frac{\partial{F}}{\partial \beta}(\alpha,\beta) \neq 0, maka secara lokal kita dapat menyelesaikan F(α,β)=0F(\alpha,\beta)=0 terhadap β\beta sebagai fungsi dari α\alpha. Dari Persamaan (1.13) diperoleh u(x,y)=f(x2+y2)\displaystyle u(x,y)=f(x^2+y^2) untuk suatu fungsi f:(0,)f:(0,\infty)\rightarrow\mathbb{R} pada daerah yang relevan.

Pada komponen awal x>0x>0, syarat u(x,x)=exu(x,x)=e^x dan z=2x2z=2x^2 memberikan f+(z)=ez/2f_+(z)=e^{\sqrt{z/2}}. Pada komponen awal x<0x<0, syarat yang sama memberikan f(z)=ez/2f_-(z)=e^{-\sqrt{z/2}}. Jadi, kedua permukaan integral lokalnya adalah u+(x,y)=e(x2+y2)/2danu(x,y)=e(x2+y2)/2,u_+(x,y)=e^{\sqrt{(x^2+y^2)/2}} \quad\text{dan}\quad u_-(x,y)=e^{-\sqrt{(x^2+y^2)/2}}, masing-masing pada lingkungan komponen awal positif dan negatif di dalam x2+y2>0x^2+y^2>0. Tidak ada satu cabang pada daerah terhubung yang memuat kedua titik awal (r,r)(r,r) dan (r,r)(-r,-r) beserta busur karakteristik di antaranya: invariansi sepanjang lingkaran akan memaksa er=ere^r=e^{-r} untuk r>0r>0.

Atribusi dan hak. Karya sumber oleh Benoit Dionne dilisensikan berdasarkan Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International. Terjemahan ini menandai perubahan dan menggunakan lisensi komponen yang sama. Ini bukan terbitan atau dukungan resmi Benoit Dionne maupun University of Ottawa.


  1. Domain parameter lengkapnya adalah D={(t,s):sI,tJs}D=\{(t,s):s\in I,\ t\in J_s\}. Mulai sekarang, kita akan sering menyatakan “untuk semua tt dan ss” dengan maksud “untuk semua (t,s)D(t,s)\in D”.↩︎

← Kembali ke Program Matematika