Ứng dụng wavelets trong nội suy dị thường độ cao

Nghiên cứu ỨNG DỤNG WAVELETS TRONG NỘI SUY DỊ THƯỜNG Đ ộ CAO LƯƠNG BẢO BÌNH Trường Đại học Bách Khoa - ĐHQG TPHCM Tóm tắt: Bài bảo giới thiệu cách tiếp cận waveỉets ứng dụng việc nội suy đại lượng trắc địa, phương pháp cỏ thể tỉnh tốn nội suy cho vị trí từ giá trị đầu vào với phân bố không gian Sau phần tảng lý thuyết hai tỉnh toán thử nghiệm sử dụng giả trị từ mơ hình EGM 2008, cho lưới đồng góc toàn cầu cho mạng lưới địa phương với vị trí phân bổ khơng đồi, nhằm minh họa xác thực tỉnh đắn phương pháp Một cách tóm tẳt, vvavelets gồm hai bước: “chia ” tín hiệu đẩu vào thành mức khác (tương ứng với độ bậc điầi hịa cầu), sau tổng hợp lại thành tín hiệu đầu (ở vị trí cần giả trị nội suv) Ẩn bên wavelets nhân tái tạo sử dung hàm sở cầu đa thức Legendre, điều giúp cho quv trình giữ tính chất điều hịa cầu cho giá trị nội suy, ưu điểm bật so với phương pháp nội suy quen thuộc khác Tính đắn độ chỉnh xác phương pháp minh chứng thơng qua hai tính tốn thử nghiệm Ở mức độ tồn cầu, lưới đồng góc 6° giá trị dị thường độ cao đến độ bậc 10 dùng để nội suy cho vị trị cho kết chênh lệch mức l ữ 12 m Ở mức độ khu vực, 5796 giá trị dị thường độ cao (độ bậc từ 32 đến 900) châu Âu “tự nội suy” với kết chênh lệch cm, với độ nhiễu ngẫu nhiên thêm vào trước Điều với độ bậc điầi hịa cao vị trí điểm đo khơng phân bố khó khăn mang tỉnh thực tế chọn cho tỉnh toán thử nghiệm thứ hai wavelets cho kết tốt Từ khóa: waveỉets, nhân tái tạo (reproducing kemel), điều hòa cầu, nội suy Đặt vấn đề Trọng lực đại lượng liên quan dị thường độ cao, dị thường trọng lực, độ lệch dây dọi thường cung cấp từ mơ hình (tồn cầu khu vực) dạng lưới điểm với kích thước lưới định, chẳng hạn 30 phút, phút, 2.5 phút, phút Để có giá trị vị trí cụ thể mà quan tâm, người dùng phải tiến hành nội suy từ giá trị mắt lưới Có nhiều phương pháp nội suy quen thuộc tuyến tính, đa thức, spline, kriging, collocation Tuy nhiên, phương pháp khơng bảo đảm việc giữ tính chất điều hịa (harmonic) đại lượng nội suy, mà nguyên tắc tính chất chung đại lượng trắc địa kể Ngày nhận bài: 1/5/2022, ngày chuyển phản biện: 5/5/2022, ngày chấp nhận phản biện: 9/5/2022, ngày chấp nhận đăng: 28/5/2022 TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 52-6/2022 Một cách tiếp cận khác cho vấn đề waveỉets.Víộc íag4p v « ÌĩflÌlÉ k fea a g trắc địã nói đến số tài liệu, chẳng hạn [1] vả P Ị Ĩ Ị f d lM f ilic ầ D quy trình tính tốn cho tốn cụ thể nội suy đại lượng b íc địa BM fcio nậy giới thiệu công thức tỉnh toán cụ thể để áp dụng nguyên tắc cùa wavelets, thông qaa nhân tảj tạo đa thức Legendre giúp bảo tồn tính hài hịa, để nội suy đại lượng trắc địa thường dùng dị thường độ cao Kèm theo hai thử nghiệm thiết kế phạm vi toàn cầu vả khu vực để kiểm chửng tính đắn đánh giá độ xác phương pháp tính tốn Cơ sở lý thuyết phương pháp Ỷ tưởng wavelets việc chia liệu đầu vào thành tín hiệu chi tiết mịn thông qua lọc thông thấp (low-pass) Áp dụng cho toán nội suy dị thường độ cao (hoặc liệu trọng lực nói chung), phương pháp gồm hai bước chính: - Phân tích: chia tách tín hiệu nguồn dừ liệu đầu vào thành nhiều cấp tín hiệu chi tiết; - Tổng hợp: kết hợp tín hiệu chi tiết để tái tạo liệu đầu Nhân tâi tạo - ♦ H ệsổổỉệ - ► 5SSU 2o - ► Tin Wệu chi tiét t Tín hiệu đầu vào Tín hiệu đằu * Tổng hợp Hình ỉ: Sơ đồ tính tốn Quy trình minh họa qua sơ đồ hình 1; đó, quan ừọng nhân tái tạo (reproducing kemel) giúp thiết lập cầu nối hai vị trí (gốc nội suy) hàm wavelet cầu lọc tín hiệu thành cấp chi tiết riêng biệt Cơ sở toán học phương pháp giải thích ngắn gọn bên Gọi Ũ.R hình càu bán kính R L 2(Cìr ) khơng gian tất hàm thực khả tích bậc F ỉìR, điều hòa cầu Yn,m(ệ) độ bậc nm tạo thành sở trực giao hoàn chỉnh L 2(Q.r) Hàm F £ (Q.r) c ó thể biểu diễn cách thơng qua chuỗi Fourier F© = £ ” với ệ I L - A ™ 1'” ” ® (1) £ ÍỊ, Fn,m hệ số Stokes tính biến đổi Fourier cầu Cho (nhiễu) trọng lực F mặt cầu Ũ.R, giá trị tiếp nối hướng lên (upward continuation) tính TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 52-6/2022 Nghiên cứu F (* )= £ " Z ” =_n Fn>m//n,m(x) (2) ữong đó, hàm Wn,mw ( ) n+1^ , m© (3) điều hịa cầu khối (solid spherical harmonics) Mối quan hệ điều hòa cầu đa thức Legendre cho [3] (4) W O - W T ) = ^ P n f t r Tl) ệ vector đon vị theo phương hướng tâm Ũ.R Phương trình (4) tảng cho việc tính wavelets cầu (xấp xỉ đến độ rtmax định) (5) với x = rệE CÌRext (iĨÌRext khơng gian bên bao gồm mặt cầu Q r), Xk = Rệk E CÌR, B„ hệ số Legendre phản ánh tính phổ tín hiệu Như vậy, tín hiệu F(x) ban đầu đại diện wavelets cầu theo phương trình (6) F ( * ) = I Ĩ =1c*B(x,xk) (6) với N số lượng vị trí X Ck hệ số đóng vai trị tương tự hệ số Stokes Fn,m cách tiếp cận điều hòa cầu Vấn đề quan trọng nhân tái tạo, giới thiệu [4] Kr^ x , x t ^ ỵ ^ 2^ ^ f n Bnp , ^ ỉk ) (7) thỏa mãn điều kiện [5] F(x) = (Krep*F )(x) (8) Viết lại tích chập (convolution) K rep * F dạng chuỗi khai triển hàm sở cầu Krep, ta thay phương trình (6) phương trình (9) sử dụng nhân tái tạo (*rep * F)(x) = ỵ N k=1 d kKrep{x,Xk) = krep(x )Td (9) với krep d vector hàm nhân tái tạo hệ số tỉ lệ Phương trình (8) cho thấy mối liên hệ hàm sờ cầu tích chập, cơng cụ lọc Phương trình (10) giới thiệu hàm tỉ lệ cầu - _ v n;'ma* 2n+l m n+1 ^ n r,T , \ $ /(* ,* * ) - , n=ũ 4TCR2 ( r ) 0j:nPn{ị ) ( 10) Cấp ỳ, với hệ số Legendre ĩj;„ = n >ô/ã iu ny cú ngha hm ti lệ hoạt động lọc thông thấp (low-pass) Từ đó, tín hiệu F;+1(x) = ( j+i * F)(x) TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN Đ SỐ 52-6/2022 ( 11) (12) (13) chứa thành phần Fj+i(x) bị thiểu Fj(x) Trong hàm tỉ lệ cầu hoạt động lọc thơng thấp hàm wavelet cầu (14) xem lọc thông dải (band-pass) xác định hệ sổ Legenđre Fy,n = Gj (x) + AFJ+1 (x) (16) bao gồm phiên mịn (ở cấp / ) , tín hiệu chi tiết cấp từ/ đến (cấp tối đa) J, phần dư AFJ+1(x) = F(x) - FJ+í(x) = (A J+1 * F )(x) (17) Tóm lại, theo tên gọi wavelets, ý tưởng phương pháp chia nhỏ tín hiệu đầu vào theo dải phổ Ngoài ra, yếu tố quan trọng khác nhân tái tạo giúp “lan truyền” tín hiệu đến vị trí cần tính tốn Và cuối cùng, tín hiệu cấp chi tiết khác tổng họp thành tín hiệu đầu Phần áp dụng sở lý thuyết vào tính tốn nội suy giá trị dị thường độ cao đế đánh giá tính hiệu độ xác cách tiếp cận Kết thử nghiệm thảo luận Đe kiểm chứng tính đắn độ xác tính tốn nội suy phương pháp, hai thử nghiệm thiết kế với ý tưởng tổng quát sau: - Thử nghiệm thiết kế đơn giản nhằm mục đích kiểm chứng minh họa tính chất lọc dải phổ wavelets Dị thường độ cao lưới đồng góc tồn cầu dùng để “nội suy” Dù khơng mang tính thực tế việc “tự nội suy” giúp loại bỏ hoàn toàn sai số liệu đảnh giá kết Dải phổ cùa liệu hạn chế mức thấp (tương ứng với độ bậc điều hòa 10) để đơn giản hóa tính tốn dễ dàng minh họa hình ảnh - Thử nghiệm giúp đánh giá độ xác phương pháp áp dụng vào tính tốn nội suy dị thường độ cao điều kiện gần với thực tế hơn: từ vị trí điểm liệu đầu vào phân bố không nội suy lưới đồng góc khu vực, giá trị dị thường độ cao tính đến độ bậc điều hòa cao bị thêm vào nhiễu giả ngẫu nhiên đóng vai trị sai sổ TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 52-6/2022 Nghiên cứu 3.1 Tái tạo dị thường độ cao độ bậc điều hịa thấp cho lưới đằng góc tồn cầu Đầu tiên, để kiểm tra tính đắn minh họa trực quan cho tiến trình wavelets, mơ mang tính đon gian hóa tiến hành sau: số liệu dị thường độ cao độ bậc điều hòa thấp phân bố lưới đồng góc tồn cầu phân tích tổng họp lại theo quy trình lý thuyết bên trên, hay nói cách khác “tự nội suy” số liệu cho vị trí cũ đê so sánh kết nhận với số liệu ban đầu, mà khác biệt sai số phương pháp DỊ thường độ cao tính từ mơ hình EGM 2008 [6], sử dụng hệ số điều hịa đến độ bậc 10, tính chưcmg trình GeoH [7] Giá trị dị thường độ cao tính cho mạng lưới đồng góc tồn cầu với kích thước lưới 6° X 6° bao gồm 30 X 60 = 1800 điểm trải dài từ 87° vĩ Bấc đến 87° vĩ Nam từ 0° kinh đến 354° kinh (6° kinh Tây) Việc chi sử dụng hệ số điều hòa bậc thấp để giảm số cấp tín hiệu chi tiết, thuận lợi cho việc minh họa hình ảnh Cụ thể với Umax = 10, áp dụng dyadic wavelets, chi cần cấp tín hiệu đủ bao phủ tồn dải phố: F2 tương úng với độ bậc từ đến 3, G2 tương ứng với độ bậc từ đến 7, G3 chứa tín hiệu với độ bậc từ đến 10 Áp dụng bước tính tốn thể hình vào số liệu này, tính hệ số ti lệ, phiên tín hiệu mịn F , cấp tín hiệu chi tiết G2 G , cuối max 10) phần dư AF khơng cần phải xét đến F2 [m] 50, -50 total [m] — 50 »♦ -50 50 •0 -50 -50 y & -m,- % ■ % L 200 20 m| -20 -50 10 g M 200 reterence [m] 50 * % 200 d:fference [m] G3 [m] 5Ơ § % —5 200 G2 [m] ạ> 50 5oj - - -50 X -12 10 õ 50 -50: -10 200 longitude -50 -2 200 longitude Hình 2: từ xuống dưói, bên trái: tín hiệu mịn, tín hiệu chi tiết cấp cấp 3; bên phải: tin hiệu tong họp đầu ra, giá trị tham khảo, khác biệt chúng TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BAN ĐỒ SỐ 52-6/2022 Nghiên cứu Hình thể tín hiệu F2, G2, G3 ỡ cột bên trái, từ xuống dưới, tín hiệu tổng (giá trị tái tạo cần tính) bên phải, ngồi cịn có giá trị ban đầu (ở bên phải) khác biệt giá trị (dưới bên phải), thể sai số tính tốn phưomg pháp Ta thấy khác biệt tín hiệu ban đầu tín hiệu tái tạo áp dụng wavelets chi mức 10'12m, hồn tồn bỏ qua mà không sợ phạm sai số Điều chứng tỏ tính đắn phương pháp quy trình tính tốn giới thiệu bên Ngồi ra, hình thể cấp tín hiệu F 2, G2, G3 minh họa rõ cho tính chất hoạt động lọc wavelets với mức chi tiết ngày cao cấp tín hiệu cao 3.2 Nội suy dị thường độ cao độ bậc điều hòa cao từ mạng lưới điểm gổc khu vực phân bố không Nếu thử nghiệm (mục 3.1) dùng liệu đơn giản, vốn thiết kế với mục đích khẳng định tính đắn phương pháp thỉ thử nghiệm này, yếu tổ mang tính thực tế xét đến Bộ liệu đầu vào dị thường độ cao 5796 vị trí phân bố khơng theo quy luật, vị trí đo dị thường trọng lực thực tế châu Âu Từ giá trị vị trí ta nội suy giá trị cho lưới đồng góc bao phủ khu vực, trải dài từ 46° vĩ Bắc đến 49.5° vĩ Bắc từ 9° kinh Đơng đến 18° kinh Đơng, kích thước lưới 0,1° X 0,1°, bao gồm 36 X 91 = 3276 điểm Giá trị nội suy so sánh với giá trị tham khảo để đánh giá tính xác phương pháp nội suy Đe có sở so sánh giá trị dị thường độ cao 5796 điểm ban đầu giá trị dị thường độ cao tham khảo 3276 điểm nội suy tính từ nguồn hệ sổ điều hịa đến độ bậc 900 mơ hình EGM 2008 Việc sử dụng đến độ bậc 900 để cân mong muốn giảm thiểu việc tính tốn q nhiều cấp tín hiệu bảo đảm tính chi tiết đủ cao Một lưu ý quan trọng khác thành phần bước sóng dài (lớn kích thước khu vực) cần phải khử trước tính tốn, cụ thể hệ số điều hịa nhỏ 32 Ngồi ra, để gần với điều kiện khó khăn thực tế hon, tập giá trị nhiễu giả ngẫu nhiên với độ lệch chuẩn cm thêm vào giá trị dị thường độ cao 5796 điểm ban đầu, đóng vai trò sai số ngẫu nhiên Áp dụng quy trình tính tốn thử nghiệm 1, cuối ta thu tín hiệu tổng giá trị nội suy 3276 điểm, thể hình bên Bên hình giá trị tham khảo 3276 điểm Các cấp tín hiệu chi tiết sổ lượng nhiều việc minh họa tính chi tiết cho cấp tín hiệu thể thử nghiệm (hình 2) Có hai lưu ý hình 3: thứ nhất, hệ trục tọa độ theo thử tự hàng (36 hàng từ xuống dưới) cột (91 cột từ trái qua phải) ô lưới; thứ hai, giá tri dị thường độ cao sau khử thành phần bước sóng dài, nên nhỏ so với giá trị đầy đủ Tuy nhiên, điều không ảnh hưởng đến việc so sánh khác biệt giá trị tinh toán nội suy giá trị tham khảo hai bị khử thành phần bước sóng dài đến độ bậc 32 TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 52-6/2022 Nghiên cửu The reterence [m] Hình 3: xng dưới: giả trị tham khảo giá trị nội suy (có áp dụng chỉnh quy hóa quasi-optimality) Hệ trục tọa độ theo thứ tự hàng cột lưới So sánh giá trị tham khảo (hình 3, trên) giá trị nội suy (hỉnh 3, dưới) cách trực quan ta thấy giống giá trị (màu sắc) kiểu mẫu Sự khác biệt thấy mắt thường khu vực nhỏ xíu sát hình Cịn số độ lệch quân phưong (RMS) 3276 điểm cm, mức độ nhiễu giả ngẫu nhiên thêm vào sổ liệu ban đầu trước đó; hay nói cách khác, sai số cuối mức độ với sai số dừ liệu, nghĩa quy trình tính tốn nội suy xác Xin nói thêm vấn đề tính tốn xuất thừ nghiệm này, đỏ hệ phương trình chuẩn (khi giải hệ số tỷ lệ) (KrepTKrev)á - KrepTf = (18) bị rơi vào tình trạng giả định yếu (ill-conditioned) mà để giải giải pháp quy hóa (regularization) thích hợp áp dụng gọi quasi-optimality [8] v ấ n đề quy hóa khơng trình bày vượt ngồi phạm vi báo Kết luận Với hai thử nghiệm thiết kế cho hai bối cành từ đơn giản đến phức tạp, từ toàn cầu đến khu vực, kết cho thấy quy trình tính tốn dựa cách tiếp cận wavelets giới thiệu báo khả thi xác, ứng dụng việc tính tốn nội suy giá trị dị thường độ cao nói riêng đại lượng trắc địa nói chung với ưu điểm bật bảo tồn tính chất điều hịa cầu đại lượng nội suy Ket từ thử nghiệm cho thấy sai số việc tính tốn mức 10'12m, hồn tồn bỏ qua ngừ cảnh trắc địa Kể thêm nhiễu giả ngẫu nhiên cm vào dừ liệu ban đầu kết nội suy đạt độ xác RMS = cm (tính cho 3276 điểm), độ lớn nhiễu Nhìn chung, cơng thức tính tốn wavelets phức tạp so với phương pháp nội suy quen thuộc khác Tuy nhiên, ưu điểm quan trọng phương pháp việc bảo toàn tính chất điều hịa cầu cho giá trị nội suy (nhờ vào việc sử dung hàm cầu đa thức Legendre), điều mà phương pháp nội suy tốn học khác khơng làm Ngồi ra, TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN Đ s ố 52-6/2022 Nghiên cứu The reterence [m] 50 100 Q u a s i-o p tím a lity 50 100 Hình 3: xuống dưới: giá trị tham khảo giá trị nội suy (cỏ áp dụng quy hỏa quasi-optimality) Hệ trục tọa độ theo thứ tự hàng cột lưới So sánh giá trị tham khảo (hình 3, trên) giá trị nội suy (hình 3, dưới) cách trực quan ta thấy giốne giá trị (màu sắc) kiếu mẫu Sự khác biệt thấy mắt thường khu vực nhỏ xíu sát hình Cịn số độ lệch quân phương (RMS) 3276 điểm cm, mức độ nhiều giả ngẫu nhiên thêm vào số liệu ban đầu trước đó; hay nói cách khác, sai sổ cuối mức độ với sai số liệu, nghĩa quy trình tính tốn nội suy xác Xin nói thêm vấn đề tính tốn xuất thử nghiệm này, hệ phưong trình chuẩn (khi giải hệ số tỷ lệ) (KrePTK r e p ) ă - K repTf = 08) bị rơi vào tình trạng giả định yếu (ill-conditioned) mà để giải giải pháp quy hóa (regularization) thích họp áp dụng gọi quasi-optimality [8] vấn đề quy hóa khơng trình bày vượt phạm vi báo Kết luận Với hai thư nghiệm thiết kế cho hai bối cảnh từ đơn giàn đến phức tạp, từ toàn cầu đến khu vực, kết cho thấy quy trình tính tốn dựa cách tiếp cận wavelets giới thiệu báo khả thi xác, ứng dụng việc tính tốn nội suy giá trị dị thường độ cao nói riêng đại lượng trắc địa nói chung với ưu điểm bật bảo tồn tính chất điều hịa cầu đại lượng nội suy Kết từ thử nghiệm cho thấy sai số việc tính tốn mức 10"12 m, hồn tồn bỏ qua ngừ cảnh trắc địa Kể thêm nhiễu giả ngầu nhiên cm vào liệu ban đầu kết nội suy đạt độ xác RMS = cm (tính cho 3276 điểm), độ lớn nhiễu Nhìn chung, cơng thức tính tốn wavelets phức tạp so với phương pháp nội suy quen thuộc khác Tuy nhiên, ưu điểm quan trọng phương pháp việc báo tồn tính chất điều hịa cầu cho giá trị nội suy (nhờ vào việc sử dung hàm cầu đa thức Legendre), điều mà phương pháp nội suy tốn học khác khơng làm Ngồi ra, TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN Đ SỐ 52-6/2022 Nghiên cứu tích họp toán chuyển đổi (giữa đại lượng khác nhau) vào nhân tái tạo phương pháp cịn áp dụng vào vấn đề xây dựng mơ hình geoid [9].Q Tài liệu tham khảo [1] Berlin Keller w (Ed) (2004) Wavelets in Geodesy and Geodynamics Walter de Gruyte [2] Eicker A (2008) Gravity ííeld refínement by radial basis iìmctions from in-situ satell data Dissertation, Bonn Ưniversity [3] Freeden w, Gervens T, Schreiner M (1998) Constructive approximation on the sphe (with applications to geomathematics) Clarendon Press, Oxford [4] Schmidt M, Fengler M, Mayer-Guerr T, Eicker A, Kusche J, Sanchez L, Han s c (200 Regional gravity modeling in terms of spherical base íunctions J Geodesy 81: 17-38 [5] Moritz H (1989) Advanced physical geodesy (2 ed.) Wichmann, Karlsruhe [6] Pavlis NK, Holmes SA, Kenyon s c , Factor JK (2008) An Earth Gravitational Mod to degree 2160: EGM2008 EGU General Assembly 2008 Vienna, Austria, April 13-18, 2008 [7] Lương Bảo Bình (2016) Tính tốn dị thường độ cao từ hệ số điều hòa độ bậc nh định dải sóng dài, Tạp chí Phát triển khoa học & công nghệ, Tập 19, K4 - 2016, 1 -1 [8] Tikhonov AN, Arsenin VY (1977) Solutions of ill-posed problems Wiley, Newyork [9] Luong BB (2011) Combined gravity Tield modeling from satellite gravity and terrestri data sources applying multi-resolution analysis Dissertation, Graz University of Technology, o Summary Apply wavelets to interpolate height anomalies Luong Bao Binh Ho Chi Minh Citv University o f Technology This paper introduces the wavelets approach applied to the interpolation of height anomalies It can be used to compute interpolateđ values for arbitrary positions from input values with the arbiữary (spatial) distribution Following the theoretical íundamental, two simulations using model EGM 2008 for a global grid and local irregular positions are presented to illustrate and authenticate the correct work o f this method.O Keywords: wavelets, reproducing kemel, spherical hannonics, interpolation TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 52-6/2022 ... trị nội suy so sánh với giá trị tham khảo để đánh giá tính xác phương pháp nội suy Đe có sở so sánh giá trị dị thường độ cao 5796 điểm ban đầu giá trị dị thường độ cao tham khảo 3276 điểm nội suy. .. F 2, G2, G3 minh họa rõ cho tính chất hoạt động lọc wavelets với mức chi tiết ngày cao cấp tín hiệu cao 3.2 Nội suy dị thường độ cao độ bậc điều hòa cao từ mạng lưới điểm gổc khu vực phân bố không... Phần áp dụng sở lý thuyết vào tính tốn nội suy giá trị dị thường độ cao đế đánh giá tính hiệu độ xác cách tiếp cận Kết thử nghiệm thảo luận Đe kiểm chứng tính đắn độ xác tính tốn nội suy phương