TÍNH ĐẠO HÀM VÀ TÍCH PHÂN CỰC HAY pdf

Chương III
TÍNH ĐẠO HÀM VÀ TÍCH PHÂN

1. TÍNH GẦN ĐÚNG ĐẠO HÀM
Trong nhiều bài toán thực tế ta cần phải tính đạo hàm của hàm số y = f(x) khi biết
giá trị của hàm này tại các mốc x
i
. Ta biết:
y
i
= f(x
i
) (3.1)
Ta có thể dùng công thức nội suy Lagrange để tính đạo hàm:
f’(x) ≈ L
n
’(x) (3.2)
với ước lượng sai số:

()
()
()
()
()
1
'
1!
n
nn
fc
d

R
xx
dx n
ω
+
⎛⎞
=
⎜⎟
⎜⎟
+
⎝⎠
(3.3)
Vì điểm c phụ thuộc x nên ước lượng (3.3) chỉ đánh giá được khi x là các mốc nội
suy x = x
i
;
Thông thường người ta xét đa thức nội suy với mốc cách đều với h = x
i+1
– x
i
.
1.1 Tính đạo hàm cấp 1
a) Đạo hàm tại các điểm biên
Khi x là điểm biên x0 hoặc xn ta dùng công thức nội suy bậc nhất với hai mốc nội
suy để tính gần đúng đạo hàm:
y’(x
0
) = (y
1
-y

0
)/h (3.4)
y’(x
n
) = (y
n
-y
n-1
)/h
Vì y
n
= y
n-1
+ y’(x
n
) h + 0(h
2
) nên sai số của ước lượng (3.4) là O(h
2
).
b) Đạo hàm tại các điểm trong
Khi x = x
i
là các điểm trong (i = 1,2, ,n-1) ta dùng công thức nội suy bậc 2 có x
i
là
điểm giữa
()
()
2

11 1
1
2
ii i
tt
y
xy ty y
−− −
−
≈+Δ+ Δ (3.5)
với x = x
i-1
+h.t
Đạo hàm (3.5) theo x ta được:
()
''' 2
11
2. 1 1

2
tx i i
t
yx yt y y
h
−−
−
⎛⎞
=≈Δ+ Δ
⎜⎟
⎝⎠

thay x = x
i
hay t = 1 vào công thức trên ta được:

() ()
'2
11 1 1
11 1 1

22
ii i i ii
yx y y y y y
hh
−− − −
⎛⎞⎛ ⎞
≈Δ + Δ =Δ + Δ −Δ
⎜⎟⎜ ⎟
⎝⎠⎝ ⎠

() ()
'
1
1
2
iii
y
xyy
h
−

≈Δ+Δ
hay

()
'
11
2
ii
i
yy
yx
h
+−
−
≈
(3.6)
với ∀i = 1,2,…,n-1.
Để tính ước lượng sai số ta có các công thức:
()
()
2
'''3
1
2
'''3
1
2
2
iii i
iii i

h
y
yhy yOh
h
y
yhy yOh
+
−
⎧
=+ + +
⎪
⎪
⎨
⎪
=− + +
⎪
⎩

Do vậy:

()
'2
11
2
ii
i
yy
y
Oh
h

+−
−
=+

hay công thức (3.6) có sai số là O(h
2
).
1.2 Đạo hàm cấp 2.
Để tính đạo hàm cấp 2 ta dùng công thức nội suy cấp 2 để tính y’’(x
i
). Đạo hàm
hai lần liên tiếp biểu thức (3.5) ta có:
() ()
'' 2
11 1
22
11
2
iiiii
yx y y y y
hh
−+ −
≈Δ = − + (3.7)
ta có các công thức sau:
()
()
23
'''34
1
23

