PHƯƠNG TRÌNH VI PHÂN TOÀN PHẦN



1. Tản mạn

Trời thì đang mưa rất là to, sáng ngủ dậy làm được gói mỳ tôm, nhai thêm hai cái bánh chocolate và làm thêm cốc nước nữa kakaka... đủ năng lượng cho một buổi sáng. Tiếng mưa vẫn rơi rả rích trên mái nhà, ngay lúc này đây, được ngồi một mình ngâm cứu một thứ gì đó thì quả thật không còn gì bằng. Làm thêm một ngụm nước, kê bàn và lật sách phương trình vi phân ra đọc. Chà ngay lúc này, cũng có một cảm nhận gì đó tí tẹo về phương trình vi phân toàn phần, tạm thời viết một bài nhỏ ở đây, có khi giờ nhớ đấy, nhưng vài hôm nữa lại quên ngay...

Đầu tiên CaolacVC sẽ trình bày tóm gọn nhất về dạng và cách giải và sau đó là đưa ra ví dụ, để khi cần thì chỉ cần mở ra nhớ và làm, cái này giống như là delta của phương trình bậc hai, nhớ và làm thôi, còn việc ngâm cứu chuyên sâu như thế nào đó là nhu cầu của mỗi người. OK! Let's Go!


2. Dạng phương trình vi phân toàn phần

Phương trình vi phân dạng

$$M(x,y)dx+N(x,y)dy=0\quad (1)$$

được gọi là phương trình vi phần toàn phần khi nó thỏa mãn điều kiện: vế trái của phương trình $(1)$ phải là vi phân toàn phần của một hàm khả vi nào đó. Tức là tồn tại một hàm $U(x,y)$ khả vi nào đó sao cho

$$dU(x,y)=M(x,y)dx+N(x,y)dy$$

Điều kiện để một phương trình vi phân dạng $(1)$ trở thành phương trình vi phân toàn phần (hay cách nhận biết phương trình vi phân toàn phần) là:

$$\displaystyle \frac{\partial M}{\partial y}=\frac{\partial N}{\partial x}$$

Ví dụ. Phương trình vi phân $(3x^2+6xy^2)dx+(6x^2y+4y^3)dy$ là phương trình vi phân toàn phần vì

$M(x,y)=(3x^2+6xy^2), N(x,y)=(6x^2y+4y^3)$

$\displaystyle \frac{\partial M}{\partial y}=\frac{\partial N}{\partial x}=12xy$


3. Cách giải

Ở đây ta chỉ nêu cách giải và áp dụng, nên mới nói là nó cũng như công thức delta mà thôi. Việc tìm hiểu thêm là nhu cầu của mỗi người.

Nếu phương trình $(1)$ là phương trình vi phân toàn phần thì tích phân tổng quát của phương trình $(1)$ là:

$$\displaystyle U(x,y)=\int_{x_0}^{x}M(x,y_0)dx+\int_{y_0}^{y}N(x,y)dy=C\quad (2.1)$$

hoặc

$$\displaystyle U(x,y)=\int_{x_0}^{x}M(x,y)dx+\int_{y_0}^{y}N(x_0,y)dy=C \quad (2.2)$$

với $(x_0,y_0)$ là một điểm điểm bất kỳ mà khi thay vào các hàm $M(x,y_0), N(x_0,y)$ xác định. Thường thì ta sẽ chọn sao cho thuận tiện trong việc tính tích phân nhất.

Cách nhớ của CaolacVC: Công thức này cực kỳ dễ nhớ nếu ta để ý một chút. Trước tiên một cách hình thức ta lấy tích phân hai vế của của phương trình $(1)$, lưu ý chỗ có đuôi vi phân là $dx$ thì cận chạy từ $x_0$ đến $x$, chỗ có phần đuôi là $dy$ thì cận chạy từ $y_0$ đến $y$.

Điều quan trọng nhất là hãy nhớ chỉ được thay $x_0$ hoặc $y_0$ vào một trong hai vị trí: phần có đuôi $dx$ thì chỉ được thay $y_0$. Phần có đuôi $dy$ thì chỉ được thay $x_0$.


4. Ví dụ áp dụng

Ví dụ luôn là thứ giải đáp mọi thắc mắc mà phần lý thuyết nếu đọc ta vẫn chưa nắm vững.

Ví dụ 1. Giải phương trình:

