Tài Liệu Học Tập: Giải Phương Trình Vi Phân Không Thuần Nhất Bằng Hàm Green

1. Giới Thiệu

Trong vật lý, chúng ta thường xuyên gặp các phương trình vi phân không thuần nhất, mô tả các hệ chịu tác động của ngoại lực hoặc nguồn nào đó. Hàm Green là một công cụ mạnh mẽ để giải quyết các phương trình này, đặc biệt khi nguồn tác động là một nguồn điểm (Dirac delta function). Tài liệu này sẽ trình bày chi tiết phương pháp sử dụng hàm Green để giải phương trình vi phân tuyến tính không thuần nhất.

2. Phương Trình Vi Phân Tuyến Tính Không Thuần Nhất

Xét phương trình vi phân tuyến tính tổng quát bậc hai:

$$L[y(x)] = f(x)$$

trong đó:

• $L$ là một toán tử vi phân tuyến tính, ví dụ: $$L = a(x) \frac{d^2}{dx^2} + b(x) \frac{d}{dx} + c(x)$$ với $a(x)$, $b(x)$, và $c(x)$ là các hàm số.
• $y(x)$ là nghiệm cần tìm.
• $f(x)$ là hàm nguồn, mô tả tác động bên ngoài.

3. Hàm Green

Hàm Green $G(x, x')$ là nghiệm của phương trình vi phân khi hàm nguồn là hàm Dirac delta:

$$L[G(x, x')] = \delta(x - x')$$

trong đó $\delta(x - x')$ là hàm Dirac delta, có tính chất:

• $\delta(x - x') = 0$ nếu $x \neq x'$
• $\int_{-\infty}^{\infty} \delta(x - x') dx = 1$

Ý nghĩa vật lý của hàm Green: $G(x, x')$ là đáp ứng của hệ tại vị trí $x$ do một nguồn điểm đặt tại $x'$.

4. Giải Phương Trình Vi Phân Bằng Hàm Green

Nếu ta tìm được hàm Green $G(x, x')$, nghiệm của phương trình vi phân không thuần nhất $L[y(x)] = f(x)$ có thể được biểu diễn dưới dạng tích phân:

$$y(x) = \int G(x, x') f(x') dx'$$

Chứng minh:

Áp dụng toán tử $L$ lên nghiệm $y(x)$:

$$L[y(x)] = L \left[ \int G(x, x') f(x') dx' \right]$$

Vì $L$ là toán tử tuyến tính, ta có thể đưa nó vào trong tích phân:

$$L[y(x)] = \int L[G(x, x')] f(x') dx'$$

Theo định nghĩa của hàm Green, $L[G(x, x')] = \delta(x - x')$, do đó:

$$L[y(x)] = \int \delta(x - x') f(x') dx'$$

Sử dụng tính chất của hàm Dirac delta:

$$L[y(x)] = f(x)$$

Vậy, nghiệm $y(x)$ thỏa mãn phương trình vi phân ban đầu.

5. Các Bước Tìm Hàm Green

1. Giải phương trình thuần nhất: Tìm hai nghiệm độc lập tuyến tính $y_1(x)$ và $y_2(x)$ của phương trình vi phân thuần nhất $L[y(x)] = 0$.

2. Xây dựng hàm Green: Hàm Green thường có dạng:

   $$G(x, x') = \begin{cases} A y_1(x) + B y_2(x), & x < x' \\ C y_1(x) + D y_2(x), & x > x' \end{cases}$$

   trong đó $A$, $B$, $C$, và $D$ là các hằng số cần xác định.

3. Điều kiện biên: Áp đặt các điều kiện biên của bài toán lên hàm Green. Điều này giúp giảm số lượng hằng số chưa biết.

4. Điều kiện liên tục và gián đoạn: Hàm Green phải thỏa mãn hai điều kiện sau tại $x = x'$:

   • Liên tục: $G(x', x'^-) = G(x', x'^+)$
   • Gián đoạn đạo hàm: $$\left. \frac{dG}{dx} \right|_{x = x'^+ } - \left. \frac{dG}{dx} \right|_{x = x'^-} = \frac{1}{a(x')}$$ trong đó $a(x)$ là hệ số của đạo hàm bậc cao nhất trong toán tử $L$.

5. Giải hệ phương trình: Từ các điều kiện biên và điều kiện liên tục/gián đoạn, ta thu được một hệ phương trình tuyến tính để xác định các hằng số $A$, $B$, $C$, và $D$.

6. Ví Dụ Minh Họa

Giải phương trình:

$$\frac{d^2y}{dx^2} = f(x), \quad 0 \le x \le 1$$

với điều kiện biên $y(0) = y(1) = 0$.

1. Tìm hàm Green: Phương trình tương ứng cho hàm Green là:

   $$\frac{d^2G(x, x')}{dx^2} = \delta(x - x')$$

   Phương trình thuần nhất $\frac{d^2G}{dx^2} = 0$ có nghiệm tổng quát $G(x) = Ax + B$. Hai nghiệm độc lập tuyến tính là $y_1(x) = x$ và $y_2(x) = 1$. Tuy nhiên, để thỏa mãn điều kiện biên y(0)=0 và y(1)=0, ta sử dụng nghiệm y1(x) = x và y2(x) = 1-x.

   Do đó, ta có:

   $$G(x, x') = \begin{cases} A x + B(1-x), & 0 \le x < x' \\ C x + D(1-x), & x' < x \le 1 \end{cases}$$

2. Áp dụng điều kiện biên:

   • $G(0, x') = 0 \Rightarrow B = 0$
   • $G(1, x') = 0 \Rightarrow C = 0$

   Vậy hàm Green có dạng:

   $$G(x, x') = \begin{cases} A x, & 0 \le x < x' \\ D(1-x), & x' < x \le 1 \end{cases}$$

3. Điều kiện liên tục: $G(x', x'^-) = G(x', x'^+)$:

   $$A x' = D (1- x')$$

4. Điều kiện gián đoạn đạo hàm:

   $$\left. \frac{dG}{dx} \right|_{x = x'^+ } - \left. \frac{dG}{dx} \right|_{x = x'^-} = 1$$ $$-D - A = 1$$

5. Giải hệ phương trình: Từ $A x' = D(1-x')$ và $-D - A = 1$, ta có:

   • $A = -(1-x')$
   • $D = -x'$

   Vậy hàm Green là:

   $$G(x, x') = \begin{cases} -x(1-x'), & 0 \le x < x' \\ -x'(1-x), & x' < x \le 1 \end{cases}$$

6. Nghiệm của phương trình:

   $$y(x) = \int_0^1 G(x, x') f(x') dx' = -\int_0^x x'(1-x) f(x') dx' - \int_x^1 x(1-x') f(x') dx'$$

7. Kết Luận

Hàm Green là một công cụ hữu ích để giải các phương trình vi phân tuyến tính không thuần nhất. Bằng cách tìm hàm Green tương ứng với một toán tử vi phân và điều kiện biên cụ thể, ta có thể biểu diễn nghiệm của phương trình dưới dạng tích phân. Phương pháp này đặc biệt hiệu quả khi nguồn tác động là nguồn điểm.