Hàm Green

Hàm Green là công cụ chuẩn để giải một phương trình vi phân tuyến tính có hạng nguồn khi các điều kiện biên hoặc điều kiện đầu đã được cố định. Nói đơn giản, nó cho biết hệ phản ứng thế nào với một nguồn điểm đơn vị đặt tại $\xi$.

Với một toán tử tuyến tính $L$, hàm Green $G(x,\xi)$ được định nghĩa bởi bài toán nguồn điểm

$$L\,G(x,\xi) = \delta(x-\xi)$$

theo biến $x$, kèm theo đúng các điều kiện biên hoặc điều kiện đầu mà nghiệm cuối cùng phải thỏa mãn. Ở đây $\delta(x-\xi)$ là delta Dirac, một phân phối biểu diễn một nguồn đơn vị tập trung tại $\xi$.

Nếu hàm Green như vậy tồn tại và bài toán là tuyến tính, thì nghiệm của

$$L\,u(x) = f(x)$$

thường có thể được viết dưới dạng

$$u(x) = \int G(x,\xi) f(\xi)\, d\xi.$$

Đó là ý tưởng chính: giải bài toán nguồn điểm một lần, rồi ghép các đáp ứng điểm đó lại để xử lý một nguồn tổng quát $f$.

Hàm Green có ý nghĩa gì

Nguồn delta là điều làm cho định nghĩa này trở nên hữu ích. Thay vì giải ngay toàn bộ phương trình, trước hết bạn hỏi hệ đáp ứng ra sao với đầu vào tập trung nhất có thể.

Vì thế, có thể hiểu $G(x,\xi)$ là “ảnh hưởng tại $x$ của một xung đơn vị tác dụng tại $\xi$”. Một nguồn phức tạp hơn $f$ khi đó được xem như sự chồng chất liên tục của các nguồn điểm. Bước này chỉ đúng vì bài toán là tuyến tính.

Vì sao điều kiện biên làm thay đổi đáp án

Hàm Green không chỉ được quyết định bởi riêng toán tử vi phân. Nó thuộc về toàn bộ bài toán, bao gồm cả các điều kiện đi kèm.

Ví dụ, toán tử $-u''$ trên $0<x<1$ với $u(0)=u(1)=0$ có một hàm Green khác với cùng toán tử đó nhưng dùng điều kiện biên Neumann. Phương trình có thể trông giống nhau, nhưng tập nghiệm được phép thay đổi, nên hạt nhân cũng thay đổi theo.

Đây là một trong những lỗi phổ biến nhất. Không có một hàm Green duy nhất cho “phương trình” nếu các điều kiện chưa được nêu rõ trong phát biểu bài toán.

Ví dụ có lời giải: bài toán Dirichlet cho $-u''(x)=f(x)$

Xét bài toán giá trị biên

$$-u''(x) = f(x), \qquad 0 < x < 1,$$

với

$$u(0)=0, \qquad u(1)=0.$$

Với bài toán cụ thể này, hàm Green là

$$G(x,\xi) = \begin{cases} x(1-\xi), & x \le \xi, \\ \xi(1-x), & x \ge \xi. \end{cases}$$

Nó liên tục theo $x$, triệt tiêu tại $x=0$ và $x=1$, và đạo hàm bậc nhất của nó có độ nhảy đúng tại $x=\xi$. Chính độ nhảy đó tạo ra nguồn delta trong phương trình $-G''(x,\xi)=\delta(x-\xi)$.

Khi đã biết hạt nhân này, mọi hạng nguồn của cùng bài toán giá trị biên đều dùng chung một công thức:

$$u(x)=\int_0^1 G(x,\xi) f(\xi)\, d\xi.$$

Chọn một hạng nguồn đơn giản:

$$f(x)=1.$$

Khi đó

$$u(x)=\int_0^1 G(x,\xi)\, d\xi = \int_0^x \xi(1-x)\, d\xi + \int_x^1 x(1-\xi)\, d\xi.$$

Tính từng phần:

$$\int_0^x \xi(1-x)\, d\xi = (1-x)\frac{x^2}{2},$$ $$\int_x^1 x(1-\xi)\, d\xi = x\frac{(1-x)^2}{2}.$$

Cộng lại, ta được

$$u(x)=\frac{x(1-x)}{2}.$$

Kiểm tra nhanh:

$$u'(x)=\frac{1-2x}{2}, \qquad u''(x)=-1,$$

nên

$$-u''(x)=1=f(x),$$

và các giá trị biên là

$$u(0)=0, \qquad u(1)=0.$$

Ví dụ này cho thấy rất rõ lợi ích của phương pháp. Bạn không phải dựng lại cách giải cho mỗi $f$ mới; bạn giữ nguyên cùng một hạt nhân và chỉ thay đổi tích phân.

Những lỗi thường gặp với hàm Green

1. Xem hàm Green như thể nó hoạt động cho bài toán phi tuyến theo đúng cách đó. Bước chồng chất phụ thuộc vào tính tuyến tính.
2. Quên rằng điều kiện biên hoặc điều kiện đầu là một phần của định nghĩa.
3. Cho rằng hạt nhân lúc nào cũng phải đối xứng. Tính đối xứng thường cần thêm cấu trúc, chẳng hạn một thiết lập tự liên hợp phù hợp.
4. Nhầm lẫn giữa hàm Green và nghiệm cơ bản. Chúng liên hệ chặt chẽ với nhau, nhưng không phải lúc nào cũng là cùng một đối tượng.
5. Dùng công thức tích phân mà không kiểm tra xem hàm Green có tồn tại cho bài toán cụ thể đó hay không.

Hàm Green được dùng ở đâu

Hàm Green xuất hiện trong phương trình vi phân thường và phương trình đạo hàm riêng khi một bài toán tuyến tính có hạng nguồn. Chúng rất phổ biến trong tĩnh điện học, khuếch tán, bài toán sóng, cơ học lượng tử và đàn hồi học.

Trong các thiết lập tuyến tính bất biến theo thời gian, hàm Green cũng liên hệ chặt chẽ với đáp ứng xung. Cách gọi có thể khác nhau giữa các lĩnh vực, nhưng ý tưởng thì giống nhau: hiểu đáp ứng với một đầu vào tập trung, rồi từ đó dựng đáp ứng với đầu vào tổng quát.

Cách nhớ nhanh

Nếu bạn quên định nghĩa hình thức, hãy giữ phiên bản này:

$$\text{Hàm Green} = \text{đáp ứng với một nguồn điểm}.$$

Khi đó, nghiệm đầy đủ được tạo ra bằng cách cộng các đáp ứng nguồn điểm trên toàn miền, và đó chính là điều phép tích phân thực hiện.

Thử một hạng nguồn tương tự

Giữ nguyên các điều kiện biên, nhưng thay $f(x)=1$ bằng $f(x)=x$. Hạt nhân vẫn giữ nguyên, và chỉ tích phân cuối cùng thay đổi:

$$u(x)=\int_0^1 G(x,\xi)\,\xi\, d\xi.$$

Đây là bước tiếp theo rất tốt vì nó tách bạch hai phần thay đổi một cách rõ ràng: hàm Green thuộc về thiết lập bài toán, còn hạng nguồn có thể thay đổi theo từng trường hợp. Nếu muốn đi xa hơn, hãy thử phiên bản của riêng bạn với một nguồn khác và kiểm tra tích phân thu được với các điều kiện biên.