Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (97.86 KB, 21 trang )
<span class='text_page_counter'>(1)</span><div class='page_container' data-page=1>
Cho hàm y = f(x) và bảng số
x x<sub>o</sub> x<sub>1</sub> x<sub>2</sub> . . . x<sub>n</sub>
y y<sub>o</sub> y<sub>1</sub> y<sub>2</sub> . . . y<sub>n</sub>
Để tính gần đúng đạo hàm, ta xấp xỉ hàm
bằng đa thức nội suy Lagrange L<sub>n</sub>(x)
Ta coù / /
/ / / /
( ) ( )
( ) ( )
<i>n</i>
<i>n</i>
<i>f</i> <i>x</i> <i>L x</i>
<i>f</i> <i>x</i> <i>L</i> <i>x</i>
x x<sub>0</sub> x<sub>1</sub>
y y<sub>0</sub> y<sub>1</sub>
h = x<sub>1</sub>- x<sub>0</sub>
y<sub>0</sub> = f(x<sub>0</sub>)
y<sub>1</sub> = f(x<sub>1</sub>) = f(x<sub>0</sub>+h)
Đa thức nội suy Lagrange
0
1
0 1
0 1 1 0
0 1
1 0
( )
( )
( )
( ) ( )
( ) ( )
<i>n</i>
<i>x</i> <i>x</i>
<i>x</i> <i>x</i>
<i>L x</i> <i>y</i> <i>y</i>
<i>x</i> <i>x</i> <i>x</i> <i>x</i>
<i>x</i> <i>x</i> <i>x</i> <i>x</i>
<i>y</i> <i>y</i>
<i>h</i> <i>h</i>
Do đó với mọi x [x<sub>0</sub>, x<sub>1</sub>] ta có
1 0 ( 0 ) ( )0
'( ) <i>y</i> <i>y</i> <i>f x</i> <i>h</i> <i>f x</i>
<i>f x</i>
<i>h</i> <i>h</i>
0 0
0
( ) ( )
'( ) <i>f x</i> <i>h</i> <i>f x</i>
<i>f x</i>
<i>h</i>
1 0
1
'( ) <i>y</i> <i>y</i>
<i>f x</i>
<i>h</i>
Thay x<sub>1</sub> baèng x<sub>0</sub>
0 0
0
( ) ( )
'( ) <i>f x</i> <i>f x</i> <i>h</i>
<i>f x</i>
<i>h</i>
0 1
2
2 max | "( ) |<sub>[</sub> <sub>, ]</sub>
2 <i>x</i> <i>x x</i>
<i>M h</i>
<i>với M</i> <i>f</i> <i>x</i>
f’(1.8) và sai số với h = 0.1, 0.01, 0.001
Ta coù <i>f</i> '(1.8) <i>f</i> (1.8 <i>h</i>) <i>f</i> (1.8)
<i>h</i>
Sai soá
2
1
"( )
<i>f</i> <i>x</i>
<i>x</i>
<sub>2</sub> max | "( ) | 1<sub>2</sub>
1.8
<i>M</i> <i>f</i> <i>x</i>
2
2(1.8)
<i>h</i>
h f’(1.8)
0.1 0.540672212 0.016
0.01 0.554018037 0.16x10-2
0.001 0.555401292 0.16x10-3
x x<sub>0</sub> x<sub>1</sub> x<sub>2 </sub>
y y<sub>0</sub> y<sub>1</sub> y<sub>2 </sub>
h = x<sub>2</sub> - x<sub>1</sub> = x<sub>1 </sub>- x<sub>0</sub>
y<sub>0</sub> = f(x<sub>0</sub>)
Đa thức nội suy Lagrange
0 2 0 1
1 2
0 1 2
0 1 0 2 1 0 1 2 2 0 2 1
0 1 0 2 1 2
2 1 0
2 2 2
( )( ) ( )( )
( )( )
( )
( )( ) ( )( ) ( )( )
( )( ) ( )( ) ( )( )
2 2
<i>n</i>
<i>x x</i> <i>x x</i> <i>x x</i> <i>x x</i>
<i>x x x x</i>
<i>L x</i> <i>y</i> <i>y</i> <i>y</i>
<i>x</i> <i>x x</i> <i>x</i> <i>x</i> <i>x</i> <i>x</i> <i>x</i> <i>x</i> <i>x</i> <i>x</i> <i>x</i>
<i>x x</i> <i>x x</i> <i>x x</i> <i>x x</i> <i>x x x x</i>
<i>y</i> <i>y</i> <i>y</i>
<i>h</i> <i>h</i> <i>h</i>
Do đó với mọi x [x<sub>0</sub>, x<sub>2</sub>] ta có
0 1 2
2 1 2 0 0 1
2 2 2
( ) ( ) ( )
'( ) ( 2 ) ( ) ( 2 )
2 2 2
<i>x x</i> <i>x x</i> <i>x x</i>
<i>f x</i> <i>y</i> <i>y</i> <i>y</i> <i>y</i> <i>y</i> <i>y</i>
<i>h</i> <i>h</i> <i>h</i>
2 1 0
2
( 2 )
"( ) <i>y</i> <i>y</i> <i>y</i>
<i>f</i> <i>x</i>
<i>h</i>
Suy ra đạo hàm cấp 1
0 1 2
0
2 0
1
0 1 2
2
( 3 4 )
'( )
2
( )
'( )
2
( 4 3 )
'( )
2
<i>y</i> <i>y</i> <i>y</i>
<i>f</i> <i>x</i>
<i>h</i>
<i>y</i> <i>y</i>
<i>f</i> <i>x</i>
<i>h</i>
<i>y</i> <i>y</i> <i>y</i>
<i>f</i> <i>x</i>
<i>h</i>
Công thức thứ 1 gọi là công thức sai phân tiến
0 0 0
0
3 ( ) 4 ( ) ( 2 )
'( )
2
<i>f x</i> <i>f x</i> <i>h</i> <i>f x</i> <i>h</i>
<i>f</i> <i>x</i>
<i>h</i>
Công thức thứ 2 gọi là công thức sai phân hướng
tâm thường viết dưới dạng (thay x<sub>1</sub> = x<sub>0</sub>)
0 0
0
( ) ( )
'( )
2
<i>f x</i> <i>h</i> <i>f x</i> <i>h</i>
<i>f</i> <i>x</i>
<i>h</i>
Công thức thứ 3 gọi là công thức sai phân lùi
thường viết dưới dạng (thay x<sub>2</sub> = x<sub>0</sub>)
0 0 0
0
( 2 ) 4 ( ) 3 ( )
'( )
2
<i>f x</i> <i>h</i> <i>f x</i> <i>h</i> <i>f x</i>
<i>f</i> <i>x</i>
<i>h</i>
0 2
2
3
3 max | "'( ) |<sub>[</sub> <sub>,</sub> <sub>]</sub>
6 <i>x</i> <i>x x</i>
<i>M h</i>
<i>với M</i> <i>f</i> <i>x</i>
đạo hàm cấp 2
2 1 0
1 2
( 2 )
''( ) <i>y</i> <i>y</i> <i>y</i>
<i>f</i> <i>x</i>
<i>h</i>
Thay x<sub>1</sub> = x<sub>0 </sub>ta được
0 0 0
0 2
( ) 2 ( ) ( )
''( ) <i>f x</i> <i>h</i> <i>f x</i> <i>f x</i> <i>h</i>
<i>f</i> <i>x</i>
<i>h</i>
0 2
2
(4)
4
4 max |<sub>[</sub> <sub>,</sub> <sub>]</sub> ( ) |
12 <i>x</i> <i>x x</i>
<i>M h</i>
<i>với M</i> <i>f</i> <i>x</i>
a. Dùng cơng thức sai phân hướng tâm, tính xấp
xỉ f’(3) với h = 0.1, 0.01, 0.001
b. Tính xấp xỉ f”(3) với h = 0.1, 0.01, 0.001
giaûi
(3 ) (3 )
'(3)
2
<i>f</i> <i>h</i> <i>f</i> <i>h</i>
<i>f</i>
<i>h</i>
h f’(3)
2
(3 ) 2 (3) (3 )
''(3) <i>f</i> <i>h</i> <i>f</i> <i>f</i> <i>h</i>
<i>f</i>
<i>h</i>
h f’’(3)
Cho hàm f(x) xác định và khả tích trên [a,b].
Ta cần tính gần đúng tích phân :
( )
<i>b</i>
<i>a</i>
<i>I</i>
Ta phân hoạch đoạn [a,b] thành n đoạn bằng
nhau với bước h = (b-a)/n
x<sub>o</sub>= a, x<sub>1</sub> = x<sub>0</sub> +h, ... , x<sub>k</sub> = x<sub>0</sub> + kh, ... , x<sub>n</sub> = b
Đa thức Lagrange trong TH các điểm cách đều
0
( 1)
( ) ( 1)...( )
!( )!( )
<i>n k</i>
<i>n</i>
<i>n</i> <i>k</i>
<i>k</i>
<i>L x</i> <i>q q</i> <i>q</i> <i>n</i> <i>y</i>
<i>k n</i> <i>k</i> <i>q</i> <i>k</i>
<i>x</i> <i>a</i>
<i>với q</i>
<i>h</i>
( 1) ( 1)...( )
( )
!( )!( )
( 1) ( 1)...( ) ( )
!( )!( )
<i>b</i> <i><sub>n</sub></i> <i>b</i> <i>n k</i>
<i>n</i> <i>k</i>
<i>k</i>
<i>a</i> <i>a</i>
<i>n</i> <i><sub>n k</sub></i>
<i>n</i>
<i>k</i>
<i>k</i>
<i>q q</i> <i>q n</i>
<i>I</i> <i>L x dx</i> <i>y dx</i>
<i>k n k</i> <i>q k</i>
<i>q q</i> <i>q n b a</i>
<i>dq y</i>
<i>k n k</i> <i>q k</i> <i>n</i>
Công thức trên gọi là công thức Newton-cotes,
các hệ số H<sub>k</sub> gọi là các hệ số cotes.
Heä số cotes có các tính chất sau :
0
1
0,
<i>n</i>
<i>k</i>
<i>k</i>
<i>n k</i> <i>k</i>
<i>H</i>
<i>H</i> <i>H</i> <i>k</i> <i>n</i>
0
* ( )
<i>n</i>
<i>k k</i>
<i>k</i>
<i>I</i> <i>I</i> <i>b a</i> <i>H y</i>
0
( 1) ( 1)...( )
!( )! ( )
<i>n</i>
<i>n k</i>
<i>k</i>
<i>q q</i> <i>q n</i>
<i>với H</i> <i>dq</i>
<i>n k n k</i> <i>q k</i>
( 1) ( 2)
1 <sub>[ , ]</sub> 2 <sub>[ , ]</sub>
| ( 1)...( ) |
( 1)!
| *|
| ( 1)...( ) |
( 2)!
max | ( ) | max | ( ) |
<i>n</i> <i>n</i>
<i>n</i>
<i>n</i> <i>n</i>
<i>n</i>
<i>n</i> <i>n</i>
<i>n</i> <i><sub>x a b</sub></i> <i>n</i> <i><sub>x a b</sub></i>
<i>M h</i>
<i>q q</i> <i>q n dq với n lẻ</i>
<i>n</i>
<i>I I</i>
<i>M h</i>
<i>q q</i> <i>q n dq với n chẵn</i>
<i>n</i>
<i>M</i> <i>f</i> <i>x vaø M</i> <i>f</i> <i>x</i>
Xét n = 1, ta coù h= b-a
I (b-a)(H<sub>o</sub>y<sub>o</sub> + H<sub>1</sub>y<sub>1</sub>)
1
0
0
1
( 1)
2
<i>H</i>
2
<i>H</i> <i>H</i>
Vaäy ( ) <sub>(</sub> <sub>0</sub> <sub>1</sub><sub>)</sub> ( ) <sub>( ( )</sub> <sub>( ))</sub>
2 2
<i>b a</i> <i>b a</i>
<i>I</i> <i>y</i> <i>y</i> <i>f a</i> <i>f b</i>
3 1 3
2 2
0
| ( 1) |
2! 12
<i>M h</i> <i>M h</i>
<i>q q</i> <i>dq</i>
Ta phân hoạch đoạn [a,b] thành n đoạn bằng
nhau [x<sub>0</sub>, x<sub>1</sub>], [x<sub>1</sub>, x<sub>2</sub>], ... , [x<sub>n-1</sub>, x<sub>n</sub>].
Ta coù
1 2
0 1 1
( ) ( ) ... ( )
<i>n</i>
<i>n</i>
<i>x</i>
<i>x</i> <i>x</i>
<i>x</i> <i>x</i> <i>x</i>
<i>I</i> <i>f x dx</i> <i>f x dx</i> <i>f x dx</i>
1 0 2 1 1
0 1 1 2 1
0 1 1 2 1
( ) ( ) ( )
( ) ( ) ... ( )
2 2 2
( ) ( ) ... ( )
2 2 2
<i>n</i> <i>n</i>
<i>n</i> <i>n</i>
<i>n</i> <i>n</i>
<i>x x</i> <i>x</i> <i>x</i> <i>x</i> <i>x</i>
<i>y</i> <i>y</i> <i>y y</i> <i>y</i> <i>y</i>
<i>h</i> <i>h</i> <i>h</i>
<i>y</i> <i>y</i> <i>y y</i> <i>y</i> <i>y</i>
Vậy
0 1 1
3 2
2 <sub>(</sub> <sub>)</sub> 2
12 12
<i>M h</i> <i>M h</i>
<i>n</i> <i>b a</i>
Xét n = 2, ta coù h = (b-a)/2
I (b-a)(H<sub>o</sub>y<sub>o</sub> + H<sub>1</sub>y<sub>1</sub>+H<sub>2</sub>y<sub>2</sub>)
1
0
0
1 1
( 1)( 2)
4 6
<i>H</i>
2 0
1
6
<i>H</i> <i>H</i>
0 1 2 1
2
1
3
<i>H</i> <i>H</i> <i>H</i> <i>H</i>
Vaäy ( ) ( <sub>0</sub> 4 <sub>1</sub> <sub>2</sub>)
6
<i>b a</i>
4 1 5
2
4 4
0
| ( 1)( 2) |
4! 90
<i>M h</i> <i>M h</i>
<i>q q</i> <i>q</i> <i>dq</i>
Ta phân hoạch đoạn [a,b] thành n đoạn bằng
nhau [x<sub>0</sub>, x<sub>1</sub>], [x<sub>1</sub>, x<sub>2</sub>], ... , [x<sub>n-1</sub>, x<sub>n</sub>].