Tài liệu Phương pháp giảm cơ sở giải phương trình elliptic bức tuyến tính phụ thuộc tham số

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC PHẠM THỊ THÙY NHUNG PHƯƠNG PHÁP GIẢM CƠ SỞ GIẢI PHƯƠNG TRÌNH ELLIPTIC BỨC TUYẾN TÍNH PHỤ THUỘC THAM SỐ LUẬN VĂN THẠC SĨ TOÁN HỌC Thái Nguyên - 2015 ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC PHẠM THỊ THÙY NHUNG PHƯƠNG PHÁP GIẢM CƠ SỞ GIẢI PHƯƠNG TRÌNH ELLIPTIC BỨC TUYẾN TÍNH PHỤ THUỘC THAM SỐ Chuyên ngành: Toán ứng dụng Mã số: 60 46 01 12 LUẬN VĂN THẠC SĨ TOÁN HỌC Người hướng dẫn khoa học: TS. NGUYỄN THANH SƠN Thái Nguyên - 2015 i Mục lục Tóm tắt nội dung iii Lời cảm ơn iv Mở đầu 1 1 0.1 Lý do chọn đề tài . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 0.2 Mục đích nghiên cứu . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 0.3 Nhiệm vụ nghiên cứu . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 0.4 Đối tượng và phạm vi nghiên cứu . . . . . . . . . . . . . . . . . . . . 4 0.5 Phương pháp nghiên cứu . . . . . . . . . . . . . . . . . . . . . . . . . 4 Kiến thức chuẩn bị 5 1.1 Không gian tích vô hướng . . . . . . . . . . . . . . . . . . . . . . . . . 5 1.2 Dạng tuyến tính và dạng song tuyến tính . . . . . . . . . . . . . . . . 6 1.2.1 Dạng tuyến tính . . . . . . . . . . . . . . . . . . . . . . . . . . 6 1.2.2 Dạng song tuyến tính . . . . . . . . . . . . . . . . . . . . . . . 7 Dạng tuyến tính và dạng song tuyến tính phụ thuộc tham số . . . . . . 8 1.3.1 Sự phụ thuộc affine vào tham số . . . . . . . . . . . . . . . . 9 1.3.2 Tính bức phụ thuộc tham số . . . . . . . . . . . . . . . . . . . 9 1.3 2 Phương pháp giảm cơ sở 11 2.1 Không gian hàm và dạng yếu của phương trình elliptic . . . . . . . . . 11 2.1.1 11 Dạng yếu phụ thuộc tham số . . . . . . . . . . . . . . . . . . . ii 2.1.2 Tích vô hướng và chuẩn khác . . . . . . . . . . . . . . . . . . 12 Không gian và cơ sở xấp xỉ hữu hạn chiều . . . . . . . . . . . . . . . . 13 2.2.1 Rời rạc hóa bài toán . . . . . . . . . . . . . . . . . . . . . . . 13 2.2.2 Phép chiếu Galerkin . . . . . . . . . . . . . . . . . . . . . . . 14 2.2.3 Phương trình đại số . . . . . . . . . . . . . . . . . . . . . . . . 15 Giảm cơ sở . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17 2.3.1 Không gian và cơ sở . . . . . . . . . . . . . . . . . . . . . . . 17 2.3.2 Phép chiếu lên không gian số chiều nhỏ . . . . . . . . . . . . 18 2.3.3 Phương trình đại số . . . . . . . . . . . . . . . . . . . . . . . . 19 2.4 Thủ tục tính toán Online- Offline . . . . . . . . . . . . . . . . . . . . 20 2.5 Thuật toán Greedy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 2.6 Ước lượng sai số hậu nghiệm . . . . . . . . . . . . . . . . . . . . . . . 23 2.6.1 Cận dưới của hằng số liên tục . . . . . . . . . . . . . . . . . . 24 2.6.2 Cận trên của hằng số liên tục . . . . . . . . . . . . . . . . . . 25 2.6.3 Những kiến thức cần thiết khác . . . . . . . . . . . . . . . . . 27 2.6.4 Ước lượng sai số . . . . . . . . . . . . . . . . . . . . . . . . . 28 2.6.5 Cận của sai số tương ứng với chuẩn trên X . . . . . . . . . . . 30 2.2 2.3 3 Ví dụ số 32 Kết luận 34 Tài liệu tham khảo 35 iii TÓM TẮT NỘI DUNG Luận văn được viết có nội dung về phương pháp giảm cơ sở cho phương trình elliptic phụ thuộc tham số. Chúng tôi trình bày đầy đủ các nguyên liệu để người đọc có thể hiểu một cách chi tiết việc xây dựng phương pháp, cụ thể là phép chiếu Galerkin, thuật toán Greedy và cận dưới của hằng số bức phụ thuộc tham số sử dụng phép cách tiếp cận min −θ. iv Lời cảm ơn Luận văn này được thực hiện và hoàn thành tại Trường Đại học Khoa học - Đại học Thái Nguyên. Tác giả xin chân thành cảm ơn các thầy cô giáo đã tận tình giảng dạy, bồi dưỡng kiến thức trong suốt quá trình học tập, nghiên cứu và rèn luyện tại trường. Tác giả xin tỏ lòng biết ơn sâu sắc tới TS. Nguyễn Thanh Sơn đã tận tình hướng dẫn trong suốt quá trình viết luận văn. Xin chân thành cảm ơn. Thái Nguyên, 2015 Phạm Thị Thùy Nhung Học viên Cao học Toán K7A, Trường ĐH Khoa học - ĐH Thái Nguyên 1 Mở đầu 0.1 Lý do chọn đề tài Để bắt đầu, ta hãy xét Bài toán Dirichlet cho phương trình Poisson −∆y = f, trong Ω, (1) y = 0, trên Γ, trong đó Ω là miền bị chặn trong Rk , k ≥ 2 có biên Γ thỏa mãn điều kiện Lipschitz và f ∈ L2 (Ω). Với dữ liệu trên Bài toán (1) không thể có nghiệm cổ điển y ∈ C 2 (Ω)∩C(Ω̄) mà thay vào đó, người ta xét nghiệm suy rộng của Bài toán (1). Nhắc lại rằng, L2 (Ω) là một không gian Hilbert với tích vô hướng Z hf, gi = f (x)g(x)dx và H01 (Ω) = {g ∈ L2 (Ω) : ∃ L2 (Ω) ∂g ∈ L2Ω và g|Γ = 0} ∂xi là không gian hàm suy rộng trên Ω. Bài toán (1) được đưa về dạng đẳng thức biến phân như sau: Tìm y ∈ H01 (Ω) Z Z 5y 5 vdx = f (x)v(x)dx, ∀y ∈ H. Ω Ω Bây giờ ta định nghĩa dạng tuyến tính f : H01 (Ω) → R trong đó, f (v) với v ∈ H01 (Ω), f (v) = R Ω f (x)v(x)dx và dạng song tuyến tính a : H01 (Ω) × H01 (Ω) → R (2) 2 với a(w, v) = R Ω 5y. 5 vdx . Theo đó Bài toán (2) có thể được viết lại như sau a(y, v) = f (v), ∀v ∈ H01 (Ω). Lưu ý rằng, H01 (Ω) là không gian con của không gian H 1 (Ω) với tích vô hướng Z hw, vi = (wv + 5w 5 v)dx Ω và chuẩn tương ứng Z ||v||H 1 (Ω) = 2 2 1/2 (v + 5v )dx . (3) Ω Theo Bất đẳng thức Friedrichs Z Z 2 v d(x) ≤ C(Ω) Ω 5v 2 dx, Ω với C(Ω) là một hằng số chỉ phụ thuộc vào Ω, ta suy ra a(v, v) ≥ 1 ||v||2H 1 (Ω) . 1 + C(Ω) (4) Bất đẳng thức (4) được gọi là tính chất bức. Thêm vào đó, từ Bất đẳng thức CauchySchwarz ta suy ra |a(w, v)| ≤ ||w||H 1 (Ω) ||v||H 1 (Ω) , (5) đây được gọi là tính liên tục của a. Trong khi xét Bài toán (1) là ta đã ngầm định rằng phương trình đó là không phụ thuộc tham số. Trong thực tế, khi sử dụng Bài toán (1) như là một mô hình toán học cho hiện tượng truyền nhiệt dừng, người ta có thể muốn giữ lại hằng số dẫn nhiệt như là một tham số. Việc này cho phép sử dụng mô hình cho nhiều chất dẫn nhiệt có hằng số dẫn nhiệt khác nhau. Tương tự như thế, hàm nguồn f cũng có thể phụ thuộc vào tham số. Tổng quát, người ta có thể xét bài toán phụ thuộc tham số a(y, v; µ) = f (v; µ), ∀v ∈ V, µ ∈ D ⊂ Rk . (6) Và cuối cùng, trong rất nhiều trường hợp, người ta không quan tâm đến toàn bộ trạng thái y mà chỉ là một phần thông tin của trạng thái được biểu diễn bởi s(µ) = `(y; µ), (7) 3 trong đó, g cũng là một dạng tuyến tính phụ thuộc tham số. Bài toán đặt ra ở đây là: với mỗi µ ∈ D, tính s(µ). Đối với những bài toán thực tế cùng yêu cầu cao về tính chính xác, cỡ của bài toán thường rất lớn, hàng nghìn cho đến hàng triệu. Với đặc thù cần phải tính nghiệm, hoặc thông tin liên quan đến nghiệm nhiều lần (many-query context), máy tính sẽ phải mất nhiều thời gian để giải mô hình cỡ lớn với nhiều giá trị khác nhau của tham số. Yêu cầu đặt ra là tìm được một thuật toán để sao cho, việc tính toán nghiệm nhanh nhưng vẫn đảm bảo kiểm soát được sai số. Do vậy, phương pháp giảm cơ sở là đặc biệt cần thiết để xử lý những bài toán dạng này. Với lí do trên, chúng tôi đã chọn "Phương pháp giảm cơ sở giải phương trình elliptic bức tuyến tính phụ thuộc tham số" làm đề tài cho luận văn thạc sĩ. Luận văn gồm 3 chương. • Chương 1: Kiến thức chuẩn bị Trong chương này, chúng tôi sẽ nhắc lại những lý thuyết cần thiết cho việc trình bày phương pháp ở chương sau. Chúng bao gồm phép chiếu Galerkin, rời rạc hóa phương trình, tính bị chặn, tính bức, tính phụ thuộc affine . . . • Chương 2: Phương pháp giảm cơ sở Trong chương này, trước tiên chúng tôi sẽ trình bày sơ lược cách rời rạc hóa bài toán (6) - (7) bằng phương pháp phần tử hữu hạn. Sau đó, phép chiếu Galerkin được sử dụng để giảm kích cỡ của cơ sở cũng đồng thời giảm kích cỡ của bài toán. Cuối cùng, chúng tôi tập trung trình bày các "nguyên liệu" cũng như kết quả phục vụ cho việc xây dựng, cơ sở giảm đó. 4 • Chương 3: Ví dụ số Chương này đưa ra 2 ví dụ minh họa được lập trình trên MATLAP với dữ liệu được lấy từ mô hình thực tế trong khoa học kĩ thuật. 0.2 Mục đích nghiên cứu Nghiên cứu lý thuyết phương pháp giảm cơ sở và áp dụng nó để giải phương trình elliptic bức tuyến tính phụ thuộc tham số. 0.3 Nhiệm vụ nghiên cứu Luận văn tập trung làm rõ một số vấn đề sau đây, trình bày ý tưởng của phương pháp giảm cơ sở giải phương trình elliptic bức tuyến tính phụ thuộc tham số, các khái niệm và tính chất liên quan đến phương pháp, nội dung phương pháp và cuối cùng là áp dụng phương pháp này cho một số ví dụ thực tế. 0.4 Đối tượng và phạm vi nghiên cứu • Đối tượng nghiên cứu: Phương pháp giảm cơ sở • Phạm vi nghiên cứu: Bài toán elliptic bức tuyến tính phụ thuộc affine vào tham số. 0.5 Phương pháp nghiên cứu • Đọc và nghiên cứu một số tài liệu liên quan như sách, báo, tạp chí, luận văn thạc sĩ, tiến sĩ. • Sử dụng nhiều kiến thức của đại số tuyến tính ứng dụng, giải tích hàm. • Kiểm chứng các kết quả lý thuyết bằng ví dụ số lập trình trên MATLAP. 5 Chương 1 Kiến thức chuẩn bị Trong chương này, chúng tôi sẽ nhắc lại những lý thuyết cần thiết cho việc trình bày phương pháp ở chương sau. Chúng bao gồm phép chiếu Galerkin, rời rạc hóa phương trình, tính bị chặn, tính bức, tính phụ thuộc affine . . . 1.1 Không gian tích vô hướng Định nghĩa 1.1. Cho Z là một không gian tuyến tính trên R. Khi đó một ánh xạ w ∈ Z, v ∈ Z → (w, v)Z ∈ R thỏa mãn với mọi w, v, z ∈ Z, α ∈ R thì • (αw + v, z)Z = α(w, z)Z + (v, z)Z và (z, αw + v)Z = α(z, w)Z + (z, v)Z , • (w, v)Z = (v, w)Z , • (w, w)Z ≥ 0. (w, w)Z = 0 ⇔ w = 0, được gọi là tích vô hướng trên Z . Định nghĩa 1.2. Cho một tích vô hướng trên Z . Khi đó, với hàm số ||·||: Z → R+ xác định bởi ||x|| = p hx, xi là một chuẩn trên Z (chuẩn sinh bởi tích vô hướng). Sau đây là kết quả kinh điển liên quan đến tích vô hướng, Bất đẳng thức CauchySchwarz |(w, v)Z | ≤ ||w||Z ||v||Z , ∀w, v ∈ Z. 6 Định nghĩa 1.3. Một không gian tuyến tính Z có tích vô hướng với chuẩn được trang bị thì không gian đó là không gian định chuẩn. Nếu không gian đó là đầy đủ thì được gọi là không gian Hilbert. 1.2 1.2.1 Dạng tuyến tính và dạng song tuyến tính Dạng tuyến tính Định nghĩa 1.4. Một ánh xạ g : Z → R được gọi là một ánh xạ tuyến tính hay dạng tuyến tính nếu, cho bất kì α ∈ R, w, v ∈ Z, g(αw + v) = αg(w) + g(v). Dạng tuyến tính g là bị chặn hay liên tục trên Z nếu |g(v)| ≤ C||v||Z , ∀v ∈ Z, với C là hằng số không phụ thuộc vào v . Định nghĩa 1.5. Cho Z là không gian định chuẩn. Chúng tôi định nghĩa không gian đối ngẫu của Z , kí hiệu Z 0 là tập hợp tất cả các phiếm hàm tuyến tính liên tục trên Z . Z’ với chuẩn kgkZ 0 = sup v∈Z g(v) 0 0 , ∀g ∈ Z , kvkZ trở thành một không gian định chuẩn, chuẩn đó được gọi là chuẩn đối ngẫu. Định lý 1.1. (Định lí Riesz) Với bất kì g ∈ Z 0 , tồn tại duy nhất wg ∈ Z sao cho (wg , v)Z = g(v), ∀v ∈ Z và kgkZ 0 = kwg kZ . 7 1.2.2 Dạng song tuyến tính Định nghĩa 1.6. Ánh xạ b : Z1 × Z2 → R là một dạng song tuyến tính nếu với bất kì α ∈ R, w, v ∈ Z1 , z ∈ Z2 , b(αw + v, z) = αb(w, z) + b(v, z) và với bất kì α ∈ R, z ∈ Z1 , w, v ∈ Z2 , b(z, αw + v) = αb(z, w) + b(z, v) một cách ngắn gọn, một dạng song tuyến tính là tuyến tính theo từng biến. Trong phần còn lại của mục này ta xét trường hợp riêng Z1 = Z2 = Z . Sau đây, chúng tôi liệt kê một số tính chất hay dùng của dạng song tuyến tính. Dạng song tuyến tính b : Z × Z → R là đối xứng nếu với bất kì w, v ∈ Z , b(w, v) = b(v, w). Dạng song tuyến tính b : Z × Z → R là phản đối xứng nếu với bất kì w, v ∈ Z , b(w, v) = −b(v, w) . Chúng tôi định nghĩa phần đối xứng và phản đối xứng của dạng song tuyến tính tổng quát b : Z × Z → R như sau: bs (w, v) = 1/2(b(w, v) + b(v, w)), ∀w, v ∈ Z và bS (w, v) = 1/2(b(w, v) − b(v, w)), ∀w, v ∈ Z, tương ứng. Một dạng song tuyến tính b : Z × Z → R là xác định dương nếu, cho bất kì v ∈ Z, b(v, v)(= bs (v, v) ≥ 0) với dấu bằng xảy ra chỉ khi v = 0. Một dạng song tuyến tính b : Z × Z → R là dạng song tuyến tính nửa xác định dương nếu với bất kì v ∈ Z, b(v, v) ≥ 0. Dạng song tuyến tính b : Z × Z → R là bức trên Z nếu với α là hằng số bức α ≡ inf w∈Z b(w, w) kwk2Z > 0. (1.1) Chú ý rằng tính bức chỉ liên quan đến phần đối xứng của b. Theo đó, chúng ta thay thế b trong (1.1) bởi bs . Một dạng song tuyến tính Z × Z → R là liên tục trên Z nếu với γ là hằng số liên tục γ ≡ sup sup w∈Z v∈Z b(w, w) < +∞. ||w||Z ||v||Z 8 1.3 Dạng tuyến tính và dạng song tuyến tính phụ thuộc tham số Đầu tiên chúng tôi giới thiệu một tham số trên miền đóng, bị chặn D ⊂ Rd . Ta kí hiệu các véc tơ tham số trên miền D là µ = (µ1 , ..., µd ). Ta nói g : Z × D → R là dạng tuyến tính phụ thuộc tham số nếu, với mỗi µ ∈ D, g(.; µ) : Z → R là một dạng tuyến tính. Ta nói dạng tuyến tính phụ thuộc tham số g là bị chặn (hoặc liên tục) nếu với mọi µ ∈ D, g(.; µ) ∈ Z 0 . Chú ý rằng chuẩn đối ngẫu của một dạng tuyến tính phụ thuộc tham số g , kg(.; µ)k cũng đúng với hàm hữu hạn phụ thuộc vào µ trên D. Tương tự, b : Z × Z × D → R là một dạng song tuyến tính phụ thuộc tham số nếu với mọi µ ∈ D, b(., .; µ) : Z × Z → R là một dạng song tuyến tính. Ta nói một dạng song tuyến tính phụ thuộc tham số b : Z × Z → R là đối xứng nếu b(w, v, µ) = b(v, w, µ), ∀w, v ∈ Z, ∀µ ∈ D. Chúng tôi định nghĩa phần đối xứng của dạng song tuyến tính phụ thuộc tham số b : Z × Z × D → R như sau bS (w, v; µ) ≡ 1/2(b(w, v; µ) + b(v, w; µ)), ∀w, v ∈ Z, ∀µ ∈ D. Dạng song tuyến tính phụ thuộc tham số b: Z × Z × D → R là bức trên Z nếu α(µ) ≡ inf w∈Z b(w, w; µ) > 0, ∀µ ∈ D. kwkZ kvkZ (1.2) Chúng tôi định nghĩa (0<) α0 ≡ minµ∈D α(µ). Dạng song tuyến tính phụ thuộc tham số b : Z × Z × D → R là liên tục trên Z nếu γ(µ) ≡ sup sup w∈Z v∈Z b(w, v; µ) < +∞, ∀µ ∈ D. kwkZ kvkZ Ta định nghĩa γ0 = maxµ∈D γ(µ)(< ∞). (1.3) 9 1.3.1 Sự phụ thuộc affine vào tham số Dạng tuyến tính phụ thuộc tham số bị chặn g : Z × D → R là affine theo tham số hay phụ thuộc affine vào tham số µ g(v; µ) = Qg X θgq (µ)g q (v), ∀v ∈ Z, q=1 cho một số Qg hữu hạn; ở đây θgq : D → R, 1 ≤ q ≤ Qg , là hàm phụ thuộc tham số, và g q (v) : Z → R, 1 ≤ q ≤ Qg , là dạng tuyến tính bị chặn độc lập với tham số µ. Cũng như vậy, chúng ta nói rằng dạng song tuyến tính phụ thuộc tham số b : Z × Z × D → R là affine theo tham số hay phụ thuộc affine vào tham số µ nếu b(w, v, µ) = Qb X θbq (µ)bq (w, v), ∀w ∈ Z1 , v ∈ Z2 , q=1 cho một số hữu hạn Qb ; ở đây θbq : D → R, 1 ≤ q ≤ Qb là hàm phụ thuộc tham số, và bq (w, v) : Z × Z → R, 1 ≤ q ≤ Qb là dạng song tuyến tính liên tục độc lập với tham số. 1.3.2 Tính bức phụ thuộc tham số Chúng ta nói một dạng tuyến tính phụ thuộc affine vào tham số b : Z ×Z ×D → R b(w, v; µ) = Qb X θbq (µ)bq (w, v) q=1 là bức tham số nếu c ≡ bS (phần đối xứng của b) thừa nhận một khai triển affine c(w, v; µ) = Qc X θcq (µ)cq (w, v), ∀w, v ∈ Z, q=1 thỏa mãn hai điều kiện: θcq (µ) > 0, ∀µ ∈ D, 1 ≤ q ≤ Qc (1.4) cq (v, v) ≥ 0, ∀v ∈ Z, 1 ≤ q ≤ Qc . (1.5) và 10 Đương nhiên, nếu bản thân hàm b đã đối xứng thì ta không cần phải dùng đến phần tuyến tính khi định nghĩa tính chất. Tức là, θbq (µ) và bq (w, v) thỏa mãn (1.4), (1.5). 11 Chương 2 Phương pháp giảm cơ sở Trong chương này, trước tiên chúng tôi sẽ trình bày sơ lược cách rời rạc hóa bài toán (6) - (7) bằng phương pháp phần tử hữu hạn. Sau đó, phép chiếu Galerkin được sử dụng để giảm kích cỡ của cơ sở cũng đồng thời giảm kích cỡ của bài toán. Cuối cùng, chúng tôi tập trung trình bày các "nguyên liệu" cũng như kết quả phục vụ cho việc xây dựng, cơ sở giảm đó. 2.1 2.1.1 Không gian hàm và dạng yếu của phương trình elliptic Dạng yếu phụ thuộc tham số Từ đây trở đi, chúng tôi sử dụng X e để kí hiệu không gian hàm mà bài toán viết dưới dạng đẳng thức biến phân được đưa ra và các kí hiệu khác có thêm chữ e để chỉ có liên quan hoặc ở trên không gian đó. Nhắc lại rằng bài toán được phát biểu như sau: cho µ ∈ D, ta tìm ue (µ) ∈ X e sao cho a(ue (µ), v; µ) = f (v; µ), ∀v ∈ X e , (2.1) se (µ) = `(ue (µ); µ). (2.2) và ước lượng Ở đây se : D → R là đầu ra tượng trưng cho quan sát nào đó của trạng thái chứa những thông tin của trạng thái mà ta quan tâm và ` là phiếm hàm phụ thuộc vào µ. 12 Ta giả sử rằng các thành phần của bài toán (2.1) - (2.2) là phụ thuộc affine vào tham số, tức là `(v; µ) = Q` X θ`q (µ)`q (v), ∀v ∈ X e , ∀µ ∈ D, q=1 f (v; µ) = Qf X θfq (µ)f q (v), ∀v ∈ X e , ∀µ ∈ D, (2.3) q=1 và a(w, v; µ) = Qa X θaq (µ)aq (w, v), ∀w, v ∈ X e , ∀µ ∈ D, (2.4) q=1 với Q` , Qf , Qa là các số tương đối nhỏ. Ta ngầm giả định rằng θ`q với 1 ≤ q ≤ Ql , θfq với 1 ≤ q ≤ Qf , và θaq với 1 ≤ q ≤ Qa là các biểu thức đại số đơn giản có thể dễ dàng tính được với độ phức tạp tính toán là O(1). Giả sử thêm rằng bài toán (2.1) - (2.2) có các tính chất sau • Tính phù hợp: `(.; µ) = f (.; µ), ∀µ ∈ D và a là đối xứng. • Tính bức tham số. 2.1.2 Tích vô hướng và chuẩn khác Nhằm mục đích xây dựng các ước lượng về sau ta cần một số chuẩn khác trên X e . Cho a là bức, chúng tôi đưa ra tích vô hướng năng lượng và năng lượng chuẩn (((w, v)))µ = a(w, v; µ), ∀w, v ∈ X e , |||w||| ≡ p a(w, w; µ), ∀w ∈ X e , (2.5) (2.6) tương ứng; chú ý rằng những đại lượng này là phụ thuộc vào tham số. Nhờ có giả định về tính bức và liên tục (2.5) tạo thành một định nghĩa tích vô hướng và (2.6) là chuẩn tương đương với chuẩn trong H 1 (Ω) (3). Chúng tôi có thể xác định tích vô hướng và 13 chuẩn trên X e . Đặc biệt, chúng ta sẽ chọn một tích vô hướng năng lượng và chuẩn phụ thuộc vào một giá trị tham số xác định µ̄ ∈ D (w, v)X e ≡ (((w, v)))µ̄ (= a(w, v; µ̄)), ∀w, v ∈ X e , || w ||X e ≡ |||w|||w̄ (= p a(w, w; µ)), ∀w ∈ X e . Nhắc lại rằng (1.2), chúng ta làm quen với hằng số bức của a trên X e như sau αe (µ) ≡ inf e w∈X α(w, w; µ) > 0, ∀µ ∈ D. kwk2X e (2.7) Tương tự, từ (1.3), hằng số liên tục của a trên X e γ e (µ) ≡ sup sup w∈X e v∈X e 2.2 2.2.1 α(w, w; µ) < +∞, µ ∈ D. kwkX e kvkX e (2.8) Không gian và cơ sở xấp xỉ hữu hạn chiều Rời rạc hóa bài toán Bài toán (2.1) - (2.2) trên không gian X e cần được rời rạc hóa để giải trên máy tính. Ta sử dụng phương pháp phần tử hữu hạn (FEM) để xấp xỉ X e bởi không gian hữu hạn chiều X N . Hãy hình dung, chẳng hạn trong không gian hai chiều và X e là H01 (Ω) và X N là không gian con sinh bởi các phần tử p1 (k) trong đó k là một phép xấp xỉ trên miền Ω và p1 (k) là tập tất cả các đa thức hai biến có bậc mỗi biến không vượt quá 1. Ta kí hiệu ϕN k , 1 ≤ k ≤ N là cơ sở của phương pháp phần tử hữu hạn (còn gọi là các phần tử hữu hạn). Khi đó, ta xấp xỉ mỗi phần tử thuộc X e bởi một tổ hợp tuyến tính v(x) = N X vi ϕ N k i=1 vi là các hằng số. Ta chọn tích vô hướng và chuẩn trên X N bằng cách đơn giản là kế thừa các tích vô hướng và chuẩn trên không gian X e (w, v)X N ≡ (w, v)X e ≡ a(w, v; µ̄), ∀w, v ∈ X N , 14 và ||w||X N ≡ ||w|| Xe ≡ q a(w, w; ¯), ∀w ∈ X N . Chú ý rằng định nghĩa các tích vô hướng và các chuẩn đưa ra là độc lập với N . Từ (1.2) ta giới thiệu hằng số bức của a trên X N như sau αN (µ) ∈ inf w∈X N a(w, w; µ) . ||w||2X N Tương tự, từ (1.3) ta giới thiệu hằng số liên tục của a trên X N γ N (µ) ≡ sup sup w∈X N v∈X N a(w, v; µ) . ||w||X N ||v||X N Ta đòi hỏi họ các không gian con đúng X N đáp ứng điều kiện xấp xỉ max inf ||u(µ) − w||X e → 0 khi N → ∞. µ∈D w∈X N (2.9) Trong (2.9), với bất kì ε > 0 luôn tồn tại N sao cho sai số trong đó thích hợp nhất với u(µ) ∈ X N là nhỏ hơn hoặc bằng ε, ∀µ ∈ D. 2.2.2 Phép chiếu Galerkin Bài toán của chúng ta như sau. Cho µ ∈ D, tìm uN (µ) ∈ X N sao cho a(uN (µ), v; µ) = f (v; µ), ∀v ∈ X N , (2.10) sN (µ) = f (uN ; µ). (2.11) và ước lượng Đây chính là hệ quả của việc áp dụng phép chiếu Galerkin Tiêu chuẩn lên bài toán (2.1) - (2.2). Người ta đã chỉ ra rằng • ||f ||(X N )0 ≤ ||f ||(X e )0 , • αN (µ) ≥ αe (µ), ∀µ ∈ D và vì thế a bức trên X N , • γ N (µ) ≤ γ e (µ), ∀µ ∈ D và vì thế a liên tục trên X N ,