Tìm hàm số cực tiểu hoặc cực đại các đại lượng tích phân, như đường ngắn nhất hay thời gian ít nhất.
Trực giácĐường nào có "chi phí" nhỏ nhất?
Giải tích thông thường tìm con số đầu vào làm cực tiểu một con số đầu ra. Giải tích biến phân đặt một câu hỏi lạ hơn: trong tất cả các hàm số có thể thỏa mãn điều kiện biên nào đó, hàm số nào làm cực tiểu một tích phân xây dựng từ nó? Một tia sáng đi qua thủy tinh bẻ cong theo hướng làm cực tiểu thời gian đi (nguyên lý Fermat); một sợi dây xích treo tự do có hình dạng làm cực tiểu thế năng; một màng xà phòng căng trên khung dây làm cực tiểu diện tích bề mặt. Trong mọi trường hợp, ẩn số không phải một con số mà là cả một đường cong.
Widget đường cong bậc ba tương tác minh họa một đường đi cạnh tranh trong bài toán biến phân.
Đường cong bậc ba y=ax3+bx2+cx+d này biểu diễn một đường đi ứng viên đang cạnh tranh để làm cực trị một phiếm hàm; kéo các thanh trượt để xem hình dạng đường (và do đó tích phân xây dựng từ nó) thay đổi ra sao khi rời khỏi cực trị thật.
Phổ thôngĐường ngắn nhất và diện tích bao lớn nhất
Định nghĩa: Phiếm hàm
Phiếm hàm là một quy tắc gán một con số duy nhất cho cả một hàm số, thay vì cho một con số. Giải tích biến phân nghiên cứu các phiếm hàm được xây dựng dưới dạng tích phân của một hàm số cùng đạo hàm của nó, và tìm hàm số làm tích phân này nhỏ nhất (hoặc lớn nhất) có thể.
J[y]=∫x0x1L(x,y,y′)dx
Ở đây x0 và x1 là hai đầu mút cố định, y(x) là đường cong cần tìm với giá trị cho trước tại hai đầu mút, y′(x) là đạo hàm của nó, và L (hàm Lagrange) là một hàm ba biến cho trước mã hóa đại lượng đang được tích lũy — độ dài cung, thời gian đi, năng lượng, v.v.
Cho y(x) là hàm khả vi liên tục đến cấp hai trên [x0,x1] với giá trị đầu mút cố định y(x0)=y0, y(x1)=y1. Nếu y làm cực trị phiếm hàm J[y]=∫x0x1L(x,y,y′)dx trong số mọi đường cong như vậy, thì y thỏa mãn ∂y∂L−dxd∂y′∂L=0 với mọi x∈(x0,x1).
Vì sao đúng?
Điều này giống với việc cho đạo hàm thông thường bằng không tại một cực tiểu, nhưng ở đây 'hướng' mà ta nhiễu loạn không phải một con số mà là cả một biến phân khả dĩ của đường cong. Yêu cầu độ biến thiên bậc nhất triệt tiêu với mọi biến phân như vậy buộc phải có sự cân bằng tại từng điểm giữa sự phụ thuộc trực tiếp vào y của L và sự phụ thuộc vào y′, được mã hóa chính xác bởi phương trình này.
Chứng minh
Cố định một hàm trơn tùy ý η(x) với η(x0)=η(x1)=0, và xét họ một tham số các đường cạnh tranh yϵ(x)=y(x)+ϵη(x), tất cả đều có cùng giá trị đầu mút với y. Đặt ϕ(ϵ)=J[yϵ]; vì y được giả sử làm cực trị J, hàm số thông thường ϕ có điểm tới hạn tại ϵ=0, nên ϕ′(0)=0.
Lấy đạo hàm dưới dấu tích phân, ϕ′(ϵ)=∫x0x1[∂y∂Lη+∂y′∂Lη′]dx, trong đó các đạo hàm riêng của L được tính dọc theo yϵ. Cho ϵ=0 ta được ∫x0x1[∂y∂Lη+∂y′∂Lη′]dx=0.
Lấy tích phân từng phần số hạng thứ hai: ∫x0x1∂y′∂Lη′dx=[∂y′∂Lη]x0x1−∫x0x1dxd∂y′∂Lηdx. Số hạng biên triệt tiêu vì η(x0)=η(x1)=0, còn lại ∫x0x1[∂y∂L−dxd∂y′∂L]ηdx=0.
Tích phân này triệt tiêu với mọiη khả dĩ. Theo bổ đề cơ bản của giải tích biến phân — nếu một hàm liên tục lấy tích phân bằng không với mọi hàm thử như vậy thì bản thân hàm đó phải đồng nhất bằng không — đại lượng trong ngoặc triệt tiêu tại mọi x, đó chính xác là ∂y∂L−dxd∂y′∂L=0.
Nếu hàm Lagrange L=L(y,y′) không phụ thuộc tường minh vào x, thì dọc theo mọi cực trị y(x), đại lượng L−y′∂y′∂L=C là hằng số.
Vì sao đúng?
Đây là một đại lượng bảo toàn, hoàn toàn tương tự bảo toàn năng lượng: khi 'luật chơi' L không đổi khi ta trượt dọc theo x, một tổ hợp cụ thể của L và y′ được giữ cố định. Nó cho một phương trình cấp một thay cho phương trình Euler–Lagrange cấp hai, vốn thường dễ giải trực tiếp hơn nhiều.
Chứng minh
Đặt H(x)=L(y,y′)−y′∂y′∂L tính dọc theo một cực trị y(x), và lấy đạo hàm theo x bằng quy tắc dây chuyền: dxdH=∂y∂Ly′+∂y′∂Ly′′−y′′∂y′∂L−y′dxd∂y′∂L.
Hai số hạng chứa y′′ triệt tiêu chính xác, còn lại dxdH=y′(∂y∂L−dxd∂y′∂L).
Nhưng y(x) là một cực trị, nên theo phương trình Euler–Lagrange thừa số trong ngoặc đồng nhất bằng không dọc theo y. Do đó dxdH=0 trên toàn khoảng, nghĩa là H(x)=L−y′∂y′∂L là hằng số, chứng minh khẳng định.
Đại họcỨng dụng thực tiễn và Ví dụ minh họa
Giải tích biến phân là nền tảng của đường trắc địa trong thuyết tương đối rộng và lập kế hoạch chuyển động robot (đường ngắn nhất/chi phí ít nhất trên không gian cong), nguyên lý Fermat trong quang học và thiết kế thấu kính, màng xà phòng diện tích cực tiểu và thiết kế vỏ mái kiến trúc, và các bài toán điều khiển tối ưu trong kinh tế học và kỹ thuật hàng không (quỹ đạo tiêu tốn nhiên liệu ít nhất). Hai ví dụ dưới đây tự tay suy ra hai đường cong biến phân nổi tiếng nhất.
Ví dụ: Bài toán đường tốc hành (brachistochrone)
Một hạt trượt không ma sát dưới tác dụng trọng lực dọc theo một sợi dây từ gốc (0,0) xuống một điểm thấp hơn (x1,y1) (lấy y dương hướng xuống). Trong số mọi hình dạng dây nối hai điểm này, hình dạng nào giúp hạt đến nơi trong thời gian ít nhất?
Lời giải
Theo bảo toàn năng lượng, xuất phát từ trạng thái nghỉ tại y=0, tốc độ của hạt ở độ cao y thỏa mãn v=2gy. Vì tốc độ là độ dài cung trên một đơn vị thời gian, tổng thời gian đi là T[y]=∫0x12gy1+(y′)2dx, dùng ds=1+(y′)2dx.
Biểu thức dưới dấu tích phân L=2gy1+(y′)2 phụ thuộc vào y và y′ nhưng không phụ thuộc tường minh vào x, nên đồng nhất thức Beltrami L−y′∂y′∂L=C áp dụng trực tiếp. Tính ∂y′∂L=2gy1+(y′)2y′ và rút gọn L−y′∂y′∂L cho 2gy1+(y′)21=const, viết lại thành y[1+(y′)2]=k với hằng số k nào đó.
Phương trình cấp một này được giải bằng phép thế y′=cot(θ/2), sau khi lấy tích phân cho họ tham số x=a(θ−sinθ),y=a(1−cosθ) với a=k/2 — một đường cycloid, đường mà một điểm trên vành bánh xe lăn vạch ra.
Vậy đường nhanh nhất không phải đường thẳng: đường cycloid ban đầu đổ dốc mạnh hơn đường thẳng, đánh đổi quãng đường dài hơn để lấy tốc độ lớn hơn sớm hơn, và sự đánh đổi này chính là điều mà bộ máy Euler–Lagrange/Beltrami làm rõ một cách chính xác.
Ví dụ: Bài toán đẳng chu: diện tích lớn nhất với chu vi cố định
Trong số mọi đường cong đóng đơn trơn trong mặt phẳng có chu vi cố định cho trước L0, đường nào bao một diện tích lớn nhất?
Lời giải
Đây là một bài toán biến phân có ràng buộc: cực đại diện tích A=21∮(xdy−ydx) với ràng buộc chu vi ∮(x′)2+(y′)2dt=L0 cố định, trong đó đường cong được tham số hóa bởi t.
Giống như tối ưu hóa có ràng buộc trong giải tích thông thường, đưa vào một nhân tử Lagrange λ và làm cực trị phiếm hàm kết hợp A−λ⋅(perimeter) một cách tự do (không ràng buộc), bằng cách áp dụng phương trình Euler–Lagrange riêng cho các thành phần x(t) và y(t).
Thực hiện biến phân này cho thấy độ cong của một đường cực trị phải là hằng số dọc theo toàn bộ chiều dài — một đường cong có độ cong hằng số trong mặt phẳng chính là một đường tròn. Vậy trong số mọi đường cong đóng có chu vi L0, chỉ có đường tròn bán kính r=L0/(2π) mới có thể là cực trị.
Vì một đường tròn bán kính đó thật sự bao diện tích A=πr2=L02/(4π), và không đường cong đóng nào cùng chu vi có thể bao diện tích lớn hơn (đây là bất đẳng thức đẳng chu), đường tròn được xác nhận là cực đại — khớp với quan sát đời thường rằng hình tròn 'dùng' đường biên của nó hiệu quả nhất.
Trong phương trình Euler–Lagrange, y(x) phải thỏa mãn điều kiện gì để phương trình đó đúng?
Đường tốc hành — đường đi xuống nhanh nhất dưới tác dụng trọng lực — có hình dạng gì?
Với phiếm hàm J[y]=∫01(y′)2dx và y(0)=0, y(1)=1, phương trình Euler–Lagrange chọn ra đường cong nào?
Trong số mọi đường cong đóng đơn có chu vi cố định, hình dạng nào làm cực đại diện tích bao?