'''33
1
26
26
iii i i
iii i i
hh
yyhy y yOh
hh
yyhy y yOh
+
−
⎧
=+ + + +
⎪
⎪
⎨
⎪
=− + − +
⎪
⎩

từ đó ta có:

()
()
'' 2
11
2
1

2
iiii
y
yy yOy
h
+−
−+ =+
Vậy sai số có bậc O(h
2
).
Chú ý:
 Chúng ta đã có công thức tính đạo hàm cấp 1 và cấp 2 tại các mốc nội suy.
Để tính đạo hàm tại các điểm không là mốc ta lại áp dụng phương pháp nội
suy Lagrange.
 Sai số khi tính đạo hàm ngoài sai số của công thức còn phải tính đến sai số
làm tròn, và các bước nội suy h phải đủ nhỏ.
Ví dụ: Hàm y = f(x) được cho tại các mốc sẽ có đạo hàm cấp 1 và cấp 2 tại các
mốc này được tính và cho trong bảng sau:
i x
i
y
i
y’
i
y
i
’’
0
1
2

3
4
5
1,0
1,1
1,2
1,3
1,4
1,4
1,266
1,326
1,393
1,469
1,553
1,647
0,6
0,635
0,715
0,8
0,89
0,94

0,7
0,9
0,8
1,0
II.TÍNH GẦN ĐÚNG TÍCH PHÂN XÁC ĐỊNH
2.1 Công thức hình thang
Giả sử chúng ta biết giá trị của hàm y = f(x) tại các mốc cách đều x
i

trên đoạn
[a,b]. Hãy lập công thức tính tích phân hàm f(x) trên [a,b] qua các giá trị tại mốc.
Chia [a,b] thành n phần bằng nhau. Khí đó ta có:
h= (b-a)/n; x
0
=a; x
n
=b; x
i
= a+i.h; y
i
= f(x
i
); (3.9)
Công thức hình thang dựa trên ý tưởng sau.Trên mỗi đoạn [x
i
, x
i+1
] ta thay diện
tích hình thang cong bởi diện tích hình thang tương ứng. Điều đó có nghĩa là:

()
1
1
2
i
i
x
ii
x

y
y
f
xdx h
+
+
+
=
∫
(3.10)








x
n
= b x
0
= a x
i
x
i+1

y
i

y
i+1


y=f(x)

Lấy tổng trên các đoạn [x
i
, x
i+1
] (i=0 ; n-1) ta có:

()
1
1
0
.
2
b
n
ii
i
a
y
y
f
xdx h
−
+

=
+
≈
∑
∫

hay

() ()
01 1
2 2
2
b
nn
a
ba
f
xdx y y y y
n
−
−
≈++++
∫
(3.11)
Ứớc lượng sai số:
Thực chất của công thức (3.11) là thay hàm f(x) trên đoạn Δx
i
bởi công thức nội
suy bậc nhất của nó trên đoạn này. Với i = 0 ta có:

() ()()
10
00
10
yy
f
xy xx Rx
xx
−
=+ − +
−


()
()
()() ()()
''
2
01 01
22
fc
M
R
x xxxx xxxx=−−≤−−
với
()
''
2
max
M

fx= ; với mọi x∈[a,b]
Vậy sai số của tích phân trên đoạn Δx
0
là
() ()()
11
00
3
10
22
01
.
22 12
xx
xx
yy
M
Mh
fxdx h xx xxdx
+
−≤−−=
∫∫

trên n đoạn ta có sai số toàn phần:

()
()
2
3
2

2
1
.
12 12
M
bah
Mh
Rn n
−
≤=
(3.12)
Ví dụ: Tính
2
1
0
x
edx
−
∫
. Ta lập bảng giá trị của hàm
i x
i
x
2
i
y
i

0
1

2
3
4
5
6
0,0
0,1
0,2
0,3
0,4
0,5
0,6
0,00
0,01
0,04
0,09
0,16
0,25
0,36
1,0000
0,9900
0,9608
0,9139
0,8521
0,7788
0,6977
7
8
9
10

0,7
0,8
0,9
1,0
0,49
0,64
0,81
1,00
0,6126
0,5273
0,4449
0,3679
Với
()
()
2
'' 2
22 1
x
y
xxe
−
=−
()
9
010
1
1
7, 462
2

i
i
yy y
=
++ =
∑

Vậy

2
1
0
0,7462
x
edx
−
≈
∫

()
''
f
x
đạt max tại x = 0 là M
2
= 2. Vậy

()
2
10

2. 0,1
0,002
12
R ≤<
Nên sau khi làm tròn ta có:

2
1
0
0,746
x
edx
−
≈
∫

2.2 Công thức Simson (Công thức parabol)
a) Xây dựng công thức
Chia đoạn [a,b] thành 2n đoạn bằng nhau, khi đó h=(b-a)/2n; Trên mỗi đoạn [x
2i
,
x
2(i+1)
] thay hàm f(x) bởi công thức nội suy bậc hai và diện tích hình thang cong
giới hạn bởi ham f(x) bởi diện tích hình thang cong giới hạn bởi parabol nội suy.

x
0
=a
x
2i

x
2i+1
x
2i+2
y = f(x)
Ta có:

()
()
2
22 2
1
2
ii i
tt
f
xyty y
−
≈+Δ+ Δ

với

2i
x
x
t
h
−
=

nên

() ()
22
2
22122
4
3
i
i
x
iii
x
h
fxdx y y y
+
++
≈++
∫
(3.13)
Lấy tổng theo i = 0, , n - 1 ta được:

() ()
01234 212
4 2 4 2 4
6
b
nn
a
ba
f
xdx y y y y y y y
n
−
−
≈ ++++++ +
∫
(3.14)
b). Ước lượng sai số
Người ta đã chứng minh công thức ước lượng sai số như sau:

()
()
4
24
2
180
ba
R
nMh
−
=

(3.15)
trong đó

()
()
4
4
ax f
M
mx= với x∈ [a,b]
Ví dụ. tính
2
1
0
x
edx
∫
.Chia đoạn [0,1] thành 10 phần bằng nhau. Khi đó ta có 2n =
10. Các giá trị của hàm
2
x
ye= cho trong bảng sau:
2
x
ye
=

i
x
i

2
i
x

i = 0 và i = 10 i lẻ i chẵn
0
1
2
3
4
5
6
7
8
9
10
0,0
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1,0
0,00
0,01

0,04
0,09
0,16
0,25
0,36
0,49
0,64
0,81
1,00
1,0000









2,7189

1,0101

1,0942

1,2840

1,6329

2,2479

1,0408

1,1725

1,4333

1,8965


Đạo hàm 4 lần liên tiếp ta được:

()
()
2
4
42
44 12 3
x
yxxe=++

Hàm này đạt giá trị cực đại tại x = 1 và M
2
= 76.e
Vậy:

4
2

76
.0,1 0,000115 0,00012
180
e
R
≤≈ <

()
2
1
0
1
3, 7188 4.7, 2685 2.5,4441 1,46268 1, 4627
30
x
edx≈++=≈
∫