(Luận văn) điều kiện cần cực trị của bài toán biến phân

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC HOÀNG THỊ QUỲNH NHƯ ĐIỀU KIỆN CẦN CỰC TRỊ lu an CỦA BÀI TOÁN BIẾN PHÂN n va gh tn to p ie Chuyên ngành: Toán ứng dụng d oa nl w Mã số: 60.46.01.12 an lu LUẬN VĂN THẠC SĨ TOÁN HỌC nf va z at nh oi lm ul NGƯỜI HƯỚNG DẪN KHOA HỌC PGS TS TẠ DUY PHƯỢNG z m co l gm @ an Lu n va Thái Nguyên - 2017 ac th si Mục lục lu an n va chuẩn 5 7 gh tn to Chương Kiến thức chuẩn bị 1.1 Khơng gian tuyến tính định chuẩn 1.1.1 Khái niệm khơng gian tuyến tính 1.1.2 Khái niệm khơng gian tuyến tính định 1.2 Phép tính vi phân 1.2.1 Dưới vi phân hàm lồi 1.2.2 Đạo hàm Gâteaux 1.2.3 Đạo hàm Fréchet p ie Chương Điều kiện cần cực trị cho tốn tối ưu khơng gian vơ hạn chiều 13 2.1 Định lí Fermat 14 2.1.1 Bài tốn trơn khơng có ràng buộc 14 2.1.2 Bài toán lồi khơng có ràng buộc 16 2.2 Qui tắc nhân tử Lagrange 19 d oa nl w nf va an lu z at nh oi lm ul Chương Điều kiện cần cho toán biến phân 30 3.1 Bài toán biến phân sở 33 3.2 Phương trình Euler 34 3.3 Bổ đề Du Bois-Reymond toán Bolza 42 3.4 Ví dụ Hilbert 49 3.5 Điều kiện Weierstrass 51 3.6 Điều kiện Legendre 53 3.7 Điều kiên Jacobi 55 3.8 Bài toán đẳng chu 59 3.9 Bài toán điều khiển tối ưu nguyên lý cực đại Pontriagin 62 3.9.1 Dẫn tới toán điều khiển tối ưu 62 3.9.2 Bài toán điều khiển tối ưu 63 z m co l gm @ an Lu n va ac th si 66 Tài liệu tham khảo lu an n va p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu n va ac th si Lời nói đầu Điều kiện cần cực trị cho toán tối ưu Fermat phát biểu cách 300 năm Lí thuyết điều kiện cần cực trị khơng gian hữu hạn chiều cho tốn có hạn chế phát triển qua nhiều thời kì nhà toán học Lagrange, Euler, Kuhn, Tucker, lu an n va gh tn to Năm 1696, Johann I Bernoulli phát biểu toán đường đoản thời (brachistochrone): Cho trước hai điểm A B mặt phẳng thẳng đứng Hãy xác định đường AM B để tác động lực trọng trường vật thể M chuyển động từ A đến B thời gian ngắn p ie Vấn đề bắt nguồn từ thực nghiệm G Galilei: Nếu cho hai viên bi giống lăn dây cung cung trịn viên bi lăn cung trịn đến điểm cuối nhanh (mặc dù đường dài hơn, tốc độ lớn hơn) d oa nl w nf va an lu Giả sử y(x) hàm mô tả đường cong chuyển động viên bi hệ trục tọa độ (x, y) với y(x0 ) = 0, y(x) ≤ x ≥ x0 Điểm A có tọa độ A(x0 , 0) điểm B có tọa độ B(x1 , y1 ) cho trước lm ul z at nh oi Giả thiết lực ma sát không đáng kể, theo định luật rơi Galilei, ta p có vận tốc viên bi −2gy(x), g ≈ 9, 8m/s2 gia tốc p rơi tự Quãng đường sau thời gian dt ds = + y 02 (x)dx p p p ds + y (x)dx + y (x) Vì v(t) = nên −2gy(x) = hay dt = p dx dt dt −2gy(x) Zx1 p + y (x) p Thời gian từ A(x0 , 0) đến B(x1 , y1 ) là: T = dx −2gy(x) z m co l gm @ x0 an Lu Bài toán trở thành: Trong số tất quỹ đạo (đường cong) y(x) nối hai điểm A(x0 , 0) đến B(x1 , y1 ) cho trước, tìm đường cong n va ac th si Zx1 p + y (x) p làm cực tiểu phiếm hàm dx −2gy(x) x0 Đây tốn cực trị có ràng buộc: Zx1 p + y (x) p T = dx → inf −2gy(x) (1) y(x0 ) = 0, y(x1 ) = y1 (2) x0 Bài tốn (1)–(2) tốn tối ưu khơng gian vô hạn chiều (không gian tất đường cong trơn nối hai điểm cho trước) Sau Bernoulli, toán tổng quát thành toán biến phân: lu Zt1 an L(t, x(t), x(t))dt ˙ → inf va J(x(.)) = n t0 to gh tn (x(t0 ), x(t1 )) ∈ Γ, p ie [t0 , t1 ] ⊂ R cho trước, Γ ⊂ Rn × Rn , L hàm liên tục miền R × Rn × Rn d oa nl w Mục đích luận văn trình bày khái niệm nghiệm yếu nghiệm mạnh địa phương tồn cục tốn biến phân, điều kiện cực trị cấp cấp hai toán biến phân an lu nf va Tác giả xin bày tỏ lòng biết ơn sâu sắc tới PGS.TS Tạ Duy Phượng tận tình hướng dẫn giúp đỡ em suốt tình học tập nghiên cứu Em xin bày tỏ lòng biết ơn chân thành tới Giáo sư, Phó giáo sư, Tiến sĩ, quý thầy cô giáo giảng dạy Đại học Khoa học, Đại học Thái Nguyên Viện Toán học, Viện Hàn lâm Khoa học Công nghệ Việt Nam, mang đến cho em nhiều kiến thức bổ ích nghiên cứu khoa học Đồng thời, xin gửi lời cảm ơn tới gia đình bạn đồng mơn giúp đỡ động viên thời gian học tập q trình hồn thành luận văn z at nh oi lm ul z l gm @ Thái Nguyên, tháng 05 năm 2017 co m Tác giả an Lu Hoàng Thị Quỳnh Như n va ac th si Chương Kiến thức chuẩn bị 1.1 Khơng gian tuyến tính định chuẩn 1.1.1 Khái niệm khơng gian tuyến tính lu an Định nghĩa 1.1.1 Tập X 6= ∅ gồm đối tượng gọi khơng gian tuyến tính trường số thực R, đó: n va gh tn to (I) Có qui tắc cho ứng với hai phần tử x, y thuộc X phần tử z thuộc X gọi “tổng” x y, ký hiệu z = x + y; p ie (II) Có qui tắc cho ứng với phần tử α ∈ R phần tử x ∈ X phần tử p thuộc X gọi tích α với x, ký hiệu p = αx oa nl w (III) Các qui tắc cho (I) (II) phải thỏa mãn tám tiên đề sau đây: d (1) ∀x, y ∈ X : x + y = y + x (tính giao hốn); lu nf va an (2) ∀x, y, z ∈ X : (x + y) + z = x + (y + z) (tính kết hợp); (3) ∃θ (phần tử 0) cho ∀x ∈ X : θ + x = x + θ = x; lm ul (4) ∀x ∈ X : ∃x0 (phần tử đối) cho: x + x0 = x0 + x = θ; z at nh oi (5) ∀x ∈ X : 1x = x; (1 ∈ R); (6) ∀α ∈ R, ∀x, y ∈ X : α(x + y) = αx + αy; z (7) ∀α, β ∈ R, ∀x ∈ X : (αβ)x = α(βx); gm @ (8) ∀α, β ∈ R, ∀x ∈ X : (α + β)x = αx + βx l m co Ví dụ 1.1.2 Khơng gian hàm liên tục từ [a, b] vào R, kí hiệu C[a, b] khơng gian tuyến tính Tương tự, khơng gian hàm khả vi liên tục [a, b], kí hiệu C1 [a, b] không gian tuyến tính an Lu n va ac th si Thật vậy, với định nghĩa (i) z = x + y ⇔ z(t) = x(t) + y(t) ∀t ∈ [a, b]; (ii) z = λ.x ⇔ z(t) = λx(t) ∀t ∈ [a, b] tám tiên đề thỏa mãn Vậy không gian hàm liên tục hàm khả vi liên tục [a, b] không gian tuyến tính 1.1.2 Khái niệm khơng gian tuyến tính định chuẩn Khơng gian tuyến tính X R gọi khơng gian (tuyến tính) định chuẩn, có qui tắc cho ứng với phần tử x ∈ X số thực không âm gọi chuẩn (hoặc độ dài) x, ký hiệu kxk, thỏa mãn tính chất sau đây: lu an (i) ∀x ∈ X : kxk ≥ kxk = ⇔ x = (tính khơng âm); n va (ii) ∀x ∈ X, ∀λ ∈ R : kλxk = |λ|kxk (tính đồng nhất); gh tn to (iii) ∀x, y ∈ X : kx + yk ≤ kxk + kyk (Bất đẳng thức tam giác) p ie Ví dụ 1.1.3 Khơng gian C[a, b] hàm liên tục ϕ : R → R không gian C1 [a, b] khơng gian tuyến tính định chuẩn với chuẩn tương ứng kxkC[a,b] = max |x(t)| oa nl w a≤t≤b d kxkC1n [a,b] = max{k x kC[a,b] , k x˙ kC[a,b] a≤t≤b a≤t≤b nf va an lu = max{max |x(t)|, max |x(t)|} ˙ Khi ba tiên đề chuẩn thỏa mãn lm ul  x1 (t)     n Ví dụ 1.1.4 Giả sử: x(·) : [a, b] → R , tức là: x(t) =   ∈   xn (t)   x˙ (t)     Rn ∀t ∈ [a, b], ta định nghĩa x(t) ˙ =    x˙ n (t) Kí hiệu: C n [a, b] C1n [a, b] khơng gian tuyến tính [a, b] với định nghĩa thơng thường phép tốn cộng véctơ nhân số  z at nh oi z m co l gm @ an Lu n va ac th si với véctơ Ta định nghĩa: kxkC n [a,b] = max {kxi kC[a,b] } 1≤i≤n kxkC1n [a,b] = max {kxi kC1 [a,b] } 1≤i≤n Khi C n [a, b] C1n [a, b] khơng gian tuyến tính định chuẩn 1.2 Phép tính vi phân 1.2.1 Dưới vi phân hàm lồi Cho f hàm lồi thật (hàm lồi thường) X Tập lu ∂f (x) := {x∗ ∈ X ∗ | f (z) − f (x) ≥ hx∗ , z − xi ∀z ∈ X} (1.1) an n va gọi vi phân f x p ie gh tn to Ví dụ 1.2.1 (Dưới vi phân hàm chỉ) Với x ∈ A ∂δ(x|A) khác rỗng chứa Từ định nghĩa ta suy (1.2) N (x|A) := {x∗ ∈ X ∗ |hx∗ , z − xi ≤ ∀z ∈ A} (1.3) d oa nl w ∂δ(x | A) = N (x | A), nf va an lu nón pháp tuyến (normal cone) tập A điểm x lm ul 1.2.2 Đạo hàm Gâteaux z at nh oi Cho X Y hai không gian tơ pơ tuyến tính, V lân cận x ∈ X F : X → Y Nếu z δF (x, h) := lim t−1 (F (x + th) − F (x)) gm @ t→0 m co l tồn với h ∈ X ánh xạ h → δF (x, h) gọi biến phân bậc F x Nếu tồn tốn tử tuyến tính liên tục Λ : X → Y cho an Lu Λh = δF (x, h) ∀h ∈ X n va ac th si Λ gọi đạo hàm Gâteaux, ký hiệu FG0 (x) Khi ta nói: F khả vi Gâteaux x Điều xảy tồn tốn tử tuyến tính liên tục Λ : X → Y cho F (x + th) = F (x) + tΛh + o(t) ∀h ∈ X Ví dụ 1.2.2 Cho r ϕ tọa độ cực x ∈ R2 f (x) = r cos 3ϕ Ta có δf (0.h) = f (h) Vì δf (0.h) khơng tuyến tính nên f khơng khả vi Gâteaux ∈ R2 Trong Giải tích lồi, vi phân đóng vai trị đạo hàm Nếu hàm lồi khả vi Gâteaux điểm vi phân điểm có phần tử đạo hàm Gâteaux lu an 1.2.3 Đạo hàm Fréchet n va gh tn to Nếu X Y không gian Banach, F : X → Y khả vi Fréchet x tồn toán tử tuyến tính liên tục Λ : X → Y cho F (x + h) = F (x) + Λh + r(h) kr(h)kY = khkX →0 khkX lim với p ie d oa nl w Khi Γ đạo hàm Fréchet, ký hiệu FF0 (x) hay F (x) Ánh xạ F gọi qui x khả vi Fréchet x ImF (x) = Y Kí hiệu L(X, Y ) khơng gian tốn tử tuyến tính liên tục từ X Y, trang bị chuẩn nf va an lu kΛk = sup kΛxkY lm ul kxkX =1 z at nh oi Nếu F : X → Y khả vi Fréchet điểm tập mở V ánh xạ x → F (x) liên tục V (hay x0 ∈ V ) theo tơ pơ L(X, Y ) ta nói F khả vi liên tục V (hay x0 ,) hay F thuộc vào lớp C1 Nếu F phiếm hàm F (x) = x gọi điểm dừng z l gm @ m co Ví dụ 1.2.3 (Đạo hàm Fréchet ánh xạ afin) Một ánh xạ A : X → Y từ khơng gian tuyến tính X vào khơng gian tuyến tính Y có dạng A(x) = Λx + a, an Lu n va ac th si với a ∈ X Λ ánh xạ tuyến tính từ X vào Y , gọi ánh xạ afin Nếu X Y không gian Banach Λ liên tục A khả vi Fréchet khắp nơi A0F (x) = Λ Mệnh đề 1.2.4 a) Nếu F khả vi Fréchet x F liên tục khả vi Gâteaux FG0 (x) = FF0 (x) b) Nếu F khả vi Gâteaux x F biến phân bậc tồn δF (x, h) = FG0 (x)h Chứng minh: Xem [4], p 35  Ví dụ 1.2.5 Hàm lu  1 x = (x )2 x 6= 0, 2 f (x1 , x2 ) = 0 trường hợp lại an n va ie gh tn to khả vi Gâteaux (0, 0) ∈ R2 khơng liên tục nên theo Mệnh đề 1.2.4 khơng thể khả vi Fréchet Ta có p Định lý 1.2.6 (Định lí giá trị trung bình) Cho X Y không gian tô pô tuyến tính, U tập mở X, ánh xạ F : U → Y khả vi Gâteaux điểm đoạn nối [x, x + h] ⊂ U Khi ta có: d oa nl w nf va an lu a) Nếu ánh xạ z → FG0 (z) ánh xạ liên tục [x, x + h] vào Y Z F (x + h) − F (x) = FG0 (x + th)h dt z at nh oi lm ul b) Nếu X Y khơng gian Banach z kF (x + h) − F (x)k ≤ sup kFG0 (x + th)k · khk gm @ 0≤t≤1 co l với Λ ∈ L(X, Y ) : 0≤t≤1 m kF (x + h) − F (x) − Λhk ≤ sup kFG0 (x + th) − Λk · khk an Lu n va ac th si 52 Theo điều kiện Weierstrass (3.22) L(t, x∗ (t), ξ) − L(t, x∗ (t), x˙ ∗ (t), ξ))Lx˙ (t, x∗ (t), x˙ ∗ (t), ξ)) ≥ (3.23) Có nghĩa là: đồ thị hàm L(t, x∗ (t), ) : R → R nằm phía đường tiếp tuyến ξ = x˙ ∗ (t) Người ta dễ đặt câu hỏi: Để đạt Rt giá trị cực tiểu tích phân t01 L(t.x∗ (t), x˙ ∗ (t))dt, không chọn x˙ ∗ (t) điểm cực tiểu L(t, x∗ (t), ) với t ∈ [t0 , t1 ]? Điều có khơng, khơng thể thỏa mãn điều kiện x∗ (t0 ) = x0 x∗ (t1 ) = x1 Sở dĩ phải đòi hỏi (3.5.1) với ξ ∈ R vì: ˙ có khơng gian C[0, 1], kx(·)kC n nhỏ |x(t)| thể lớn tùy ý lu an n va gh tn to Dễ thấy rằng: điều kiện Weierstrass (3.22) hiển nhiên thỏa mãn hàm L tựa qui, tức L(t.x, ·) lồi (theo biến thứ ba) với t x p ie Quay lại ví dụ 3.2.4, ta có x∗ (t) = t, x˙ ∗ (t) = w ε(t, x∗ (t), x˙ ∗ (t), ξ) = ξ − − (ξ − 1)3 d oa nl = (ξ − 1)2 (ξ + 2) < ξ < −2 nf va an lu Như vậy, nghiệm x∗ (t) = t không thỏa mãn điều kiện Weierstrass, nên khơng thể nghiệm tối ưu địa phương mạnh toán (3.7) z at nh oi lm ul Để chứng minh định lí K Weierstrass sử dụng loại hàm nhiễu đặc biệt:   / [τ, τ + ]  0 cho t ∈ h(t, λ, ) := λ cho t = τ + λ, ξ ∈ R, (3.24)    afin đoạn [τ, τ + λ] [τ + λ, τ + ], z l gm @ m co gọi biến phân Weierstrass Hàm có phần đồ thị (khác 0) hình tam giác, mà độ dài đáy  tiến tới tam giác tiến dần đến triệt tiêu (để hàm bị nhiễu x∗ (·) + h(., λ, ) tiến dần đến x∗ (·) theo chuẩn k.kC n ), độ nghiêng x˙ = ξ cạnh đoạn [τ, τ + λ] an Lu n va ac th si 53 không thay đổi Như vậy, tam giác nhiễu giống kim Vì vậy, người ta cịn gọi hàm nhiễu biến phân hình kim Với τ thuộc (t0 , t1 )  > 0, khảo sát Z t1 ˙ λ, ))dt ϕ λ := L(t, x∗ (t) + h(t, λ, ), x˙ ∗ (t) + h(t, t0 ta có lim ϕ0 (+0) = L(τ, x∗ (τ ), x˙ ∗ (τ ), ξ)−L(τ, x∗ (τ ), x˙ ∗ (τ ))−ξLx˙ (τ, x∗ (τ ), x˙ ∗ (τ )) →0 Do ϕ (+0) ≥ (với  > 0) nên đẳng thức vừa cho ta kết luận Định lí 3.5.1 (xem [1], trang 54) lu 3.6 Điều kiện Legendre an n va Khác với hai điều kiện xét, điều kiện Legendre điều kiện cần bậc hai Cụ thể, ta có: p ie gh tn to Định lý 3.6.1 (Điều kiện Legendre cho nghiệm cực tiểu yếu) Giả thiết L khả vi liên tục ba lần miền U ⊂ R3 , với t ∈ [t0 , t1 ], nl w (t, x∗ (t), x˙ ∗ (t)) ∈ U d oa hàm x∗ (·) khả vi liên tục hai lần [t0 , t1 ] Nếu x∗ (·) nghiệm tối ưu địa phương yếu toán (3.3) an lu A(t) := Lx˙ x˙ (t, x∗ (t), x˙ ∗ (t)) ≥ ∀t ∈ [t0 , t1 ] nf va (3.25) z at nh oi lm ul (3.25) gọi điều kiện Legendre Nó nói lên quĩ đạo tối ưu qua điểm mà L không lõm theo biến thứ ba z Bây ta xét toán gần giống với (3.7) Ví dụ 3.2.4, thay đổi chút điều kiện biên @ F(x(·)) = x˙ (t)dt → inf; x(0) = 0, x(1) = γ, (3.26) an Lu m Z co l gm Ví dụ 3.6.2 (Nghiệm Phương trình Euler khơng thỏa mãn điều kiện Legendre) n va ac th si 54 Tương tự Ví dụ 3.2.4, Phương trình Euler điều kiện biên 6x(t)ă x(t) = 0, x(0) = 0, x(1) = γ, cho nghiệm x∗ (t) = γt Khi L = x˙ dẫn đến A(t) = Lx˙ x˙ (t, x∗ (t), x˙ ∗ (t)) = 6x˙ ∗ (t) = 6γ Rõ ràng, điều kiện Legendre thỏa mãn γ ≥ Do đó, theo Định lí 3.2.2 3.6.1, tốn (3.26) khơng có nghiệm tối ưu địa phương yếu nên γ < Chứng minh Định lí (3.6.1): Kí hiệu lu an B(t) = Lxx |x∗ (t) , C(t) = Lxx ˙ |x∗ (t) (3.27) va n Ta có biến phân bậc hai F(·) x∗ (·) (3.28) p ie gh tn to δ F(x∗ (·), x(·)) Z t1 = (A(t)x˙ (t) + B(t)x2 (t) + 2C(t)x(t)x(t))dt ˙ t0 Z t1 d (A(t)x˙ (t) + (B(t) − C(t))x2 (t)dt = dt t0 w d oa nl (3.29) nf va an lu cho x(·) ∈ L0 , L0 = {x(·) ∈ C1 [t0 , t1 ] | x(t0 ) = x(t1 ) = 0} Vì x∗ (·) nghiệm tối ưu địa phương tốn (3.3) nên theo Định lí 3.2 (điều kiện cần bậc cao) ta có lm ul δ F(x∗ (·), x(·)) ≥ 0, ∀x(·) ∈ L0 (3.30) z at nh oi Giả sử điều kiện (3.25) không thỏa mãn giả thiết khả vi L x∗ (·) kéo theo tồn điểm τ ∈ (t0 , t1 ) cho z gm @ A(τ ) < m co l Để điều dẫn đến mâu thuẫn với (3.30), ta sử dụng phép biến an Lu n va ac th si 55 phân sau: √  λ− h(t, γ, τ ) = 0 |t−τ | √ λ cho |t − τ | ≤ λ2 , cho |t − τ | ≥ λ2 Biến phân có nét đặc biệt √ ≤ h(t, λ, τ ) ≤ λ , √ ˙ λ, τ )| ≤ λ √1 = , |h(t, λ, τ )h(t, λ Z h˙ (t, λ, τ )dt = |t−τ |≤λ/2 lu an Vì vậy, thay vào (3.28) ta có va λ2 δ F(x∗ (·), h(., λ, τ )) ≤ max A(t) + max |B(t)| + max |C(t)|λ |t−τ |≤λ/2 |t−τ |≤λ/2 |t−τ |≤λ/2 n tn to ie gh → A(τ ) < λ → p Như vậy, λ đủ nhỏ δ F(x∗ (·), h(., λ, τ )) < Có điều h(., λ, τ ) khơng khả vi liên tục nên khơng thuộc L0 Tuy nhiên, ˜ λ, τ ) ∈ L0 thỏa mãn cần làm trơn h(., λ, τ ) ta thu h(., ˜ λ, τ )) < Điều mâu thuẫn với (3.30) Do đó, x∗ (·) δ F(x∗ (·), h(., phải thỏa mãn điều kiện (3.25)  d oa nl w an lu Điều kiên Jacobi nf va 3.7 z at nh oi lm ul Trong phần trên, ta thường dùng cách hàm biến phân khác khoảng nhỏ để gây nhiễu cục Vì vậy, điều kiện cần thu đưa đảm bảo mang tính chất cục Nghĩa là: khoảng đủ nhỏ nghiệm thỏa mãn điều kiện cần đủ tốt Song, sống hàng ngày, nhiều tốt diện hẹp gộp lại với chưa cho thành phẩm tối ưu diện rộng Ví dụ sau minh họa điều z co l gm @ m Với hai điểm A B cho trước bề mặt cầu, ta phải tìm đường ngắn chạy mặt cầu nối hai điểm Khi hai điểm A B đủ gần đường nối ngắn phải đoạn an Lu n va ac th si 56 nằm giữa, thuộc giao mặt cầu mặt phẳng qua tâm cầu hai điểm Rõ ràng, điều kiện cần cực tốt, đến mức A B đủ gần cho nghiệm nhất, nghiệm tối ưu Nhưng A B nằm xa nhau, gần đối xứng qua tâm, ta bắt đầu lúng túng với khái niệm “nằm giữa”, lẽ có hai đoạn cong khác mặt phẳng qua tâm “nằm giữa” hai điểm Trừ trường hợp A B thực đối xứng qua tâm, lại hai “ứng cử viên” mà điều kiện cần đưa đường nối ngắn lu Làm để khắc phục tình trạng trên? Ta thêm địi hỏi sau vào điều kiện cần: Ngoài điểm B ra, đoạn nối cần tìm khơng có điểm đối xứng với A qua tâm cầu! Điều kiện giúp loại bớt trường hợp xấu an n va p ie gh tn to Với kinh nghiệm này, ta hiểu phần xuất điều kiện Jacobi, điều kiện cần mang tính chất tồn cục quan trọng cho nghiệm cực tiểu địa phương Hãy để ý rằng: Biến phân dùng để chứng minh điều kiện khác gần khắp khoảng khảo sát [t0 , t1 ] oa nl w d Đặt K(x(·)) = δ F(x∗ (·), x(·)), (3.29)-(3.30) cho ta Z t1 d K(x(·)) = (A(t)x˙ (t) + B(t) − C(t)x2 (t))dt ≥ ∀x(·) ∈ L0 dt t0 (3.31) nf va an lu lm ul z at nh oi Nếu tồn biến phân x(·) ∈ L0 không tầm thường (;tức x(·) đồng 0), cho giá trị cực tiểu K(·), nghĩa gây ảnh hưởng tới F(x∗ (·)), x(·) phải thỏa mãn phương trình Euler sau:   d d − (A(t)x) ˙ + B(t) − C(t) x = (3.32) dt dt z co l gm @ m Người ta gọi Phương trình Jacobi Với giả thiết tính khả vi Định lí 3.6.1, (3.32) có an Lu n va ac th si 57 dạng ˙ x˙ + (B(t) C(t))x A(t)ă x A(t) = Nếu thêm vào điều kiện Legendre ngặt A(t) = Lx˙ x˙ (t, x∗ (t), x˙ ∗ (t)) > ∀t ∈ [t0 , t1 ] (3.33) thỏa mãn ta đưa (3.32) dạng sau xă P (t)x Q(t)x = 0, ú ˙ ˙ P (t) = −A−1 (t)A(t), Q(t) = −A−1 (t)(B(t) − C(t)) (3.34) Gọi Φ(., t0 ) nghiệm phương trình (3.34) với điều kiện ban đầu lu an Φ(t0 , t0 ) = 0, ˙ , t0 ) = Φ(t n va gh tn to (ta biết tốn Cauchy có nghiệm) Một điểm τ > t0 gọi điểm liên hợp t0 Φ(τ, t0 ) = p ie Định lý 3.7.1 (Điều kiện cần Jacobi) Giả thiết rằng: L khả vi liên tục ba lần miền U ⊂ R3 , với t ∈ [t0 , t1 ], d oa nl w (t, x∗ (t), x˙ ∗ (t)) ∈ U nf va an lu hàm L khả vi liên tục hai lần [t0 , t1 ] thỏa mãn điều kiện Legendre ngặt (3.33) Nếu x∗ (·) nghiệm tối ưu địa phương yếu tốn (3.3) khơng thể tồn điểm liên hợp t0 khoảng (t0 , t1 ) z at nh oi lm ul z Có thể so sánh vai trò điểm liên hợp t0 điểm đối xứng qua A tâm cầu (trong ví dụ kể trên) Cho nên, điều kiện “khơng thể tồn điểm liên hợp t0 khoảng (t0 , t1 )” giống đòi hỏi “đoạn nối cần tìm khơng có điểm đối xứng với A qua tâm cầu (ngoài điểm B ra)” co l gm @ m Chứng minh (tóm tắt): Giả sử ngược lại: tồn τ ∈ (t0 , t1 ) với an Lu Φ(τ, t0 ) = (3.35) n va ac th si 58 Do phương trình (3.34) có nghiệm nên ˙ Φ(τ, t0 ) 6= (3.36) Chọn  Φ(t, t ) h(t) = 0 cho t ∈ [t0 , τ ], cho t ∈ [τ, t1 ] lu Thực tích phân phần sử dụng (3.32) (3.35), ta nhận Z τ ˙ K(h(·)) = (A(t)Φ˙ (t, t0 ) + (B(t) − C(t))Φ (t, t0 ))dt t0  Z τ d ˙ t0 )] + (B(t) − C(t))Φ(t, ˙ t0 ) Φ(t, t0 )dt = − [A(t)Φ(t, dt t0 an va n = tn to p ie gh Để mâu thuẫn, dùng biến phân sau hàm h(·) :    h(t) cho t0 ≤ t ≤ τ − λ, h(t, λ, τ, ) = afin [τ − λ, τ + ],    cho t ≥ τ + , d oa nl w lu nf va an với λ  thông số nguyên dương đủ nhỏ Với ϕ(λ) = K(h(., λ, τ, )), z at nh oi lm ul ϕ0 (+0) = −A(τ )Φ˙ (τ, t0 ) z (xem chi tiết [4], tiếng Anh: P 133) Do điều kiện Legendre ngặt (3.33) (3.36), ta có ϕ0 (+0) < 0, nghĩa tồn λ0 > 0 > đủ nhỏ để K(h(., λ0 , τ, 0 )) < Bây ta cần “làm trơn” h(., λ0 , τ, 0 ) để nhận h(·) ∈ L0 với K(h(·)) < Điều mâu thuẫn với (3.31)  m co l gm @ an Lu n va ac th si 59 Ví dụ 3.7.2 T Z (x˙ (t) − x2 (t))dt → inf; F(x(·)) = x(0) = x(T ) = 0 Phương trỡnh Euler cho bi toỏn ny l xă + x = Trước hết, xét trường hợp T = 2π Tất hàm có dạng xc (t) = c sin t thỏa mãn phương trình điều kiện biên x(0) = x(2π) = Tuy nhiên, không hàm số nghiệm tối ưu địa phương yếu, chúng, phương trình (3.34) tng ng u l xă + x = 0, x(0) = 0, x(0) ˙ =1 lu cho nghiệm Φ(t, 0) = sin t Do Φ(π, 0) = nên π điểm liên hợp 0, lại nằm khoảng (0, 2π) Nghĩa xc (·) không thỏa mãn điều kiện cần Jacobi an n va p ie gh tn to Trong trường hợp T < π, hàm x∗ (t) ≡ thỏa mãn phương trình Euler điều kiện biên x(0) = x(T ) = Nó thỏa mãn điều kiện cần Jacobi, điểm liên hợp lớn T Trên thực tế, tối ưu địa phương yếu, Z T Z T F(x(·)) = (x˙ (t) − x2 (t))dt = (x(t) ˙ − x(t) cot t)2 dt d oa nl w 0 lu 3.8 nf va an T < π x(·) ∈ L0 Bài toán đẳng chu lm ul z at nh oi Phép tính biến phân khơng dừng lại tốn sở đơn giản (3.2) hay (3.3) Trong phần này, ta xét lớp rộng hơn, gọi toán đẳng chu: Z t1 F(x(·)) = f0 (t, x, x)dt ˙ → inf; t0 Z t1 fj (t, x, x)dt ˙ = αj , j = 1, , m, (3.37) z m hi : Rn → Rsi , i = 0, an Lu fj : R × Rn → R, j = 1, , m, co h0 (x(t0 )) = h1 (x(t1 )) = 0, l gm @ t0 (3.38) n va ac th si 60 “Đẳng chu” (isoperimetric) có nghĩa “chu vi không thay đổi” Nhưng chu vi nào? Thực toán đẳng chu cổ điển, điều kiện “chu vi khơng thay đổi” thể dạng: tích phân biểu diễn chu vi số cho trước Về sau, người ta dùng tên “Đẳng chu” cho tốn mở rộng hay có dạng tương tự Sử dụng qui tắc nhân tử Lagrange (Định lí 3.4), ta chứng minh định lí sau: Định lý 3.8.1 Giả thiết hàm xuất toán khả vi liên tục Nếu x∗ (·) nghiệm tối ưu địa phương yếu (3.37) tồn nhân tử Lagrange λj ∈ R, ≤ j ≤ m, li ∈ Rsi , i = 0, 1, cho chúng không triệt tiêu hàm Lagrange lu an L(t, x, x) ˙ = m X λj fj (t, x, x) ˙ n va j=0 p ie gh tn to thỏa mãn phương trình Euler   d − Lx˙ + Lx |x∗ (t) = dt nl w với điều kiện biên T d oa Lx|x ˙ ∗ (t0 ) = h (x∗ (t0 ))l0 , T Lx|x ˙ ∗ (t1 ) = h (x∗ (t1 )) lu nf va an Ví dụ 3.8.2 Xác định hình dáng sợi dây đồng treo lên hai đầu lm ul z at nh oi Bài tốn có ý nghĩa thực tế: Khi căng dây điện, người ta cần biết dáng điệu để xác định độ võng Điều thực cần thiết cho việc đảm bảo an toàn đường điện z m co l gm @ Ta biểu diễn hình dáng sợi dây hệ trục tọa độ (t, x) Để đơn giản, xét trường hợp hai đầu treo cao nhau, tức hai điểm (−t0 , 0) (t0 , 0) Ta có s  2 p ∆x ∆s = ∆2 x + ∆2 t = ∆t + ∆t an Lu n va ac th si 61 Suy ds = p X Z s= + x˙ (t)dt Độ dài sợi t0 ds = p + x˙ dt Z t0 gρx p + x˙ dt (g ≈ 9.8m/s2 ) −t0 −t0 Gọi l độ dài cho trước sợi dây, ta có tốn biến phân sau đây: Z t0 p x + x˙ dt → inf; −t Z t00 p + x˙ dt = l, x(−t0 ) = x(t0 ) = −t0 lu Bây ta áp dụng Định lí 3.8.1 Hàm Lagrange tương ứng p p p 2 L = λ0 x + x˙ + λ1 + x˙ = + x˙ (λ0 x + λ1 ) an n va Vì L khơng phụ thuộc vào t, theo Hệ 3.2.3, phương trình Euler có tích phân lượng p ie gh tn to d oa nl w H(t) = Lx˙ (x∗ (t), x˙ ∗ (t))x˙ ∗ (t) − L(x∗ (t), x˙ ∗ (t)) p = x˙ (1 + x˙ )−1/2 (λ0 x + λ1 ) − + x˙ (λ0 x + λ1 )   p x˙ = (λ0 x + λ1 ) √ − + x˙ 2 + x˙ (−1) = (λ0 x + λ1 ) √ = const, + x˙ nf va an lu tức lm ul λ0 x + λ1 = C + x˙ z at nh oi Nếu λ0 = p λ1 + x˙ = C  2 λ1 ⇒1 + x˙ = C  2 λ1 ⇒x˙ = − = const C p z m co l gm @ an Lu ⇒x˙ = const n va ac th si 62 Điều kiện biên kéo theo x ≡ l = 2t0 Nếu l > 2t0 ta đặt λ0 = λ1 = λ Sau giải phương trình ta có p x + λ = C + x˙   t +D ⇒ x + λ = C cosh C Do điều kiện biên đối xứng nên D = Độ dài sợi dây x = C cosh(t/C), t0 ≤ t ≤ t0 , 2C sinh(t0 /C) Gọi C0 nghiệm phương trình C sinh(t0 /C) = ¯ = C0 cosh(t0 /C0 ), ta nhận dạng sợi dây l/2 λ ¯ x(t) = C0 cosh(t/C0 ) + λ lu an 3.9 n va Bài toán điều khiển tối ưu nguyên lý cực đại Pontriagin Dẫn tới toán điều khiển tối ưu tn to 3.9.1 p ie gh Một bước phát triển toán cực trị Bài toán điều khiển tối ưu Phần trình bày tốn điều khiển tối ưu phát triển biến phân cổ điển nl w d oa Xét tập hợp hàm liên tục tuyệt đối nf va an lu  AC([0, T ], Rn ) := x(.) : [0, T ] → Rn |x(.)liên tục tuyệt đối cho hàm Z J(x) = z at nh oi lm ul Bài tốn đặt là: Tìm hàm x(t) ∈ AC([0, T ], Rn ), với điều kiện cho trước x(0) = x0 , x(T ) = x1 , T L(t, x(t), x(t))dt, ˙ z @ l gm L : R × Rn × Rn → R hàm liên tục, đạt giá trị nhỏ m co Nếu ta đặt x(t) ˙ = u(t), t ∈ [0, T ] coi u(t) ( hàm đo x(t) liên tục tuyệt đối) biến điều khiển, ta có tốn điều khiển tối ưu: Trong số tất hàm đo u(t) [0, T ], an Lu n va ac th si 63 tìm quỹ đạo tương ứng x(t) hệ  x(t) ˙ = u(t), t ∈ [0, T ] x(0) = x0 , x(T ) = x1 cho Z J(x, u) = T L(t, x(t), u(t))dt, lu đạt cực tiểu Như vậy, tốn biến phân cổ điển nói đưa tốn tìm cực tiểu hàm mục tiêu J(u) tập điều khiển chấp nhận u(.) đưa trạng thái x0 x1 phương trình điều khiển x(t) ˙ = u(t) Ví dụ cho thấy cần thiết nghiên cứu toán điều khiển tối ưu mơ tả hệ phương trình vi phân tổng quát an n va Bài toán điều khiển tối ưu tn to 3.9.2 p ie gh Trong mục này, ta xét toán điều khiển tối ưu với khoảng thời gian [t0 , t1 ] cố định trạng thái kết thúc x(t1 ) hoàn toàn tự Cụ thể là: Z t1 L(t, x(t), u(t))dt → inf, (3.39) F(x, u) := oa nl w t0 d x(t) ˙ = f (t, x(t), u(t)) ∈ Rn , u(t) ∈ U ⊂ Rr , (3.41) x(t0 ) = x0 (3.42) nf va an lu (3.40) lm ul z at nh oi L f liên tục theo biến khả vi liên tục theo x So với mơ hình tổng qt xem [1] trang 70 ta có ϕ0 (t, x) = x − x0 , ϕ1 (t, x) ≡ z @ điều khiển tối ưu liên tục khúc m co l gm Để chứng minh Nguyên lí cực đại Pontryagin cách đơn giản, ta sử dụng giả thiết quan trọng là: an Lu Cùng với giả thiết x(t1 ) tự do, ta sử dụng biến phân hình kim, n va ac th si 64 làm chứng minh Điều kiện Weierstrass cho nghiệm cực tiểu địa phương mạnh tốn biến phân (Định lí 3.5.1) Định lý 3.9.1 (Ngun lí cực đại Pontryagin) Cho x∗ (·), u∗ (·) q trình tối ưu tốn (3.39)-(3.42) điều khiển tối ưu u∗ (·) liên tục khúc Lúc tồn hàm vector p(·) cho: a) hàm p(·) thỏa mãn phương trình liên hợp p(t) ˙ = −Hx (t, x∗ (t), u∗ (t), p(t), 1) = −fxT (t, x∗ (t), u∗ (t))p(t) + Lx (t, x∗ (t), u∗ (t)) (3.43) lu an điều kiện hoành va n p(t1 ) = 0; (3.44) p ie gh tn to b) điều kiện cực đại H(t, x∗ (t), u∗ (t), p(t), 1) = sup H(t, x∗ (t), u.p(t), 1) (3.45) w u∈U oa nl thỏa mãn hầu khắp [t0 , t1 ] d Nhắc lại an lu nf va H(t, x, u, p, 1) = (p| f (t, x, u)) − L(t, x, u) lm ul z at nh oi Kết luận trường hợp đặc biệt mơ hình tổng quát xem [1] trang 66 Ở không cần nhắc đến hai vector l0 vàl1 Đặc biệt, λ0 6= 0, giả sử λ0 = p(·) phải nghiệm phương trình vi phân z p(t1 ) = 0, gm @ p(t) ˙ = −fxT (t, x∗ (t), u∗ (t))p(t), m co l tức p(t) ≡ 0, mâu thuẫn với điều kiện nhân tử Lagrange khơng đồng thời triệt tiêu, Vì vậy, đặt λ0 = an Lu Chứng minh:Xem [1] trang 75-78 n va ac th si 65 Kết luận Luận văn trình bày điều kiện cần toán cực trị, điện kiện cần Fermat cho toán tối ưu hàm trơn đến tốn cực trị có ràng buộc, cho hàm mục tiêu trơn lồi không gian hữu hạn chiều Đặc biệt, luận văn trình bày điều kiện cần tối ưu cho toán biến phân (bài toán tối ưu không gian vô hạn chiều) lu an n va p ie gh tn to Trong luận văn trình bày chứng minh chi tiết điều kiện cần tối ưu cho toán biến phân sở toán biến phân Bolza, với chứng minh chi tiết Bổ đề Lagrange Bổ đề Du Bois-Reymond d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu n va ac th si 66 Tài liệu tham khảo Tiếng Việt [1] Hoàng Xuân Phú (1997), Bài giảng cao học Lý thuyết tốn cực trị, Hà Nội lu [2] Hồng Tụy (2005), Hàm thực giải tích hàm, NXB ĐHQG Hà Nội an n va Tiếng Anh ie gh tn to [3] V M Alekseev, V M Tikhomirov, S V Formin (1979), Optimal Control, tiếng Nga: Nauka, Moscow p [4] A D Ioffe, V M Tikhomirov (1974) Theory of Extremal Problems, tiếng Nga: Nauka, Moscow; tiếng Anh (1979), NorthHolland, Amsterdam oa nl w d [5] L S Pontryagin, V G Boltyasliĩ, R V Gamkrelidze, E F Mishchenko (1961), The mathematical Theory of Optimal Processes, tiếng Anh (1961), Fizmatgiz, Moscow; tiếng Nga (1962), Interscience, New York nf va an lu lm ul z at nh oi [6] Mike Mesterton-Gibbons (2009), STUDENT MATHEMATICAL LIBRARY Volume 50 A Primer on the Calculus of Variations and Optimal Control Theory z [7] Richard Vinter (2000), Optimal Control, Birkhauser, Boston, Basel, Berlin m co l gm @ an Lu n va ac th si