Luận văn thạc sỹ toán học "Giải gần đúng phương trình phi tuyến và phương trình vi phân trên máy tính điện tử luận văn thạc sĩ toán học"
Trang 1ĐẠI HỌC THÁI NGUYÊN
TRƯỜNG ĐẠI HỌC SƯ PHẠM
TRẦN THỊ HOÀN
GIẢI GẦN ĐÚNG PHƯƠNG TRÌNH PHI TUYẾN
VÀ PHƯƠNG TRÌNH VI PHÂN TRÊN MÁY TÍNH ĐIỆN TỬ
LUẬN VĂN THẠC SĨ TOÁN HỌC
THÁI NGUYÊN - 200
Trang 2ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC SƯ PHẠM
LUẬN VĂN THẠC SĨ TOÁN HỌC
Người hướng dẫn khoa học: TS Tạ Duy Phượng
THÁI NGUYÊN - 2007
Trang 31
MỤC LỤC
Trang
Lời nói đầu 2-3
Chương 1 Giải gần đúng phương trình phi tuyến trên máy tính điện tử……… …… ………… ………4
Chương 2 Giải gần đúng nghiệm của bài toán Cauchy cho phương trình vi phân thường trên máy tính điện tử …48
Đ1 Phương pháp giải gần đúng bài toán Cauchy cho phương trình vi phân thường……….….……… 48
Đ2 Phương pháp Euler ………… ……… …… ….…52 Đ3 Phương pháp Runge-Kutta ………… ……… ….…57 Đ4 Giải bài toán Cauchy cho phương trình vi phân trên máy tính điện tử
………… ……….……… ……… 64
Kết luận 82 Tài liệu tham khảo 83
Trang 4LỜI NÓI ĐẦU
Các bài toán thực tế (trong thiên văn, đo đạc ruộng đất,…) dẫn đến việc cần phải giải các phương trình phi tuyến (phương trình đại số hoặc phương trình vi phân), tuy nhiên, các phương trình này thường phức tạp, do đó nói chung khó có thể giải được (đưa được về các phương trình cơ bản) bằng các biến đổi đại số Hơn nữa, vì các công thức nghiệm (của phương trình phi tuyến hoặc phương trình vi phân) thường phức tạp, cồng kềnh, nên cho dù có công thức nghiệm, việc khảo sát các tính chất nghiệm qua công thức cũng vẫn gặp phải rất nhiều khó khăn Vì vậy, ngay từ thời Archimedes, các phương pháp giải gần đúng đã được xây dựng Nhiều phương pháp (phương pháp Newton-Raphson giải gần đúng phương trình phi tuyến, phương pháp Euler và phương pháp Runge-Kutta giải phương trình vi phân) đã trở thành kinh điển và được sử dụng rộng rãi trong thực tế
Với sự phát triển của công cụ tin học, các phương pháp giải gần đúng lại càng có ý nghĩa thực tế lớn Để giải một phương trình bằng tay trên giấy, có khi phải mất hàng ngày với những sai sót dễ xảy ra, thì với máy tính điện tử, thậm chí với máy tính điện tử bỏ túi, chỉ cần vài phút Tuy nhiên, việc thực hiện các tính toán toán học trên máy một cách dễ dàng càng đòi hỏi người sử dụng có hiểu biết sâu sắc hơn về lí thuyết toán học Mặt khác, nhiều vấn đề lí thuyết (sự hội tụ, tốc độ hội tụ, độ chính xác, độ phức tạp tính toán,…) sẽ được soi sáng hơn trong thực hành tính toán cụ thể Vì vậy, việc sử dụng thành thạo công cụ tính toán là cần thiết cho mọi học sinh, sinh viên Công cụ tính toán sẽ hỗ trợ đắc lực cho việc tiếp thu các kiến thức lí thuyết, giảng dạy lí thuyết gắn với thực hành tính toán, sẽ giúp học sinh, sinh viên không chỉ tiếp thu tốt hơn các kiến thức khoa học, mà còn tiếp cận tốt hơn với các phương pháp và công cụ tính toán hiện đại
Nói chung, trong các trường phổ thông và đại học hiện nay, việc gắn giảng dạy lí thuyết với tính toán thực hành còn chưa được đẩy mạnh Điều này hoàn toàn không phải vì thiếu công cụ tính toán, mà có lẽ là vì việc phổ biến cách sử dụng các công cụ tính toán còn ít được quan tâm
Với mục đích minh họa khả năng sử dụng máy tính điện tử trong dạy và học
môn Giải tích số, chúng tôi chọn đề tài luận văn Giải gần đúng phương trình phi
Trang 53
tuyến và phương trình vi phân trên máy tính điện tử Luận văn gồm hai chương:
Chương 1 trình bày ngắn gọn các phương pháp giải gần đúng phương trình phi tuyến và đặc biệt, minh họa và so sánh các phương pháp giải gần đúng phương trình
thông qua các thao tác thực hành cụ thể trên máy tính điện tử khoa học Casio fx-570 ES Chương 2 trình bày phương pháp Euler, phương pháp Euler cải tiến và phương
pháp Runge-Kutta giải phương trình vi phân thường Các phương pháp này được so
sánh và minh họa qua thực hành tính toán trên máy tính Casio fx-570 ES và trên chương trình Maple
Có thể coi các qui trình và chương trình trong luận văn là các chương trình mẫu để giải bất kì phương trình phi tuyến hoặc phương trình vi phân nào (chỉ cần khai báo lại phương trình cần giải) Điều này đã được chúng tôi thực hiện trên rất nhiều phương trình cụ thể
Tác giả xin chân thành cám ơn TS Tạ Duy Phượng (Viện Toán học), người Thầy đã hướng dẫn tác giả hoàn thành luận văn này Xin được cảm ơn Trường Đại học Sư phạm (Đại học Thái Nguyên), nơi tác giả đã hoàn thành chương trình cao học dưới sự giảng dạy nhiệt tình của các Thầy Xin được cám ơn Phòng Giáo dục Phổ Yên (Thái Nguyên), nơi tác giả công tác, đã tạo mọi điều kiện thuận lợi để tác giả hoàn thành khóa học và luận văn Cuối cùng, xin được cám ơn Gia đình đã động viên, giúp đỡ và chia xẻ những khó khăn với tác giả trong thời gain học tập
Thái Nguyên, 20.9.2007
Trần Thị Hoàn
Trang 6CHƯƠNG I
GIẢI GẦN ĐÚNG PHƯƠNG TRÌNH PHI TUYẾN TRÊN MÁY TÍNH ĐIỆN TỬ
Đ1 GIẢI GẦN ĐÚNG PHƯƠNG TRÌNH f x( ) 0
Phương trình f x ( ) 0 thường gặp nhiều trong thực tế Tuy nhiên, ngoài một số lớp phương trình đơn giản như phương trình bậc nhất, phương trình bậc hai, phương trình bậc ba và bậc bốn là các phương trình có công thức nghiệm biểu diễn qua các hệ số, và một vài lớp phương trình được giải nhờ các kĩ thuật của đại số (phân tích ra thừa số, đặt ẩn phụ,…) để đưa về các phương trình bậc nhất hoặc bậc hai, hầu hết các phương trình phi tuyến là không giải được chính xác (không có công thức biểu diễn nghiệm qua các hệ số của phương trình), vì vậy người ta thường tìm cách tìm nghiệm gần đúng của phương trình Và ngay cả khi biết công thức nghiệm, do tính phức tạp của công thức, giá trị sử dụng của công thức nhiều khi cũng không cao Thí dụ, ngay cả với lớp phương trình đơn giản là phương trình đa thức bậc ba ax3bx2cx d 0, mặc dù có công thức Cardano để giải, nhưng vì công thức này chứa nhiều căn thức khá cồng kềnh (xem, thí dụ:
Eric W Weisstein: CRS Concise Encyclopedia of Mathematics, CRS Press, New
York, 1999, mục Cubic Equation, trang 362-365),
nên thực chất chúng ta cũng chỉ có thể tìm được nghiệm gần đúng Hơn nữa, đa số các phương trình, thậm chí những phương trình rất đơn giản về mặt hình thức nhưng lại xuất phát từ các bài toán thực tế, thí dụ, phương trình xcosx không có công thức biểu diễn nghiệm thông qua các phép toán cơ bản (cộng, trừ, nhân, chia, khai căn, lũy thừa), nói cách khác, không giải được hoặc rất khó giải bằng các phép biến đổi đại số, nhưng có thể giải gần đúng đến độ chính xác bất kì rất dễ dàng nhờ phép lặp xn1cosxn, nhất là trên máy tính điện tử bỏ túi (chỉ cần bấm liên tiếp một phím )
Những phương trình xuất hiện trong các bài toán thực tế (thí dụ, khi đo đạc,…) nói chung có thông tin đầu vào (thể hiện trên các hệ số, trong công thức) chỉ
Trang 75
là gần đúng (sai số trong đo đạc, đánh giá, tính toán sơ bộ, ) Vì vậy việc tìm nghiệm chính xác cũng không có ý nghĩa thực tế lớn, trong khi đó với các phương pháp giải gần đúng phương trình, ta thường có công thức đánh giá độ chính xác của nghiệm gần đúng và có thể tìm nghiệm đến độ chính xác bất kì cho trước, nên phương pháp giải gần đúng phương trình có ý nghĩa rất quan trọng trong giải quyết các bài toán thực tế
Các phương pháp giải chính xác phương trình chỉ mang tính đơn lẻ (cho từng lớp phương trình), còn các phương pháp giải gần đúng phương trình mang tính phổ dụng: một phương pháp có thể dùng để giải cho những lớp phương trình rất rộng, thí dụ, chỉ đòi hỏi hàm số là liên tục chẳng hạn, vì vậy khả năng ứng dụng của giải gần đúng là rất cao
Giải gần đúng phương trình liên quan đến nhiều vấn đề quan trọng khác của toán học Thí dụ, theo điều kiện cần cực trị (Định lí Fermat), điểm x0 là điểm cực trị (địa phương) của hàm số y F x ( ) thì nó phải là điểm dừng, tức là
'( )'( ) 0
y xF x Như vậy, để tìm điểm cực trị, trước tiên ta phải giải phương trình y'F x'( ) : f x( ) 0 để tìm điểm dừng (điểm được nghi ngờ là điểm cực trị) Trong thực tế để tìm nghiệm tối ưu, ta thường đi tìm các điểm dừng (nghi ngờ là cực trị) nhờ giải gần đúng phương trình y'F x'( ) : f x( ) 0
Bởi vì một trong những thế mạnh của máy tính điện tử là khả năng lặp lại một công việc với tốc độ cao, mà giải gần đúng phương trình thực chất là việc thực hiện một dãy các bước lặp, nên nhờ máy tính mà việc giải gần đúng phương trình trở nên đơn giản, nhanh chóng và thuận tiện Không những thế, máy tính còn cho phép, thông qua lập trình, mô phỏng quá trình thực hiện bước lặp giải phương trình, bởi vậy nó là công cụ tốt trợ giúp học sinh và sinh viên tiếp thu các kiến thức toán học nói chung, các phương pháp giải gần đúng phương trình nói riêng Do đó thực hành giải gần đúng trên máy tính điện tử có một ý nghĩa nhất định trong giảng dạy và học tập bộ môn toán trong các trường phổ thông và đại học
Trong chương này, để giải gần đúng phương trình, chúng ta luôn giả thiết rằng, f x( ) là một hàm xác định và liên tục trên một đoạn nào đó của đường thẳng
Trang 8thực Nhiều khi điều kiện này đã là đủ để xây dựng phương pháp giải gần đúng Trong một số phương pháp, ta sẽ giả thiết rằng f x( ) khả vi đến cấp cần thiết (có đạo hàm cấp một hoặc có đạo hàm cấp hai)
Nếu f x( ) 0 thì điểm x được gọi là nghiệm hoặc không điểm của
phương trình f x( ) 0 Ta cũng giả thiết rằng các nghiệm là cô lập, tức là tồn tại một lân cận của điểm x không chứa các nghiệm khác của phương trình Khoảng lân cận (chứa x ) này được gọi là khoảng cách li của nghiệm x
Các bước giải gần đúng phương trình
Giải gần đúng phương trình f x( ) 0 được tiến hành theo hai bước:
Bước 1 Tìm khoảng chứa nghiệm
Một phương trình nói chung có nhiều nghiệm Ta cần tìm khoảng chứa nghiệm, tức là khoảng ( , )a b trong đó phương trình có nghiệm (có duy nhất nghiệm), bằng một trong các tiêu chuẩn sau
Định lí 1 (Bolzano-Cauchy) Nếu hàm f x( ) liên tục trên đoạn a b, và thỏa mãn điều kiện f a f b( ) ( ) 0 thì phương trình f x( ) 0 có ít nhất một nghiệm trong
khoảng ( , )a b
Ý nghĩa hình học của Định lí này khá rõ ràng: Đồ thị của một hàm số liên tục là một đường cong liên tục (liền nét), khi chuyển từ điểm A a f a( , ( )) sang điểm
( , ( ))
B b f b nằm ở hai phía khác nhau của trục hoành, đường cong này phải cắt trục
hoành tại ít nhất một điểm (có thể tại nhiều điểm)
Thí dụ, hàm số yf x( )x33x1 có f( 2) 3; f( 1) 1 ; f(0) 1 và
(2) 1
f nên phương trình x3 3x 1 0 có ba nghiệm phân biệt trong các khoảng ( 3, 1) ; ( 1,0) và (0,2)
Trang 97
Định lí 2 (Hệ quả của Định lí 1) Giả sử f x( ) là một hàm liên tục và đơn điệu chặt trên đoạn a b, Khi ấy nếu f a f b( ) ( ) 0 thì phương trình f x( ) 0 có duy nhất một nghiệm trong khoảng ( , )a b
Ý nghĩa hình học của Định lí này là: Đồ thị của một hàm số liên tục tăng chặt (giảm chặt) là một đường cong liên tục (liền nét) luôn đi lên (đi xuống) Khi di chuyển từ điểm A a f a( , ( )) sang điểm B b f b( , ( )) nằm ở hai phía khác nhau của trục hoành thì đồ thị phải cắt và chỉ cắt trục hoành một lần (Hình vẽ)
Hai định lí trên chỉ đòi hỏi tính liên tục mà không đòi hỏi tính khả vi (tồn tại đạo hàm) của f x( ) Nếu f x( ) có đạo hàm thì có thể dùng tiêu chuẩn dưới đây
Trang 10Định lí 3 (Hệ quả của Định lí 2) Giả sử hàm số f x( ) có đạo hàm f x( ) và đạo hàm f x( ) của nó không đổi dấu (luôn dương hoặc luôn âm) trên đoạn a b, Khi ấy nếu f a f b( ) ( ) 0 thì phương trình f x( ) 0 có duy nhất một nghiệm trong khoảng ( , )a b
Từ ba định lí trên, ta đi đến hai phương pháp tìm khoảng cách li nghiệm của phương trình f x( ) 0 (khoảng chứa duy nhất một nghiệm): phương pháp hình học và phương pháp giải tích
Phương pháp giải tích
Giả sử ta phải tìm nghiệm của phương trình f x( ) 0 trong khoảng ( , )a b Ta đi tính giá trị f a( ), f b( ) và các giá trị f x( )i của hàm số tại một số điểm
( , )
xa b , i1,2, ,n Nếu hàm f x( ) đơn điệu chặt trên khoảng x xi, i1 và điều kiện f x f x( ) (ii1) 0 được thỏa mãn thì x xi, i1 là một khoảng cách li nghiệm của phương trình f x( ) 0 Nếu thông tin về hàm f x( ) quá ít thì ta thường dùng quy trình chia đoạn thẳng (chia khoảng ( , )a b thành 2, 4, 8,…phần) và thử điều kiện f x f x( ) (ii1) 0 để tìm khoảng cách li nghiệm
Một đa thức bậc n có không quá n nghiệm Vì vậy phương trình đa thức có không quá n khoảng cách li nghiệm
Khi hàm f x( ) đủ tốt (có đạo hàm, có dạng cụ thể, ), ta có thể khảo sát đồ thị để chia trục số thành các khoảng đổi dấu của đạo hàm (khoảng đồng biến và nghịch biến của hàm số) và xác định khoảng cách li nghiệm
Phương pháp hình học
Trong trường hợp đồ thị hàm số tương đối dễ vẽ, ta có thể vẽ phác đồ thị để tìm khoảng cách li nghiệm hoặc giá trị thô của nghiệm như là giao điểm (gần đúng) của đồ thị với trục hoành Cũng có thể dùng các máy tính đồ họa (máy tính có khả
năng vẽ hình như Casio Algebra fx-2.0 Plus hoặc Sharp EL-9650) hoặc các phần
Trang 119
mềm tính toán (Maple, Matlab,…) để vẽ đồ thị Sau đó, nhờ tính toán, ta “tinh
chỉnh” để đi đến khoảng cách li nghiệm chính xác hơn
Bước 2 Giải gần đúng phương trình
Có bốn phương pháp cơ bản giải gần đúng phương trình: phương pháp chia đôi, phương pháp lặp, phương pháp dây cung và phương pháp tiếp tuyến (phương
pháp Newton-Raphson) Nhằm làm cơ sở lí thuyết cho các tính toán trong Đ3,
trong Đ2 chúng tôi sẽ vắn tắt trình bày nội dung của các phương pháp này, chủ yếu
là dựa vào các giáo trình Giải tích số [1] - [6]
Đ2 CÁC PHƯƠNG PHÁP TÌM NGHIỆM GẦN ĐÚNG CỦA PHƯƠNG TRÌNH ( ) 0f x
1 Phương pháp chia đôi
Nội dung của phương pháp chia đôi rất đơn giản: Giả sử f x( ) là một hàm
liên tục trên đoạn a b, và f a f b( ) ( ) 0 Khi ấy theo Định lí Bolzano-Cauchy, phương trình f x( ) 0 có ít nhất một nghiệm trong khoảng ( , )a b
Chia đôi đoạn a b, và tính ( )2
a b
Nếu ( ) 02
a b
2 a b
x là một nghiệm của phương trình f x( ) 0
Nếu ( ) 02
a b
a hoặc ( , )2
a bb Gọi khoảng mới (khoảng nhỏ) chứa nghiệm là ( , )a b1 1
Lại chia đôi khoảng ( , )a b1 1 và tính giá trị tại điểm giữa 11
Tiếp tục mãi quá trình này ta đi đến:
Trang 12Hoặc tại bước thứ n nào đó ta có ( ) 02
b a
Sự hội tụ của phương pháp chia đôi
Dãy an là dãy đơn điệu tăng, bị chặn trên bởi b, dãy bn là đơn điệu
giảm và bị chặn dưới bởi a nên cả hai dãy đều có giới hạn
Do
Nếu chọn nghiệm gần đúng là x b n thì
2
Trang 1311
Như vậy, sau bước thứ n, nên chọn nghiệm gần đúng là
Nếu tại mỗi bước n ta đều chọn
Do đó khi tính toán (trên máy tính bỏ túi với màn hình hiển thị được 10 chữ số chẳng hạn), ta có thể dừng tính toán khi xn1xn xn1 đúng đến số thập phân cần thiết (thí dụ, ta có thể dừng tính toán khi được nghiệm chính xác đến 10 chữ số, tức là 1010)
2 Phương pháp lặp
Giả sử ( , )a b là khoảng cách li nghiệm của phương trình f x( ) 0 Giải phương trình f x( ) 0 bằng phương pháp lặp gồm các bước sau:
Bước 1 Đưa phương trình f x( ) 0 về phương trình tương đương x g x ( )
Bước 2 Chọn x0( , )a b làm nghiệm gần đúng đầu tiên
Bước 3 Thay x x 0 vào vế phải của phương trình x g x ( ) ta được nghiệm gần đúng thứ nhất x1g x( )0 Lại thay x1g x( )0 vào vế phải của phương trình
( )
x g x ta được nghiệm gần đúng thứ hai x2g x( )1 Lặp lại quá trình trên, ta nhận được dãy các nghiệm gần đúng
1( )0
xg x , x2 g x( )1 , x3g x( )2 , x4 g x( )3 , ,xn g x( n1),
Trang 14Nếu dãy các nghiệm gần đúng xn , n1,2, hội tụ, nghĩa là tồn tại lim
x g x tương đương với phương trình f x( ) 0 trên đoạn a b, Nếu g x( ) và
Trang 16(1 )lg
xxqxx ta có kết luận: nếu dãy xn hội tụ
thì khi n đủ lớn hai nghiệm gần đúng xn và xn1 xấp xỉ bằng nhau Vì vậy khi sử dụng máy tính ta thường dừng quá trình lặp khi các kết quả liên tiếp xn1, xn,
1
Trang 17Khi ấy với mỗi x0 a b; , dãy xn xây dựng theo phương pháp lặp xn g x( n1)sẽ hội tụ tới điểm bất động (tức là x g x( )) duy nhất x trong khoảng ( , )a b của
ánh xạ g
3 Phương pháp dây cung
Giả sử ( , )a b là khoảng cách li nghiệm Ta thay cung của đường cong
( )
yf x trên đoạn [ , ]a b bằng dây trương cung ấy và coi giao điểm của dây cung (đường thẳng) với trục hoành là nghiệm xấp xỉ của phương trình f x( ) 0
Để xây dựng dãy xấp xỉ xn, ta xét hai trường hợp:
f(b) x1
x
Trang 18Thay khoảng ( , )a b bằng khoảng ( , )x b1 , ta đi đến nghiệm 11
( )()( )( )
f b b a
Nghiệm x bây giờ nằm trong khoảng ( , )a x1
Thay ( , )a b bằng khoảng ( , )a x1 , ta đi đến nghiệm xấp xỉ
( )()( )( )
x
Trang 1917 Ta có thể tổng kết thành một công thức như sau:
( )()( )( )
nên hội tụ Hơn nữa, chuyển qua giới hạn trong công thức
( )()( )( )
f x x dxx
Trang 20nnnnn
Trang 2119
4 Phương pháp tiếp tuyến (Phương pháp Newton-Raphson)
Giả sử ( , )a b là khoảng cách li nghiệm Ta thay cung của đường cong
Trang 22nghiệm xấp xỉ: 1 ( )'( )
f a
f a
Nghiệm x bây giờ nằm trong khoảng ( , )x b1 (xem hình 5)
Thay khoảng ( , )a b bằng khoảng ( , )x b1 , ta đi đến nghiệm 1
( )'( )
nnn
Trang 2321
Thật vậy, chuyển qua giới hạn trong biểu thức 1 ( )'( )
'( )
f xx x
f x Suy ra f x( ) 0 Do ( , )a b là khoảng cách li nghiệm nên x x
chính là nghiệm ban đầu
Đánh giá sai số
Giả sử 0m1f x'( ) và f x''( )M1 Khi ấy ta có đánh giá sai số:
max( ) ,[ , ]( )n
nên f x( n1)f x'( n1)(xn xn1) 0 Thay vào đẳng thức trên ta được:
21
Trang 24Từ công thức trên và công thức ( )
Như vậy, tốc độ hội tụ của phương pháp tiếp tuyến là bậc hai
Đ3 TÌM NGHIỆM GẦN ĐÚNG CỦA PHƯƠNG TRÌNH
( ) 0
Như ta thấy trong Đ2, cả bốn phương pháp giải gần đúng phương trình
Do đã được thiết kế sẵn cho phép thực hiện các thao tác tính toán liên tiếp,
việc thực hiện dãy lặp trên các máy tính điện tử khoa học (Casio fx-500 MS, Casio fx-570 MS, Sharp 506 WM, Casio fx-500 ES, Casio fx-570 ES) là khá đơn giản và
thuận tiện Trong Đ này, chúng tôi trình bày cách sử dụng các loại máy này cho mục đích giải gần đúng phương trình theo các phương pháp đã trình bày ở mục trên Thực hành giải gần đúng phương trình trên máy tính điện tử khoa học cho phép cảm nhận rõ hơn các vấn đề của giải tích số (sự hội tụ, tốc độ hội tụ,…) của từng phương pháp
Để tiện trình bày, chúng tôi chọn máy tính điện tử khoa học Casio 570 ES, là
loại máy có nhiều ưu điểm trong các tính năng giải toán và được sử dụng tương đối phổ biến hiện nay trong các trường phổ thông và đại học
Máy tính điện tử khoa học Casio 570 ES có một số phím rất tiện dùng trong
tính toán Nó được thiết kế để có thể tính toán đại số (tính toán theo công thức) kêt hợp với những ô nhớ (9 ô nhớ) Thiết kế này đặc biệt thích hợp cho thực hiện dãy lặp, do đó đặc biệt thuận tiện cho việc giải gần đúng phương trình
Bài 1 Giải phương trình đại số bậc cao x9 x 10 0
Trang 2523
Đây tuy chỉ là một phương trình đa thức, tuy nhiên bậc của nó khá cao nên khó có thể giải được bằng các kĩ thuật của đại số (đặt ẩn phụ, nhóm số hạng,…) để đưa về phương trình bậc thấp hơn
Đặt yf x( )x9 x10
Do y' 9x8 1 0 với mọi x nên hàm số đồng biến trên toàn trục số Ta dễ dàng tính được f(1) 8 và f(2) 504 0 Do đó, phương trình x9 x 10 0 có duy nhất một nghiệm trong khoảng (1;2) Để so sánh, ta sẽ tìm nghiệm gần đúng của phương trình này theo cả bốn phương pháp đã trình bày trong Đ2
Phương pháp chia đôi khoảng chứa nghiệm
Đưa giá trị x1 vào ô nhớ A : 1 SHIFT STO A Đưa giá trị x2 vào ô nhớ B : 2 SHIFT STO B
Chú ý: Từ nay về sau, để cho tiện, các phím số được viêt như là các số, còn các
phím chữ trên màn hình được để trong các ô vuông Thí du, phím số 2 ta vẫn viết là số 2, còn phím ô A ta viết là A
Khai báo công thức f x( )x9 x10:
ALPHA Xx^ 9 ALPHA X 10 (3.1) Để tính giá trị của biểu thức tại một điểm nào đó ta bấm phím CALC (Calculate-tính): CALC
Máy hỏi (hiện trên màn hình): X?
Khai báo 32
c (đang ở trong ổ C ): ALPHA C (15331512 )
Chú ý: Từ nay về sau, để tiện trong trình bày, ta ghi ngay đáp số của kết quả tính
toán hiển thị trên màn hình sau phím và để ở trong ngoặc Thí dụ, sau khi khai
Trang 26báo công thức (3.1) và bấm phím CALC , máy hỏi X?, ta khai báo ALPHA C
(tức là giá trị của đối số x cần tìm giá trị của hàm số f x( )x9 x10 chính bằng giá trị trong ổ C ) Sau khi bấm phím , máy hiện ngay trên màn hình đáp số 15331
512 , tức là 315331( )
f và ta viết (15331512 )
a ba c Ô B được giải phóng (giá trị cũ
x trong ô nhớ B không cần dùng nữa)
Chú ý: Vì tại mỗi bước i ta chỉ cần nhớ ba giá trị a bi, i hoặc
c nên ta cũng chỉ cần sử dụng ba ô nhớ A , B , C là đủ
Lại chia đôi khoảng 1 1 3( , ) ( , ) (1; )
B (vừa được giải phóng):
Sử dụng phím để đi về dòng công thức (3.1) và dùng phím CALC để tính giá trị của hàm số tại c11,25:
Bấm phím CALC (Calculate-tính): CALCMáy hỏi: X?
Khai báo 111
Trang 2725
Chú ý: Để tiện trình bày, ta thường viết dấu bằng thay cho dấu gần đúng , thí
( ) 1,2994194034
f , mặc dù chính xác hơn, phải viết là:
( ) 1,2994194034
f (chính xác đến 10 chữ số) Chứng tỏ nghiệm của phương trình nằm trong khoảng
5 3( , ) ( , ) ( ; )
Sử dụng phím để đi về dòng công thức (3.1) và dùng phím CALC để tính giá trị của hàm số tại 2 11
4 8
a ba c Ô nhớ B được giải phóng (giá trị
cũ 2 32
b trong ô nhớ B không cần dùng nữa)
Trang 28Lại chia đôi khoảng 3 3 2 2 5 11( , ) ( , ) ( ; )
4 8
5 1121
Sử dụng phím để đi về dòng công thức (3.1) và dùng phím CALC để tính giá trị của hàm số tại 3 21
4 16
a bc c Ô nhớ C được giải phóng (giá trị
cũ 118
c trong ô nhớ C không cần dùng nữa)
Lại chia đôi khoảng 4 4 1 3 5 21( , ) ( , ) ( ;)
4 16
5 2141
Sử dụng phím để đi về dòng công thức (3.1) và dùng phím CALC để tính giá trị của hàm số tại 4 41
1,2812532
Trang 2927 Bấm phím CALC (Calculate-tính): CALCMáy hỏi: X?
Khai báo 4 41
4 32
a bc c Ô nhớ C được giải phóng (giá trị
cũ 3 2116
c trong ô nhớ C không cần dùng nữa)
Tiếp tục quá trình này, ta đi đến bảng sau, trong đó mỗi dòng cho ta các giá trị cn, ( )n
f c và giá trị f a( )n hoặc f b( )n đã tính được ở dòng trên, trái dấu với f c( )n Vì các giá trị f a( )n hoặc f b( )n đã được tính đến 10 chữ số ở dòng trên, và chỉ cần dấu của chúng, nên để tiết kiệm chỗ, ta chỉ viết giá trị khác 0 đầu tiên của nó sau dấu phảy Mặt khác, nhìn vào mỗi dòng, ta cũng biết ô nào (không có) được giải phóng, nhường chỗ để đưa giá trị mới của cnvào
Trang 30bình: Ta thấy phương pháp chia đôi tuy đơn giản nhưng nói chung hội tụ chậm,
thao tác trên máy khá phức tạp
Phương pháp lặp
Trang 31Sau 6 lần bấm phím, ta đi đến dãy nghiệm xấp xỉ:
1,5; 1,268436614; 1,27223043; 1,272168998; 1,272169993; 1,272169977; 1,272169977;
Kết quả nghiệm gần đúng là: 1,272169977 (chính xác đến 10 chữ số thập phân)
Phương pháp dây cung
Xét f x( )x9 x10 Ta có f x'( ) 9x8 1 0 và f x( ) 72x7 0với mọi x trong khoảng (1;2) Theo công thức dây cung
( )()( )( )
nnnn
Trang 32Khai báo công thức
Phương pháp tiếp tuyến
Theo công thức tiếp tuyến 1 ( )'( )
Khai báo giá trị ban đầu: 2
Khai báo công thức
Ans ( Ans x^ 9 Ans 10 ) ( 9 Ans x^ 8 1 )
Tính các giá trị xn1 bằng cách bấm liên tiếp phím Ta được dãy các giá trị xấp xỉ:
1,781344902; 1,592631378; 1,438653785; 1,331291548; 1,281547657; 1,272435766; 1,272170196; 1,27216977; 1,27216977
Sau chín lần lặp ta đã đi đến đáp số
Kết luận: Cả bốn phương pháp (chia đôi đoạn chứa nghiệm, lặp, dây cung và tiếp
tuyến) đều cho nghiệm gần đúng tới 10 chữ số thập phân Để tìm được nghiệm gần
Trang 3331
đúng đến 10 chữ số thập phân theo phương pháp chia đôi cần 30 bước lặp, và thao tác phức tạp trong mỗi bước Để tìm nghiệm gần đúng đến 10 chữ số thập phân theo phương pháp lặp chỉ cần 6 bước Để tìm nghiệm gần đúng đến 10 chữ số thập phân theo phương pháp dây cung phải cần đến 300 bước lặp Để tìm nghiệm gần đúng đến 10 chữ số thập phân theo phương pháp tiếp tuyến chỉ cần 9 bước lặp Thao tác lặp theo ba phương pháp sau rất đơn giản: sau khi khai báo giá trị ban đầu và công thức lặp, ta chỉ cần liên tiếp bấm phím cho đến khi được giá trị không đổi (điểm bất động), là nghiệm gần đúng đên 10 chữ số thập phân
Bài 2 Giải phương trình chứa căn thức 3x28 x 5 0
Đặt y f x( ) 3 x28 x5 Không cần máy ta cũng có thể dễ dàng tính được f(0) 5 và f(5) 3.5 2 5 5 10 2 5 0 8 8 Do đó, phương trình
3x2 x 5 0 có ít nhất một nghiệm trong khoảng (0;5) Để xem dáng điệu
hình học của đồ thị, ta có thể nhờ Maple để vẽ đồ thị trong khoảng (0;5) nhờ lệnh
plot (vẽ đồ thị):
[> plot(3*x-2*x^(1/8)-5,x=0 5);
Nhìn vào đồ thị hàm số y f x( ) 3 x28 x5 ta thấy, đồ thị cắt trục hoành tại
một điểm trong khoảng (2;3) Tuy nhiên, để chính xác hơn ta cần tính giá trị của hàm số tại các điểm x 2 và x3
Ta có thể tính các giá trị của hàm số đã cho tại điểm x2 nhờ lệnh eval (evaluate - tính giá trị của hàm số tại x2) và lệnh evalf(%) (tính giá trị của biểu thức trên):
Trang 34[> eval(3*x-2*x^(1/8)-5,x=2);
1 2 2(1 8/ )
[> evalf(%);
-1.181015466Cũng có thế kết hợp hai lệnh trên vào một lệnh: > evalf(eval(3*x-2*x^(1/8)-5,x=2));
4 2 3(1 8/ )
Tất nhiên, ta cũng có thể dùng máy tính khoa học để tính các giá trị trên Thí dụ, sử
dụng Casio fx-570 ES để tính giá trị của hàm số 3x28 x 5 tại x2 như sau: Khai báo hàm số 3x28 x 5:
Trang 3533
3 ALPHA X 2 ALPHA Xx 1 8 5 (3.2) Tính giá trị hàm số:
Bấm phím CALC (Calculate-tính) Máy hỏi: X?
Khai báo x 2: 2
Máy trả lời (đáp số hiện trên màn hình), màn hình có dạng:
Khai báo x 3: 3
Máy trả lời (đáp số hiện trên màn hình):
Phương pháp chia đôi khoảng chứa nghiệm
Đưa giá trị x2 vào ô nhớ A : 2 SHIFT STO A Đưa giá trị x3 vào ô nhớ B : 3 SHIFT STO B
Sử dụng phím để đi về dòng công thức (3.2) và dùng phím CALC để tính giá trị của hàm số tại 5
c :
Trang 36Bấm phím CALC (Calculate-tính) Máy hỏi: X?
Khai báo 52
Sử dụng phím để đi về dòng công thức (3.2) và dùng phím CALC để tính giá trị của hàm sô tại 1 9
Bấm phím CALC (Calculate-tính): Máy hỏi: X?
Khai báo 1 94
c (đang ở trong ổ B ): CALC
f bf Chứng tỏ nghiệm của phương trình nằm trong
khoảng 2 2 1 1 9 5( ; ) ( , ) ( ; )
4 2
a ba c Ô nhớ A được giải phóng (giá trị cũ a2
trong ô nhớ A không cần dùng nữa)
Trang 37Sử dụng phím để đi về dòng công thức (3.2) và dùng phím CALC để tính giá trị của hàm số tại 2 19
f bf Chứng tỏ nghiệm của phương trình nằm trong
( ; ) ( ; ) (2,375; )2
a bc b Ô nhớ B được giải phóng (giá trị cũ 1 94
c
trong ô nhớ B không cần dùng nữa)
( ; ) ( ; ) (2,375; )2
19 539
Sử dụng phím để đi về dòng công thức (3.2) và dùng phím CALC để tính giá trị của hàm số tại 3 39
Bấm phím CALC (Calculate-tính): CALC
Trang 38Máy hỏi: X?Khai báo 3 39
c (đang ở trong ổ B ):
ALPHA B (0,076879549)
Như vậy, ta có ( )3 (39) 0,076879549 016
( )( )(2,375)0.103373306 08
trình nằm trong khoảng ( ; ) ( ; ) (4 4 2 3 19 39; )8 16
a b c c Ô nhớ C được giải phóng (giá trị cũ 5
c trong ô nhớ C không cần dùng nữa) Lại chia đôi khoảng ( ; ) ( ; ) (4 4 2 3 19 39; )
8 16
a b c c và gửi
19 3977
Sử dụng phím để đi về dòng công thức (3.2) và dùng phím CALC để tính giá trị của hàm số tại 4 77
2,4062532
Trang 3937 trình nằm trong khoảng ( ; ) ( ; ) (4 4 2 3 19 39; )
8 16
a b c c Ô nhớ C được giải phóng (giá trị cũ 5
c trong ô nhớ C không cần dùng nữa) Ta có bảng kết quả sau:
Trang 40Dãy lặp có dạng 1 28 53
Chọn giá trị ban đầu x0 2,5: 2.5
Khai báo ( )x0 : ( 2 Ans x 1 8 5 ) 3 Tính giá trị của ( )x0 : Bấm phím
Tính các giá trị tiếp theo: Bấm phím Ta được dãy xấp xỉ:
2.5; 2.414235595; 2.410980682; 2.410855171; 2.410850329; 2.410850142; 2.410850135; 2.410850134; 2.410850134; …
Sau 9 bước ta đi đến nghiệm gần đúng đến 10 chữ số