›
›

Giáo trình Giải tích 3 - Huỳnh Thế Phùng

pdf 40 trang ngocly 2760
Download
Bạn đang xem 20 trang mẫu của tài liệu "Giáo trình Giải tích 3 - Huỳnh Thế Phùng", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pdfgiao_trinh_giai_tich_3_huynh_the_phung.pdf

Nội dung text: Giáo trình Giải tích 3 - Huỳnh Thế Phùng

  1. GIÁO TRÌNH GIẢI TÍCH III Huỳnh Thế Phùng, Khoa Toán, ĐHKH Huế Ngày 26 tháng 9 năm 2006
  2. 1 Mục lục Chương 1. Phép Tính Vi phân Hàm nhiều biến 3 1.1. Giới hạn và Liên tục . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 1.1.1. Hàm nhiều biến . . . . . . . . . . . . . . . . . . . . . . . . . . 3 1.1.2. Giới hạn . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 1.1.3. Sự liên tục . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 1.2. Đạo hàm và Vi phân . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 1.2.1. Đạo hàm riêng . . . . . . . . . . . . . . . . . . . . . . . . . . 5 1.2.2. Đạo hàm theo hướng . . . . . . . . . . . . . . . . . . . . . . . 6 1.2.3. Vi phân . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 1.2.4. Đạo hàm hàm số hợp và tính bất biến của vi phân . . . . . . 8 1.2.5. Đạo hàm hàm ẩn . . . . . . . . . . . . . . . . . . . . . . . . . 9 1.3. Đạo hàm cấp cao và Công thức Taylor . . . . . . . . . . . . . . . . . 11 1.3.1. Đạo hàm cấp cao . . . . . . . . . . . . . . . . . . . . . . . . . 11 1.3.2. Vi phân cấp cao . . . . . . . . . . . . . . . . . . . . . . . . . . 12 1.3.3. Công thức Taylor . . . . . . . . . . . . . . . . . . . . . . . . . 13 1.4. Cực trị . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 1.4.1. Điều kiện cần . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 1.4.2. Điều kiện đủ . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 1.4.3. Cực trị có điều kiện . . . . . . . . . . . . . . . . . . . . . . . . 16 1.5. Thực hành tính toán trên Maple . . . . . . . . . . . . . . . . . . . . . 17 1.5.1. Giới hạn và đồ thị hàm nhiều biến . . . . . . . . . . . . . . . 17 1.5.2. Tính đạo hàm . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 1.5.3. Khai triển Taylor . . . . . . . . . . . . . . . . . . . . . . . . . 21 1.6. Bài tập . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 Chương 2 Ứng dụng của phép tính vi phân trong hình học 24 2.1. Các hệ toạ độ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24 2.1.1. Hệ toạ độ cực . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
  3. 2 2.1.2. Hệ toạ độ trụ . . . . . . . . . . . . . . . . . . . . . . . . . . . 25 2.1.3. Hệ toạ độ cầu . . . . . . . . . . . . . . . . . . . . . . . . . . . 25 2.2. Hàm vectơ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 2.2.1. Khái niệm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 2.2.2. Giới hạn - Liên tục - Đạo hàm . . . . . . . . . . . . . . . . . . 27 2.3. Các đối tượng liên quan đến đường cong . . . . . . . . . . . . . . . . 28 2.3.1. Tiếp tuyến và pháp tuyến của đường cong phẳng . . . . . . . 28 2.3.2. Tiếp tuyến và pháp diện của đường cong trong không gian . . 29 2.3.3. Độ cong . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29 2.3.4. Hình bao của họ đường cong . . . . . . . . . . . . . . . . . . . 32 2.4. Mặt cong . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32 2.4.1. Khái niệm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32 2.4.2. Tiếp diện và pháp tuyến của mặt cong . . . . . . . . . . . . . 33 2.5. Thực hành tính toán . . . . . . . . . . . . . . . . . . . . . . . . . . . 35 2.5.1. Vẽ đường cong trong mặt phẳng . . . . . . . . . . . . . . . . . 35 2.5.2. Vẽ mặt cong trong không gian . . . . . . . . . . . . . . . . . . 36 2.5.3. Vận động đồ thị . . . . . . . . . . . . . . . . . . . . . . . . . . 37 2.6. Bài tập . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
  4. Chương 1. PHÉP TÍNH VI PHÂN HÀM NHIỀU BIẾN 1.1. Giới hạn và Liên tục 1.1.1. Hàm nhiều biến Cho E là một tập con khác rỗng của Rn. Một ánh xạ f từ E vào R được gọi là một hàm nhiều biến (cụ thể là n biến) xác định trên E: f :E −→ R; x = (x1, ··· , xn) ∈E −→ f(x) = f(x1, ··· , xn) ∈ R. Khi n = 1, f trở thành hàm một biến thực, khi n = 2, 3 ta có hàm hai, ba biến mà thường được viết đơn giản là f(x, y), f(x, y, z) với x, y, z ∈ R. Tập E được gọi là miền xác định của f. Thông thường hàm f được cho dưới dạng công thức còn miền xác định được hiểu là tập hợp các điểm x làm cho f(x) có nghĩa. Chẳng hạn hàm hai biến f(x, y) = ln((x2 + y2)x) có miền xác định là tập E = {(x, y) ∈ R2 | x > 0}. Tương tự đồ thị hàm một biến, đồ thị của hàm n biến f là tập hợp con của Rn+1 mà được định nghĩa như sau: Gr(f) = {(x, f(x)) | x ∈ E}. Bây giờ cho f và g là các hàm nhiều biến trên E và λ là một số thực, ta ký hiệu λf, f ± g, fg, f/g, f ∨ g, f ∧ g là các hàm mới được xác định bởi, ∀x ∈ E : (λf)(x) := λf(x); (f ± g)(x) := f(x) ± g(x); (fg)(x) := f(x)g(x); µ ¶ f f(x) (x) := , (g(x) 6= 0; g g(x)
  5. 4 (f ∨ g)(x) := max{f(x), g(x)}; (f ∧ g)(x) := min{f(x), g(x)}. Ta nói f g và f ≥ g được định nghĩa hoàn toàn tương tự. 1.1.2. Giới hạn Cho f là hàm xác định trên E và x0 ∈ E. Một số thực L được gọi là giới hạn của hàm f tại x0 nếu ∀² > 0, ∃δ > 0, ∀x ∈ E : 0 0, ∀x ∈ E : 0 M; x→x0 lim f(x) = −∞ ⇐⇒ ∀m ∈ R, ∃δ > 0, ∀x ∈ E : 0 < d(x, x0) < δ ⇒ f(x) < m. x→x0 Ví dụ 1.2. 1 lim = +∞. (x,y)→(0,0) x2 + y2 Định lý sau đây được chứng minh tương tự đối với hàm một biến: Định lý 1.2. Giả sử lim f(x) = L ∈ R, lim g(x) = M ∈ R và λ ∈ R. Lúc đó, x→x0 x→x0 a) limx→x0 (f ± g)(x) = L ± M; b) limx→x0 (λf)(x) = λL; c) limx→x (fg)(x) = LM; 0 µ ¶ f L d) Nếu M 6= 0 thì lim (x) = ; x→x0 g M e) Nếu f ≤ g thì L ≤ M. Các phát biểu a)-c) được hiểu là vế trái tồn tại và bằng vế phải mỗi khi vế phải có nghĩa.
  6. 5 1.1.3. Sự liên tục Cho hàm f xác định trên tập E ⊂ Rn và x0 ∈ E. Ta nói f liên tục tại x0 nếu giới hạn của f tại x0 tồn tại và bằng f(x0): lim f(x) = f(x0). x→x0 Nếu f liên tục tại mọi điểm x ∈ E, ta nói f liên tục trên E. Định lý 1.3. Hàm f liên tục tại điểm x0 ∈ E nếu, và chỉ nếu, với mọi dãy vectơ (xk) ⊂ E hội tụ về x0, dãy số (f(xk)) hội tụ về f(x0). 0 Định lý 1.4. Cho hàm n biến f liên tục tại điểm x và n hàm m biến ϕj(u) liên 0 m 0 0 tục tại điểm u ∈ R . Ngoài ra, ϕj(u ) = xj với mọi 1 ≤ j ≤ n. Lúc đó hàm hợp F (u) := f(ϕ1(u), ϕ2(u), ··· , ϕn(u)) là hàm m biến liên tục tại u0. Hệ quả 1.1. Cho f và g là hai hàm xác định trên E, liên tục tại x0 ∈ E và λ là một số thực. Lúc đó, các hàm λf, f ± g, fg đều liên tục tại x0. Hơn nữa, nếu 0 f g(x ) 6= 0 thì hàm g cũng liên tục tại điểm đó. Định lý 1.5. Cho E là tập đóng và bị chặn trong Rn và f là hàm liên tục trên E. Lúc đó ∗ ∗ a) Tồn tại hai điểm x∗, x ∈ E sao cho f(x∗) ≤ f(x) ≤ f(x ) với mọi x ∈ E. b) f liên tục đều trên E, tức là ∀² > 0, ∃δ > 0, ∀x, x0 ∈ E : d(x, x0) < δ ⇒ |f(x) − f(x0)| < ². 1.2. Đạo hàm và Vi phân 1.2.1. Đạo hàm riêng Để đơn giản, trước tiên ta sẽ xét trường hợp hàm hai biến. Cho f : E ⊂ R2 → R và (x0, y0) ∈ Int(E). Lúc đó, tồn tại số dương ² sao cho với mọi số gia ∆x ∈ (−², ²) ta có (x0 + ∆x, y0) ∈ E. Ta sẽ gọi biểu thức sau ∆xf := f(x0 + ∆x, y0) − f(x0, y0) ∆xf là số gia của hàm f tương ứng với số gia ∆x. Nếu tồn tại giới hạn của ∆x khi ∆x → 0 thì giới hạn này được gọi là đạo hàm riêng của hàm f theo biến x tại điểm 0 ∂f (x0, y0) và được ký hiệu là fx(x0, y0) hay ∂x (x0, y0). Vậy 0 ∂f f(x0 + ∆x, y0) − f(x0, y0) fx(x0, y0) = (x0, y0) := lim . ∂x ∆x→0 ∆x
  7. 6 Tương tự, ta có định nghĩa đạo hàm riêng của f theo biến y tại (x0, y0): 0 ∂f f(x0, y0 + ∆y) − f(x0, y0) fy(x0, y0) = (x0, y0) := lim ∂y ∆y→0 ∆y 0 0 0 và ngay cả với hàm n biến f(x1, x2, ··· , xn) tại một điểm x = (x1, ··· , xn). Chẳng hạn, ∂f f(x0 + ∆x , x0, ··· , x0 ) − f(x0, x0, ··· , x0 ) (x0) := lim 1 1 2 n 1 2 n . ∂x1 ∆x1→0 ∆x1 Nếu tại điểm x0 ∈ E đạo hàm riêng của f theo n biến đều tồn tại thì ta gọi vectơ µ ¶ ∂f ∂f ∂f ∇f(x0) := (x0), (x0), ··· , (x0) ∂x1 ∂x2 ∂xn là građiên của f tại x0. Có khi người ta cũng ký hiệu vectơ này là gradf(x0). Trong thực tế, để tính đạo hàm riêng của một hàm f theo biến xi ta chỉ việc xem f như là hàm một biến xi còn các biến khác là hằng số. y 2 Ví dụ 1.3. Với f(x, y) = x và g(x, y, z) = x y sin(x + z) ta có µ ¶ y 1 ∇f(x, y) = − , ; x2 x ¡ ¢ ∇g(x, y, z) = 2xy sin(x + z) + x2y cos(x + z), x2 sin(x + z), x2y cos(x + z) . 1.2.2. Đạo hàm theo hướng Cho f là một hàm xác định trong một lân cận của điểm x0 ∈ Rn và v ∈ Rn là một vectơ khác không. Lúc đó, nếu giới hạn sau tồn tại ta gọi nó là đạo hàm của hàm f tại x0 theo hướng v: ∂f f(x0 + tv) − f(x0) f 0(x0; v) = (x0) := lim . ∂v t→0+ t Có thể kiểm chứng được rằng, nếu đạo hàm riêng theo biến x1 của f tồn tại thì với e1 = (1, 0, ··· , 0) ta có ∂f ∂f ∂f ∂f (x0) = (x0); (x0) = − (x0). ∂e1 ∂x1 ∂(−e1) ∂x1 Ngược lại, nếu tồn tại đạo hàm của f theo các hướng ±e1 có giá trị đối nhau thì đạo hàm riêng của f theo biến x1 cũng tồn tại. Các bạn tự phát biểu và chứng minh các khẳng định tương tự đối với e2, ··· , en. Chú ý. Một hàm có các đạo hàm riêng, thậm chí có đạo hàm theo mọi hướng, tại một điểm có thể không liên tục tại điểm đó. Chẳng hạn, trong Ví dụ 1.1, nếu ta định nghĩa thêm g(0, 0) = 0 thì có thể kiểm chứng được hàm g xác định trên R2, có 0 0 các đạo hàm riêng gx, gy nhưng g không liên tục tại (0, 0).
  8. 7 1.2.3. Vi phân Cho hàm y = f(x) xác định trong một lân cận V của điểm x0. Với các số gia 0 ∆xi đủ bé sao cho x + ∆x ∈ V , với ∆x = (∆x1, ··· , ∆xn), ta có số gia của hàm số là ∆f = f(x0 + ∆x) − f(x0). Nếu ∆f có thể biểu diễn dưới dạng Xn 0 ∆f = Ai∆xi + α(∆x)k∆xk; x + ∆x ∈ V, i=1 trong đó Ai, 1 ≤ i ≤ n, là các hằng số còn α là hàm n biến sao cho lim α(∆x) = 0, ∆x→0 thì f được gọi là khả vi tại điểm x0 và biểu thức Xn dy = df := Ai∆xi i=1 được gọi là vi phân của hàm f tại điểm x0 (tương ứng với vectơ gia ∆x). Mệnh đề 1.6. Nếu f khả vi tại x0 thì f liên tục tại điểm đó. Mệnh đề 1.7. Nếu f khả vi tại x0 thì f có các đạo hàm riêng tại điểm đó và Xn ∂f df = h∇f(x0), ∆xi = (x0)∆x . (1.2) ∂x i i=1 i Hơn nữa, f có đạo hàm theo mọi hướng tại x0 và ∂f (x0) = h∇f(x0), vi; ∀v ∈ Rn. ∂v Vì một hàm có các đạo hàm riêng tại một điểm có thể không liên tục tại điểm đó nên cũng không khả vi tại điểm đó. Tuy nhiên ta có kết quả sau Định lý 1.8. Nếu f có các đạo hàm riêng trong một lân cận của x0 và các đạo hàm này liên tục tại x0, thì f khả vi tại điểm dó. Nếu gi là hàm chiếu xuống tọa độ thứ i: gi(x1, ··· , xn) = xi thì ta sẽ ký hiệu dxi := dgi. Mặt khác, gi khả vi tại mọi điểm và dgi = ∆xi. Vậy, dxi = ∆xi. Do đó công thức (1.2) có thế viết lại: Xn ∂f df = dx . (1.3) ∂x i i=1 i
  9. 8 1.2.4. Đạo hàm hàm số hợp và tính bất biến của vi phân n Cho y = f(x1, x2, ··· , xn) là hàm xác định trên tập mở G ⊂ R và xi = ϕi(t), 1 ≤ i ≤ n, là n hàm số thực xác định trên khoảng (a, b) sao cho (ϕ1(t), ϕ2(t), ··· , ϕn(t)) ∈ G; ∀t ∈ (a, b). Lúc đó, ta có hàm hợp t −→ y = f(ϕ1(t), ϕ2(t), ··· , ϕn(t)) =: g(t) từ (a, b) vào R. 0 Định lý 1.9. Nếu các hàm ϕi khả vi tại t0 ∈ (a, b) còn hàm f khả vi tại x = (ϕ1(t0), ··· , ϕn(t0)), thì g cũng khả vi tại t0 và Xn ∂f g0(t ) = (x0)ϕ0 (t ). 0 ∂x i 0 i=1 i Nếu các hàm ϕi khả vi trên (a, b) và f khả vi trên G, thì g cũng khả vi trên (a, b) và Xn ∂f g0(t) = (ϕ (t), ··· , ϕ (t))ϕ0 (t). ∂x 1 n i i=1 i Từ định lý trên ta thường viết dy Xn ∂y dx = i , dt ∂x dt i=1 i hay Xn ∂y dy = dx . (1.4) ∂x i i=1 i n Bây giờ giả sử y = f(x1, ··· , xn) là hàm khả vi trên tập mở G ⊂ R và m xi = ϕi(u) = ϕi(u1, ··· , um), 1 ≤ i ≤ n, là các hàm khả vi trên một tập mở E ⊂ R sao cho (ϕ1(u), ··· , ϕn(u)) ∈ G với mọi u ∈ E. Lúc đó ta có hàm hợp g : E → R là một hàm m biến, xác định bởi g(u) = f(ϕ1(u), ··· , ϕn(u)); u ∈ E. Bằng cách sử dụng Định lý 1.9 và xem g là hàm theo một biến uj ta có ∂g Xn ∂f ∂ϕ = i ; 1 ≤ j ≤ m. ∂u ∂x ∂u j i=1 i j Từ đó, ta nhận được vi phân của hàm g: ( ) ( ) Xm ∂g Xm Xn ∂f ∂ϕ Xn ∂f Xm ∂ϕ dy = du = i du = i du . ∂u j ∂x ∂u j ∂x ∂u j j=1 j j=1 i=1 i j i=1 i j=1 j
  10. 9 Lại sử dụng Định lý 1.9 cho các hàm xi = ϕi(u) ta được Xn ∂y dy = dx . (1.5) ∂x i i=1 i Đối chiếu (1.3), (1.4) và (1.5) ta thấy dạng vi phân của y không hề thay đổi cho m dù xi là các biến độc lập, hàm của một biến t ∈ R hay là hàm của m biến u ∈ R . Ta nói dạng vi phân bậc nhất có tính bất biến. Vận dụng các kết quả trên một cách thích hợp ta có các công thức tính vi phân sau Hệ quả 1.2. Cho u và v là các hàm nhiều biến khả vi trên miền chung E ⊂ Rn. Lúc đó, trên miền này ta có a) d(u ± v) = du ± dv; b) d(λu) = λdu, λ ∈ R; c) d(u.v) = udv + vdu; ³u´ vdu − udv d) d = , v 6= 0; v v2 1.2.5. Đạo hàm hàm ẩn Cho F (x, y), x ∈ Rn, y ∈ R là một hàm n + 1 biến, xác định trong một tập mở G ⊂ Rn+1. Xét phương trình F (x, y) = 0. (1.6) Nếu tồn tại hàm n biến y = f(x); x ∈ E ⊂ Rn sao cho F (x, f(x)) = 0; ∀x ∈ E, thì f được gọi là hàm ẩn xác định bởi phương trình (1.6). 0 0 Định lý 1.10. Giả sử hàm hai biến F liên tục cùng với các đạo hàm Fx, Fy trong 2 0 một lân cận của điểm (x0, y0) ∈ R . Ngoài ra, F (x0, y0) = 0; Fy(x0, y0) 6= 0. Lúc đó a) Tồn tại duy nhất một hàm y = f(x) thoả mãn f(x0) = y0 và F (x, f(x)) = 0 với mọi x thuộc vào một lân cận ∆ = [x0 − δ, x0 + δ] của x0, b) f liên tục, có đạo hàm liên tục trên ∆ và 0 0 Fx(x, f(x)) f (x) = − 0 , ∀x ∈ ∆. Fy(x, f(x)) Định lý 1.11. Giả sử hàm n + 1 biến F (x1, ··· , xn, y) liên tục cùng với các đạo 0 0 0 0 0 n+1 hàm riêng Fx1 , ··· ,Fxn , Fy trong một lân cận của điểm (x , y ) ∈ R . Ngoài ra, 0 0 0 0 0 F (x , y ) = 0; Fy(x , y ) 6= 0. Lúc đó,
  11. 10 a) Tồn tại duy nhất một hàm y = f(x) thoả mãn f(x0) = y0 và F (x, f(x)) = 0 0 0 0 0 0 với mọi x thuộc vào một lân cận ∆ = [x1 − δ, x1 + δ] × · · · × [xn − δ, xn + δ] của x , b) f liên tục, có các đạo hàm riêng liên tục trên ∆ và F 0 (x, f(x)) f 0 (x) = − xi , ∀x ∈ ∆, 1 ≤ i ≤ n. xi 0 Fy(x, f(x)) n m Bây giờ cho Fi(x, y), x ∈ R , y ∈ R , 1 ≤ i ≤ m, là m hàm n + m biến, xác định trong một tập mở G ⊂ Rn+m. Xét hệ phương trình Fi(x, y) = 0; 1 ≤ i ≤ m. (1.7) n Nếu tồn tại m hàm n biến yi = fi(x); x ∈ E ⊂ R , 1 ≤ i ≤ m sao cho Fi(x, f1(x), ··· , fm(x)) = 0; ∀x ∈ E, 1 ≤ i ≤ m, thì {fi | 1 ≤ i ≤ m} được gọi là hệ hàm ẩn xác định bởi hệ phương trình (1.7). Nếu tồn tại các đạo hàm riêng của các hàm Fi theo các biến yj thì định thức sau được gọi là Định thức Jacobi của hệ hàm Fi đối với các biến yj:   ∂F1 (x, y) ··· ∂F1 (x, y) ∂y1 ∂ym  ∂F2 (x, y) ··· ∂F2 (x, y)  ∂y1 ∂ym  DJy(x, y) := det  . . .  .  . .  ∂Fm (x, y) ··· ∂Fm (x, y) ∂y1 ∂ym Định lý 1.12. Giả sử các hàm Fi(x1, ··· , xn, y1, ··· , ym) liên tục cùng với các đạo 0 0 n+m hàm riêng ∂Fi/∂yj, 1 ≤ i, j ≤ m, trong một lân cận của điểm (x , y ) ∈ R . 0 0 0 0 Ngoài ra, F (x , y ) = 0 và DJy(x , y ) 6= 0. Lúc đó, 0 0 a) Tồn tại duy nhất hệ hàm yi = fi(x), 1 ≤ i ≤ m, thoả mãn fi(x ) = yi và 0 Fi(x, f1(x), ··· , fm(x)) = 0 với mọi x thuộc vào một lân cận ∆ của điểm x , b) Các fi liên tục, có các đạo hàm riêng liên tục trên ∆. Hơn nữa, nếu đặt yi = fi(x1, ··· , xn), thì với mọi i ∈ {1, ··· , m}, j ∈ {1, ··· , n} ta có ∂f DJ (x, y) i (x) = − (y1,··· ,yi−1,xj ,yi+1,··· ,ym) . ∂xj DJy(x, y) Hệ quả 1.3 (Đạo hàm hàm ngược). Giả sử F : D ⊂ Rm → Rm sao cho các hàm thành phần Fi(y) liên tục cùng với các đạo hàm riêng ∂Fi/∂yj, 1 ≤ i, j ≤ m, trên 0 0 tập mở D 3 y . Ngoài ra, ma trận Jacobian JF (y ) = (∂Fi/∂yj) không suy biến. Lúc đó tồn tại một lân cận U của y0 và một lân cận V của z0 = F (y0) và một ánh xạ F −1 : V → U, có các hàm thành phần khả vi liên tục, thoả mãn a) ∀y ∈ U, ∀z ∈ V : z = F (y) ⇔ y = F −1(z) b) ∀z ∈ V : J(F −1)(z) = JF (y)−1, với y = F −1(z).
  12. 11 1.3. Đạo hàm cấp cao và Công thức Taylor 1.3.1. Đạo hàm cấp cao Để đơn giản trước tiên ta xét hàm hai biến. Cho z = f(x, y) là hàm xác định 2 0 0 trên tập mở G ⊂ R , có các đạo hàm riêng fx(x, y), fy(x, y) trên G. Đây cũng là các hàm hai biến. Nếu các hàm này cũng có các đạo hàm riêng thì các đạo hàm đó được gọi là đạo hàm riêng cấp 2 của f. Nói chung f có 4 đạo hàm riêng cấp 2: µ ¶ 2 ∂ ∂f ∂ f 00 00 =: =f 2 =z 2 ; ∂x ∂x ∂x2 x x µ ¶ ∂ ∂f ∂2f =: =f 00 =z00 ; ∂x ∂y ∂x∂y yx yx µ ¶ ∂ ∂f ∂2f =: =f 00 =z00 ; ∂y ∂x ∂y∂x xy xy µ ¶ 2 ∂ ∂f ∂ f 00 00 =: =f 2 =z 2 . ∂y ∂y ∂y2 y y Tương tự, ta có các khái niệm đạo hàm cấp cao hơn và của những hàm nhiều biến hơn. Chẳng hạn, với hàm u = f(x, y, z) ta có đạo hàm riêng cấp 4: µ µ µ ¶¶¶ ∂4f ∂ ∂ ∂ ∂f u(4) = := . xyzx ∂x∂z∂y∂x ∂x ∂z ∂y ∂x 00 00 000 000 000 Nói chung, fxy 6= fyx, fx2y 6= fxyx 6= fyx2 . Tuy nhiên những đạo hàm hỗn hợp này sẽ trùng nhau trong trường hợp chúng liên tục. Điều đó được thể hiện trong định lý sau Định lý 1.13. Giả sử z = f(x, y) là hàm xác định trên tập mở G, có các đạo hàm 00 00 riêng cấp hai hỗn hợp fxy, fyx. Nếu các đạo hàm này liên tục tại điểm (x0, y0) ∈ G, thì 00 00 fxy(x0, y0) = fyx(x0, y0). Chứng minh. Đặt ∆ := f(x0 + h, y0 + k) − f(x0 + h, y0) − f(x0, y0 + k) + f(x0, y0), g1(t) := f(t, y0 + k) − f(t, y0), với h và k lần lượt là số gia của x và y. Sử dụng Định 0 lý Lagrange cho g1 rồi cho fx ta tìm được các số δ, θ ∈ (0, 1) (phụ thuộc vào h, k) thoả mãn 0 ∆ = g1(x0 + h) − g1(x0) = g1(x0 + δh)h 0 0 = [fx(x0 + δh, y0 + k) − fx(x0 + δh, y0)]h 00 = fxy(x0 + δh, y0 + θk)hk.
  13. 12 Tương tự, nếu đặt g2(s) := f(x0 + h, s) − f(x0, s), và áp dụng Định lý Lagrange lần 0 lượt cho g2 rồi cho fy, ta cũng tìm được các số α, β ∈ (0, 1) thoả mãn 00 ∆ = fyx(x0 + αh, y0 + βk)hk. Từ đó: 00 00 fyx(x0 + αh, y0 + βk) = fyx(x0 + αh, y0 + βk). Cho h, k → 0 ta nhận được điều phải chứng minh. Định lý này cũng được mở rộng không mấy khó khăn cho các trường hợp đạo hàm cấp cao hơn, hoặc với hàm nhiều biến hơn với điều kiện các đạo hàm hỗn hợp đó liên tục. Chẳng hạn với hàm u = x3 sin(y + z2), các bạn có thể kiểm tra các đạo (4) (4) (4) (4) (4) 2 hàm ux2yz, uxyxz, uxyzx, uyxzx, uyzx2 , đều bằng nhau và bằng −12xz sin(y + z ). 1.3.2. Vi phân cấp cao Để đơn giản, trước tiên ta xét hàm hai biến. Cho z = f(x, y) là hàm xácđịnh và khả vi trên tập mở G ⊂ R2. Vi phân của f tại mỗi điểm (x, y) ∈ G là ∂f ∂f df(x, y) = (x, y)∆x + (x, y)∆y. ∂x ∂y Như vậy, df là một hàm hai biến trên G. Nếu df cũng khả vi thì vi phân của nó sẽ được gọi là vi phân cấp hai của f. Lúc đó, µ ¶ ∂ ∂ ∂f ∂f ∂2f ∂2f df(x, y) = (x, y)∆x + (x, y)∆y = (x, y)∆x + (x, y)∆y; ∂x ∂x ∂x ∂y ∂x2 ∂x∂y µ ¶ ∂ ∂ ∂f ∂f ∂2f ∂2f df(x, y) = (x, y)∆x + (x, y)∆y = (x, y)∆x + (x, y)∆y. ∂y ∂y ∂x ∂y ∂y∂x ∂y2 Tóm lại, vi phân cấp hai của f tương ứng với cặp số gia (∆x, ∆y) là ∂ ∂ d2f(x, y) := d(df)(x, y) = df(x, y).∆x + df(x, y).∆y ∂x ∂y µ ¶ µ ¶ ∂2f ∂2f ∂2f ∂2f = (x, y)∆x + (x, y)∆y ∆x + (x, y)∆x + (x, y)∆y ∆y ∂x2 ∂x∂y ∂y∂x ∂y2 ∂2f ∂2f ∂2f ∂2f = (x, y)∆x2 + (x, y)∆x∆y + (x, y)∆x∆y + (x, y)∆y2. ∂x2 ∂x∂y ∂y∂x ∂y2 Nếu các đạo hàm riêng cấp hai hỗn hợp liên tục thì theo Định lý 2.15 vi phân cấp hai của f có thể viết gọn hơn: ∂2f ∂2f ∂2f d2f(x, y) = (x, y)∆x2 + 2 (x, y)∆x∆y + (x, y)∆y2, ∂x2 ∂x∂y ∂y2
  14. 13 mà, để đơn giản người ta viết lại một cách hình thức như sau µ ¶ ∂ ∂ 2 d2f(x, y) = ∆x + ∆y f(x, y). ∂x ∂y Tương tự, ta cũng có định nghĩa của các vi phân cấp cao hơn. Hơn nữa, bằng quy nạp ta có thể chứng minh được mệnh đề sau Định lý 1.14. Nếu hàm hai biến f(x, y) có các đạo hàm riêng đến cấp m liên tục trên tập mở G ⊂ R2 thì f khả vi cấp m trên G và ta có vi phân cấp m của f tương ứng với cặp số gia (∆x, ∆y) là µ ¶ ∂ ∂ m Xm ∂mf dmf(x, y) = ∆x + ∆y f(x, y) := Ck ∆xk∆ym−k. ∂x ∂y m ∂xk∂ym−k k=0 Bằng một lược đồ tương tự ta nhận được khái niệm vi phân cấp cao của hàm nhiều biến cũng như công thức tính của nó. Cụ thể ta có mệnh đề Định lý 1.15. Nếu hàm nhiều biến f(x1, ··· , xn) có các đạo hàm riêng đến cấp m liên tục trên tập mở G ⊂ Rn thì f khả vi cấp m trên G và vi phân cấp m của f tương ứng với vectơ gia ∆x = (∆x1, ··· , ∆xn) là Ã !m Xn ∂ dmf(x , ··· , x ) = ∆x f(x , ··· , x ), (x , ··· , x ) ∈ G. 1 n ∂x i 1 n 1 n i=1 i 1.3.3. Công thức Taylor Định lý 1.16. Giả sử y = f(x) là một hàm có các đạo hàm riêng liên tục đến cấp m trên tập mở G ⊂ Rn, x0 là một điểm trong G và ∆x là vectơ sao cho đoạn thẳng [x0, x0 + ∆x] nằm gọn trong G. Lúc đó, tồn tại số θ ∈ (0, 1) sao cho mX−1 1 1 f(x0 + ∆x) = f(x0) + dkf(x0) + dmf(x0 + θ∆x), (1.8) k! m! k=1 trong đó, dkf(x) ký hiệu vi phân cấp k của f tại x tương ứng với vectơ gia ∆x. Chứng minh. Đặt F là ham một biến F (t) = f(x0 + t∆x). Khai triển MacLaurin hàm này đến cấp m ta có mX−1 F (k)(0) F (m)(θ) F (t) = F (0) + tk + tm. (1.9) k! m! k=1
  15. 14 Chú ý rằng Xn ∂f F 0(t) = (x0 + t∆x).∆x = df(x0 + t∆x), ∂x i i=1 i Xn ∂2f F 00(t) = (x0 + t∆x).∆x ∆x = d2f(x0 + t∆x), ··· ∂x ∂x i j i,j=1 i j F (m)(t) = dmf(x0 + t∆x). Thay vào (1.9) với t = 1 ta được điều phải chứng minh. (1.8) được gọi là Công thức Taylor đến cấp m của hàm f tại điểm x0, tương ứng với vectơ gia ∆x. Hệ quả 1.4 (Định lý giá trị trung bình). Giả sử f là hàm có các đạo hàm riêng liên tục trên tập mở U ⊂ Rn và a, b ∈ G là hai điểm phân biệt sao cho [a, b] ⊂ G. Lúc đó tồn tại điểm c ∈ (a, b) thoả mãn f(b) − f(a) = h∇f(c), b − ai. 1.4. Cực trị 1.4.1. Điều kiện cần Cho f là hàm nhiều biến xác định trên G ⊂ Rn. Ta nói hàm f đạt cực tiểu (cực đại) địa phương tại điểm x0 ∈ G nếu tồn tại số dương δ sao cho ¡ ¢ f(x0) ≤ f(x) f(x0) ≥ f(x) ; ∀x ∈ B(x0, δ) ∩ G. Trong cả hai trường hợp ta nói f đạt cực trị địa phương tại x0. Định lý 1.17. Nếu f đạt cực trị địa phương tại một điểm trong x0 của G, tại đó tồn tại các đạo hàm riêng của f, thì các đạo hàm này phải bằng 0. Tức là ∇f(x0) = 0. Một điểm tại đó gradiên của f bằng không được gọi là điểm dừng của f. Định lý 1.17 cho thấy mọi điểm cực trị của f đều là điểm dừng. Tuy vậy điều ngược lại nói chung không còn đúng. Chẳng hạn, hàm f(x, y) = x2 − y2 có ∇f(x, y) = (2x, −2y) với (x, y) ∈ R2, vì vậy hàm này có một điểm dừng là (0, 0) nhưng đó không phải là điểm cực trị. Thật vậy, trong một lân cận bé tuỳ ý của (0, 0) ta luôn tìm được hai điểm tại đó hàm f có một giá trị bé hơn f(0, 0) và một giá trị lớn hơn f(0, 0).
  16. 15 1.4.2. Điều kiện đủ Trước khi phát biểu điều kiện đủ cực trị ta nhắc lại Công thức Taylor đến cấp hai của f tại một điểm x0: 1 f(x0 + ∆x) = f(x0) + df(x0) + d2f(x0 + θ∆x), θ ∈ (0, 1). 2 Nếu f có các đạo hàm riêng đến cấp hai liên tục thì theo Định lý 1.15 Ã !2 Xn ∂f 1 Xn ∂ f(x0 + ∆x) = f(x0) + (x0)∆x + ∆x f(x0 + θ∆x) ∂x i 2 ∂x i i=1 i i=1 i Xn ∂f 1 Xn Xn ∂2f = f(x0) + (x0)∆x + (x0 + θ∆x)∆x ∆x . ∂x i 2 ∂x ∂x i j i=1 i i=1 j=1 i j Tóm lại, 1 f(x0 + ∆x) = f(x0) + h∇f(x0), ∆xi + ∆xT ∇2f(x0 + θ∆x)∆x, (1.10) 2 trong đó ∆xT là vectơ chuyển vị của ∆x còn ∇2f(x) ký hiệu ma trận Hessian của f tại một điểm x. Đó là ma trận vuông cấp n × n mà phần tử ở hàng i cột j chính 2 là ∂ f (x). Đây là một ma trận đối xứng khi các đạo hàm riêng cấp hai liên tục. ∂xi∂xj Nếu x0 là điểm dừng thì ∇f(x0) = 0, nên (1.10) được viết lại như sau 1 f(x0 + ∆x) − f(x0) = ∆xT ∇2f(x0 + θ∆x)∆x, 2 Nhắc lại rằng một ma trận A vuông cấp n × n được gọi là xác định dương (nửa xác định dương, xác định âm, nửa xác định âm) nếu uT Au > 0 (uT Au ≥ 0) (uT Au < 0) (uT Au ≤ 0); ∀u ∈ Rn \{0}. A được gọi là không xác định dấu nếu tồn tại hai vectơ u, v ∈ Rn sao cho uT Au < 0 < vT Av. Định lý 1.18. Gỉa sử f có đạo hàm đến cấp hai liên tục trên một tập mở G, nhận điểm x0 ∈ G làm điểm dừng. Lúc đó nếu ∇2f(x) nửa xác định dương (nửa xác định âm) trong một lân cận của x0, thì x0 là điểm cực tiểu (cực đại) địa phương. Bây giờ ta nhắc lại một kết quả quen biết trong đại số tuyến tính. Cho A = (aij) là một ma trận thực vuông đối xứng cấp n. Ta đặt   a11 a12 ··· a1k   a21 a22 ··· a2k ∆k(A) := det  . . . .  , 1 ≤ k ≤ n.  . . .  ak1 ak2 ··· akk
  17. 16 Định lý 1.19. A là ma trận xác định dương (âm) khi và chỉ khi k ∆k(A) > 0 ((−1) ∆k(A) > 0); 1 ≤ k ≤ n. Định lý 1.20. Gỉa sử f có đạo hàm đến cấp hai liên tục trên một tập mở G, nhận điểm x0 ∈ G làm điểm dừng. Lúc đó a) Nếu ∇2f(x0) xác định dương thì x0 là điểm cực tiểu. b) Nếu ∇2f(x0) xác định âm thì x0 là điểm cực đại. c) Nếu ∇2f(x0) không xác định dấu thì x0 không phải là điểm cực trị. Khi f là hàm hai biến ta có Ã ! ∂2f ∂2f 2 ∂x2 (x, y) ∂x∂y (x, y) ∇ f(x, y) = ∂2f ∂2f . ∂y∂x (x, y) ∂y2 (x, y) Do đó, bằng cách đặt ∂2f ∂2f ∂2f A := (x , y ); B := (x , y ); C := (x , y ) ∂x2 0 0 ∂x∂y 0 0 ∂y2 0 0 2 2 2 ta có ∆1(∇ f(x0, y0)) = A và ∆2(∇ f(x0, y0)) = AC − B =: D. Từ Định lý 1.20 ta có hệ quả sau Hệ quả 1.5. Gỉa sử f(x, y) nhận (x0, y0) là điểm dừng. Ngoài ra, các đạo hàm riêng đến cấp hai tồn tại và liên tục tại (x0, y0). Lúc đó a) Nếu D > 0 và A > 0 thì (x0, y0) là điểm cực tiểu. b) Nếu D > 0 và A < 0 thì (x0, y0) là điểm cực đại. c) Nếu D < 0 thì (x0, y0) không phải là điểm cực trị. Như vậy, nếu D = 0 thì ta vẫn chưa xác định được (x0, y0) có phải là điểm cực trị hay không. 1.4.3. Cực trị có điều kiện Ta xét bài toán cực trị  min f(x1, x2, ··· , xn),  g (x , x , ··· , x ) = 0, (P) 1 1 2 n ···  gm(x1, x2, ··· , xn) = 0, n ở đó m ≤ n. Tập hợp X := {x ∈ R | gi(x) = 0; 1 ≤ i ≤ m} được gọi là tập chấp nhận được của bài toán và mỗi điểm x ∈ X được gọi là một điểm chấp nhận được.
  18. 17 Một điểm x¯ ∈ X được gọi là nghiệm (địa phương) của bài toán (P) nếu tồn tại lân cận U của x¯ sao cho f(¯x) ≤ f(x); ∀x ∈ U ∩ X. Định lý sau cho ta một điều kiện cần của cực trị có điều kiện: Định lý 1.21. Giả sử x¯ là một nghiệm của bài toán (P) tại đó các hàm f và gi khả vi liên tục. Hơn nữa, ma trận Jacobi   ∂g1 (¯x) ∂g1 (¯x) ··· ∂g1 (¯x) ∂x1 ∂x2 ∂xn  ∂g2 (¯x) ∂g2 (¯x) ··· ∂g2 (¯x)  ∂x1 ∂x2 ∂xn   . . . .   . . .  ∂gm (¯x) ∂gm (¯x) ··· ∂gm (¯x) ∂x1 ∂x2 ∂xn ¯ ¯ ¯ có hạng bằng m. Lúc đó tồn tại các số thực λ1, λ2 ··· , λm sao cho Xm ¯ ∇f(¯x) − λi∇gi(¯x) = 0. i=1 ¯ ¯ Các số λ1, , λm trong định lý trên được gọi là các nhân tử Lagrange của bài toán (P) đối với điểm cực trị x¯. 1.5. Thực hành tính toán trên Maple 1.5.1. Giới hạn và đồ thị hàm nhiều biến a) Định nghĩa hàm nhiều biến số. Để đơn giản ta chỉ xét hàm hai biến. Cú pháp: [> f:=(x, y)− > (biểu thức hàm theo x, y); Ví dụ: [> f:= (x, y)− > 3*x∧2*sin(x*y); f := (x, y) → 3x2sin(xy) b) Tính giới hạn của hàm nhiều biến. Chẳng hạn, ta xét hàm hai biến, trường hợp tổng quát được viết tương tự. Cú pháp: [> limit(f(x, y), {x=a, y=b}); (Limit sẽ cho công thức hình thức) Chú ý rằng nếu viết limit(limit(f(x,y), x=a), y=b) thì máy sẽ tính giới hạn lặp, trước tiên theo x và sau đó theo y. Có khi hàm không tồn tại giới hạn tại (a, b) nhưng vẫn tồn tại các giới hạn lặp. Ví dụ:
  19. 18 [> limit(x*y/(x∧2+y∧2), {y=0,x=0}); µ ¶ xy limit , {y = 0, x = 0} x2 + y2 Tức là máy không tính nổi giới hạn này (vì không tồn tại). Tuy nhiên: [> limit(limit(x*y/(x∧2+y∧2), y=0), x=0}); 0 Muốn vẽ đồ thị hàm nhiều biến ta cần khởi động các gói lệnh plottools, plots. c) Vẽ đồ thị hàm z = f(x, y). Cú pháp: [> plot3d(f(x,y), x=a b, y=c d); Lúc đó, đồ thị là một mặt trong không gian Oxyz với miền xác định là hình chữ nhật [a, b] × [c, d]. Ví dụ: [> f:= (x, y) − > x∧2 + y∧2: [> with(plots): with(plottools): [> plot3d(f(x,y), x=-3 3, y=-2 2); Hình 1.1: Đồ thị hàm z = x2 + y2 Nếu vẽ nhiều mặt trên cùng một không gian toạ độ thì ta viết Cú pháp: [> plot3d({f(x,y), g(x,y), }, x=a b, y=c d); Ví dụ: (xem Hình 1.2) [> plot3d({x∧2+y∧2, sqrt(1-x∧2-y∧2)},y=0 sqrt(1-x∧2),x=-1 1); d) Vẽ mặt được cho dưới dạng tham số. Giả sử mặt S được cho bởi hệ  x = x(u, v), y = y(u, v), (u, v) ∈ [a, b] × [c, d].  z = z(u, v),
  20. 19 p Hình 1.2: Đồ thị các hàm z = x2 + y2 và z = 1 − x2 − y2 Để vẽ mặt S ta dùng lệnh (chú ý đừng nhầm lẫn với lệnh vẽ nhiều mặt cùng lúc) Cú pháp: [> plot3d([x(u,v), y(u,v), z(u,v)], u=a b, v=c d); e) Vẽ mặt được cho bởi phương trình ẩn dạng F (x, y, z) = 0. Cú pháp: [> implicitplot3d(F(x,y,z)=0,x=a b, y=c d, z=e f); x2 y2 2 Ví dụ: Để vẽ mặt 4 + 9 − z = 1 trong hình hộp [−5, 5] × [−6, 6] × [−1 1], ta viết [> implicitplot3d(x∧2/4+y∧2/9- z∧2-1 =0, x=-5 5, y=-6 6, z=-1 1); Hình 1.3: Đồ thị hàm ẩn x2/4 + y2/9 − z2 − 1 = 0 f) Vẽ các đường mức của một hàm hai biến. Cú pháp: [> contourplot(f(x,y), x=a b, y=c d); Lúc đó, máy sẽ vẽ trên mặt phẳng Oxy các đường cong dạng f(x, y) = α, với các α khác nhau.
  21. 20 1.5.2. Tính đạo hàm Việc tính đạo hàm của hàm một biến hoặc đạo hàm riêng của hàm nhiều biến được thực hiện tương tự nhau. a) Tính đạo hàm riêng cấp một. Cú pháp: [> diff(hàm số, tên biến riêng); (dùng Diff thì cho công thức hình thức) Ví dụ: [> diff(sqrt(1+yx∧2), x); yx √ 1 + x2 [> f:=(x,y)− >x*sin(x*y): [> diff(f(x,y), y); x2 cos(xy) [> Diff(f(x,y), x); ∂ x sin(xy) ∂x b) Tính đạo hàm cấp cao. Cú pháp: [> diff(f(x,y), x$k, y$m); Ví dụ: (với hàm f như trên) [> diff(f(x, y ),x$2, y$3); 8x3 sin(xy)y − 12x2 cos(xy) + x4 cos(xy)y2 c) Tính građiên (chẳng hạn, hàm ba biến): Cú pháp: [> grad(f(x,y,z), [x,y,z]); d) Tính Hessian. Cú pháp: [> hessian(f(x,y,z), [x,y,z]); Ví dụ: (với f ở trên) [> grad(f(x,y), [x,y]); [sin(xy) + x cos(xy)y, x2 cos(xy)] [> hessian(f(x,y), [x,y]); · ¸ 2 cos(xy)y − x sin(xy)y2 2 cos(xy)x − x2 sin(xy)y 2 cos(xy)x − x2 sin(xy)y −x3 sin(xy)
  22. 21 1.5.3. Khai triển Taylor Khai triển Taylor hàm f(x, y) tại (a, b) đến cấp n: Cú pháp: [> mtaylor(f(x,y), [x=a, y=b], n); Nếu khai triển tại (0, 0) ta chỉ cần viết [x, y] thay cho [x = a, y = b]; nếu không khai báo n thì ngầm định n = 6. Ví dụ: [> mtaylor(sin(x + y∧3), [x, y ], 8); 1 1 1 1 1 1 x + y3 − x3 − y3x2 + x5 − x7 − y6x + y3x4 6 2 120 5040 2 24 [> mtaylor(sin(x + y∧3), [x, y ]); 1 1 1 x + y3 − x3 − y3x2 + x5 6 2 120 1.6. Bài tập 1.1. Cho hàm hai biến ( x2 + y2, x 6= 0; f(x, y) = 0, x = 0. 0 0 Cho biết hàm này liên tục tại những điểm nào? Tính fx(0, 0), fy(0, 0). 1.2. Tìm, hoặc chứng minh không tồn tại, các giới hạn lim(lim f(x, y)); lim(lim f(x, y)); lim f(x, y)) x→0 y→0 y→0 x→0 (x,y)→(0,0) với |x + y| f(x, y) = . x2 + y2 − xy 1.3. Tìm, hoặc chứng minh không tồn tại, các giới hạn lim(lim f(x, y)); lim(lim f(x, y)); lim f(x, y)) x→0 y→0 y→0 x→0 (x,y)→(0,0) với x2 + y2 f(x, y) = . x2 + y2 + (x − y)2
  23. 22 1.4. Cho hàm f : R × R → R xác định bởi ( 3 x y , (x, y) 6= (0, 0), f(x, y) := x6+y2 0, (x, y) = (0, 0). 0 Chứng minh rằng hàm f gián đoạn tại điểm (0, 0) nhưng các đạo hàm riêng fx(0, 0), 0 fy(0, 0) tồn tại. 1.5. Khảo sát tính liên tục và sự tồn tại các đạo hàm riêng cấp một của hàm hai biến ( (x2 + y2) sin √ 1 , (x, y) 6= (0, 0); f(x, y) = x2+y2 0, (x, y) = (0, 0). 1.6. Phác hoạ đồ thị của các hàm hai biến sau. Khảo sát tính liên tục và sự tồn tại các đạo hàm riêng cấp một của chúng (p 1 − x2 − y2, x2 + y2 0 với mọi x ∈ S(0, 1). Chứng minh rằng phương trình f(x) = 0 có vô số nghiệm. 1.10. Cho f : Rn → R, với mỗi x ∈ Rn và v ∈ Rn ta định nghĩa hàm một biến thực: gx,v(t) := f(x + tv); t ∈ R. n a) Chứng minh f lồi ⇔ gx,v lồi với mọi x, v ∈ R . b) Từ đó suy ra nếu f khả vi cấp hai trên Rn thì f lồi khi và chỉ khi ∇2f(x) nửa xác định dương tại mọi điểm x ∈ Rn. 1.11. Dùng vi phân cấp một của hàm hai biến để tính gần đúng biểu thức sau ln(e0.15+3 sin(0.01) + cos(0.01) − 1).
  24. 23 µ ¶ 1.02 1.12. Tính gần đúng biểu thức tan , biết π ≈ 0.785. 0.96 4 π 1.13. Tính gần đúng biểu thức arcsin (0.49 + ln(1.02)), biết ≈ 0.523. 6 1.14. Viết biểu thức vi phân toàn phần của hàm z = 5x3 sin(x2 + y)ey. 1.15. Viết biểu thức vi phân toàn phần của hàm u = 5 sin(x + y2z) cos(xy) phụ theo dt, biết x = 3t2 + 1, y = ln t, z = et. 1.16. Viết biểu thức vi phân toàn phần của hàm z = 3x2yex phụ theo dr và dϕ, biết x = r cos ϕ, y = r sin ϕ. 1.17. Cho hàm 2 biến ( 2 2 xy x −y , x2 + y2 6= 0, f(x, y) = x2+y2 0, x2 + y2 = 0. 00 00 Chứng minh rằng fxy(0, 0) 6= fyx(0, 0). 1.18. Khảo sát sự tồn tại hàm ẩn dạng y = f(x) hoặc x = g(y) của phương trình 3xex+y−cos(π ln(y+1)) + (x2 − 4) sin y − 3 = 0 trong lân cận điểm M(1, 0). Nếu tồn tại, hãy xác định đạo hàm của hàm ẩn tại điểm tương ứng. 1.19. Cho hàm 3 biến f(x, y, z) = 5x2 cos(y + z2). Xác định ma trận Hessian ∇2f(1, π, 0) và tính định thức của nó. 1.20. Cho hàm 2 biến f(x, y) = 3(x + 1)ey+x2 . Hãy tính vi phân df(1, 0) và vi phân cấp hai d2f(1, 0) tương ứng với vectơ gia (∆x, ∆y). 1.21. Khảo sát sự tồn tại hàm ẩn y = f(x) của các phương trình sau a) 5x2ey + 3y cos x = 0 tại (0, 0); π π b) sin x + cos y − 1 = 0 tại ( 2 , 2 ). Tính đạo hàm hàm ẩn tại các điểm tương ứng. 1.22. Khảo sát cực trị của các hàm hai biến thực: z = x3 + y3 − 3xy; z = 2x4 + y4 − x2 − 2y2; z = x3y + xy3 − xy; z = xy2 − x2y; z = 2x2 + 3y2 − x4; z = xy2 + x2y; x + y + 1 1 z = ; z = ; x2 + y2 + 1 1 + x2 + y2 p z = 1 − x2 + y2; z = (x2 + y2)e−x2−y2 ; p z = xy3 + x3y; z = xy 1 − 4x2 − 9y2. 1.23. Chứng minh hàm số sau có vô số điểm dừng, trong đó có một điểm cực đại z = (x2 + y2)2 − 2x2 − 2y2.
  25. Chương 2 ỨNG DỤNG CỦA PHÉP TÍNH VI PHÂN TRONG HÌNH HỌC 2.1. Các hệ toạ độ 2.1.1. Hệ toạ độ cực Để biểu diễn một điểm trong mặt phẳng ta thường sử dụng hệ toạ độ Đêcac vuông góc. Nghĩa là người ta sử dụng hai trục toạ độ vuông góc Ox, Oy cắt nhau tại một điểm O. Mỗi điểm M trong mặt phẳng được xác định hoàn toàn bởi cặp số thực (x, y) ∈ R2; trong đó x = OP , y = OQ với P và Q lần lượt là hình chiếu vuông góc của M lên Ox và Oy. Nếu bây giờ ta đặc trưng mỗi điểm M trong mặt phẳng bởi cặp giá trị (r, ϕ), trong đó r là độ dài đoạn OM và ϕ là góc định hướng lập bởi trục Ox và vectơ −−→ OM, thì cặp số này sẽ được gọi là toạ độ cực của điểm M. Rõ rằng, mỗi cặp (r, ϕ) ∈ G = R+ × R xác định duy nhất một điểm trên mặt phẳng. Tuy vậy, tương ứng với mỗi điểm M trên mặt phẳng thì có thể có nhiều cặp toạ độ cực khác nhau cho nó với các giá trị ϕ sai khác một bội của 2π. Đặc biệt để biểu diễn điểm O ta có thể dùng cặp (0, ϕ) với ϕ tuỳ ý. Để việc biểu diễn là duy nhất (đối với các điểm M 6= O) người ta thường thu hẹp phạm vi của ϕ, chẳng hạn chỉ cho phép (r, ϕ) ∈ G = [0, +∞) × [0, 2π) hoặc (r, ϕ) ∈ G = [0, +∞) × [−π, π) Mối liên hệ giữa toạ độ cực và toạ độ Đêcac được cho bởi hệ sau ( x = r cos ϕ, y = r sin ϕ. Lúc đó, phương trình các đường trong mặt phẳng f(x, y) = 0
  26. 25 cũng có thể được cho dưới dạng g(r, ϕ) = 0 hay ngược lại. Chẳng hạn đường thẳng x + y = 1, đường tròn đơn vị x2 + y2 = 1 và đường parabol y = x2 có thể được viết lại lần lượt là r(cos ϕ + sin ϕ) = 1, r = 1 và r sin ϕ − r2 cos2 ϕ = 0. Ngoài ra có một số đường mà cách biểu diễn theo hệ toạ độ này cho chúng ta một ý nghĩa hình học sáng sủa hơn hẳn dùng toạ độ khác. Chẳng hạn đường thẳng y = 1 hoàn toàn không nên viết r sin ϕ = 1. Hay ngược lại đường xoáy trôn ốc r = ϕ rất khó tìm được một cách biểu diễn trong toạ độ Đêcac. 2.1.2. Hệ toạ độ trụ Tương tự trong mặt phẳng, trong không gian ta thường sử dụng hệ toạ độ Đêcac vuông góc. Nghĩa là người ta sử dụng ba trục toạ độ vuông góc Ox, Oy, Oz cắt nhau tại một điểm O. Mỗi điểm M trong không gian được xác định hoàn toàn bởi bộ ba số thực (x, y, z) ∈ R3; trong đó x = OP , y = OQ, z = OR với P , Q và R lần lượt là hình chiếu vuông góc của M lên Ox, Oy và Oz. Nếu gọi N là hình chiếu của M lên mặt phẳng Oxy và ta đặc trưng cho điểm M trong không gian bởi bộ ba (r, ϕ, z) trong đó z = OR, r là độ dài đoạn ON và ϕ −−→ là góc định hướng lập bởi trục Ox và vectơ ON, thì bộ ba này sẽ được gọi là toạ độ trụ của điểm M. Nói cách khác, ta thay cặp (x, y) trong toạ độ Đêcac bởi toạ độ cực tương ứng (r, ϕ) và giữ nguyên toạ độ thứ ba: z. Cách xác định toạ độ như vậy cũng không duy nhất. Nếu thu hẹp phạm vi toạ độ là (r, ϕ, z) ∈ [0, +∞) × [0, 2π) × R thì việc biểu diễn là duy nhất đối với các điểm không thuộc trục Oz. Mối liên hệ giữa toạ độ trụ và toạ độ Đêcac được cho bởi hệ sau  x = r cos ϕ, y = r sin ϕ,  z = z. Lúc đó, phương trình r2 + z2 = 1 xác định mặt cầu đơn vị còn phương trình z = r xác định một mặt nón hướng về chiều dương của trục Oz. 2.1.3. Hệ toạ độ cầu Trong không gian người ta còn sử dụng hệ toạ độ cầu để biểu diễn các điểm. Cụ thể, ứng với một điểm M trong không gian có hình chiếu vuông góc lên mặt phẳng Oxy là N, ta gọi toạ độ cầu của M là bộ ba (ρ, ϕ, θ) ∈ [0, +∞) × [0, 2π) × [−π, π] −−→ trong đó ρ là độ dài đoạn OM, ϕ là góc lập bởi Ox với vectơ ON còn θ là góc lập −−→ bởi Oz và vectơ OM.
  27. 26 Mối liên hệ giữa toạ độ cầu và toạ độ Đêcac được cho bởi hệ sau  x = ρ sin θ cos ϕ, y = ρ sin θ sin ϕ,  z = ρ cos θ. Dễ thấy rằng, mặt cầu đơn vị bây giờ có phương trình là ρ = 1 còn mặt nón π r = z trong toạ độ trụ bây giờ có phương trình trong toạ độ cầu là θ = 4 . Trong khi đó, mặt parabol z = x2 + y2 có phương trình cos θ = ρ sin2 θ. 2.2. Hàm vectơ 2.2.1. Khái niệm Cho I là một khoảng trong R. Lúc đó, một ánh xạ f :I −→ Rn; t ∈I −→ f(t) ∈ Rn được gọi là một hàm vectơ (n chiều) biến số t xác định trong I. Với mỗi giá trị t ∈ I, f(t) là một n−vectơ với n thành phần: f(t) = (f1(t), ··· , fn(t)). Vì vậy việc cho một hàm vectơ f(t) tương đương với việc cho n hàm số thực fi(t), 1 ≤ i ≤ n, gọi là các hàm thành phần của f. Tập hợp C := {f(t) | t ∈ I} được gọi là một đường cong trong Rn. Khi n = 2 ta nhận được một đường cong phẳng và khi n = 3 ta có đường cong trong không gian. Thông thường, đường cong phẳng được cho bởi hàm (x(t), y(t)) hay được viết dưới dạng hệ ( x = x(t), t ∈ I. y = y(t), Tương tự, đường cong trong không gian được xác định bởi  x = x(t), y = y(t), t ∈ I.  z = z(t), Chẳng hạn, hệ  x = cos t, y = sin t, t ∈ R  z = t, cho ta đường xoắn dạng lò xo bao quanh mặt trụ x2 + y2 = 1.
  28. 27 2.2.2. Giới hạn - Liên tục - Đạo hàm n n Cho hàm vectơ f : I → R và t0 ∈ I. Ta nói vectơ v ∈ R là giới hạn của hàm f tại t0 và ký hiệu v = lim f(t) t→t0 nếu ∀² > 0, ∃δ > 0, ∀t ∈ I : 0 < |t − t0| < δ ⇒ kf(t) − vk < ². Sử dụng Bổ đề 3.3 trong giáo trình Giải tích 2 dễ chứng minh được rằng v = lim f(t) ⇐⇒ vi = lim fi(t); 1 ≤ i ≤ n. (2.1) t→t0 t→t0 Hàm f được gọi là liên tục tại t0 nếu f(t0) = lim f(t), t→t0 được gọi là liên tục nếu nó liên tục tại mọi điểm thuộc khoảng I. Từ (2.1) ta có ngay kết quả sau Mệnh đề 2.1. Hàm vectơ f liên tục trên I (tại t0) khi và chỉ khi các hàm thành phần fi cũng liên tục trên I (tại t0). Đạo hàm của hàm vectơ f tại điểm t0 được định nghĩa là giới hạn sau nếu nó tồn tại: 0 f(t0 + ∆t) − f(t0) f (t0) := lim . ∆t→0 ∆t Mệnh đề 2.2. Hàm vectơ f có đạo hàm tại t0 khi và chỉ khi các hàm thành phần fi cũng có đạo hàm tại điểm đó. Hơn nữa, ta có 0 0 0 f (t0) = (f1(t0), ··· , fn(t0)). Mệnh đề 2.3. Cho f và g là các hàm vectơ n chiều, có đạo hàm tại t0 và λ là một số thực. Lúc đó, các hàm vectơ f ± g, λf và hàm số hf, gi cũng có đạo hàm. Hơn nữa 0 0 0 a) (f ± g) (t0) = f (t0) ± g (t0), 0 0 b) (λf) (t0) = λf (t0), 0 0 0 c) hf, gi (t0) = hf (t0), g(t0)i + hf(t0), g (t0)i. Hệ quả 2.1. Nếu f là một hàm vectơ có đạo hàm sao cho kf(t)k là hằng số, thì tại mọi điểm ta có f 0(t)⊥f(t). Nếu f là một hàm vectơ có đạo hàm tại mọi điểm trên I, thì f 0 cũng là một hàm vectơ. Nếu f 0 cũng có đạo hàm thì ta gọi đạo hàm này là đạo hàm cấp hai của f và ký hiệu là f 00. Tương tự, ta có thể định nghĩa các đạo hàm cấp cao hơn của f. Sử dụng Mệnh đề 2.2 nhiều lần ta cũng thu được kết quả sau
  29. 28 Mệnh đề 2.4. Hàm vectơ f có đạo hàm đến cấp k tại t0 khi và chỉ khi các hàm thành phần fi cũng có đạo hàm cấp k tại điểm đó. Hơn nữa, ta có (k) (k) (k) f (t0) = (f1 (t0), ··· , fn (t0)). Một đường cong được gọi là liên tục, khả vi, khả vi liên tục v.v. nếu hàm vectơ f tương ứng có các tính chất đó. 2.3. Các đối tượng liên quan đến đường cong 2.3.1. Tiếp tuyến và pháp tuyến của đường cong phẳng Cho đường cong phẳng C xác định bởi hệ ( x = x(t), t ∈ [a, b]. (2.2) y = y(t), Với t0 ∈ (a, b), M0 = f(t0) = (x(t0), y(t0)) là một điểm trên đường cong. Giả thiết 0 thêm rằng f (t0) tồn tại và khác không. Với số gia ∆t đủ bé điểm Mt = f(t0 +∆t) = f(t0+∆t)−f(t0) (x(t0 + ∆t), y(t0 + ∆t)) nằm trên đường cong C. Lúc đó ∆t là vectơ chỉ phương của cát tuyến M0Mt. Khi ∆t → 0 ta có Mt → M0 và cát tuyến M0Mt tiến 0 dần về tiếp tuyến của đường cong tại M0. Do đó, nếu f (t0) tồn tại và khác không thì đó cũng chính là vectơ chỉ phương của tiếp tuyến đường cong C tại M0. Đường thẳng đi qua M0 và vuông góc với tiếp tuyến đường cong tại đó được gọi là pháp tuyến của đường cong. Rõ ràng phương trình của tiếp tuyến và pháp tuyến đường cong tại M0 lần lượt được cho bởi ( x = x(t ) + tx0(t ), 0 0 t ∈ R 0 y = y(t0) + ty (t0), và 0 0 x (t0)(x − x(t0)) + y (t0)(y − y(t0)) = 0. Bây giờ giả sử đường cong C được cho bởi phương trình F (x, y) = 0, 2 với F : R → R là hàm hai biến, có các đạo hàm riêng tại điểm M0(x0, y0) ∈ C. M0 được gọi là điểm kỳ dị nếu ∇F (x0, y0) = 0. Trường hợp ngược lại ta nói M0 là điểm chính quy. Lúc đó, theo Định lý ?? tồn tại một hàm ẩn y = y(x) (hay x = x(y)) với y(x0) = y0, F (x, y(x)) = 0 với mọi x và 0 0 0 Fx(x0, y0) + Fy(x0, y0).y (x0) = 0.
  30. 29 Mặt khác, chú ý rằng C có thể được biểu diễn bởi hệ phương trình tham số trong lân cận M0 là ( x = x, y = y(x). 0 Nên (1, y (x0)) cũng là vectơ chỉ phương của tiếp tuyến của C tại M0. Do đó, 0 0 (Fx(x0, y0),Fy(x0, y0)) là vectơ chỉ phương của pháp tuyến đường cong tại M0 mà ta gọi là vectơ pháp của C. Lúc này, phương trình đường tiếp tuyến của C tại M0 là 0 0 (x − x0).Fx(x0, y0) + (y − y0).Fy(x0, y0) = 0 2 2 0 Ví dụ 2.1. Đường tròn đơn vị có phương trình x + y = 1. Ở đây, Fx = 2x và 0 −→ Fy = 2y nên tại mỗi điểm M0(x0, y0) trên đường tròn có thể chọn vectơ pháp là n = (x0, y0). Vậy phương trình tiếp tuyến đường tròn tại M0 là x0(x−x0)+y0(y−y0) = 0, hay x0x + y0y = 1. 2.3.2. Tiếp tuyến và pháp diện của đường cong trong không gian Giả sử C là đường cong khả vi trong không gian được xác định bởi hệ  x = x(t), y = y(t), t ∈ [a, b]. (2.3)  z = z(t), Bằng lập luận tương tự trên đường cong phẳng ta thấy, với mỗi t0 ∈ (a, b), v := 0 0 0 (x (t0), y (t0), z (t0)) (nếu khác vectơ không) là vectơ chỉ phương của đường thẳng tiếp tuyến với C tại M0((x(t0), y(t0), z(t0)). Do đó phương trình tham số của tiếp tuyến này là  0 x = x(t0) + tx (t0), y = y(t ) + ty0(t ), t ∈ R.  0 0  0 z = z(t0) + tz (t0), Mặt phẳng chứa M0, vuông góc với v được gọi là pháp diện của đường cong. Dĩ nhiên, phương trình của pháp diện của C tại M0 là 0 0 0 (x − x(t0))x (t0) + (y − y(t0))y (t0) + (z − z(t0))z (t0) = 0. 2.3.3. Độ cong a. Trường hợp đường cong phẳng Cho C là đường cong phẳng khả vi hai lần, xác định bởi hệ (2.2). Lấy điểm M(x(t), y(t)) ∈ C. Tập hợp {(x(τ), y(τ)) | τ ∈ [a, t]} chính là cung của đường cong
  31. 30 có hai mút là A(x(a), y(a)) và M. Trong chương trình Giải tích I ta đã biết độ dài cung này được tính bởi Z _ t p s(t) := AM = x0(τ)2 + y0(τ)2dτ a và do đó, vi phân cung là p ds(t) = x0(t)2 + y0(t)2dt. Bây giờ cố định một điểm M0(x(t0), y(t0)) với t0 ∈ (a, b). Với mỗi ∆t đủ bé t = t0 + ∆t ∈ (a, b), ta ký hiệu ∆α là độ lớn của góc lập bởi hai vectơ tiếp tuyến của _ đường cong tại M0 và Mt(x(t), y(t)), ∆s là độ dài cung M0Mt. Độ cong của C tại M0 là giới hạn sau nếu nó tồn tại ∆α C(M0) := lim . (2.4) ∆t→0 ∆s 1 Nếu C(M0) 6= 0, đại lượng ρ = được gọi là bán kính độ cong. C(M0) Ví dụ 2.2. Nếu C là một đoạn thẳng thì ∆α luôn luôn bằng không. Do đó, độ cong tại mọi điểm của nó bằng không. Nếu C là đường tròn tâm I, bán kính R > 0 thì với ∆α bé, ∆α cũng chính là góc lập bởi hai bán kính IM0 và IMt. Do đó, độ dài cung ∆s bằng R∆α. Suy ra ∆α 1 C(M0) = lim = . ∆t→0 R∆α R Lúc này, ρ = R hay bán kính độ cong của đường tròn trùng với bán kính của nó. Ở đây ta thấy khái niệm độ cong khá hợp lý ở chỗ là độ cong của đường tròn càng bé khi bán kính của nó càng lớn. Bổ đề 2.1. −→ −→ a) Cho hai vectơ m = (m1, m2), n = (n1, n2) lập thành một góc nhọn. Lúc đó, |m n − m n | sin(−→\m, −→n ) = 1 2 2 1 . (2.5) k−→mk.k−→n k −→ −→ b) Nếu m = (m1, m2, m3) và n = (n1, n2, n3), thì (2.5) được thay bằng s¯ ¯ ¯ ¯ ¯ ¯ ¯m m ¯2 ¯m m ¯2 ¯m m ¯2 ¯ 1 2¯ + ¯ 2 3¯ + ¯ 3 1¯ ¯ n1 n2 ¯ ¯ n2 n3 ¯ ¯ n3 n1 ¯ sin(−→\m, −→n ) = . k−→mk.k−→n k
  32. 31 Định lý 2.5. Giả sử C là đường cong phẳng, khả vi đến cấp hai, xác định bởi hệ (2.2) và M(x(t), y(t)), t ∈ (a, b) là một điểm trên C sao cho (x0(t), y0(t)) 6= (0, 0). Lúc đó |x0y00 − x00y0| C(M) = . (2.6) 2 2 3 (x0 + y0 ) 2 Đặc biệt, nếu C được biểu diễn bởi phương trình y = f(x), x ∈ [a, b] thì với M(x, f(x)) ta có |f 00| C(M) = . 2 3 (1 + f 0 ) 2 Ví dụ 2.3. Đường cycloid: x = a(t − sin t), y = a(1 − cos t) (a > 0) là quỹ tích của một điểm P cố định trên đường tròn bán kính a khi đường tròn này lăn trên trục Ox. Độ cong của đường cycloid tại mỗi điểm M(a(t − sin t), a(1 − cos t)) là | cos t − 1| 1 C(M) = = ¯ ¯. √ 3 ¯ t ¯ 2 a(1 − cos t) 2 4a sin 2 Như vậy, độ cong chỉ xác định tại các điểm M với t 6= 2kπ. Nếu đường cong được cho trong toạ độ cực: r = f(ϕ), thì bằng cách viết x = f(ϕ). cos ϕ, y = f(ϕ). sin ϕ, xem như đó là phương trình tham số của đường cong, từ (2.6) ta có |r2 + 2r02 − rr00| C(M) = . 2 3 (r2 + r0 ) 2 b. Trường hợp đường cong trong không gian Cho đường cong C khả vi đến cấp hai trong không gian, được xác định bởi hệ (2.3). Tương tự như trong mặt phẳng ta có công thức vi phân cung p ds(t) = x0(t)2 + y0(t)2 + z0(t)2dt. _ Với mỗi t0 ∈ (a, b) và ∆t đủ bé, ta vẫn ký hiệu ∆s là độ dài cung M0Mt còn ∆α là số đo góc lập bởi các vectơ tiếp tuyến với đường cong tại M0 và Mt, trong đó M0 và Mt lần lượt là các điểm trên đường cong tương ứng với tham số t0 và t0 + ∆t. Độ cong của C tại M0 cũng được định nghĩa bởi công thức (2.4). Định lý sau cho chúng ta công thức tính độ cong của đường trong không gian. Định lý 2.6. Giả sử M(x(t), y(t), z(t)), t ∈ (a, b) là một điểm trên đường cong C sao cho (x0(t), y0(t), z0(t)) 6= (0, 0, 0). Lúc đó s¯ ¯ ¯ ¯ ¯ ¯ ¯x0 y0 ¯2 ¯y0 z0 ¯2 ¯z0 x0 ¯2 ¯ ¯ + ¯ ¯ + ¯ ¯ ¯x00 y00¯ ¯y00 z00¯ ¯z00 x00¯ C(M) = . 2 2 2 3 (x0 + y0 + z0 ) 2
  33. 32 2.3.4. Hình bao của họ đường cong Cho một họ đường cong phẳng C(λ) phụ thuộc tham số λ. Nếu mọi đường cong C(λ) đều tiếp xúc với một đường cong L và ngược lại, tại mỗi điểm M ∈ L đều tồn tại một đường cong C(λ) của họ tiếp xúc với L tại M, thì L được gọi là hình bao của họ C(λ). Ví dụ 2.4. * Phương trình (x − λ)2 + y2 = R2, trong đó R là số cố định, biểu diễn một họ đường tròn bán kính R có tâm (λ,0) chạy trên trục Ox. Hình bao của họ này là hai đường thẳng y = ±R. * Phương trình x cos λ + y sin λ − 1 = 0 biểu diễn một họ đường thẳng mà khoảng cách từ gốc O đến đường thẳng ấy bằng 1. Hình bao của họ này là đường tròn đơn vị. Định lý 2.7. Cho họ đường cong F (x, y, λ) = 0 phụ thuộc tham số λ. Nếu các đường của họ ấy không chứa điểm kỳ dị, thì phương trình của hình bao L của chúng được xác định bởi hệ ( F (x, y, λ) = 0, 0 Fλ(x, y, λ) = 0. 2.4. Mặt cong 2.4.1. Khái niệm Cho D là một miền trong R2 và một ánh xạ f :D −→ R3; (u, v) ∈D −→ f(u, v) ∈ R3. Với mỗi điểm (u, v) ∈ D, f(u, v) là một vectơ trong không gian với 3 thành phần: f(u, v) = (x(u, v), y(u, v), z(u, v)). Vì vậy việc cho ánh xạ f tương đương với việc cho 3 hàm số thực x(u, v), y(u, v) và z(u, v). Tập hợp S := {f(u, v) | (u, v) ∈ D} được gọi là một mặt cong trong không gian. Thông thường, mặt cong được cho bởi hệ sau mà ta gọi là phương trình tham số của mặt.  x = x(u, v), y = y(u, v), (u, v) ∈ D. (2.7)  z = z(u, v),
  34. 33 Nếu từ phương trình tham số ta giải được u, v theo x, y và từ đó đưa được về phương trình z = g(x, y) (2.8) (hoặc x = g(y, z), y = g(z, x)) thì (2.8) được gọi là dạng hiển của mặt cong S. Nếu ta loại được hai biến u, v từ (2.7) để nhận được phương trình F (x, y, z) = 0, (2.9) trong đó F là hàm ba biến, thì phương trình này được gọi là dạng ẩn của mặt. Mặt được gọi là liên tục nếu các hàm thành phần liên tục; được gọi là trơn nếu các hàm thành phần khả vi liên tục và ma trận Jacobi của f: µ ¶ 0 0 0 xu yu zu Jf := 0 0 0 xv yv zv có hạng bằng 2. Sử dụng Định lý hàm ẩn ta có thể chứng minh được trong trường hợp mặt cong trơn thì tại lân cận của mỗi điểm của mặt, luôn tồn tại một biểu diễn dưới dạng hiển và một biểu diễn dưới dạng ẩn của mặt. Lúc này, các hàm g, F tương ứng cũng khả vi liên tục. Hơn nữa ∇F 6= 0 tại mọi điểm của mặt. Ví dụ 2.5. Mặt cầu tâm O bán kính R > 0 có phương trình tham số  x = R sin θ cos ϕ, y = R sin θ sin ϕ, (θ, ϕ) ∈ [0, π] × [0, 2π).  z = R cos θ, Hoặc dưới dạng phương trình ẩn: x2 + y2 + z2 − R2 = 0. 2.4.2. Tiếp diện và pháp tuyến của mặt cong Giả sử S là một mặt cong trơn trong không gian và M0 là một điểm của mặt. Ta nói một mặt phẳng (P ) đi qua M0 là mặt phẳng tiếp xúc hay tiếp diện của S tại M0 nếu d(M, P ) lim = 0, (2.10) S M M M→M0 0 trong đó d(M, P ) là khoảng cách từ M đến mặt phẳng (P ). Ở đây, (2.10) được hiểu là d(M, P ) ∀² > 0, ∃δ > 0, ∀M ∈ S : 0 < M0M < δ ⇒ < ². M0M Đường thẳng đi qua M0 và vuông góc với tiếp diện của S tại M0 được gọi là pháp tuyến của mặt tại điểm đó.
  35. 34 Định lý 2.8. Giả sử S là mặt cong trơn được xác định bởi hệ (2.7) và M0 ∈ S là điểm ứng với cặp tham số (u0, v0). Lúc đó tiếp diện của S tại M0 là mặt phẳng đi qua M0 và nhận hai vectơ sau làm vectơ chỉ phương 0 0 0 0 fu(u0, v0) = (xu(u0, v0), yu(u0, v0), zu(u0, v0)); 0 0 0 0 fv(u0, v0) = (xv(u0, v0), yv(u0, v0), zv(u0, v0)). Từ định lý này suy ra phương trình tham số của tiếp diện tại M0 của S là  0 0 x = x(u0, v0) + u.xu(u0, v0) + v.xv(u0, v0), y = y(u , v ) + u.y0 (u , v ) + v.y0 (u , v ), (u, v) ∈ R2.  0 0 u 0 0 v 0 0  0 0 z = z(u0, v0) + u.zu(u0, v0) + v.zv(u0, v0), Đặc biệt, nếu mặt được cho dưới dạng hiển (2.8) thì phương trình của tiếp diện tại M0(x0, y0, g(x0, y0)) là 0 0 2 z = (x − x0)gx(x0, y0) + (y − y0)gy(x0, y0) + g(x0, y0), (x, y) ∈ R . −→ 0 0 Lúc đó ta có vectơ pháp của mặt: n = (gx(x0, y0), gy(x0, y0), −1). Định lý 2.9. Giả sử S là mặt cong trơn được xác định bởi phương trình ẩn (2.9) và M0(x0, y0, z0) ∈ S. Lúc đó tiếp diện của S tại M0 là mặt phẳng đi qua M0 và nhận ∇F (M0) làm vectơ pháp. Lúc này phương trình của tiếp diện là 0 0 0 (x − x0)Fx(x0, y0, z0) + (y − y0)Fy(x0, y0, z0) + (z − z0)Fz(x0, y0, z0) = 0 (2.11) và pháp tuyến có phương trình x − x0 y − y0 z − z0 0 = 0 = 0 . Fx(x0, y0, z0) Fy(x0, y0, z0) Fz(x0, y0, z0) Ví dụ 2.6. * Viết phương trình tiếp diện và pháp tuyến của mặt parabol z = x2 + y2 tại M(1, 2, 5). Ta có phương trình mặt ở dạng hiển với g(x, y) = x2 + y2. Tại M ta có −→n = (2, 4, −1). Nên phương trình tiếp diện là z = 2(x − 1) + 4(y − 2) + 5 hay 2x + 4y − z − 5 = 0 và phương trình pháp tuyến là x − 1 y − 2 z − 5 = = − . 2 4 1 * Viết phương trình tiếp diện của mặt Êlip x2+2y2+5z2 = 8 tại điểm M(1, 1, 1). Đây là phương trình dạng ẩn với F (x, y, z) = x2 +2y2 +5z2 −8. Áp dụng (2.11) ta có phương trình tiếp diện là 2(x − 1) + 4(y − 1) + 10(z − 1) = 0 hay x + 2y + 5z − 8 = 0.
  36. 35 2.5. Thực hành tính toán Phần thực hành trong chương này chủ yếu ta sẽ nghiên cứu cách vẽ các đường cong và mặt cong trong mặt phẳng và trong không gian, bao gồm cả các trường hợp đã xét trong Chương 1. Để thực hiện các lệnh trong mục này, nói chung, ta cần khởi động các gói công cụ plots, plottools. 2.5.1. Vẽ đường cong trong mặt phẳng a. Dùng toạ độ Đê-các. Trong toạ độ Đê-các, một đường cong phẳng (C) thường được biểu diễn như là đồ thị của hàm một biến f nào đó: (C): y = f(x), x ∈ [a, b], hoặc được biểu diễn dưới dạng tham số: ( x = x(t), (C): t ∈ [a, b], y = y(t), hoặc dưới dạng một phương trình ẩn (C): F (x, y) = 0. Để vẽ đường cong trong trường hợp thứ nhất ta dùng lệnh [> plot(f(x), x=a b); trong trường hợp thứ hai ta dùng lệnh [> plot([x(t), y(t), t=a b]); và trong trường hợp thứ ba: [> implicitplot(F(x,y)=0, x=a b, y=c d); b. Dùng toạ độ cực. Trong toạ độ cực, một đường cong phẳng thường có hai cách biểu diễn (C): r = f(ϕ); ϕ ∈ [a, b] hoặc ( r = r(t), (C): t ∈ [a, b]. ϕ = ϕ(t), Để vẽ đường cong trong trường hợp thứ nhất ta dùng lệnh [> polarplot(f(phi), phi=a b); và trong trường hợp thứ hai ta dùng lệnh [> polarplot([r(t), phi(t), t=a b]);
  37. 36 2.5.2. Vẽ mặt cong trong không gian a. Dùng toạ độ Đê-các. Trong toạ độ Đê-các, một mặt cong phẳng (S) cũng thường được biểu diễn như là đồ thị của hàm hai biến f nào đó: (S): z = f(x, y), (x, y) ∈ [a, b] × [c, d], hoặc được biểu diễn dưới dạng tham số:  x = x(s, t), (S): y = y(s, t), (s, t) ∈ [a, b] × [c, d],  z = z(s, t), hoặc bởi một phương trình ẩn: (S): F (x, y, z) = 0. Để vẽ mặt cong trong trường hợp thứ nhất ta dùng lệnh [> plot3d(f(x,y), x=a b, y=c d); trong trường hợp thứ hai ta dùng lệnh [> plot3d([x(s,t), y(s,t), z(s,t)], s=a b, t=c d); và trong trường hợp thứ ba ta dùng lệnh [> implicitplot3d(F(x,y,z)=0, x=a b, y=c d, z=e f); b. Dùng toạ độ trụ. Trong toạ độ trụ, một mặt cong phẳng (S) thường được biểu diễn bởi một trong hai cách: (S): r = f(ϕ, z), (ϕ, z) ∈ [a, b] × [c, d], hoặc  r = r(s, t), (S): ϕ = ϕ(s, t), (s, t) ∈ [a, b] × [c, d].  z = z(s, t), Để vẽ mặt cong trong trường hợp thứ nhất ta dùng lệnh [> cylinderplot(f(phi, z), phi=a b, z=c d); và trong trường hợp thứ hai ta dùng lệnh [> cylinderplot([r(s,t), phi(s,t), z(s,t)], s=a b, t=c d); c. Dùng toạ độ cầu. Trong toạ độ cầu, một mặt cong phẳng (S) thường được biểu diễn bởi một trong hai cách: (S): ρ = f(ϕ, θ), (ϕ, θ) ∈ [a, b] × [c, d],
  38. 37 hoặc  ρ = ρ(s, t), (S): ϕ = ϕ(s, t), (s, t) ∈ [a, b] × [c, d].  θ = θ(s, t), Để vẽ mặt cong trong trường hợp thứ nhất ta dùng lệnh [> sphereplot(f(phi, theta), phi=a b, theta=c d); và trong trường hợp thứ hai ta dùng lệnh [> sphereplot([rho(s,t), phi(s,t), theta(s,t)], s=a b, t=c d); Ví dụ: Để vẽ mặt cầu đơn vị ta có thể dùng một trong các lệnh sau (xem Hình 3.1). [> plot3d([sin(s)*cos(t),sin(s)*sin(t),cos(s)],s=0 Pi,t=0 2*Pi); [> cylinderplot(sqrt(1-z∧2), phi=0 2*Pi, z=-1 1); [> cylinderplot([sin(s), phi, cos(s)], phi=0 2*Pi, s=0 Pi); [> sphereplot(1, phi=0 2*Pi, theta=0 Pi); Hình 2.1: Mặt cầu đơn vị 2.5.3. Vận động đồ thị Vận động đồ thị là sự biến thiên của đồ thị theo tham số. Điều đó có nghĩa là ta cho một họ đường cong (Ct) hay mặt cong (St) phụ thuộc vào một tham số t. Sau đó vẽ tất cả các đường/mặt này ứng với các giá trị t khác nhau. Họ đường cong, mặt cong có thể biểu diễn dưới các dạng khác nhau và theo các hệ toạ độ khác nhau. Ở đây, chúng ta chỉ xét họ được viết dưới dạng đơn giản: (Ct): y = f(x, t), x ∈ [a, b], t ∈ [t1, t2] (St): z = f(x, y, t), (x, y) ∈ [a, b] × [c, d], t ∈ [t1, t2].
  39. 38 Để vận động họ (Ct) ta dùng lệnh [> animate(f(x, t), x=a b, t=t1 t2); và vận động họ (St) bằng lệnh [> animate3d(f(x, y, t), x=a b, y=c d, t=t1 t2); Chú ý là khi thực hiện lệnh này ta thấy đồ thị chưa vận động bởi vì máy chỉ vẽ một đường/mặt ứng với một tham số cụ thể nào đó. Nếu đưa con trỏ chuột vào vùng đồ thị và kích trái chuột thì một bảng lệnh điiêù hành sẽ hiện ra ngay dưới thanh công cụ, gồm các ký hiệu play, continuous, stop quen thuộc. Nếu bạn muốn đồ thị vận động liên tục theo các tham số thì nhấn continuous/play, sau đó muốn dừng thì nhấn stop. 2.6. Bài tập 2.1. Viết phương trình tiếp diện và phương trình pháp tuyến của mặt y = 3z2−xz+1 tại các điểm A(x = 0, y = 1, z = 0) và B(x = 2, y = 2, z = 1). 2.2. Viết phương trình tiếp diện và phương trình pháp tuyến của mặt x2+2y2+3z2 = 1 tại các điểm A(1, 0, 0) và B( √1 , √1 , 1 ). 3 6 3 2.3. Viết phương trình tiếp tuyến và phương trình pháp diện của đường cong x = cos t, y = sin t, z = t; t ∈ R π tại các điểm A(1, 0, 0) và B(0, 1, 2 ). 2.4. Viết phương trình tiếp diện và phương trình pháp tuyến của mặt z = 2x2 + 3y2 + 1 tại các điểm A(1, 0, 3) và B(−1, −1, 6). 2.5. Viết phương trình tiếp tuyến và phương trình pháp diện của đường cong sau tại các điểm (−1, −1, 2) và (0, 0, 0): ( z = x2 + y2, x + y + z = 0. 2.6. Viết phương trình của tiếp tuyến và pháp diện của đường cong ( x2 + y2 + z2 − 1 = 0, x − y = 0, √ √ √ 3 3 3 tại điểm M( 3 , 3 , 3 ).
  40. 39 2.7. Viết phương trình của tiếp tuyến và pháp diện của đường cong  x = 2 sin t, y = 5 cos(t2 + t),  z = 1 + 3 sin t, tại điểm M(0, 5, 1). π 2.8. Xác định độ cong của các đường cong sau tại điểm M(0, 1, 2 ): π a) x = cos t, y = sin t, z = t; b) x = cos t, y = sin t, z = 2 . 2.9. Xác định độ cong của đường cong phẳng r = ϕ tại điểm (π, π) (tức là điểm (−π, 0) theo toạ độ Đê-các). 2.10. Hai mặt cong trơn (S) và (T ) được gọi là tiếp xúc nhau tại điểm M0 ∈ (S)∩(T ) nếu tiếp diện tại M0 của (S) và của (T ) trùng nhau. Giả sử (S) và (T ) lần lượt được biểu diễn bởi các phương trình ẩn f(x, y, z) = 0 và g(x, y, z) = 0. Chứng minh rằng (S) và (T ) tiếp xúc nhau tại M0(x0, y0, z0) khi và chỉ khi ( f(x0, y0, z0) = g(x0, y0, z0) = 0, ∇f(x0, y0, z0) = ∇g(x0, y0, z0). 2.11. Chứng minh rằng họ đường cong 2 3 2 6 (Cλ): x − 2λ x + y − 2λy + λ = 0 có hình bao chứa đường thẳng y = 0.