Xấp xỉ biến đổi laplace ngược và công thức cầu phương nội suy

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TPHCM LÊ DUY THỨC Chuyên Ngành Mã Số : : Toán Giải Tích 604601 LUẬN VĂN THẠC SĨ TOÁN HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: TS NGUYỄN CAM Thành Phố Hồ Chí Minh – Năm 2006 LỜI CẢM ƠN Tôi xin bày tỏ lòng biết ơn sâu sắc đến TS Nguyễn Cam, người Thầy dạy dỗ, dìu dắt từ năm đầu đại học Xin chân thành cảm ơn PGS.TS Lê Hoàn Hóa, PGS.TS Nguyễn Bích Huy, TS Nguyễn Thành Long, người Thầy quan tâm, giúp đỡ truyền đạt cho kiến thức tảng thời gian học đại học cao học Xin cảm ơn Thầy hội đồng chấm luận văn cho nhận xét quý báu, Thầy- Cô truyền đạt kiến thức học phần Cảm ơn BGH Trường PTTH Mạc Đónh Chi TPHCM, đồng nghiệp tạo điều kiện, động viên để hoàn thành khoá học Cảm ơn gia đình, bạn bè người thân hỗ trợ, giúp đỡ nhiều mặt Xin cảm ơn bạn Thúy Trang, University of Western Australia, động viên cung cấp nhiều tài liệu bổ ích trình làm luận văn LỜI MỞ ĐẦU Phép biến đổi Laplace có nhiều áp dụng quan trọng khoa học kỹ thuật Bài toán khôi phục hàm gốc từ hàm ảnh phép biến đổi Laplace nhiều nhà Toán học quan tâm khảo cứu đến có nhiều phương pháp đưa Trong luận văn này, khảo sát số phương pháp tính xấp xỉ biến đổi Laplace ngược thông qua công thức cầu phương nội suy Trong chứng minh hội tụ công thức nội suy, tính ổn định nghiệm xấp xỉ thu được, minh hoạ việc giải số máy tính thông qua ví dụ cụ thể Luận văn chia làm chương sau : Chương : Trình bày kiến thức chuẩn bị cho việc tính toán tích phân Mellin Chương : Khảo sát số phương pháp tính tích phân Mellin công thức cầu phương nội suy Sau định lý hội tụ trình nội suy tính ổn định nghiệm xấp xỉ Chương : Đưa công thức cầu phương nội suy với độ xác cao Chương : Xây dựng công thức tính toán cho công thức cầu phương nội suy với hệ số cân Cuối ví dụ giải số máy tính CHƯƠNG I : MỘT SỐ KIẾN THỨC CHUẨN BỊ 1.1 Định nghóa Cho f (t ) với t  , f (t ) khả tích đoạn [a,b], (0  a  b) với số phức p    i  , ta định  nghóa F ( p )  e  pt (1) f (t )dt F(p) gọi biến đổi Laplace f (t ) 1.2 Định lý Nếu F(p) xác định p0  0  i0 F(p) xác định p    i  thoaû Re( p  p0 )    0  Chứng minh t Đặt (t )   f (u )e  p0 u du (t  0)  Vì F ( p0 )   e p0 t f (t )dt xaùc định nên lim (t ) tồn Suy tồn số Q cho: t  (t )  Q, t  Xeùt Re( p  p0 )    0  c (với c >0), a>0, b>0 Ta có : ab  a f (t )e  pt a b dt   e  ( p  p0 )t d (t ) a =e  ( p  p0 )t (t ) a b a ab   a ( p0  p )e  ( p  p0 )t (t )dt (tích phân phần)  (a  b).e  ( p  p0 )( a  b )  (a ).e ( p  p0 ) a a b  p  p0  e ( p  p0 )t (t ) dt a  Q e ( p  p0 )( a  b ) Qe ab ( p  p0 ) a  Q p  p0  e ( 0 )t dt a  Q.e c(a  b)  Qe ca a b  Q p  p0  e ct dt a  Q.e  ca  Q.e ca e ct  Q p  p0 c a b a Q p  p0 (e ca  e c ( a  b ) ) c Q  2Qe ca  p  p0 e ca c p  p0  ca  Q(2  ).e   c  2Qe ca  (khi a đủ lớn)  Vậy theo điều kiện Cauchy e  pt f (t )dt hội tụ nên F(p) xác định p với Re p  Re p0 1.3 Định lý Cho F(p) xác định p0  0  i0 F(p) hàm quy nửa mặt phẳng Rep> 0 Chứng minh Xét p    i  với   0 Lấy miền D đóng bị chặn chứa nửa mặt phẳng Rep > 0 , p  D Khi  Re( p  p0 )  2c  p  p0  M , p  D Tồn c >0, M >0, thỏa :  n Xét dãy Fn ( p )   e  pt f (t )dt , p  D, n  1,2,3 Ta coù Fn hàm giải tích D vì: n  ( p  h )t Fn ( p  h)  Fn ( p ) n  ( p  h )t e  e  pt  pt f (t )dt  lim  (e  e ) f (t )dt  lim  lim h h h0 h0 h h0 n  lim  e  pt h0 e  ht  f (t )t dt ht n    e  pt f (t )tdt Với m, n > 0, m < n ta coù : n Fn ( p)  Fm ( p)   e  pt m f (t )dt   e  pt n f (t )dt   e  pt f (t )dt m Theo chứng minh định lý mục 1.2 ta coù: n  pt  e f (t )dt  Q(2  m M  cm )e c Fn ( p)  Fm ( p)  , p  D , m đủ lớn Vậy Fn(p) hội tụ F(p) Fn(p) giải tích, F ( p ) giải tích D Hay F(p) quy nửa mặt phẳng Re p  0 1.4 Nhận xét  Đặt E     R /  e  pt f (t )dt hội tụ  ( P    i ) Và   inf E + Nếu    : E    F ( p ) không xác định p + Nếu    : E  R  F ( p ) xác định với p + Nếu   R : - Với    F(p) không xác định - Với    F(p) xác định F(p) hàm qui 1.5 Định nghóa f(t) gọi hàm gốc thỏa: 1) f(t) xác định với t  R , f(t) = với t  2) f (t ) khả tích đoạn hữu hạn  3) F ( p )   e pt f (t )dt xác định p Lúc ta gọi F(p) hàm ảnh biến đổi Laplace f 1.6 Định lý Cho M  0,   R cho : f (t )  Me t , t  ø, ta có :    (  định nghóa 1.4) Chứng minh Đặt Re p   giả sử    thì:   e  pt f (t ) dt    f (t ) et dt    M et et dt   M  e(  )t dt    F ( p ) xác định với Re p     Vaäy    1.7 Định lý Cho F(p) xác định p0  0  i0 lim F ( p)  p  ( nửa mặt phẳng Re p   ) Chứng minh Vì F(p) xác định p0 nên F(p) xác định p có Re p   Xét A>0, ta có :  F ( p)   e  pt A f (t )dt   e  pt  f (t )dt  e A  pt f (t )dt  M1  M  với M  e  pt  f (t )dt  e A     p0 t f (t ).e ( p0  p )t  dt  A  e  p0 t f (t ) e(0  )t dt A e  p0 t f (t ) dt (vì e(0  )t  1) A Vì F ( p0 ) xác định nên với A đủ lớn M   Ta coù f (t ) khả tích đoạn [0,A] nên tồn hàm g(t) khả vi liên tục thỏa : A  f (t )  g (t ) e0 t dt   A Do M1   e  pt A f (t )dt   e  pt A [ f (t )  g (t )]dt   e  pt g (t )dt  M  M 0 A A Trong M   et f (t )  g (t ) dt   e A M4   e  pt g (t )dt   g (t )e  pt p 0 t f (t )  g (t ) dt   A A A  pt   e g '(t )dt p0 1 = g (0)  g ( A).e  pA   e  pt g '(t )dt p p p0 suy M  p  với p  A  0 A   e 0 t g '(t ) dt  g (0)  g ( A).e     L (L phụ thuộc A) p L M   F ( p )  3 Vaäy lim F ( p)  p   1.8 Định nghóa Cho hàm số g(t) xác định R ta gọi g biểu diễn tích phân Fourier với t ta coù : 1  i t   g (t  0)  g (t  0)   e  g (t )ei x dxd  (2) 2   Phương trình (2) gọi công thức Fourier 1.9 Định lý Cho hàm gốc f(t) cho g (t )  f (t ).e ct (với c  R ) thỏa mãn :  i)  g (t ) dt hội tụ  it  ii) g (t )  e  g ( x )e i x dxd   2   iii) f (t )   f (t  0)  f (t  0) c  i F ( p ).e pt dp Thì f (t )   2i c  i (coâng thức Mellin) Chứng minh Cho f(t) hàm gốc c  R thoả f (t ).e ct khả tích tuyệt đối  0,  Đặt g(t)= f (t ).e ct Giả sử g(t) thoả (2) để đơn giản cách ghi, ta viết : g (t )   g (t  0)  g (t  0) g (t )   it  e  g ( x )e i x dxd   2   Công thức (2) trở thành :  (i c )t  f (t )  e g ( x)e ix dxd    2  Suy  Với F ( p )   f (t )e  pt  dt xác định p  c  i   f (t ) e (3)  pt dt   f (t ) e ct dt nên f (t ).e ct khả tích tuyệt đối đường thẳng c  i  (với      ) Ta coù :  F ( p )  F (c  i )   f (t ).e  (c  i )t dt   (i c )t  f (t )  e f ( x )e  (c  i )t dxd    2   Vaø (3) cho ta :  (i c )t F (c  i)d  = e 2  Treân đường thẳng p = c  i  dp  id  , nên ta lại có : c  i f (t )  F ( p ).e pt dp  2i c  i (4) Coâng thức (4) gọi công thức Mellin 1.10 Định lý  Xét phương trình e  pt f (t )dt  F ( p ) (*) Trong F(p) cho trước f(t) hàm gốc phải tìm, toán không chỉnh theo nghóa Hadamard Chứng minh Với F(p) hàm giải tích (*) vô ngiệm Bây ta xét f(t) ứng với hàm ảnh F(p) f1(t) ứng với hàm ảnh F1(p) cho:  n  f ( t ),  t   n2 f1 (t )    f (t ), t  n d ( f , f1 )  f  f1  n (chọn n lớn)  Trong khi: F1 ( p )  F ( p )   e  pt n2 ( f1 (t )  f (t ))dt   n.e  pt dt ( x  t ) Pn ( x) Pn (t )  Bn [ Pn 1 ( x) Pn (t )  Pn ( x ) Pn 1 (t )]  Dn [ Pn 1 ( x) Pn (t )  Pn ( x ) Pn 1 (t )] Làm tương tự cho Pn 1  x  , Pn   x  , , P1  x  , P0  x  ta hệ phương trình sau : ( x  t ) Pn 1 ( x) Pn 1 (t )  Bn 1[ Pn ( x) Pn 1 (t )  Pn 1 ( x ) Pn (t )]  Dn 1[ Pn  ( x) Pn 1 (t )  Pn 1 ( x ) Pn  (t )] ( x  t ) Pn  ( x ) Pn  (t )  Bn  [ Pn ( x) Pn  (t )  Pn  ( x) Pn (t )]  Dn  [ Pn  ( x) Pn  (t )  Pn  ( x) Pn  (t )] ………………………………………………………………………………………………………………… ( x  t ) P1 ( x ) P1 (t )  B1[ P2 ( x ) P1 (t )  P1 ( x) P2 (t )]  D1[ P0 ( x) P1 (t )  P1 ( x ) P0 (t )] ( x  t ) P0 ( x) P0 (t )  B0 [ P1 ( x) P0 (t )  P0 ( x ) P1 (t )] Nhaân hai vế phương trình cho : 1, Dn 1 D D D2 Dn 1 Dn  D1 , , n 1 n  , Bn  Bn  Bn  B1 Bn  Bn  B0 Ta hệ phương trình tương đương : ( x  t ) Pn 1 ( x) Pn 1 (t )  Bn 1[ Pn ( x) Pn 1 (t )  Pn 1 ( x ) Pn (t )]  Dn 1[ Pn  ( x) Pn 1 (t )  Pn 1 ( x ) Pn  (t )] Dn 1 ( x  t ) Pn  ( x ) Pn  (t ) Bn   Dn 1[ Pn 1 ( x) Pn  (t )  Pn  ( x ) Pn 1 (t )  Dn 1 Dn  [ Pn  ( x) Pn  (t )  Pn  ( x) Pn  (t )] Bn  …………………………………………………………………………………………………………………………… Dn 1 Dn  D2 ( x  t ) P1 ( x) P1 (t ) Bn  Bn  B1  Dn 1 Dn  D2 [ P2 ( x) P1 (t )  P1 ( x ) P2 (t )] Bn  Bn  B2  Dn 1 Dn  D1 [ P0 ( x ) P1 (t )  P1 ( x) P0 (t )] Bn  Bn  B1 Dn 1 Dn  D2 D D D2 ( x  t ) P0 ( x) P0 (t )  n 1 n  [ P1 ( x) P0 (t )  P0 ( x ) P1 (t )] Bn  Bn  B0 Bn  Bn  B1 Cộng phương trình theo vế ta : n 1 ( x  t )   m pm ( x) Pm (t )  Bn 1[ Pn ( x ) Pn 1 (t )  Pn 1 ( x) Pn (t )] (3.3.5) m0 Trong : m  Trong (3.3.5) đặt t  xk  cho x  xk  Dn 1 Dn  D1 m , m  0,1, , n  2;  n 1  Bn  Bn  Bm xk nghiệm đa thức Pn ( x) Chia hai vế (3.3.5) pk 1  ta có : p pk n 1  m Pm ( x) Pm ( xk )  Bn 1 m0 Pn ( x ) Pn 1 ( xk ) x  xk Nhaân hai vế cho e p p  s sau lấy tích phân ta : P ( x) Pn 1 ( xk )   i  p  s n 1   i p  s ( ) ( ) e p  P x P x dp  e p Bn 1 n dp  m m m k   2i  i i x x   m0 k   i Ta coù :   i p  s Pm ( xk ) e p P ( )dp  0, m  m 2i  i p (do điều kiện trực giao (3.1.7))   i p  s 1 Pm ( xk ) e p Pm ( )dp  , m0  p 2i   i (s ) Và Do ta coù :   i 0 1  Bn 1 pn 1 ( xk ) e p P s  (s ) 2i  i Pn ( ) p dp 1 (  ) p pk Chia hai vế cho P 'n ( xk ) ta : Pn ( ) 1 p 0  Bn 1 pn 1 ( xk ) e p P s dp  1 2i  i  ( s) P 'n ( xk ) (  ) P 'n ( ) p pk pk  i  Vậy Ta cần phải tính Ak  0 Bn 1( s ) Pn 1 ( xk ) P 'n ( xk ) 0 : Bn 1 0 D D D  n 1 n  Bn 1 Bn 1Bn  B1 B0 Ta coù : Dn 1 (n  1)(2n  s  2)(2 n  s  3)  Bn 1 (2n  s  3)(2n  s  4)(n  s  2)  (n  1)(2n  s  2) ; (2n  s  4)(n  s  2) B0  s Do : 0 (1)n 1 (n  1)!(2n  s  2)(2n  s  4) (2  s ) s  Bn 1 (2 n  s  4) (2  s ) s(n  s  2)(n  s  3) (1  s ) s (1)n 1 (n  1)!(2n  s  2)  (n  s  2)(n  s  3) (1  s) s (1)n 1 (n  1)!(2n  s  2)( s)  (n  s  1) Thay (3.3.7) vào (3.3.6) ta có : (3.3.7) (1) n 1 (n  1)!(2n  s  2) Ak  (n  s  1) Pn ( xk ) P 'n ( xk ) (3.3.8) Để tính điểm pk công thức cầu phương, nói chúng nghiệm đa thức Pn( s ) ( x) Để tính pk ta phải xác định hệ số Pn( s ) ( x) sử dụng phương pháp để tìm nghiệm chúng Điều gặp hai khó khăn : thứ với giá trị tham số s ta phải xác định hệ số Pn( s ) ( x) thông qua công thức truy hồi Thứ hai việc tính toán nghiệm đa thức Pn( s ) ( x) nhờ công thức Newton làm giảm độ xác n có giá trị lớn Để tránh điều ta sử dụng phương pháp khác để tìm điểm công thức cầu phương Với nghiệm xk đa thức Pn( s ) ( x) ta xây dựng hệ phương trình đại số chứa xk tham số s Xét phương trình vi phân có nghiệm đa thức Pn( s ) ( x) : x Pn( s) '' ( x)  ( sx  1) Pn( s ) ' ( x)  n(n  s  1) Pn( s ) ( x)  (3.3.9) Ta có nhận xét xk2 Pn( s ) ' ( xk )  Thaät vaäy từ (3.2.1) số hạng tự Pn( s ) ( x) khác nên xk  Hơn nữa, từ (3.3.9) Pn( s ) ' ( xk )  Pn( s ) '' ( xk )  Đạo hàm phương trình (3.3.9) n lần ta có : (4) Pn( s ) ''' ( xk )  0, Pn( s ) ( xk )  0, , Pn( s ) Điều vô lý Pn( s ) (n) (n) ( xk )  ( xk )  Trong (3.3.9) lấy x  xk chia hai vế cho xk2 Pn( s ) ' ( xk )  ta : Pn( s ) ' ( xk ) Pn( s ) '' ( xk )  s  0 xk x (k ) Ta vieát Pn( s ) ( x) dạng : Pns ( x) n  a ( x  x1 ( x  x2 ) ( x  xn )  a  ( x  xm ) m 1 : (3.3.10) n n Pn( s ) ' ( x)  a   ( x  xm ) i 1 m 1 mi n Pn( s ) ' ( xk )  a  ( xk  xm ) m 1 m k Pn( s ) '' ( x) n n  a n   ( x  xm ) j 1 i 1 m 1 i  j m  i, j n Pn( s ) '' ( xk )  2a  n  ( xk  xm ) j 1 m 1 j  k m  j, k Suy : n n 2a  Pn( s ) '' ( xk ) Pn( s ) ' ( xk )  ( xk  xm ) j 1 m 1 j  k m  j, k n  a  ( xk  xm ) n   (x j 1 jk k  xj ) m 1 mk Phương trình (3.3.10) viết lại sau : n  (x j 1 jk Thay xk k s   0  x j ) xk xk (3.3.11) ta : pk (3.3.11) n pk p j  (p j 1 jk n   (p j 1 jk  Pk ) j 2pj j  Pk ) n   2(1  j 1 jk  spk  pk2   s  pk  pk )  s  pk  p j  pk n  2(n  1)   (p j 1 jk n   (p j 1 j k k pk  s  pk  0, j  Pk ) 2n  s  1 0 pk  Pj ) (3.3.12) Từ (3.3.12) thay k  1, 2, , n ta có hệ phương trình xác định điểm pk cho công thức cầu phương CHƯƠNG CÁC PHƯƠNG PHÁP TÍNH BIẾN ĐỔI LAPLACE NGƯC THÔNG QUA CÔNG THỨC CẦU PHƯƠNG VỚI HỆ SỐ CÂN BẰNG 4.1 Xây Dựng Công Thức Tính Toán Giả sử toán tính biến đổi Laplace ngược quy toán tính tích phân : c  i p  s e p ( p )dp 2i c i (4.1.1) Ta xây dựng công thức cầu phương với hệ số cân : n c  i p  s e p ( p )dp  Cn  ( pk ) 2i c i k 1 (4.1.2) Các đại lượng chưa biết (4.1.2) số Cn pk ( k  1,2 , n ) Ta chọn chúng cho (4.1.2) xác theo đa thức bậc n theo biến xác theo haøm ( p )  pk Yêu cầu tương đương với công thức (4.1.2) p (k  0,1, , n) Nhân tử Cn xác định từ điều kiện mà (4.1.2) xác với ( p )  :   i p  s e p dp  nCn 2i i neân : 1  i p  s Cn  e p dp   n 2i   i n( s) Với giá trị Cn vừa tìm (4.1.2) trở thành : (4.1.3)   i p  s n e p ( p )dp   ( pk ) 2i i n ( s) k 1 Chuyển từ biến p sang biến x  (4.1.4) 1 đặt xk  ta có hệ phương trình sau : pk p n( s )   i p  s 1 n( s )  x1  x2  xn  e p p dp   2i  i  ( s  1)  n( s )   i p  s 2 n (s )  x12  x22  xn2  e p p dp   i i ( s  2)   n( s )   i p  s  n n( s )  n n n x1  x2  xn  e p p dp   i i  (s  n)  (4.1.5) Từ hệ ta nhận giá trị xk nên nhận pk Nhưng hệ phi tuyến nên gặp khó khăn việc tìm nghiệm Chúng cố gắng tìm phương pháp khác để tính xk Ta giới thiệu hàm n ( x) mà nghiệm số xk : n ( x )  ( x  x1 )( x  x2 ) ( x  xn ) Khai trieån đa thức theo biến x : n ( x )  x n  A1 x n 1  A2 x n    An 1 x  An Để ý hệ số Ak hàm đối xứng nghiệm Cộng phương trình (4.1.5) theo vế ta có : n Si  n( s )  xik  (i  s) (i  1, 2, n) (4.1.6) k 1 Từ tính chất đa thức ta có quan hệ sau gọi quan hệ Newton : S1  A1  S  A1S1  A2     S3  A1S2  A2 S1  A3     S n  A1S n 1  A2 S n    nAn   (4.1.7) Thay (4.1.6) vào (4.1.7) ta : 1   A1    ( S  1) n (s )  1   A1  A2   ( s  2)  ( s  1) n( s)  1   A1  A2  A3   ( s  3)  (s  2)  ( s  1) n( s )    1   A1  A2    ( s  n) ( s  n  1) ( s  n  2)  n  An  An 1   n( s ) ( s  1)  (4.1.8) Các giá trị Ak tìm từ hệ Vì ma trận hệ số có dạng tam giác nên có định thức : k! nk [ ( s)]k neân Ak   n k [(s )]k  , : k!  n( s )  (s  1) ( s  2) n( s ) ( s  1) n ( s) ( s  k  2)  ( s  k  1)  ( s  k  3) ( s  k  2)  ( s  k  4)  ( s  k  3) 0 0 ( s  1)  (s  2) ( s  3) k 1 n ( s) ( s  1) ( s  k  1) ( s  k ) n s Với  k  n A1   Tìm tất Ak ta xây dựng đa thức n ( x) Sau tìm nghiệm xk nó, ta xác định tất điểm pk  cho công thức cầu phương xk (4.1.4) 4.2 Một Ví Dụ Về Lời Giải Số Trong phần này, sử dụng MATLAB để xây dựng chương trình tính xấp xỉ f(t) công thức cầu phương nội suy với mốc nội suy cách (1.2.1) Chương trình cho phép tính hệ số đa thức xấp xỉ hàm f(t) (hàm gốc) biết F(p) (hàm ảnh) số khoảng nội suy Sau vẽ đồ thị đa thức xấp xỉ Để minh họa, chương trình sau lấy F ( p )  f (t )  e t mà hàm gốc cho baûng p 1 %tinh xap xi f(t) clc disp('chuong trinh tinh xap xi f(t)') disp('nhap vao so khoang noi suy'); n = input(' so khoang noi suy n =' ); %khoi phuc da thuc tu so lk(1/p) % vong lap tinh akj Z=zeros(1,(n+1)); for k=0:n %tinh k giai thua if (k==0) kgt=1; else kgt =1; for l=1:k kgt=kgt*l; end end % tinh (n-k) giai thua if ((n-k)==0) nkgt=1; else nkgt =1; for m=1:(n-k) nkgt=nkgt*m; end end A=[1:(n+1)] A(:,(k+1))=[] B=poly(A) E=(((-1)^(n-k))*(k+1)^n)/((kgt)*(nkgt)) M=B*E X=((k+1)/(k+2))*M Z=Z+X end heso=[]; for j=0:n % tinh j giai thua if (j==0) jgt=1; else jgt =1; for w=1:j jgt=jgt*w; end end a=(1/jgt)*Z((j+1)) heso=[heso,a] end hsdt=heso(:,(n+1):-1:1) t=[0:0.1:n+1]; y=exp(-t); plot(t,polyval(hsdt,t),t,y,'r-') % ( p  1)  ( p  1) /( p  2) axis equal legend('hamxapxi','y=exp(-t)') Đồ thị đa thức xấp xỉ hàm gốc sau cho thấy hội tụ phương pháp : KẾT LUẬN Cho đến nay, toán tính biến đổi Laplace ngược toán mở Ngoài phương pháp nhận nhiều quan tâm nghiên cứu sử dụng tính giải tích hàm ảnh, khai triển hàm gốc dựa vào đa thức trực giao…thì phương pháp số để tính tích phân Mellin có nhiều áp dụng thực tiễn Trong chương thu công thức xấp xỉ hàm gốc sau: f (t )  n  c  i pt  s  n    e p l  p  r p dp ( ) ( )  Ak (t )( pk ) k n  k   p 2i c i k 0  k 0    n Ak (t ) = :  j 0 ak j t s  j 1 ( s  j ) n vaø f (t )   Ak ( pk ) k 1 : (1) n 1 (n  1)!(2n  s  2) Ak  (n  s  1) Pn ( xk ) P 'n ( xk ) Cuối ví dụ cụ thể giải số chương TÀI LIỆU THAM KHẢO Tiếng Anh Krylov V.I and Skoblya N.S (1985), A handbook of Methods of Approximate Fourier Transfomation and Inversion of The Laplace Transformation, Mir publisher, Moscow Krylov V.I (1962), Approximate Calculation of Integrals, the Macmillan Company, New York Krylov V.I., Skoblya N.S., Handbook of Numerical Inversion of Laplace Transforms, IPST Press, Jerusalem Smirnov V.I., and Lebedev, N.A (1964), A Constructive Theory of functions of a Complex Variable, Nauka, Moscow Sveshnikov A.G., Tikhonov N.A (1978), The Theory of Functions of a Complex Variable, Mir publisher, Moscow Tiếng Việt 1.Đặng Đình ng, Trần Lưu Cường, Huỳnh Bá Lân, Nguyễn văn Nhân (2001), Biến Đổi Tích Phân, NXB Giáo Dục ... văn này, khảo sát số phương pháp tính xấp xỉ biến đổi Laplace ngược thông qua công thức cầu phương nội suy Trong chứng minh hội tụ công thức nội suy, tính ổn định nghiệm xấp xỉ thu được, minh hoạ... i ( 3.1.1 ) Trong chương ta xây dựng công thức cầu phương nội suy với đa thức nội suy cấp n -1 theo 1 Độ xác công thức cầu phương nội suy này, với mốc nội suy p p  A  2 cho trước nửa mặt phẳng... nội suy tính ổn định nghiệm xấp xỉ Chương : Đưa công thức cầu phương nội suy với độ xác cao Chương : Xây dựng công thức tính toán cho công thức cầu phương nội suy với hệ số cân Cuối ví dụ giải