Luận văn Phương pháp lặp Banach cho bài toán bất đẳng thức biến phân

0đại học thái nguyên trường đại học khoa học ---------------------------- Phan Thế Nghĩa phương pháp lặp Banach cho BàI TOáN BấT ĐẳNG THứC BIếN PHÂN luận văn thạc sĩ toán học ứng dụng Thái Nguyên-2009 Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 0đại học thái nguyên trường đại học khoa học ---------------------------- Phan Thế Nghĩa phương pháp lặp Banach cho BàI TOáN BấT ĐẳNG THứC BIếN PHÂN Chuyên ngành: Toán ứng dụng Mã số: 60.46.36 luận văn thạc sĩ toán học ứng dụng NGUờI HướNG DẫN KHOA HọC: TS. PHạM NGọC ANH Thái Nguyên-2009 Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 01 Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 2Mục lục Trang phụ bìa Mục lục 2 Lời cảm ơn 3 Một số ký hiệu và chữ viết tắt 4 Lời nói đầu 5 Chương 1. Bài toán Bất đẳng thức biến phân 1.1. Một số khái niệm cơ bản 7 1.2. Phát biểu bài toán và ví dụ 10 1.3. Sự tồn tại nghiệm của bài toán VI 18 Chương 2. Phương pháp lặp Banach giải bài toán (VI) đơn điệu mạnh 2.1. Tính không giãn của ánh xạ nghiệm 23 2.2. Mô tả thuật toán và sự hội tụ 27 Chương 3. Phương pháp lặp Banach giải bài toán đồng bức 3.1. Tính không giãn của ánh xạ nghiệm 30 3.2. Mô tả thuật toán và sự hội tụ 35 3.3. Kết quả tính toán thử nghiệm 3.3.1. Mô hình cân bằng bán độc quyền 38 3.3.2. Kết quả tính toán thử nghiệm 43 Tài liệu tham khảo 44 Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 3Lời cảm ơn Bản luận văn này được hoàn thành tại trường Đại học Khoa học-Đại học Thái Nguyên dưới sự hướng dẫn của TS. Phạm Ngọc Anh. Tác giả xin bày tỏ lòng kính trọng và biết ơn sâu sắc tới thầy về sự tận tình hướng dẫn trong suốt thời gian tác giả làm luận văn. Trong quá trình học tập và làm luận văn, thông qua các bài giảng và xêmina, tác giả thường xuyên nhận được sự quan tâm giúp đỡ và đóng góp những ý kiến quý báu của PGS. TS. Lê Thị Thanh Nhàn, TS. Nguyễn Thị Thu Thủy và các thầy các cô trong trường Đại học Khoa học-Đại học Thái Nguyên. Từ đáy lòng mình, tác giả xin bày tỏ lòng biết ơn sâu sắc đến các thầy các cô. Tác giả xin bày tỏ lòng biết ơn tới các thầy, các cô khoa Khoa học Cơ bản, Ban Chấp Hành Đoàn trường Cao đẳng Công nghiệp Thái Nguyên đã đã tạo điều kiện giúp đỡ tác giả trong thời gian làm cao học. Xin chân thành cảm ơn anh chị em học viên cao học và bạn bè đồng nghiệp gần xa đã trao đổi, động viên và khích lệ tác giả trong quá trình học tập, nghiên cứu và làm luận văn. Luận văn sẽ không hoàn thành được nếu không có sự thông cảm, giúp đỡ của những người thân trong gia đình tác giả. Đây là món quà tinh thần, tác giả xin kính tặng gia đình thân yêu của mình với tấm lòng biết ơn chân thành và sâu sắc. Tác giả Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 4Một số ký hiệu và chữ viết tắt R n không gian Euclide n-chiều |β| trị tuyệt đối của số thực β x := y x được định nghĩa bằng y ∀x với mọi x ∃x tồn tại x I ánh xạ đồng nhất A ⊂ B tập A là tập con thực sự của tập B A ⊆ B tập A là tập con của tập B A ∪B A hợp với B A ∩B A giao với B AìB tích Đề-các của hai tập A và B convD bao lồi của tập D argmin{f(x) | x ∈ C} tập các điểm cực tiểu của hàm f trên C AT ma trận chuyển vị của ma trận A xk → x dãy {xk} hội tụ mạnh tới x V I bài toán bất đẳng thức biến phân Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 5Lời nói đầu Theo Harker và Pang, bài toán bất đẳng thức biến phân được giới thiệu lần đầu tiên vào năm 1966 bởi Hartman và Stampacchia. Những nghiên cứu đầu tiên về bất đẳng thức biến phân liên quan tới việc giải các bài toán biến phân, bài toán điều khiển tối ưu và các bài toán biên có dạng của phương trình đạo hàm riêng. Bài toán biến phân trong không gian vô hạn chiều và các ứng dụng của nó được giới thiệu trong cuốn sách "An introduction to variational inequalities and their application" của Kinderlehrer và Stampacchia xuất bản năm 1980 và trong cuốn sách "Variational and quasivariational inequalities: Application to free boundary problems" của Baiocchi và Capelo xuất bản năm 1984. Năm 1979 Michael J. Smith đưa ra bài toán cân bằng mạng giao thông và năm 1980 Defermos chỉ ra rằng: Điểm cần bằng của bài toán này là nghiệm của bài toán bất đẳng thức biến phân. Từ đó bài toán bất đẳng thức biến phân được phát triển và trở thành một công cụ hữu hiệu để nghiên cứu và giải các bài toán cân bằng trong kinh tế tài chính, vận tải, lý thuyết trò chơi và nhiều bài toán khác (xem [7]). Bài toán bất đẳng thức biến phân có quan hệ mật thiết với các bài toán tối ưu khác. Bài toán bù phi tuyến, xuất hiện vào năm 1964 trong luận án tiến sĩ của Cottle, là một trường hợp đặc biệt của bài toán bất đẳng thức biến phân (xem [5]). Gần đây, bài toán bất đẳng thức biến phân cũng là một đề tài được nhiều người quan tâm nghiên cứu vì vai trò của nó trong lý thuyết toán học và trong các ứng dụng thực tế (xem [5, 7]). Một trong các hướng nghiên cứu quan trọng của bài toán bất đẳng thức biến phân là việc xây dựng các phương pháp giải. Thông thường các phương pháp giải được chia thành các loại sau: Loại thứ nhất là các phương pháp chuyển bài toán về hệ phương trình và dùng các phương pháp thông dụng như phương pháp Newton, phương pháp điểm trong để giải hệ phương trình này. Loại thứ hai là phương pháp có tính chất kiểu đơn điệu. Điển hình của phương pháp này là các phương pháp gradient sau này được tổng quát bởi Cohen thành nguyên Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 6lý bài toán phụ (xem [5]), phương pháp điểm gần kề của Rockafellar (xem [3]), phương pháp hiệu chỉnh Tikhonov (xem [5]),.... Các phương pháp này khá hiệu quả, dễ thực thi trên máy tính nhưng các điều kiện hội tụ chỉ được đảm bảo dưới các giả thiết khác nhau về tính chất đơn điệu. Loại thứ ba là các phương pháp được dựa trên kỹ thuật hàm chắn (xem [5]). Nội dung chính của phương pháp này là chuyển bài toán bất đẳng thức biến phân về cực tiểu của hàm chắn và sau đó sử dụng kỹ thuật tối ưu trơn hoặc không trơn để tìm cực tiểu của hàm chắn. Phương pháp này có thể giải được các bài toán với các giả thiết rất nhẹ. Tuy nhiên, tốc độ hội tụ của thuật toán được đề xuất là chậm (xem [5]). Loại thứ tư là các phương pháp dựa trên cách tiếp cận điểm bất động. Nội dung chính của phương pháp này là chuyển bài toán bất đẳng thức biến phân về tìm điểm bất động của ánh xạ nghiệm. Luận văn này trình bày phương pháp giải bài toán bất đẳng thức biến phân thông qua tìm điểm bất động của ánh xạ nghiệm được viết trong bài báo "P. N. Anh, L. D. Muu, V. H. Nguyen and J. J. Strodiot (2005), On the contrac- tion and nonexpansiveness properties of the marginal mapping in generalized variational inequalities involving cocoercive operators, in: Generalized Con- vexity and Generalized Monotonicity and Applications. Eds A. Eberhard, N. Hadjisavvas and D. T. Luc, Springer, pp. 89-111". Ngoài lời nói đầu và phần tài liệu tham khảo, luận văn được chia làm ba chương. Chương 1 có tiêu đề là "bài toán bất đẳng thức biến phân". Chương này nhắc lại các kiến thức cơ bản về bài toán bất đẳng thức biến phân, các ví dụ, các kiến thức liên quan và các ứng dụng của bài toán bất đẳng thức biến phân. Chương 2 gồm hai phần cơ bản: Phần thứ nhất trình bày mối quan hệ giữa nghiệm của bài toán bất đẳng thức biến phân và ánh xạ nghiệm. Phần thứ hai chỉ ra ánh xạ nghiệm là co khi hàm giá là đơn điệu mạnh và Lipschitz. Chương 3 trình bày phương pháp lặp Banach cho ánh xạ đồng bức và một vài tính toán ứng dụng của thuật toán được đề xuất. Khi đó, ánh xạ nghiệm chỉ là không giãn và việc tìm điểm bất động của ánh xạ không giãn được tìm theo kiểu điểm bất động của Nadler. Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 7CHƯƠNG I BàI TOáN BấT ĐẳNG THứC BIếN PHÂN 1.1. Một số khái niệm cơ bản Cho hai véc tơ x := (x1, x2, ..., xn) T , y := (y1, y2, ..., yn) T ∈ Rn 〈x, y〉 = n∑ i=1 xiyi được gọi là tích vô hướng của hai véc tơ x và y. Chuẩn Euclide và khoảng cách được xác định tương ứng bởi ||x|| := √ 〈x, x〉, d(x, y) := ||x− y||. Ta nhắc lại một số kiến thức cơ bản của giải tích lồi sẽ được dùng cho các chương tiếp theo. Định nghĩa 1.1. • Tập con C ⊂ Rn được gọi là tập lồi, nếu λx + (1− λ)y ∈ C ∀x, y ∈ C, λ ∈ (0, 1). • Tập con C ⊂ Rn được gọi là nón, nếu λx ∈ C ∀x ∈ C, λ ≥ 0. • Cho C ⊂ Rn là một tập lồi và x ∈ C , nón pháp tuyến ngoài của C tại x, ký hiệu NC(x), được xác định bởi công thức NC(x) := {w ∈ R n : 〈w, y − x〉 ≤ 0 ∀y ∈ C}. Cho C ⊂ Rn là một tập lồi, ánh xạ f : C → Rn. Khi đó, Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 8Định nghĩa 1.2. • Miền hữu hiệu của f , ký hiệu dom f , được xác định bởi domf := {x ∈ Rn : f(x) < +∞}. • f được gọi là chính thường, nếu domf 6= ∅, f(x) > −∞ ∀x ∈ C. • f được gọi là hàm lồi trên C , nếu f(λx1 + (1− λ)x2) ≤ λf(x1) + (1− λ)f(x2) ∀x1, x2 ∈ C, λ ∈ [0, 1]. • f được gọi là hàm lồi chặt trên C , nếu f(λx1 + (1− λ)x2) < λf(x1) + (1− λ)f(x2) ∀x1 6= x2 ∈ C, λ ∈ (0, 1). • f được gọi là hàm lồi mạnh với hệ số β > 0 trên C , nếu ∀x1 6= x2 ∈ C, λ ∈ (0, 1), ta có f(λx1 + (1− λ)x2) < λf(x1) + (1− λ)f(x2)− λ(1− λ)β||x1 − x2||2. Bây giờ ta giả sử rằng f là một hàm lồi trên tập lồi C trong không gian Rn. Khi đó, véc tơ w ∈ Rn được gọi là dưới gradient của hàm f tại x ∈ C , nếu f(y)− f(x) ≥ 〈w, y − x〉 ∀y ∈ C. Tập tất cả các dưới gradient của hàm f tại x được gọi là dưới vi phân của f , ký hiệu ∂f(x), hay ∂f(x) := {w ∈ Rn : f(y)− f(x) ≥ 〈w, y − x〉 ∀y ∈ C}. Khi đó, f được gọi là khả dưới vi phân trên C , nếu ∂f(x) 6= ∅ ∀x ∈ C. Ví dụ 1.1. Cho C là một tập lồi khác rỗng của không gian Rn. Xét hàm chỉ trên tập C δ(x) :=   0 nếu x ∈ C, +∞ nếu x /∈ C. Khi đó ∂δC(x) = NC(x). Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 9Thật vậy, nếu x ∈ C thì δC(x) = 0 và ∂δC(x) = {w ∈ R n : δC(y) ≥ 〈w, y − x〉 ∀y ∈ C}. Hay ∂δC(x) = {w ∈ R n : 0 ≥ 〈w, y − x〉 ∀y ∈ C} = NC(x). Ví dụ 1.2. Trong không gian R n cho hàm chuẩn f(x) := ||x|| x ∈ Rn. Khi đó, ∂f(x) :=   {w ∈ Rn : ||w|| = 1, 〈w, x〉 = ||x||} nếu x 6= 0, B¯(0, 1) nếu x = 0, trong đó B¯(0, 1) là hình cầu đóng, tâm tại 0 và bán kính 1. Thật vậy, ta xét các trường hợp sau: Trường hợp 1. Với x 6= 0, ta cần chứng minh ∂f(x) = {w ∈ Rn : ||w|| = 1, 〈w, x〉 = ||x||}. Nếu w thỏa mãn ||w|| = 1, 〈w, x〉 = ||x|| thì 〈w, x〉 ≤ ||w||.||x|| = ||x||. Do đó 〈w, x− y〉 ≤ ||x|| − ||y||. Hay w ∈ ∂f(x). Ngược lại, nếu w ∈ ∂f(x), thì −||x|| = ||0|| − ||x|| ≥ 〈w, 0− x〉 = −〈w, x〉, ||x|| = ||2x|| − ||x|| ≥ 〈w, 2x− x〉 = 〈w, x〉 suy ra ||x|| = 〈w, x〉. (∗) Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 10 Mặt khác ||λz + x|| − ||x|| ≥ 〈w, λz + x− x〉 = 〈w, λz〉 ∀λ > 0, z ∈ Rn. Suy ra ||z + x λ || − 1 λ ||x|| ≥ 〈w, x〉. Cho λ → ∞, ta nhận được ||z|| ≥ 〈w, z〉 ∀z ∈ Rn. Do vậy ||w|| ≤ 1. Hơn nữa nếu ||w|| < 1 thì với mọi z ∈ Rn, ||z|| = 1 ta có |〈w, z〉| < 1. Khi đó, thay z = x||x|| ta có |〈w, z〉| = |〈w, x ||x|| 〉| < 1. Do đó 〈w, x〉 < ||x||. Điều này mâu thuẫn với (∗). Vậy ||w|| = 1. Trường hợp 2. Với x = 0. Ta có ∂f(x) = {w ∈ Rn : 〈w, y〉 ≤ ||y|| ∀y} = {w ∈ Rn : ||w|| ≤ 1} = B¯(0, 1). 1.2. Phát biểu bài toán và ví dụ "Bài toán bất đẳng thức biến phân" là một trong những bài toán được quan tâm nhiều trong toán học nói chung và đặc biệt trong ngành tối ưu tính toán nói riêng. Luận văn này sẽ trình bày một phương pháp giải bài toán bất đẳng thức biến phân trong không gian hữu hạn chiều. Chương này bao gồm việc nhắc lại các kiến thức cơ bản nhất về bài toán bất đẳng thức biến phân sẽ được sử dụng cho các chương sau. Bài toán bất đẳng thức biến phân trong không gian hữu hạn chiều có thể được phát biểu như sau: Cho C là một tập con lồi, đóng khác rỗng của không gian Euclidean n-chiều R n , F : C → Rn là ánh xạ liên tục. Bài toán bất đẳng thức biến phân Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 11 (viết tắt là: VI) là bài toán tìm điểm x∗ ∈ C , sao cho: 〈F (x∗), x− x∗〉 ≥ 0 ∀x ∈ C. (1.1) Tập nghiệm của VI ký hiệu là S∗. Định nghĩa 1.3. Cho C là tập lồi, đóng trong Rn, và cho F : C → Rn là một ánh xạ. Khi đó, F được gọi là: (a) đơn điệu trên C , nếu: 〈F (u)− F (v) , u− v〉 ≥ 0 ∀u, v ∈ C (b) đơn điệu ngặt trên C , nếu: 〈F (u)− F (v) , u− v〉 > 0 ∀u, v ∈ C, u 6= v. (c) đơn điệu mạnh trên C với hằng số τ > 0 (viết tắt là:τ -đơn điệu mạnh) nếu: 〈F (u)− F (v) , u− v〉 ≥ τ ‖u− v‖2 ∀u, v ∈ C. (d) đồng bức với mô đun δ (viết tắt là:δ-đồng bức) trên C nếu tồn tại một số δ > 0 sao cho: 〈F (u)− F (v), u− v〉 ≥ δ||F (u)− F (v)||2 u, v ∈ C. Ta nhắc lại kết quả tương đương sau: Nhận xét 1.1. Cho C là một tập lồi và F : C → Rn là một ánh xạ khả vi liên tục trên tập mở chứa C . Khi đó, i) F đơn điệu trên C khi và chỉ khi ∇F (x) là nửa xác định dương trên C hay 〈y,∇F (x)y〉 ≥ 0 ∀y ∈ C. ii) F đơn điệu chặt trên C khi và chỉ khi ∇F (x) là xác định dương trên C hay 〈y,∇F (x)y〉 > 0 ∀y ∈ C, y 6= 0. iii) F đơn điệu mạnh trên C khi và chỉ khi ∇F (x) là xác định dương đều trên C hay tồn tại β > 0 sao cho 〈y,∇F (x)y〉 > β||y||2 ∀y ∈ C, y 6= 0. Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 12 Các ví dụ dưới đây cho ta thấy được cách tiếp cận của bài toán bất đẳng thức biến phân Ví dụ 1.3. Cho f(x) là một hàm thực khả vi trên tập mở chứa C = [a, b]. Tìm x0 ∈ C sao cho f(x0) = min x∈C f(x). x0 ∈ [a, b], suy ra có 3 trường hợp xảy ra: TH1: Nếu x0 ∈ (a, b), theo định lý Fermat, ta có f ′(x0) = 0. TH2: Nếu x0 = a, f ′(x0) = lim x→x0+ f(x)−f(x0) x−x0 ≥ 0. TH2: Nếu x0 = b, f ′(x0) = lim x→x0 − f(x)−f(x0) x−x0 ≤ 0. Kết hợp lại, ta có thể viết: x0 là nghiệm của bài toán f ′(x0).(x− x0) ≥ 0 ∀x ∈ C. Như vậy x0 là nghiệm của bài toán bất đẳng thức biến phân VI với F = f ′ trên C = [a, b]. Bây giờ, ta xét ví dụ tổng quát hơn Ví dụ 1.4. Cho f(x) là một hàm thực khả vi trên tập mở chứa C ⊆ IRn. Tìm x0 ∈ C sao cho f(x0) = min x∈C f(x). Mệnh đề 1.1. Nếu x0 là nghiệm của bài toán trên, thì x0 là nghiệm của bài toán bất đẳng thức biến phân VI với F (x) := ∇f(x). Chứng minh. Với mọi y ∈ C , do C lồi nên (1− t)x0 + ty ∈ C ∀t ∈ [0, 1]. Đặt ϕ(t) := f(x0 + t(y − x0)). Giả thiết cho x0 là nghiệm hay t = 0 là nghiệm của ϕ(t) trên [0, 1]. Theo Ví dụ 1.3, ta có ϕ′(t0).(t− t0) ≥ 0 ∀t ∈ [0, 1]. Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 13 Hay 〈∇f(x0), x− x0〉 ≥ 0 ∀x ∈ C. 2 Mệnh đề 1.2. Cho f là hàm lồi khả vi trên tập lồi C ⊆ Rn. Khi đó, x0 ∈ C là nghiệm của bài toán min x∈C f(x) khi và chỉ khi x0 là nghiệm của bài toán bất đẳng thức VI với F (x) := ∇f(x). Chứng minh. Điều kiện cần được suy ra từ Mệnh đề 1.1. Do f là hàm lồi trên C , nên f(x)− f(x0) ≥ 〈∇f(x0), x− x0〉 ∀x ∈ C. Giả thiết cho 〈∇f(x0), x− x0〉 ≥ 0 ∀x ∈ C. Do đó f(x) ≥ f(x0) ∀x ∈ C. Hay x0 là nghiệm của bài toán min x∈C f(x). 2 Ví dụ 1.5. (Bài toán bù, ký hiệu CP) Cho C = Rn+ và F : C → R n . Bài toán được đặt ra là: Tìm điểm x0 ∈ C sao cho F (x0) ∈ C, 〈F (x0), x0〉 = 0. Mệnh đề 1.3. x0 ∈ C = Rn+ là nghiệm của bài toán bù CP khi và chỉ khi x 0 là nghiệm của bài toán VI hay 〈F (x0), x− x0〉 ≥ 0 ∀x ∈ C. Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 14 Chứng minh. (⇒) Giả sử x0 là nghiệm của bài toán bù CP hay F (x0) ∈ C, 〈F (x0), x0〉 = 0. Khi đó 〈F (x0), x− x0〉 = 〈F (x0), x〉 − 〈F (x0), x0〉 = 〈F (x0), x〉 ≥ 0 ∀x ∈ C. (⇐) Giả sử x0 là nghiệm của bài toán bất đẳng thức biến phân VI hay x0 ∈ C : 〈F (x0), x− x0〉 ≥ 0 ∀x ∈ C. Gọi ei = (0, 0, ..., 0, 1, 0, ...0) T (1 ở vị trí thứ i). Khi đó, x1 = x0 + ei ∈ C . Thay x1 vào bất đẳng thức biến phân, ta có 〈F (x0), x1 − x0〉 ≥ 0. Hay 〈F (x0), ei〉 ≥ 0 ∀i = 1, 2, ...n. Vậy F (x0) ∈ C . Từ 0 ∈ C và 〈F (x0), x− x0〉 ≥ 0 ∀x ∈ C. suy ra −〈F (x0), x0〉 ≥ 0. Do đó 〈F (x0), x0〉 = 0. 2 Dưới đây ta xét hai ví dụ thực tế của bài toán VI. Ví dụ 1.6. Bài toán cân bằng mạng giao thông Xét một mạng giao thông được cho bởi một mạng luồng hữu hạn. Gọi •N : tập hợp các nút của mạng. •A: là tập hợp các cạnh (mỗi cạnh được gọi là một đoạn đường). Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 15 Giả sử O ⊆ N , D ⊆ N sao cho O ∩ D = ∅. Mỗi phần tử của O được gọi là điểm nguồn, còn mỗi phần tử của D được gọi là điểm đích. Mỗi điểm nguồn và điểm đích được nối với nhau bởi một tập hợp liên tiếp các cạnh (được gọi là một tuyến đường). Ký hiệu: •f ia là mật độ giao thông của phương tiện i trên đoạn đường a ∈ A. Đặt f là véc tơ có các thành phần là f ia với i ∈ I và a ∈ A (I là tập hợp các phương tiện giao thông. •cia là chi phí khi sử dụng phương tiện giao thông i trên đoạn đường A. Đặt c là véc tơ có các thành phần là cia với i ∈ I, a ∈ A. •diw là nhu cầu sử dụng loại phương tiện i ∈ I trên tuyến đường w = (O,D) với O ∈ O, D ∈ D. Giả sử rằng chi phí giao thông phụ thuộc vào lưu lượng, tức là c = c(f) là một hàm của f . •λiw là mức độ chi phí trên tuyến đường w của phương tiện giao thông i. •xiw là mật độ giao thông của phương tiện i ∈ I trên tuyến w ∈ O ìD. Giả sử trong mạng trên, phương trình cân bằng sau được thoả mãn diw = ∑ p∈Pw xip ∀i ∈ I, w ∈ O ìD, (1.2) trong đó, Pw ký hiệu tập hợp các tuyến đường của w = (O,D) (nối điểm nguồn O và điểm đích D). Theo phương trình (2.18), thì nhu cầu sử dụng loại phương tiện i trên tuyến đường w bằng đúng tổng mật độ giao thông của phương tiện đó trên mọi tuyến đường nối điểm nguồn và điểm đích của tuyến đường đó. Khi đó ta có f ia = ∑ p∈Pw xipδap ∀i ∈ I, w ∈ O ìD, (1.3) trong đó δap :=   1 nếu a ∈ p, 0 nếu a /∈ p. Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 16 Với mỗi tuyến đường p nối một điểm nguồn và một điểm đích, đặt cip = ∑ a∈A ciaδap. (1.4) Như vậy, cip là một chi phí khi sử dụng phương tiện i trên tuyến đường p. Đặt d là véc tơ có các thành phần là diw (i ∈ I, w ∈ O ìD) và đặt f là véc tơ có các thành phần là dia (i ∈ I, a ∈ O ìD). Một cặp (d ∗, f ∗) thoả mãn các điều kiện (2.18) và (2.21) được gọi là điểm cân bằng của mạng giao thông nếu cip(f ∗) =   λiw(d ∗) khi xip > 0, > λiw(d ∗) khi xip = 0, với mỗi i ∈ I và mỗi tuyến đường p. Theo định nghĩa này, tại điểm cân bằng đối với mọi loại phương tiện giao thông và mọi tuyến đường, chi phí sẽ thấp nhất khi có lưu lượng giao thông trên tuyến đó. Trái lại, chi phí sẽ không phải thấp nhất. Đặt K = {(f, d) | ∃ x ≥ 0 sao cho (2.18) và (2.21) đúng}. Khi đó, ta có định lý sau. Định lý 1.1. Một cặp véc tơ (f ∗, d∗) ∈ K là một điểm cân bằng của mạng giao thông khi và chỉ khi nó là nghiệm của bất đẳng thức biến phân sau: Tìm (f ∗, d∗) ∈ K sao cho 〈 ( c(f ∗)), λ(d∗) ) , (f, d)−(f ∗, d∗)〉 ≥ 0 ∀(f, d) ∈ K. Ví dụ 1.7. Bài toán kinh tế bán độc quyền Giả sử có n công ty cùng sản xuất một loại sản phẩm và lợi nhuận pi của mỗi công ty i phụ thuộc vào tổng số lượng sản phẩm của tất cả các công ty σ := ∑n i=1 xi. Ký hiệu hi(xi) là chi phí của công ty i khi sản xuất ra lượng hàng hoá xi. Giả sử rằng lợi nhuận của công ty i được cho bởi fi(x1, ..., xn) = xipi( n∑ i=1 xi)− hi(xi) (i = 1, ..., n), (1.5) Số húa bởi Trung tõm Học liệu – Đại học Thỏi Nguyờn www.VNMATH.com 17 trong đó p( ∑n j=1 xj) là giá của một đơn vị sản phẩm, phụ thuộc vào tổng sản phẩm, còn hàm chi phí của mỗi công ty i chỉ phụ thuộc vào mức độ sản xuất của công ty đó. Đặt Ui ⊂ IR, (i = 1, ..., n) là tập chiến lược của công ty i. Lẽ dĩ nhiên, mỗi công ty cần xác định cho mình một mức độ sản xuất để đạt được lợi nhuận cao nhất. Tuy nhiên, trong trường hợp tổng quát, việc tất cả các công ty đều có lợi nhuận cực đại là khó có thể được. Vì vậy người ta dùng đến khái niệm cân bằng: Một điểm x∗ = (x∗1, ..., x ∗ n) ∈ U := U1 ì ... ì Un được gọi là điểm cân bằng Nash nếu fi(x ∗ 1, ..., x ∗ i−1, yi, x ∗ i+1, ..., x ∗ n) 6 fi(x ∗ 1, ..., x ∗ n) ∀yi ∈ Ui, ∀i = 1, ..., n. Trong mô hình cân bằng Cournot cổ điển, hàm chi