Các phương pháp tính tích phân gần đúng cho hàm số có số biến rất lớn và ứng dụng trong tài chính

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN NGUYỄN TIẾN ĐÀ CÁC PHƯƠNG PHÁP TÍNH TÍCH PHÂN GẦN ĐÚNG CHO HÀM SỐ CÓ SỐ BIẾN RẤT LỚN VÀ ỨNG DỤNG TRONG TÀI CHÍNH Chuyên ngành : TOÁN GIẢI TÍCH Mã số : 60 46 01 02 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: GS TSKH ĐINH DŨNG HÀ NỘI - Năm 2013 Mục lục Lời nói đầu CÁC PHƯƠNG PHÁP TÍNH TÍCH PHÂN CÓ SỐ CHIỀU RẤT LỚN 1.1 Phân rã ANOVA 1.1.1 Phân rã ANOVA cổ điển 1.1.2 Phân hoạch ANOVA có điểm neo 1.2 Phương pháp tính tích phân theo số chiều 1.2.1 Sự chặt cụt rời rạc hóa 1.3 Sai số chi phí 1.4 Xây dựng tiên nghiệm sử dụng không gian hàm có trọng 5 12 14 14 15 18 PHƯƠNG PHÁP TÍNH TÍCH PHÂN TỐI ƯU TRÊN LƯỚI THƯA 20 2.1 Lưới thưa tổng quát 20 2.2 Mối quan hệ phương pháp tính tích phân lưới thưa với phương pháp tính tích phân theo số chiều 21 2.3 Lưới thưa tối ưu không gian có trọng 23 2.4 Tỉ lệ chi phí lợi nhuận 24 2.5 Phân tích chi phí 27 2.6 Phân tích sai số 28 2.7 Phân tích sai số so với chi phí 30 MỘT SỐ ỨNG DỤNG TRONG TÀI CHÍNH VÀ 3.1 Kết số 3.1.1 Xây dựng đường ngẫu nhiên (RW) 3.1.2 Xây dựng cầu Brown (BB) 3.1.3 Xây dựng thành phần chủ yếu (PCA) 3.2 Tùy chọn kiểu Châu Á 3.2.1 Mô hình toán 3.2.2 Số chiều hiệu dụng 3.2.3 Sai số chi phí tích phân 3.3 Trái phiếu lãi suất không 3.3.1 Mô hình toán 3.3.2 Số chiều hiệu dụng KẾT QUẢ SỐ 33 33 34 34 35 36 37 38 39 40 40 42 3.3.3 Sai số chi phí tích phân 3.4 Trái vụ bảo đảm tài sản chấp 3.4.1 Mô hình toán 3.4.2 Số chiều hiệu dụng 3.4.3 Chi phí sai số tích phân 42 43 44 45 45 Kết luận 53 Tài liệu tham khảo 54 Lời nói đầu Lý thuyết xấp xỉ nhánh toán học nói chung giải tích nói riêng Hiện phát triển mạnh mẽ, thâm nhập vào hầu hết ngành toán học ngành khoa học khác hóa học, vật lí chí tài toán Đôi phải làm việc với toán mà việc tính giá trị xác gặp nhiều khó khăn số lí biểu thức toán học cồng kềnh, phức tạp phải tính tích phân có số chiều lớn nhiều lúc người ta xem thảm họa cần phải khắc phục loại trừ hay làm giảm theo nghĩa đó.Với toán việc tính gần giá trị cho sai số tính toán nhỏ cần thiết Điều đáng nói xuất nhiều tích phân có số chiều lớn nhiều mô hình bao gồm từ toán học, vật lí, hóa học đến tài Trong hầu hết trường hợp, tích phân xuất tính toán theo công cụ giải tích thông thường phương pháp số cần phải áp dụng Ở vấn đề tiên thảm họa số chiều cần phải tránh theo nghĩa Thảm họa số chiều thể chỗ chi phí để xấp xỉ tích phân với độ xác ε cho trước phụ thuộc theo hàm mũ vào số chiều toán Nó trở ngại lớn cho phương pháp số truyền thống với toán có số chiều cao Điều nói đến [8] Hơn thảm họ số chiều tìm thấy theo quan điểm định lý phức tạp lý thuyết số Ở đó, thể vài toán tích phân với thuật toán tốt chí không tránh khỏi thảm họa số chiều, toán gọi không khả thi Tuy nhiên nhiều toán ứng dụng lại xuất lớp toán nhỏ khả thi, thêm vào tồn thuật toán phá vỡ thảm họa số chiều, thuật toán ngẫu nhiên Monte Carlo (MC) thuật toán có dạng Mặc dù tốc độ hội tụ thấp phải sử dụng số lượng tương đối lớn điểm đánh giá Sau phương pháp tựa Monte Carlo (QMC) thay cho phương pháp Monte Carlo, phương pháp giành tốc độ hội tụ nhanh Sai số phương pháp đạt ε(n) = O(n−1 (log n)d ) cho hàm dấu tích phân có phương sai bị chặn Với nhiều toán tài chính, chuyên gia lý thuyết số chứng minh phương pháp tựa Monte Carlo hội tụ gần độc lập với số chiều, đồng thời nhanh xác phương pháp Monte Carlo Với hàm đủ trơn, kết tương tự tìm thấy cho phương pháp tính tích phân lưới thưa Một giải thích cho thành công phương pháp dựa phân tích phân rã phương sai gọi tắt ANOVA Trong luận văn, tác giả trình bày lại số phương pháp tính tích phân gần cho hàm số có số chiều lớn dựa theo nội dung báo: Michael Griebel, Markus Holtz, "Dimension - wise integration of high - dimensional functions with applications to finance", Journal of Complexity 26 (2010), pp 455 - 489 Đồng thời dựa theo nội dung báo tác giả đưa vài kết số quan trọng số ứng dụng tài Cụ thể là, phương pháp tính tích phân xây dựng sở chặt cụt rời rạc hóa phân rã ANOVA có điểm neo Những phương pháp thiết lập nhằm khai thác số chiều hiệu dụng thấp hàm dấu tích phân mà phương pháp lưới thưa trường hợp đặc biệt Hơn phương pháp áp dụng theo cách thích nghi theo số chiều thích nghi địa phương Hiệu chúng kiểm tra chuyên gia số học tích phân có số chiều lớn xuất phát từ tài Nội dung luận văn gồm chương Chương Các phương pháp tính tích phân có số chiều lớn Chương nhắc lại hai phân rã ANOVA, phân rã ANOVA cổ điển phân rã ANOVA có điểm neo Qua đưa khái niệm tương ứng số chiều chặt chụt số chiều hiệu dụng cho loại phân rã Sau sử dụng phân rã ANOVA có điểm neo để giới thiệu lớp phương pháp cho việc tính tích phân nhiều biến Chương Phương pháp tính tích phân tối ưu lưới thưa Chương trình bày phương pháp lưới thưa tổng quát, lưới thưa cổ điển lưới thưa có trọng Qua phân tích mối quan hệ yếu tố sai số, chi phí lợi nhuận phương pháp cho việc tính tích phân nhiều biến Chương Một số kết số ứng dụng tài Chương trình bày kết số số ứng dụng tài mô hình tùy chọn kiểu Châu Á, trái phiếu không lãi suất toán CMO ( trái vụ bảo đảm tài sản chấp) Luận văn hoàn thành hướng dẫn tận tình GS.TSKH Đinh Dũng, tác giả xin bày tỏ lòng kính trọng biết ơn sâu sắc đến thầy, thông qua luận văn tác giả xin lời cảm ơn chân thành đến thầy cô hội đồng phản biện đọc đóng góp ý kiến quý báu cho luận văn Một lần tác giả xin gửi lời cảm ơn tới gia đình, bạn bè, đồng nghiệp động viên giúp đỡ nhiều trình hoàn thành luận văn Do thời gian, kinh nghiệm lực nhiều hạn chế nên luận văn chắn không tránh khỏi thiếu sót ý muốn, tác giả mong nhận đóng góp ý kiến phê bình kịp thời thầy cô, bạn bè đồng nghiệp để luận văn hoàn thiện mặt nội dung lẫn hình thức Tác giả xin chân thành cảm ơn! Hà Nội, ngày 28 tháng 10 năm 2013 Tác giả Nguyễn Tiến Đà Chương CÁC PHƯƠNG PHÁP TÍNH TÍCH PHÂN CÓ SỐ CHIỀU RẤT LỚN 1.1 Phân rã ANOVA Trong mục giới thiệu phép phân tích phương sai cổ điển (ANOVA cổ điển) phép phân tích phương sai có điểm neo (ANOVA có điểm neo) hàm nhiều biến f Dựa phân rã này, định nghĩa khái niệm khác số chiều hiệu dụng f Để bắt đầu, cho Ωd ⊆ Rd tập xác định f giả sử dµ(x) = n ∏ dµj (x) kí hiệu tích độ đo định nghĩa j=0 tập Borel Ωd Ở x = (x1 , , xd ) µj với j = 1, 2, , d độ đo Ω Với ∫ dµj (xj ) = (1) Ω[0,1] Đồng thời kí hiệu V (d) không gian Hilbert gồm tất hàm f : Ωd → R Trong tích vô hướng xác định ∫ (f, g) = f (x)g(x)dµ(x) Ωd Với tập u ⊆ D, D = {1, , d} tập số, độ đo µ cảm sinh phép chiếu Pu : V (d) → V (|u|) ∫ cho Pu f (xu ) = f (x)dµD\u (x) Ωd−|u| Ở ta kí hiệu xu kí hiệu véc tơ có số chiều |u| gồm thành phần x mà số thuộc u dµD\u (x) := ∏ j̸∈u dµj (xj ) Phép chiếu định nghĩa phân rã f ∈ V (d) thành tổng hữu hạn f (x1 , , xd ) = f0 + d ∑ d ∑ fj1 (xj1 ) + j1 fj1 ,j2 (xj1 , xj2 ) + + fj1 , ,jd (xj1 , , xjd ) j1 k2 > Liên quan tới số chiều d = 512, nhìn thấy hình 2b, |∆k | < 10−10 k2 > Điều phản ánh số chiều chặt cụt toán Vasicek với việc xây dựng cầu Brown Những kết thể hình 3.2 biểu thị cho ví dụ đặc biệt có tập số I nhỏ đủ để chụp tất số k tương ứng với phân bố đáng kể giá trị tích phân Điều giải thích độ xác cao lại đạt với chi phí phương pháp SGH Hình 3.2: Lát mỏng hai chiều qua tập số I với ε = 10−3 cho toán Vasicek với xây dựng cầu Brownian 3.4 Trái vụ bảo đảm tài sản chấp Cuối xem xét toán bổn phận trả nợ việc vay chấp hay gọi toán CMO (CMO từ viết tắt cụm từ "Collaralized Mortgage Obligations") Bài toán mô hình diễn đạt chi tiết [10],[11] có thêm vài tài liệu tham khảo khác Nó dùng để nghiên cứu hiệu phương pháp tựa Monte Carlo (QMC) Nó xem xét [12] để chứng minh hiệu phương pháp lưới thưa (SG) 43 3.4.1 Mô hình toán Bây xem xét toán CMO [10] Cụ thể xét đảm bảo chấp tài sản với kì hạn d tháng Ta kí hiệu tk = ∆t.k với k = 1, 2, , d, ∆t kí hiệu cho chiều dài chu kì tháng Khi giá trị tổng tất tiền nợ cho g(z) := d ∑ u k mk , k=1 uk := k−1 ∏ (1 + ij )−1 j=0 yếu tố giảm giá cho tháng k, phụ thuộc vào tỉ lệ lãi suất ij Tỉ lệ lãi suất ik cho tháng k tính công thức ik := K0k eσ(z1 +z2 + +zk ) i0 Trong zj với j = 1, 2, , d số ngẫu nhiên phân phối chuẩn Ở i0 tỉ lệ lãi suất thời điểm ban đầu, σ số dương K0 := e−σ /2 cho E[ik ] = i0 Tiền nợ hàng tháng thời điểm tk cho mk := crk ((1 − ωk ) + ωk ck ), c kí hiệu tiền trả hàng tháng Ngoài ck := d−k ∑ (1 + i0 )−j , j=0 rk := k−1 ∏ (1 − ωj ), j=1 ωk :=K1 + K2 arctan(K3 ik + K4 ) Ở rk kí hiệu cho phần lại chấp, phụ thuộc vào phần trả trước ωk tháng k Hơn K1 , K2 , K3 , K4 số tỉ lệ toán trước ωk Giá trị kì vọng tổng tất tiện nợ viết ∫ tích phân d chiều P V := g(z)φd (z)dz, Rd 44 T φd (z) := e−x x/2 (2π)d/2 hàm trọng Gauss Tiếp theo tập trung vào hiệu tính toán số tích phân Theo chuyên gia số học, sử dụng tham số i0 = 0.007, c = 1, K1 = 0.01, K2 = −0.005, K3 = 10, K4 = 0.5 σ = 0.0004 Với d = 256 giành PV= 119.21588257 3.4.2 Số chiều hiệu dụng Bây sử dụng tham số [12], đồng thời xem xét trường hợp d = 256 Đầu tiên xem xét số chiều chồng chất trường hơp có điểm neo cho toán CMO Chúng ta lại giành ds với α ∈ [0.9, 0.9999] tất cách xây dựng đường, nhìn vào bảng 3.4b Số chiều chặt cụt phương pháp trình bày bảng 3.3b Nó cho thấy đường xây dựng có ảnh hưởng nhỏ vào số chiều chặt cụt trường hợp có điểm neo, nhiên lợi phương pháp BB, PCA, LT so sánh với RW không rõ ràng cho toán CMO Với α = 0.9 có dt = 123, trường hợp RW LT Số chiều chặt cụt cắt giảm tới dt = 18, dt = 13 tương ứng cho BB, PCA Tuy nhiên với yêu cầu độ xác cao hơn, ví dụ α 0.99 tất cách xậy dựng cắt giảm không đáng kể chí không Chú ý với toán CMO, số chiều hiệu dụng trường hợp cổ điển rõ ràng khác với số chiều chặt cụt trường hợp có điểm neo Ở BB, PCA, LT đưa việc cắt giảm số chiều đáng kể Thậm chí LT cắt giảm số chiều chồng chất tới trường hợp cổ điển Hình 3.3: Số chiều chặt cụt toán Vasicek CMO 3.4.3 Chi phí sai số tích phân Tiếp theo kết số tương ứng minh họa hình 3.7 Chúng ta quan sát phương pháp tựa Monte Carlo (QMC) hội tụ nhanh 45 hơn, dao động so với phương pháp Monte Carlo chuyển từ RW tới BB, PCA LT Phương pháp SGP thực tương tự phương pháp tựa Monte Carlo (QMC) trường hợp BB, PCA tồi tệ trường hợp của RW LT Phương pháp hỗn hợp CUHRE/QMC, COW, CPW CAD giành kết tốt kết hợp với PCA Trong trường hợp chúng tốt phương pháp QMC phương pháp SGP Cuối nhìn thấy đến phương pháp SGH kết hợp với BB PCA phương pháp hiệu cho toán CMO Nó đạt tốc độ hội tụ cao kết xác Với 104 điểm đánh giá hàm, phương pháp SGH giành sai số tương đối nhỏ khoảng 100 lần so với sai số tương đối phương pháp QMC 1-α RW 1e - 1 1e - 1 1e - 1e - (a) Vasicek BB PCA LT 1 1 1 1 1 problem (d=512) 1-α RW BB PCA LT 1e - 1 1 1e - 2 2 1e - 2 2 1e - 2 2 (b) CMO problem (d=256) Bảng 3.4 Số chiều chồng chất trường hợp cổ điển Nhắc lại LT không cải thiện hội tụ phương pháp tính tích phân theo số chiều so sánh với RW Sự nhận xét giải thích số chiều hiệu dụng trường hợp cổ điển LT giành kết tối ưu dt = cho toán CMO Tuy nhiên quan sát cho thấy LT không đưa việc cắt giảm số chiều chặt cụt trường hợp có điểm neo Điều nói lên hiệu phương pháp tính tích phân theo số chiều phụ thuộc vào số chiều chặt cụt trường hợp có điểm neo Chú ý phương pháp QMC hội tụ nhanh dao động với LT RW Điều cho thấy dáng điệu hội tụ phương pháp QMC liên quan tới số chiều hiệu dụng dt trường hợp cổ điển có điểm neo Cuối ý số chiều chặt cụt dt trường hợp có điểm neo giải thích ảnh hưởng cách thức xây dựng không hiệu cao phương pháp SGH dt cao cho toán Tốc độ hội tụ nhanh giải thích số 46 chiều chồng chất thấp ds độ trơn hàm dấu tích phân Hình 3.4: Dáng điệu hội tụ phương pháp khác cho toán Asian với số chiều d=16 47 Hình 3.5: Dáng điệu hội tụ phương pháp số khác cho toán Asian 100 với số chiều d=16 48 Hình 3.6: Dáng điệu hội tụ phương pháp khác cho toán Vasicek số chiều d=512 49 Hình 3.7: Dáng điệu hội tụ phương pháp khác cho toán CMO với số chiều d=256 50 Kết luận Nội dung luân văn chủ yếu dựa vào kết đề cập báo: Michael Griebel, Markus Holtz , "Dimension - wise integration of high - dimensional functions with applications to finance", Journal of Complexity 26 (2010), pp 455 - 489 Ở trình bày lại lớp phương pháp cho việc tính toán tích phân có số chiều cao, thông qua việc giới thiệu số khái niệm số chiều hiệu dụng trường hợp có điểm neo Bên cạnh giải thích phương pháp lưới thưa kết hợp chặt cụt phân rã ANOVA có điểm neo kết hợp với rời rạc hóa chuỗi cho phép làm cân sai số mô hình hóa sai số rời rạc hóa theo cách tới ưu Cũng dựa vào nội dung báo đề cập đến toán ứng dụng từ tài mà dẫn đến tích phân hàm trơn với số chiều lên tới 512 Với toán mô hình này, phương pháp lưới thưa thích nghi theo số chiều dựa công thức Gauss - Hermite hay gọi phương pháp SGH mang lại hiệu tốt Phương pháp làm lợi từ số chiều hiệu dụng thấp tích phân làm mịn lưới thưa thích nghi theo số chiều khai thác cách tối đa độ trơn của hàm dấu tích phân tránh việc biến đổi đặc biệt hình cầu đơn vị Hiển nhiên kết mở rộng theo nhiều hướng khác Ví dụ thú vị để xác định lớp hàm cho số chiều hiệu dụng trường hợp có điểm neo trùng với số chiều hiệu dụng trường hợp cổ điển Chúng ta toán giá cổ phiếu trái phiếu không lãi suất thuộc lớp hàm toán CMO Một điều thú vị làm người quan tâm nên chọn điểm neo tốt mang lại hiệu cao Những lĩnh vực khác nghiên cứu tương lai, việc cải thiện phương pháp tính tích phân theo số chiều dạng lưới thưa ví dụ cân sai số mô hình hóa rời rạc hóa Với ứng dụng từ tài phương pháp hỗn hợp CUHRE/QMC cải thiện thêm thay sử dụng phương pháp CUHRE 51 phương pháp thích nghi địa phương khác mà xử lý trực tiếp tích phân Rd cho tránh việc biến đổi đặc biệt hình cầu đơn vị Chúng ta hi vọng phương pháp tương tự chí cho kết tốt giành với phương pháp tính tích phân theo số chiều dựa lưới thưa Cuối thấy mức độ hầu hết phương pháp kết mang lại không bị hạn chế tài mà sử dụng theo lĩnh vực ứng dụng khác vật lý hóa học 52 Tài liệu tham khảo [1] F.J Hikernell (1998), "A generalized discepancy quadrature error bound", Math Comput, 67, pp 299 - 322 [2] H Rabitz, O Alis (1999), "General foundation of high - dimension model representaions", Journal of mathematical chemistry, 25, pp 197 - 233" [3] Hegland (2003), "Adaptive spare grids", ANZIAM J 44, C335 - C353 [4] G Wasilkowski , H Woz’niakowski (1999), "Weighted tensor product algorithms for linear multivariate problems", J complexity 15, pp 402 - 447 [5] Imai,Tan (2006), "A general dimension reduction technique for derivative pricing", Journal of computaion finance 10 (2), pp 129 -155 [6] J Bersten, T Espelid, A Genz (1991), "Algorithms 698: DCUHRE - An adaptive multidimensional integration routine for a vector of integrals", ACM Transactions on mathematical Software, 17, pp 452 - 456 [7] Michael Griebel, Markus Holtz (2010), "Dimension - wise integration of high dimensional functions with applications to finance", Journal of Complexity 26, pp 455 - 489 [8] M Griebel (2006), "Sparse grids and related approximation schemes for higher dimensional problems", in: L Pardo, A Pinkus, E.Suli, M Todd (Eds), Foundtaions of Computational Mathematics (FoCM05), Santander, Cambridge University Press, pp 106 -161 [9] Ninomiya, S Tezuka (1996), "Toward real - time pricing of complex financial derivatives", Applied Mathematical Finance 3, pp - 20 [10] R Caflisch, W.Morokoff, A Owen (1997), "Valuation of Mortgage backed securities using Brownian bridges to reduce effective dimention", J Comp Finance, 1, pp 27 - 46 53 [11] S Paskov, J Traub (1995), "Faster valuation of finance derivatives", J Portfolio Management 22, pp 113-120 [12] T Gerstner, M Griebel (1998), "Numerical integration using spare grid", Numerical Algorithms 18, pp 209 -232 [13] T Gerstner, M Griebel (2003), "Dimension - adaptive tensor - product quadrature", computing 71 (1), pp 65 -87 [14] T Hahn (2005), "Cuba - a library for mutidimensional numerical integration", Computer Physics Comumunications 168, pp.78 -95 [15] X.Wang, K - T Fang (2003), "The effective dimension and quasi - Monte Carlo integration", J complexity 19 (2), pp 101 -124 [16] Wang (2006), "On the effects of dimension reduction techniques on some high - dimensional problems in finance", Operator Reseach 54 (6), pp 1063 -1078 54 [...]... thuộc vào chuẩn ||f ||1,γ của hàm dưới dấu tích phân Chuẩn này có thể tăng nhanh theo số mũ so với việc tăng của số chiều d là tác nhân của bài toán có số chiều cao 32 Chương 3 MỘT SỐ ỨNG DỤNG TRONG TÀI CHÍNH VÀ KẾT QUẢ SỐ 3.1 Kết quả số Trong mục này chúng ta sử dụng một số tích phân có số biến rất lớn xuất phát từ tài chính để kiểm tra hiệu quả của phương pháp lưới thưa và phương pháp tính tích phân có. .. sử dụng phân rã ANOVA có điểm neo để định nghĩa một lớp mới các phương pháp cho việc tính toán tích phân với số chiều cao Ta kí hiệu ∫ f (x)dx (14) f (z)φd (z)dz, (15) If := [0,1]d và ∫ Iφ f := Rd trong đó φd là hàm trọng Gauss Hai loại này là điển hình trong ứng dụng cho tính tích phân có số chiều cao 1.2.1 Sự chặt cụt và rời rạc hóa Tiếp theo chúng ta phát triển lớp mới các phương pháp tính tích phân. .. sai số nhỏ nếu α gần 1 Các điểm tựa ngẫu nhiên là được phân bố đều theo số chiều lớn đến mức mà ta có thể hi vọng rằng Ifdt là được xấp xỉ tốt khi số chiều chặt cụt dt nhỏ Đồng thời đánh giá |If − Iftr | ≤ √ 1 − ασ(f ) (10) giải thích một phần thành công của phương pháp lưới thưa cho tích phân có số chiều cao với những hàm số có số chiều hiệu dụng thấp vì những phương pháp này có thể tính toán Iftr rất. .. điển là rất hữu ích để phân tích tầm quan trọng của các số hạng với số chiều khác nhau và tương tác giữa chúng nhưng nó không được sử dụng như một công cụ cho việc thiết kế để xây dựng lược đồ tích phân vì số hạng hằng trong phân hoạch cổ điển yêu cầu phải tính tích phân Phân rã ANOVA có điểm neo lại có thuận lợi là số hạng con của nó là khá rẻ để tính toán vì rằng ta chỉ cần thay việc tính tích phân. .. ( 12 ) và công thức tính tích phân sai phân ∆k := Umk − Umk−1 (26) với Um0 = 0 và mọi k ≥ 1 Bây giờ cho f : [0, 1]d → R là hàm nhiều biến Khi đó tích phân d - chiều có thể được viết thành tổng thu gọn vô hạn If = ∑ ∆k , f (27) k∈Nd ở đây k ∈ Nd là kí hiệu tích các chỉ số với kj > 0 và ∆k f := (∆k1 ⊗ ∆k2 ⊗ ∆kd )f (28) Một lớp đặc biệt các phương pháp tính tích phân cho xấp xỉ của If là dựa vào sự chặt... tính tích phân đặc biệt Qu sao cho phương pháp pháp (20) tương ứng lớp phương pháp lưới thưa tổng quát Trong mục này chúng ta sẽ xác định tập chỉ số I, tập sẽ làm cân bằng sai số mô hình hóa và sai số rời rạc hóa theo một cách tối ưu cho hàm dưới dấu tích phân được lấy từ không gian hàm Sobolev Bây giờ chúng ta sẽ xem xét công thức tính tích phân một biến Umk trong (25) được cho bởi công thức hình thang,... u∈S là số các điểm đánh giá của hàm f , ở đây nu là kí hiệu của các số điểm đánh giá hàm của Qu 1.3 Sai số và chi phí Đầu tiên chúng ta xem xét trường hợp của phương pháp tính tích phân tùy ý Qu Đầu tiên chúng ta có đánh giá sai số |If − AS f | = | ∑ u⊆D Ifu − ∑ qu | ∑ u∈S u∈S 15 |Ifu − qu | + ∑ u̸∈S |Ifu | (21) Điều này chứng tỏ rằng sai số của (20) phụ thuộc vào phương pháp tính tích phân Qu và sự... số S Ở đây số hạng thứ hai diễn tả sai số mô hình hóa theo nghĩa chặt cụt, tuy nhiên số hạng thứ nhất mô tả sai số rời rạc hóa nghĩa là nó dựa vào sai số của từng số hạng trong phân rã ANOVA có điểm neo Mục đích của chúng ta trong phần tiếp theo là làm cân bằng chi phí và độ chính xác bằng việc liên hệ chi phí của phương pháp tính tích phân Qu với tầm quan trọng của số hạng fu trong phân rã ANOVA có. .. dụ như trong phương pháp lưới thưa cổ điển người ta chọn tập chỉ số { I = k ∈ Nd : |k|1 ở đây |k|1 = d ∑ } l+d−1 , (31) kj , hoặc đối với phương pháp tích thì tập I lại có dạng j=1 { I = k ∈ Nd : |k|∞ } l , (32) ở đây |k|∞ := max {kj : j = 1, 2 , d} 2.2 Mối quan hệ giữa phương pháp tính tích phân trên lưới thưa với phương pháp tính tích phân theo số chiều Có một mối quan hệ gần gũi giữa phương pháp trên... tập con tương ứng với tập chỉ số của phương pháp lưới thưa cổ điển, thông thường thì có 2d tập con được sử dụng, với mỗi tập con ứng với số hạng con của phân rã ANOVA Hình 2.1: Các tập chỉ số tối ưu Il,γ trên mức l=7 2.5 Phân tích chi phí Trong mục tiếp theo chúng ta sử dụng n(d, l, γ) để kí hiệu số các điểm đánh giá hàm của phương pháp SGl,γ Bổ đề 2.3 (Chi phí của lưới thưa cổ điển) Cho mi ) ( l+d−2 ... TÀI CHÍNH VÀ KẾT QUẢ SỐ 3.1 Kết số Trong mục sử dụng số tích phân có số biến lớn xuất phát từ tài để kiểm tra hiệu phương pháp lưới thưa phương pháp tính tích phân có số chiều lớn khác Các toán... giải thích cho thành công phương pháp dựa phân tích phân rã phương sai gọi tắt ANOVA Trong luận văn, tác giả trình bày lại số phương pháp tính tích phân gần cho hàm số có số chiều lớn dựa theo... Chương CÁC PHƯƠNG PHÁP TÍNH TÍCH PHÂN CÓ SỐ CHIỀU RẤT LỚN 1.1 Phân rã ANOVA Trong mục giới thiệu phép phân tích phương sai cổ điển (ANOVA cổ điển) phép phân tích phương sai có điểm neo (ANOVA có