$$(3x^2+6xy^2)dx+(6x^2y+4y^3)dy \quad (3.1)$$

Trước tiên ta phải kiểm tra điều kiện để phương trình đã cho là phương trình vi phân toàn phần hay không. Nếu có thì ta mới áp dụng được công thức của bài viết này, nếu không thì bài viết này coi như vứt, chẳng giúp ích được gì. (Mách nhỏ là ở trên, lúc nãy hình như kiểm tra là phương trình vi phân toàn phần rồi hay sao đó)

$M(x,y)=(3x^2+6xy^2), N(x,y)=(6x^2y+4y^3)$

$\displaystyle \frac{\partial M}{\partial y}=\frac{\partial N}{\partial x}=12xy$

Vậy $(3.1)$ là phương trình vi phân toàn phần. Ta chọn $(x_0;y_0)=(0,0)$. Khi đó theo công thức $(2.2)$ ta được

$\displaystyle \int_{0}^{x}(3x^2+6xy^2)dx+\int_{0}^{y}4y^3dy=C$

Hay tích phân tổng quát của $(3.1)$ là

$x^3+3x^2y^2+y^4=C$


Một số bài tập có lời giải

Bài tập 1. Giải phương trình $$(y-x)dx+(y^3+x)dy=0$$

$M(x;y)=y-x$

$N(x;y)=y^3+x$

$\dfrac{\partial M}{dy}=\dfrac{\partial N}{dx}=1$ (thỏa mãn điều kiện của phương trình vi phân toàn phần)

Chọn $(x_0;y_0)=(0;0)$. Tích phân tổng quát

$\displaystyle \int_0^x (0-x)dx+\int_0^y (y^3+x)dy=C$

$\displaystyle \Rightarrow \int_0^x(-x)dx+\int_0^y (y^3+x)dy=C$

$\displaystyle \Rightarrow -\frac{x^2}{2}+\frac{y^4}{4}+xy=C$

Bài tập 2. Giải phương trình $$\left[ \left( 1+x+y \right){{e}^{x}}+{{e}^{y}} \right]dx+\left( {{e}^{x}}+x{{e}^{y}} \right)dy=0$$

$\displaystyle M(x;y)=(1+x+y){{e}^{x}}+{{e}^{y}}$

$\displaystyle N(x;y)={{e}^{x}}+x{{e}^{y}}$

$\displaystyle \frac{\partial M}{dy}=\frac{\partial N}{dx}={{e}^{x}}$

Chọn $({{x}_{0}};{{y}_{0}})=(0;0)$

Suy ra tích phân tổng quát

$\displaystyle \int\limits_{0}^{x}{\left[ (1+x+0){{e}^{x}}+{{e}^{0}} \right]dx}+\int\limits_{0}^{y}{\left( {{e}^{x}}+x{{e}^{y}} \right)dy}=C$

$\displaystyle \Leftrightarrow \int\limits_{0}^{x}{\left[ (1+x){{e}^{x}}+1 \right]dx}+\int\limits_{0}^{y}{\left( {{e}^{x}}+x{{e}^{y}} \right)dy}=C$

$\displaystyle \Leftrightarrow \left[ (1+x){{e}^{x}}-{{e}^{x}}+x \right]\left| _{0}^{x} \right.+\left( y{{e}^{x}}+x{{e}^{y}} \right)\left| _{0}^{y} \right.=C$

$\displaystyle \Leftrightarrow {{e}^{x}}(x+y)+x{{e}^{y}}=C$

Post a Comment

13 Comments

  1. This comment has been removed by the author.

    ReplyDelete
  2. Ad ơi,M(x0,y) và N(x,y0) phải đồng thời xác định ạ?

    ReplyDelete
    Replies
    1. Đúng rồi, chọn sao cho tính tích phân nó dễ là được!

      Delete
  3. (X0, y0) có cần phải thoả mãn khi thay vào pt ban đầu thì đúng = 0 k ạ?

    ReplyDelete
  4. giả sử có bài cho y>0 thì đặt (x0,y0) thế nào ạ

    ReplyDelete
  5. cảm ơn ad, bài viết rất kỹ và chi tiết

    ReplyDelete

Vui lòng đăng nhập google để bình luận
Để gõ công thức toán, hãy đặt [biểu thức toán] trong dấu $$
Ví dụ: $[biểu thức toán]$