Tốc độ thay đổi của hàm nhiều biến theo từng hướng tọa độ, kết hợp lại thành vectơ gradient.
Trực giácHai độ dốc cùng lúc trên một bề mặt cong
Hãy tưởng tượng bạn đang đứng trên một sườn đồi có độ cao cho bởi hàm hai biến z=f(x,y), trong đó x chỉ hướng đông và y chỉ hướng bắc. Khác với đường cong một biến, không có một độ dốc duy nhất dưới chân bạn: bước thẳng về hướng đông làm độ cao thay đổi theo một tốc độ, bước thẳng về hướng bắc lại thay đổi theo tốc độ khác, còn bước xiên góc sẽ pha trộn cả hai. Cố định y và chỉ cho x thay đổi tương đương với việc cắt sườn đồi bằng một mặt phẳng thẳng đứng đông-tây, đưa bài toán về đạo hàm một biến quen thuộc.
Ngọn đồi hình chuông trơn trong không gian ba chiều với đỉnh ở giữa, bị cắt bởi các mặt phẳng đứng theo hai hướng vuông góc.
Đồi có độ cao z=e−(x2+y2). Giữ nguyên y và chỉ cho x biến đổi sẽ cắt mặt bởi một mặt phẳng đứng đông-tây; hệ số góc của đường cắt đó ngay tại chỗ bạn đứng là ∂x∂f. Giữ nguyên x thay vào đó cho ∂y∂f.
Đại họcĐịnh nghĩa chặt chẽ: đạo hàm riêng, mặt phẳng tiếp xúc và gradient
Định nghĩa: Đạo hàm riêng và vectơ gradient
Cho f:U⊆Rn→R xác định trên tập mở U. Đạo hàm riêng của f theo biến xi tại a∈U, ký hiệu ∂xi∂f(a) hoặc fxi(a), là giới hạn khi giữ cố định mọi tọa độ khác xj (j=i). Ghép tất cả n đạo hàm riêng thành một vectơ duy nhất ta được vectơ gradient ∇f(a)=(∂x1∂f(a),…,∂xn∂f(a)).
Sự tồn tại đơn thuần của các đạo hàm riêng tại (x0,y0) yếu hơn tính khả vi toàn phần. Hàm f(x,y) được gọi là khả vi toàn phần tại (x0,y0) nếu số gia Δf=f(x0+h,y0+k)−f(x0,y0) biểu diễn được dưới dạng fx(x0,y0)h+fy(x0,y0)k+o(h2+k2) khi (h,k)→(0,0). Khi f khả vi toàn phần, đồ thị z=f(x,y) có mặt phẳng tiếp xúc xác định tại (x0,y0), và đạo hàm theo hướng của mọi vectơ đơn vị u=(u1,u2) với ∥u∥=1 được cho bởi tích vô hướng Duf=∇f⋅u.
Lưới mũi tên của trường gradient của sin x cos y với các đường dòng hội tụ về các đỉnh của hàm thế.
Trường gradient ∇φ=(cosxcosy,−sinxsiny) của thế φ(x,y)=sinxcosy: tại mỗi điểm mũi tên chỉ hướng tăng nhanh nhất (nơi Duf=∇f⋅u lớn nhất), và các đường dòng bám theo nó tới các cực đại của φ.
Tổng hợp các toán tử vi phân nhiều biến và đại lượng cấp hai
Đối tượng
Công thức
Ý nghĩa hình học / vật lý
Đạo hàm riêng
∂xi∂f
Hệ số góc của lát cắt dọc theo trục tọa độ xi.
Vectơ gradient
∇f=(fx1,…,fxn)
Vuông góc với mặt mức, chỉ hướng tăng nhanh nhất; chuẩn ∥∇f∥ là tốc độ tăng cực đại.
Đạo hàm theo hướng
Duf=∇f⋅u
Tốc độ thay đổi tức thời khi di chuyển theo vectơ vận tốc đơn vị u.
Ma trận Hessian và biệt thức
D=fxxfyy−fxy2
Đo độ cong địa phương để phân loại điểm dừng thành cực trị hoặc điểm yên ngựa.
Đại họcCác định lý cốt lõi: Hướng tăng nhanh nhất, đối xứng Schwarz và tiêu chuẩn Hessian
Nếu f khả vi tại a và ∇f(a)=0, thì với mọi vectơ đơn vị u (thỏa ∥u∥=1), đạo hàm theo hướng thỏa mãn Duf(a)=∇f(a)⋅u=∥∇f(a)∥cosθ, trong đó θ là góc giữa ∇f(a) và u. Hệ quả là Duf(a) đạt giá trị lớn nhất ∥∇f(a)∥ khi u=∥∇f(a)∥∇f(a), và đạt giá trị nhỏ nhất −∥∇f(a)∥ theo hướng ngược lại.
Vì sao đúng?
Tích vô hướng chiếu vectơ gradient lên hướng di chuyển đã chọn: bạn nhận toàn bộ độ lớn của gradient chỉ khi đi trùng khít với hướng của nó, và độ cao không đổi khi bạn đi vuông góc với nó dọc theo đường đồng mức.
Chứng minh
Xét hàm một biến g(t)=f(a+tu). Theo định nghĩa đạo hàm theo hướng, Duf(a)=g′(0). Vì f khả vi toàn phần tại a, ta có f(a+tu)−f(a)=∇f(a)⋅(tu)+o(∣t∣) khi t→0. Chia hai vế cho t và lấy giới hạn khi t→0 ta thu được Duf(a)=∇f(a)⋅u.
Theo công thức hình học của tích vô hướng Euclid và điều kiện vectơ đơn vị ∥u∥=1, ta có ∇f(a)⋅u=∥∇f(a)∥∥u∥cosθ=∥∇f(a)∥cosθ. Do −1≤cosθ≤1, giá trị lớn nhất ∥∇f(a)∥ đạt được duy nhất tại θ=0, tức khi u=∥∇f(a)∥∇f(a), và giá trị nhỏ nhất −∥∇f(a)∥ đạt được tại θ=π.
Nếu f(x,y) có các đạo hàm riêng cấp hai liên tục trong một lân cận của (x0,y0), thì ∂x∂y∂2f(x0,y0)=∂y∂x∂2f(x0,y0). Hơn nữa, nếu ∇f(x0,y0)=0 và D=fxx(x0,y0)fyy(x0,y0)−fxy(x0,y0)2, thì (x0,y0) là điểm cực tiểu địa phương chặt nếu D>0 và fxx>0, là điểm cực đại địa phương chặt nếu D>0 và fxx<0, và là điểm yên ngựa nếu D<0.
Vì sao đúng?
Tính đối xứng của đạo hàm riêng hỗn hợp đảm bảo ma trận Hessian đối xứng, và việc đưa dạng toàn phương của nó về tổng bình phương cho biết mặt cong lõm lên như lòng chảo, úp xuống như mái vòm hay xoắn lại như yên ngựa quanh mặt phẳng tiếp xúc nằm ngang.
Chứng minh
Với h,k nhỏ khác không, xét sai phân cấp hai Δ(h,k)=f(x0+h,y0+k)−f(x0+h,y0)−f(x0,y0+k)+f(x0,y0). Áp dụng Định lý giá trị trung bình một biến trước hết cho g(x)=f(x,y0+k)−f(x,y0) trên [x0,x0+h] rồi cho fx theo biến y ta được Δ(h,k)=hkfyx(c1,d1), trong khi áp dụng theo thứ tự ngược lại cho Δ(h,k)=hkfxy(c2,d2) tại các điểm trung gian hội tụ về (x0,y0). Cho hai biểu thức bằng nhau và lấy giới hạn (h,k)→(0,0) theo tính liên tục ta chứng minh được ∂x∂y∂2f(x0,y0)=∂y∂x∂2f(x0,y0).
Tại điểm dừng nơi ∇f(x0,y0)=0, công thức Taylor cho f(x0+h,y0+k)−f(x0,y0)=21Q(h,k)+o(h2+k2), với Q(h,k)=fxxh2+2fxyhk+fyyk2. Khi fxx=0, tách bình phương đủ viết lại Q(h,k)=fxx1[(fxxh+fxyk)2+Dk2] với D=fxxfyy−fxy2. Nếu D>0, biểu thức trong ngoặc vuông dương chặt với mọi (h,k)=(0,0), nên dấu của Q(h,k) trùng với dấu của fxx (cho cực tiểu chặt khi fxx>0 và cực đại chặt khi fxx<0). Nếu D<0, Q(h,k) nhận cả giá trị dương lẫn âm dọc theo các đường thẳng khác nhau qua gốc, tạo thành điểm yên ngựa.
Đại họcỨng dụng thực tiễn và Ví dụ minh họa
Đạo hàm riêng và gradient là động lực của khoa học và kỹ thuật hiện đại: trong vật lý, định luật dẫn nhiệt Fourier phát biểu rằng thông lượng nhiệt tỉ lệ với đối của gradient nhiệt độ −∇T; trong học máy, việc huấn luyện mạng nơ-ron cực tiểu hóa hàm mất mát L(θ) bằng cách lặp lại các bước dịch chuyển ngược hướng gradient −∇L(θ) (hạ gradient); còn trong kinh tế học và cơ học, tiêu chuẩn Hessian giúp phân biệt các trạng thái cân bằng bền (cực tiểu địa phương của thế năng) với các điểm yên ngựa không bền.
Ví dụ: Đạo hàm theo hướng và hướng tăng nhanh nhất của trường nhiệt độ
Một tấm kim loại có nhiệt độ T(x,y)=x2+3xy+y3 (tính bằng độ C) tại tọa độ (x,y) tính bằng mét. Hãy tìm vectơ gradient ∇T(1,2), đạo hàm theo hướng tại (1,2) theo hướng của vectơ v=(3,4), và tốc độ tăng nhiệt độ lớn nhất tại (1,2).
Lời giải
Trước hết, tính các đạo hàm riêng bằng cách lấy đạo hàm theo một biến và giữ nguyên biến kia: Tx(x,y)=2x+3y và Ty(x,y)=3x+3y2. Thay tọa độ (1,2) vào ta được Tx(1,2)=2(1)+3(2)=8 và Ty(1,2)=3(1)+3(2)2=15, vậy vectơ gradient là ∇T(1,2)=(8,15).
Tiếp theo, chuẩn hóa v=(3,4) để được vectơ đơn vị u=∥v∥v=(53,54) vì ∥v∥=32+42=5. Đạo hàm theo hướng là DuT(1,2)=∇T(1,2)⋅u=8⋅53+15⋅54=584=16.8. Tốc độ tăng nhiệt độ lớn nhất tại (1,2) đạt được theo hướng của ∇T(1,2) và bằng ∥∇T(1,2)∥=82+152=289=17 độ C trên mét.
Ví dụ: Phân loại điểm dừng bằng định thức Hessian
Tìm và phân loại tất cả các điểm dừng của hàm số f(x,y)=x3+y3−3xy bằng tiêu chuẩn đạo hàm cấp hai Hessian.
Lời giải
Cho cả hai đạo hàm riêng cấp một bằng không: fx(x,y)=3x2−3y=0 và fy(x,y)=3y2−3x=0. Từ phương trình thứ nhất ta có y=x2; thế vào phương trình thứ hai được 3(x2)2−3x=3x(x3−1)=0, có nghiệm thực x=0 và x=1. Do đó hai điểm dừng là (0,0) và (1,1).
Tiếp theo tính các đạo hàm riêng cấp hai: fxx=6x, fyy=6y và fxy=−3, suy ra biệt thức Hessian D(x,y)=fxxfyy−fxy2=36xy−9. Tại (0,0), ta có D(0,0)=−9<0, nên (0,0) là điểm yên ngựa. Tại (1,1), ta có D(1,1)=36−9=27>0 và fxx(1,1)=6>0, nên (1,1) là điểm cực tiểu địa phương chặt với giá trị cực tiểu f(1,1)=−1.
Mặt hình yên ngựa đi lên theo một trục và đi xuống theo trục còn lại, có điểm dừng tại gốc tọa độ.
Mặt yên ngựa mẫu z=x2−y2: gốc tọa độ là điểm dừng (fx=fy=0) nhưng không phải cực đại hay cực tiểu — mặt đi lên dọc trục x và đi xuống dọc trục y. Tiêu chuẩn Hessian nhận ra hình dạng này qua D=fxxfyy−fxy2=−4<0.
Cho hàm số f(x,y)=x3y2−4x+5y. Vectơ gradient ∇f(1,2) bằng bao nhiêu?
Nếu ∇f(x0,y0)=(6,−8), giá trị lớn nhất của đạo hàm theo hướng Duf(x0,y0) trên mọi vectơ đơn vị u là bao nhiêu?
Trong học máy, thuật toán hạ gradient cập nhật tham số θ để cực tiểu hóa hàm mất mát khả vi L(θ). Công thức cập nhật nào dịch chuyển θ theo hướng làm L giảm nhanh nhất tại lân cận địa phương?
Giả sử (x0,y0) là điểm dừng của hàm trơn f(x,y) có fxx(x0,y0)=−4, fyy(x0,y0)=−3 và fxy(x0,y0)=2. Điểm (x0,y0) thuộc loại nào?
Tài liệu tham khảo
Jerrold E. Marsden, Anthony J. Tromba (2012). Vector Calculus