phuong phap tinh

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 (469.78 KB, 72 trang )

(1)Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. GIỚI THIỆU MÔN HỌC I. GIỚI THIỆU CHUNG Phương pháp tính là một lĩnh vực của toán học chuyên nghiên cứu các phương pháp giải các bài toán (chủ yếu là gần ñúng) bằng cách dựa trên những dữ liệu số cụ thể và cho kết quả cũng dưới dạng số. Ngày nay phần lớn các công việc tính toán ñều ñược thực hiện trên máy tính. Tuy vậy thực tế chứng tỏ rằng việc áp dụng các thuật toán và phương pháp tính toán khác nhau có thể cho tốc ñộ tính toán và ñộ chính xác rất khác nhau. Lấy ví dụ ñơn giản như tính ñịnh thức của ma trận chẳng hạn, nếu tính trực tiếp theo ñịnh nghĩa thì việc tính ñịnh thức của một ma trận vuông cấp 25 cũng mất hàng triệu năm (ngay cả với máy tính hiện ñại nhất hiện nay); trong khi ñó nếu sử dụng phương pháp khử Gauss thì kết quả nhận ñược gần như tức thời. Như vậy, phương pháp tính là công cụ không thể thiếu trong các công việc cần thực hiện nhiều tính toán với tốc ñộ tính toán nhanh và ñộ chính xác cao như vật lý, ñiện tử viễn thông, công nghệ thông tin.... Phương pháp tính ñược nghiên cứu từ rất lâu và cho ñến nay những thành tựu ñạt ñược là một khối lượng kiến thức ñồ sộ ñược in trong nhiều tài liệu sách, báo... Tuy nhiên, môn học "Phương pháp tính" chỉ nhằm cung cấp những kiến thức căn bản nhất về phương pháp tính. Với lượng kiến thức này sinh viên có thể áp dụng vào giải quyết những bài toán thông thường trong thực tế và có khả năng tự tìm hiểu ñể nâng cao kiến thức cho mình khi gặp các vấn ñề phức tạp hơn. Trong phương pháp tính chúng ta thường quan tâm ñến hai vấn ñề: • Phương pháp ñể giải bài toán. • Mối liên hệ giữa lời giải số gần ñúng và lời giải ñúng, hay sai số của lời giải. II. MỤC ðÍCH Môn học phương pháp tính cung cấp cho sinh viên kiến thức căn bản nhất về một số phương pháp giải gần ñúng trên dữ liệu số với mục ñích • Tạo cơ sở ñể học tốt và nghiên cứu các nghành khoa học kỹ thuật. • Góp phần rèn luyện phương pháp suy luận khoa học, tư duy logic, phương pháp nghiên cứu thực nghiệm và xây dựng thế giới quan khoa học và tác phong khoa học cần thiết cho người kỹ sư tương lai.. 1.

(2) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. III. PHẠM VI NGHIÊN CỨU Nghiên cứu một số phương pháp cơ bản nhất của phương pháp tính, ñược ứng dụng nhiều trong thực tế như các phương pháp tính trong ñại số tuyến tính, bài toán nội suy, tìm nghiệm gần ñúng các phương trình phi tuyến, tính gần ñúng ñạo hàm và tích phân, giải gần ñúng một số dạng của phương trình vi phân... Tìm hiểu các lĩnh vực ứng dụng của các phương pháp trong thực tế. IV. PHƯƠNG PHÁP HỌC TẬP: ðể học tốt môn học này, sinh viên cần lưu ý những vấn ñề sau: 1. Kiến thức chuẩn bị: • Sinh viên phải có kiến thức cơ bản về toán học cao cấp. • Thành thạo sử dụng máy tính cầm tay (sẽ ñược giảng viên hướng dẫn trên lớp) 2. Tài liệu và dụng cụ học tập: • Giáo trình Phương pháp tính của trường ðHCN Tp HCM. • Máy tinh cầm tay (Casio 570 MS, ES hoặc Vinacal 570 MS) Nếu cần sinh viên nên tham khảo thêm: • Giải tích số. Phạm Kỳ Anh, nhà xuất bản ñại học Quốc Gia Hà Nội, 1966. • Phương pháp tính. Tạ Văn ðỉnh, Nhà xuất bản Giáo dục - 1995. • Phương Pháp tính. Dương Thuỳ Vỹ, Nhà xuất bản Khoa học và Kỹ thuật, 2001. 3. Tham gia ñầy ñủ các buổi hướng dẫn học tập: Thông qua các buổi hướng dẫn học tập, giảng viên sẽ giúp sinh viên nắm ñược nội dung tổng thể của môn học và giải ñáp thắc mắc, ñồng thời sinh viên cũng có thể trao ñổi, thảo luận với những sinh viên khác về nội dung bài học. 4. Chủ ñộng liên hệ với bạn học và giảng viên: Cách ñơn giản nhất là tham dự các diễn dàn học tập trên mạng Internet, qua ñó có thể trao ñổi trực tiếp các vấn ñề vướng mắc với giảng viên hoặc các bạn học khác ñang online. ðịa chỉ email ñể trao ñổi với giảng viên : 5. Tự ghi chép lại những ý chính: Việc ghi chép lại những ý chính là một hoạt ñộng tái hiện kiến thức, kinh nghiệm cho thấy nó giúp ích rất nhiều cho việc hình thành thói quen tự học và tư duy nghiên cứu. 6. Học ñi ñôi với hành Học lý thuyết ñến ñâu thực hành làm bài tập ngay ñến ñó ñể hiểu và nắm chắc lý thuyết.. 2.

(3) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. CHƯƠNG 1 SỐ XẤP XỈ VÀ SAI SỐ MỤC ðÍCH, YÊU CẦU Sau khi nghiên cứu chương 1, yêu cầu sinh viên: 1. Hiểu ñược sai số tuyệt ñối và sai số tương ñối. 2. Nắm ñược cách viết số xấp xỉ. 3. Nắm ñược các qui tắc tính sai số. 4. Hiểu và biết cách ñánh giá sai số tính toán và sai số phương pháp . 1.1. ðẶT VẤN ðỀ Trong kỹ thuật giá trị các thông số chúng ta tiếp cận nói chung không phải là giá trị ñúng (vì nó là kết quả của các phép ño và thí nghiệm). Như vậy chúng ta ñã sử dụng giá trị gần ñúng thay cho giá trị ñúng, việc này nẩy sinh nhiều vấn ñề phức tạp vì giá trị ñúng chỉ có một nhưng giá trị gần ñúng thì rất nhiều. ðể có cơ sở khoa học trong việc sử dụng các số gần ñúng người ta ñưa ra khái niệm sai số ñể ño ñộ chênh lệch giữa các giá trị ñúng và giá trị gần ñúng. Chú ý rằng khi sử dụng số gần ñúng thay cho một số ñúng nào ñó người ta luôn phải dùng ñồng thời hai ñại lượng ñó là : giá trị gần ñúng và sai số. Hai ñại lượng này có vai trò như nhau. 1.2. SAI SỐ TUYỆT ðỐI VÀ SAI SỐ TƯƠNG ðỐI 1.2.1. Sai số tuyệt ñối Xét ñại lượng ñúng A và ñại lượng gần ñúng của nó là a. Ta nói a xấp xỉ A và viết a ≈ A. Trị tuyệt ñối ∆ = |a-A| ñược gọi là sai số tuyệt ñối của a (khi dùng a ñể xấp xỉ A). Trong thực tế ta không biết ñược số ñúng A, do ñó nói chung sai số tuyệt ñối không tính ñược. Vì vậy ta tìm cách ước lượng sai số tuyệt ñối của a bằng số ∆a > 0 sao cho |a - A| ≤ ∆a Số dương ∆a ñược gọi là sai số tuyệt ñối giới hạn của a. Chú ý: Nếu ∆a là sai số tuyệt ñối giới hạn của a thì mọi số thực lớn hơn ∆a ñều là sai số tuyệt ñối giới hạn của a nhưng nếu sai số tuyệt ñối giới hạn quá lớn so với sai số tuyệt ñối thì nó không còn có ý nghĩa về phương diện sai số nữa. Trong những ñiều kiện cụ thể người ta cố gắng chọn ∆a là số dương bé nhất.. 3.

(4) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. 1.2.2. Sai số tương ñối Gọi ∆ là sai số tuyệt ñối của a khi dùng a ñể xấp xỉ A, khi ñó ñại lượng δ=. ∆ A. ñược gọi là sai số tương ñối của a. Tuy nhiên một lần nữa ta thấy rằng A thường không biết, vì vậy người ta ñịnh nghĩa ñại lượng δa =. ∆a a. là sai số tương ựối giới hạn của a. đôi khi người ta biểu diễn sai số tương ựối dưới dạng %. Ví dụ . với a =10, ∆a = 0.05, khi ñó ta có δ a =. 0.05 = 0.5% . 10. Vì trong thực tế chúng ta chỉ có thể thao tác với các sai số giới hạn, do ñó người ta thường gọi một cách ñơn giản ∆a là sai số tuyệt ñối, δ a là sai số tương ñối. 1.2.3. Chú thích: Sai số tuyệt ñối không nói lên ñầy ñủ "chất lượng" của một số xấp xỉ, “chất lượng” ấy còn ñược phản ánh qua sai số tương ñối. 1.3. CÁCH VIẾT SỐ XẤP XỈ 1.3.1. Chữ số có nghĩa Một số viết dưới dạng thập phân có thể gồm nhiều chữ số, nhưng ta chỉ kể các chữ số từ chữ số khác không ñầu tiên tính từ trái ñến chữ số cuối cùng khác không phía bên phải là các chữ số có nghĩa. Chẳng hạn số 2.740 có 3 chữ số có nghĩa, số 0.02078 có 4 chữ số có nghĩa. 1.3.2. Chữ số ñáng tin Mọi số thập phân ñều có dạng n. a = ± ( α n 10 n + … + α110 + α 0100 + α −110−1 + … + α − m 10− m ) = ± ∑ α s , α s ∈ {0,1,..., 9} . s =− m. Giả sử a là xấp xỉ của số A với sai số tuyệt ñối là ∆a. s. Nếu ∆a ≤ 0.5 × 10 thì ta nói rằng chữ số αs là ñáng tin (như vậy các chữ số có nghĩa bên trái αs ñều là ñáng tin). s. Nếu ∆a > 0.5 × 10 thì ta nói rằng chữ số αs là ñáng nghi (như vậy các chữ số bên phải αs ñều là ñáng nghi). Ví dụ. Cho số xấp xỉ a = 4.67329 hãy xác ñịnh các chữ số ñáng tin và các chữ số ñáng ngờ khi ∆a = 0.004726 hoặc ∆a= 0.005726. 4.

(5) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Giải -2. Ta có ∆a = 0.004726 ≤ 0.5 × 10 do ñó các chữ số ñáng tin là: 4,6,7; các chữ số ñáng ngờ là 3,2, 9. -1. Khi ∆a= 0.005726 ta có ∆a ≤ 0.5 × 10 do ñó các chữ số ñáng tin là: 4,6; các chữ số ñáng ngờ là 7, 3, 2, 9. 1.3.3. Cách viết số xấp xỉ a. Kèm theo sai số Nếu ∆a là sai số tuyệt ñối giới hạn của a khi xấp xỉ A thì ta quy ước viết: A = a ± ∆a với ý nghĩa a – ∆a ≤ A ≤ a + ∆a Hoặc A = a(1 ± δ a ) b. Mọi chữ số có nghĩa ñều ñáng tin Cách thứ hai là viết theo quy ước: mọi chữ số có nghĩa ñều ñáng tin; có nghĩa là sai số tuyệt ñối giới hạn không lớn hơn một nửa ñơn vị ở hàng cuối cùng. Ví dụ. Khi viết a = 4.67329 thì ta hiểu lúc này ∆a= 0.5 × 10. -5. 1.4. CÁC LOẠI SAI SỐ KHI XỬ LÝ BÀI TOÁN KỸ THUẬT Khi giải một bài toán phức tạp người ta thường thay bài toán ñó bằng bài toán ñơn giản hơn ñể có thể giải ñược bằng tay hoặc bằng máy. Phương pháp thay bài toán phức tạp bằng một phương pháp ñơn giản tính ñược như vậy gọi là phương pháp gần ñúng. Sai số do phương pháp gần ñúng tạo ra gọi là sai số phương pháp. Mặc dầu bài toán ñã ở dạng ñơn giản, nhưng trong quá trình giải ta thường xuyên phải làm tròn các kết quả hoặc xử dụng các số xấp xỉ , sai số tạo ra trong quá trình này gọi là sai số tính toán. Trong thực tế việc ñánh giá các loại sai số, nhất là sai số tính toán nhiều khi là bài toán rất khó thực hiện. Tóm lại khi thực hiện một bài toán bằng phương pháp gần ñúng ta thường gặp những loại sai số sau ñây: • Sai số trong việc mô hình hóa bài toán : xuất hiện do việc giả thiết bài toán ñạt ñược một số ñiều kiện lý tưởng nhằm làm giảm ñộ phức tạp của bài toán. • Sai số phương pháp : xuất hiện do việc giải bài toán bằng phương pháp gần ñúng. • Sai số của số liệu : xuất hiện do việc ño ñạc và cung cấp giá trị ñầu vào không chính xác. • Sai số tính toán : xuất hiện do làm tròn số hoăc xử dụng các số xấp xỉ trong quá trình tính 5.

(6) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. toán, quá trình tính càng nhiều thì sai số tích luỹ càng lớn. Những sai số trên ñây tổng hợp lại nhiều khi dẫn ñến những lời giải quá cách xa so với lời giải ñúng và vì vậy không thể dùng ñược. Chính vì vậy việc tìm ra những thuật toán hữu hiệu ñể giải các bài toán thực tế là ñiều rất cần thiết. 1.5. SAI SỐ TÍNH TOÁN THƯỜNG GẶP 1.5.1. Sai số quy tròn các số xấp xỉ Khi tính toán với các con số ta thường làm tròn các số theo quy ước: Nếu chữ số bỏ ñi ñầu tiên ≥ 5 thì thêm vào chữ số giữ lại cuối cùng một ñơn vị, còn nếu chữ số bỏ ñi ñầu tiên < 5 thì ñể nguyên chữ số giữ lại cuối cùng. Giả sử a là xấp xỉ của A với sai số tuyệt ñối giới hạn là ∆. Giả sử ta quy tròn a thành a' với sai số quy tròn tuyệt ñối giới hạn là θ, tức là: | a' - a| ≤ θ. Khi ñó |a' - A| = | a' - a + a -A| ≤ | a' - a| + | a -A| ≤ θ + ∆ Vậy có thể lấy θ + ∆ làm sai số tuyệt ñối giới hạn của a'. Như vậy việc quy tròn làm tăng sai số tuyệt ñối giới hạn. 1.5.2. Sai số khi tính toán trên các số xấp xỉ Bài toán Cho u = f(x1, x2,..., xn) . Biết các ñối số x1, x2,..., xn là các số xấp xỉ với các sai số tuyệt ñối tương ứng là ∆1 , ∆2 , ... ∆n và f là hàm khả vi liên tục theo các ñối số xi. Hãy xác ñịnh ∆u, δ u . Giải Theo công thức vi phân của hàm nhiều biến ta có: du = Từ ñây suy ra ∆u ≈. ∂u ∂u dx1 + … + dxn ∂x1 ∂xn. ∂u ∂u ∂u ∂u ∆x1 + … + ∆xn ≤ ∆1 + … + ∆n ∂x1 ∂xn ∂x1 ∂xn. Vì vậy có thể chọn : ∆u =. ∂u ∂u ∆1 + … + ∆n ∂x1 ∂xn. ðể tìm δ u ta dùng công thức : δ u =. ∆u u. 6.

(7) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Ví dụ. Cho hàm u = f ( x, y, z ) = x 2 y + yz. Hãy xác ñịnh giá trị hàm số u, sai số tuyệt ñối và. sai số tương ñối của u biết x = 0.983, y = 1.032(1 ± 0.05), z = 2.114 ± 0.02. ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ......................................................................................................................................... 7.

(8) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. BÀI TẬP. Trong các bài tập dưới ñây chúng ta ngầm hiểu sai số tương ñối và sai số tuyệt ñối là sai số tương ñối giới hạn và sai số tuyệt ñối giới hạn o. o. Bài 1. Khi ño 1 số góc ta ñược các giá trị : a= 21 37’3”; b=1 10’ . Hãy xác ñịnh sai số tương ñối của các số xấp xỉ ñó biết rằng sai số tuyệt ñối trong phép ño là 1”. Bài 2. Hãy xác ñịnh sai số tuyệt ñối của các số xấp xỉ sau ñây cho biết sai số tương ñối của. chúng: a) a= 13267 ; δa=0,1%. b) b=2.32; δb=0.7%. Bài 3. Hãy xác ñịnh số các chữ số ñáng tin trong các số a,b với sai số như sau:. a) a= 0,3941; ∆a=0,25.10. -2. b) b=38,2543; ∆a= 0,27.10. -2. Bài 4. Hãy xác ñịnh số những chữ số ñáng tin trong các số a với sai số tương ñối như sau:. a) a=1,8921; δa=0,1.10. -2. b) a=22,351; δa=0,1. Bài 5. Hãy qui tròn các số dưới ñây( xem là ñúng) với 3 chữ số có nghĩa ñáng tin và xác ñịnh sai số tuyệt ñối ∆ và sai số tương ñối δ của chúng:. a) a= 2,514. b) 0,16152. c) 0,01204. d) –0,0015281. Bài 6. Hãy xác ñịnh: Giá trị của các hàm số, Sai số tuyệt ñối giới hạn và Sai số tương ñối. giới hạn. Biết giá trị của các ñối số cho với mọi chữ số có nghĩa ñều ñáng tin.. a ) u = f ( x, y, z ) = tg ( x 2 y + yz ), b) u = f ( x, y, z ) = zesin( xy ) ,. x = 0.983, y = 1.032, z = 2.114. x = 0.133, y = 4.732, z = 3.015.. Bài 7. Hãy xác ñịnh: Giá trị của các hàm số, sai số tuyệt ñối và sai số tương ñối. Biết giá trị. của các ñối số cho với mọi chữ số có nghĩa ñều ñáng tin.. a ) u = x 2 sin( yz ),. x = 1.113; y = 0.102; z = 2.131.. b) u = z ln( xy ) ,. x = 0.162; y = 4.531; z = 1.91.. c) u = 2 x + 2 y ,. x = 0.085; y = 0.055; z = 2.152.. d ) u = (1 + xyz ) x ,. x = 2.918; y = 1.032; z = 2.114.. 2. Bài 8. Tính thể tích V của hình cầu và chỉ ra sai số tuyệt ñối, biết rằng ñường kính ño ñược. d = 1,112m và sai số của phép ño là 1 mm. Bài 9. Hãy xác ñịnh sai số tương ñối giới hạn và sai số tuyệt ñối giới hạn và chữ số ñáng tin. của cạnh hình vuông a. Biết rằng diện tích hình vuông là S = 16, 45cm 2 , ∆ S = 0, 01.. 8.

(9) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. CHƯƠNG 2 TÍNH GẦN ðÚNG NGHIỆM CỦA PHƯƠNG TRÌNH PHI TUYẾN MỤC ðÍCH, YÊU CẦU Sau khi học xong chương 2, yêu cầu sinh viên:. 1. Kiểm tra ñược khoảng cách ly nghiệm 2. Tìm ñược nghiệm gần ñúng và ñánh giá ñược sai số của phương trình phi tuyến. 3. Biết vận dụng các phương pháp trên vào các bài toán thực tế. 2.1. GIỚI THIỆU CHUNG 2.1.1. ðặt vấn ñề. Khi giải quyết bài toán kỹ thuật chúng ta thường gặp loại yêu cầu : Xác ñịnh thông số ñầu vào, ñể ñầu ra của một hệ thống nào ñó ñạt một mức cho trước. yêu cầu này có thể phát biểu bằng ngôn ngữ toán học như sau: Xác ñịnh giá trị x ∈ ( a, b) sao cho f ( x) = 0 ,. (2.1). Như chúng ta ñã biết việc giải phương trình (2.1) không ñơn giản (vì không có phương pháp chung) ngay cả khi f ( x) là ña thức có bậc lớn hơn 3. Trong kỹ thuật người ta có thể chấp nhận giá trị x∗ (sao cho f ( x∗ ) ≈ 0 ) thay cho nghiệm ñúng α của phương trình nhưng với ñiều kiện ñánh giá ñược sai số tuyệt ñối giữa x∗ và α. ðiều này cũng hoàn toàn hợp lý bởi thực tế ngay cả khi chúng ta xác ñịnh ñược chính xác. giá trị thông số ñầu vào thì khi qua hệ thống kết quả ñầu ra cũng chỉ ñạt ñược kết quả gần bằng với yêu cầu. Giá trị x∗ nói ở trên gọi là nghiệm gần ñúng của phương trình (2.1). Việc ñi tìm giá trị x∗ và ñánh giá sai số gọi là giải gần ñúng phương trình. Chú ý: Khi ñánh giá sai số chúng ta cần phải tính. ∆ x* = x * −α và ∆ f ( x*) = f ( x *) − f (α ) = f ( x *) .. {. }. Sai số chung của bài toán ñược tính bởi ∆ = max ∆ x*; ∆ f ( x*) . Trong bài giảng chỉ tập chung tính ∆ x* = x * −α . 2.1.2. Các bước giải gần ñúng phương trình phi tuyến. Khi giải gần ñúng nghiệm của phương trình (2.1) ta cần tuân thủ các bước sau: Bước 1: Kiểm tra (2.1) có nghiệm ñúng duy nhất trên [a,b] (hay [a,b] là khoảng cách ly) Bước 2: Dùng các thuật toán ñể tìm giá trị x∗ và ñánh giá sai số 9.

(10) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. 2.1.3. Một số ñịnh lý cần thiết trong việc thực hiện giải gần ñúng ðể thực hiện bước 1, 2 ta dùng ñịnh lý dưới ñây ðịnh lý1.. Nếu hàm số f(x) liên tục, ñơn ñiệu trên ñoạn [a,b] và f(a)f(b)<0 thì ñoạn [a,b] là một khoảng cách ly nghiệm của phương trình (2.1). ðịnh nghĩa2.. Gọi S = { x : x − x0 ≤ C} là một lân cận ñóng của x0 ∈ R , A là ánh xạ từ S vào S. Ta nói A là ánh xạ co trên S nếu tồn tại hằng số q < 1 sao cho. ∀x, y ∈ S :. A( x) − A( y ) ≤ q x − y .. ðịnh lý3.. Giả sử α là nghiệm ñúng của phương trình x = A( x) và tồn tại lân cận ñóng S của α sao cho A là ánh xạ co trên S thì α là nghiệm duy nhất của phương trình x = A( x) trên S và. α có thể thu ñược bằng cách lấy giới hạn của dãy xn +1 = A( xn ); n = 0,1,.. với x0 là một ñiểm nào ñó thuộc S. ðịnh lý4.. Với hàm f(x) liên tục và khả vi trên ñoạn [a,b], ngoài ra tồn tại m sao cho 0 < m ≤ |f'(x)| với mọi x thuộc [a,b] khi ñó ta có ñánh giá: xn − α ≤. f ( xn ) . m. 2.2. MỘT SỐ THUẬT TOÁN TÌM GIÁ TRỊ GẦN ðÚNG NGHIỆM 2.2.1. Phương pháp lặp ñơn a. Mô tả phương pháp - Giả sử phương trình (2.1) có khoảng cách ly nghiệm là [a,b]. - Biến ñổi (2.1) ñược về dạng tương ñương x = ϕ ( x) ( hàm ϕ gọi là hàm lặp) - Chọn x0 =. a+b . Tính các nghiệm xấp xỉ xn +1 theo công thức 2 xn +1 = ϕ ( xn ), n = 0,1, 2,.... - đánh giá sai số ∆ n = xn − α . đặt ∆n ≤. q = max { ϕ '( x) } Ta có: x∈[ a,b ]. q xn − xn −1 1− q 10.

(11) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. b. ðiều kiện hội tụ của phương pháp ðịnh lý. Nếu hàm ϕ ( x) có ñạo hàm ϕ '( x) và thỏa mãn: ϕ '( x) ≤ q < 1, ∀x ∈ [ a, b] thì phương pháp lặp hội tụ, tức là: lim xn = α n →∞. Chú ý.. Khi sử dụng phương pháp lặp chúng ta sẽ gặp khó khăn trong việc tìm hàm ϕ ( x) (vì phải thỏa ñiều kiện : ϕ '( x) ≤ q < 1, ∀x ∈ [ a, b] ). ðể khắc phục ñiều này ta làm theo hướng dẫn - Nếu f '( x) > 0, ∀x ∈ [a, b] ta ñặt ϕ ( x) = x −. f ( x) với M ≥ max { f '( x) } M x∈[ a ,b ]. - Nếu f '( x) < 0, ∀x ∈ [ a, b] ta ñặt ϕ ( x) = x +. f ( x) với M ≥ max { f '( x) } M x∈[ a ,b ].  f (c ) ≠ 0 - Nếu ∃c ∈ [ a, b] :  thì ta thu hẹp ñoạn [ a, b] thành [c + ε , b] hoặc  f '(c) = 0 [ a, c − ε ] với ε là hằng số dương ñủ nhỏ sao cho ñoạn thu hẹp vẫn là ñoạn cách ly nghiệm. Chú ý: Với cách ñặt như trên thì. q = max { ϕ '( x) } = 1 − x∈[ a ,b ]. Trong ñó m = min. x∈[ a ,b ]. m <1 M. { f '( x) } .. Ví dụ1. Tìm xấp xỉ nghiệm trên ñoạn [1,2] của phương trình:. x3 − x − 1 = 0. thỏa yêu cầu sai số 10-1 Giải ðặt f ( x) = x3 − x − 1 suy ra f '( x) = 3 x 2 − 1 , f ''( x) = 6 x.  f ''( x) = 0 ⇔ x = 0 ∉ (1, 2)  Ta có  f '(1) = 2 suy ra  f '(2) = 11 .  f '( x) > 0, ∀x ∈ [1, 2]   max { f '( x) } = 11  x∈[1,2]  { f '( x) } = 2 m = xmin ∈[1,2] . 11.

(12) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Bước 1: Kiểm tra [1,2] là khoảng cách ly nghiệm. f(x) liên tục trên [1,2], Ta có. f (1) f (2) = −5 < 0, f '( x) > 0, ∀x ∈ [1, 2].. . Vậy [1,2] là khoảng cách ly nghiệm.. Bước 2: Tính giá trị nghiệm và ñánh giá sai số. f ( x) x3 − x − 1 − x3 + 12 x + 1 Chọn M=11. ðặt ϕ ( x) = x − = x− = M 11 11 Ta có q = max { ϕ '( x) } = x∈[1,2]. ðặt x0 =. 9 < 1 . Vậy hàm ϕ ( x) thỏa ñiều kiện của phương pháp lặp. 11. 1+ 2 = 1.5 ta tính các giá trị x1, x2… theo công thức lặp dưới ñây 2.  x1 = ϕ ( x0 ) = 1.420455   ;  q − 1 ∆ ≤ x − x ≈ 0.36 > 10  1 1 − q 1 0.  x2 = ϕ ( x1 ) = 1.379947    ∆ ≤ q x − x ≈ 0.18 > 10−1  2 1 − q 2 1.  x3 = ϕ ( x2 ) = 1.357418    ∆ ≤ q x − x ≈ 0.1 > 10−1  3 1 − q 3 2.  x4 = ϕ ( x3 ) = 1.344351    ∆ ≤ q x − x ≈ 5.9 × 10−2 < 10−1  4 1 − q 4 3. ;. Vậy x4 = 1.34451 là nghiệm gần ñúng thỏa yêu cầu về sai số. Ví dụ2. Tìm xấp xỉ nghiệm trên ñoạn [0.5,1] của phương trình: 3x 2 − e x = 0. thỏa yêu cầu sai số 10-2. Ví dụ3. Giải gần ñúng trên [4,5] của phương trình:.  x+1  cos  π  + 0.148x − 0.9062 = 0  8  thỏa yêu cầu sai số 10-2. ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 12.

(13) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ……………………………………………………………………………………………………………………………………………………… ……………………………………………………………………………………………………………………………………………………… .................................................................................................................................................................................................................................... 13.

(14) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. 2.2.2. Phương pháp Newton-Rapson hay ( phương pháp tiếp tuyến ) a. Mô tả phương pháp. - Giả sử phương trình (2.1) có khoảng cách ly nghiệm là [a,b]. - Chọn x0 thuộc [a,b] sao cho f(x0) cùng dấu với f’’(x), ∀x ∈ (a, b) - Tính giá trị của nghiệm gần ñúng thứ n+1 theo công thức x n +1 = x n −. f (x n ) ; n=0,1,2... f '(x n ). - đánh giá sai số ∆ n +1 = x n +1 − α . ðặt. M = max. x∈[a,b]. { f ''(x) },. m = min. x∈[a,b]. { f '(x) }.. Ta có ∆ n +1 ≤. M x n +1 − x n 2m. 2. b. ðiều kiện hội tụ của phương pháp. Chú ý phương pháp Newton-Rapson cũng là một dạng của phương pháp lặp với hàm lặp là ϕ( x ) = x −. f (x) f '(x). Do vậy muốn phương pháp Newton-Rapson hội tụ thì hàm ϕ ( x ) phải thỏa ñiều kiện. ϕ '( x) =. f ( x) f ''( x). [ f '( x)]2. < 1; ∀x ∈ [ a, b]. Việc kiểm tra ñiều kiện trên khá vất vả nên trong thực hành người ta thường sử dụng ñiều kiện ñủ dưới ñây ðịnh lý.. Giả sử hàm f(x) có f’(x) khác không trên ñoạn [a,b] và f''(x) không ñổi dấu trong (a,b). Nếu x0 , xn ñược chọn như trong mục a) thì phương pháp Newton-Rapson hội tụ, tức là: lim xn = α n →∞. 14.

(15) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Ví dụ1. Tìm xấp xỉ nghiệm trên ñoạn [-3,-2] của phương trình: x3 + 3 x 2 − 1 = 0 thỏa yêu cầu sai số 10-3. Giải ðặt f ( x ) = x3 + 3 x 2 − 1 suy ra f '( x) = 3 x 2 + 6 x , f ''( x ) = 6 x + 6 , f (3) ( x ) = 6 Xét f’(x).  f ''( x) = 0 ⇔ x = −1 ∉ (−3, −2)  Ta có  f '(−3) = 9  f '(−2) = 0  Vì f '( −2) = 0 nên ta phải thu hẹp [-3,-2] thành [-3,-2.5].  f ''( x ) = 0 ⇔ x = −1 ∉ (−3, −2.5)  suy ra Khi ñó  f '( −3) = 9  f '(−2.5) = 3.75 .  f '( x ) > 0, ∀x ∈ [ −3, −2.5]    m = min { f '( x) } = 3.75 x∈[ −3, −2.5] . Xét f’’(x).  f (3) ( x) = 6 > 0  Ta có  f ''( −3) = −12 suy ra  f ''( −2.5) = −9 .  f ''( x) < 0, ∀x ∈ [−3, −2.5]    M = max { f ''( x) } = 12 x∈[ −3, −2.5] . Bước 1: Kiểm tra [-3,-2.5] là khoảng cách ly nghiệm. Ta có. f(x) liên tục trên [-3,-2.5], f (−3) f ( −2.5) = −1 × 2.215 < 0,. Vậy [-3,-2.5] là khoảng cách ly nghiệm.. f '( x) > 0, ∀x ∈ [ −3, −2.5]. Bước 2: Tính giá trị nghiệm và ñánh giá sai số. - ðặt x0 = −3 (Vì f ( −3) cùng dấu với f ''( x), ∀x ∈ [ −3, −2.5] ) - Tính các giá trị x1, x2… theo công thức lặp. f ( x0 )   x1 = x0 − f '( x ) = −2.888889 0  ;   ∆1 ≤ M x1 − x0 2 ≈ 0.0198 > 10−3  2m. f ( x1 )  x = x − = −2.879452 2 1  f '( x1 )    ∆ 2 ≤ M x2 − x1 2 ≈ 1.43 × 10−4 < 10−3  2m. Vậy x2 = −2.879452 là nghiệm gần ñúng thỏa yêu cầu sai số. 15.

(16) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. x 2. Ví dụ2. Cho phương trình: 2 ln x − + 1 = 0 . Tìm xấp xỉ nghiệm trên ñoạn [0.2,1] sau 4 lần. lặp. đánh giá sai số khi nhận giá trị nghiệm ở lần lặp thứ tư là xấp xỉ nghiệm. Ví dụ3. Giải gần ñúng trên [0,1] của phương trình: x 2 − cosπx = 0 thỏa yêu cầu sai số 10-4. ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 16.

(17) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 17.

(18) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. BÀI TẬP Bài 1. Dùng một trong hai phương pháp (Lặp hoặc Newton-Rapson) tìm nghiệm gần ñúng của phương trình dưới ñây thỏa yêu cầu sai số 10-4. 1) x3 − x − 1 = 0;. x0 ∈ [1;2]. 7) 3x 2 − e x = 0. ;. 2) x 4 − 3 x 2 − 3 = 0;. x0 ∈ [1;2]. 8) x − cos x = 0. ; x 0 ∈ [0;1]. 3) x 4 − 2 x3 − 4 = 0;. x0 ∈ [ 2;3]. 9) x+ ln x − 5 = 0. ; x 0 ∈ [3;5]. x0 ∈ [ 0.2;1]. 10). 4). x − tgx = 0;. x 5) π +0,5sin = x; 2 6) x − 2. −x. = 0;. x 0 ∈ [0;1]. x − x3 + 1 = 0. 11) esin x − x 4 + 3 = 0. x0 ∈ [ 0;2π ]. 12) tg. x0 ∈ [0.3;1]. x 2x −1 −e − x 2 + 10 = 0 2. x 0 ∈ [1;2]. ; ;. x 0 ∈ [1; 2]. ;. x 0 ∈ [3;4]. Bài 2. Tự tìm khoảng cách ly nghiệm và giải bằng phương pháp lặp ñơn hoặc tiếp tuyến. 1) ( x − 1) 2 =. ex 2. 4) 4 x − 5ln x = 5. 2) x =ln(x+1). 5) x = x − sin( x + lg x − 2). 3) sinx+cosx=4x. 6) x 2 = e x + 2. 2. Bài 3. Tìm nghiệm dương nhỏ nhất của các phương trình : x. 1) e 2 − 2x = 0. 2) 1,8x 2 − sin(10x) = 0. 3) 2 x − 4x = 0. 18.

(19) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. CHƯƠNG 3 GIẢI GẦN ðÚNG HỆ PHƯƠNG TRÌNH TUYẾN TÍNH MỤC ðÍCH, YÊU CẦU: Sau khi nghiên cứu chương 3, yêu cầu sinh viên:. 1. Nắm ñược các xu hướng xử lý các bài toán ñại số tuyến 2. Hiểu và thực hiện ñược các phương pháp tìm nghiệm nghiệm xấp xỉ của hệ phương trình tuyến tính. 3. Biết cách ñánh giá sai số của từng phương pháp 3.1. GIỚI THIỆU CHUNG 3.1.1. ðặt vấn ñề. Khi xác ñịnh giá trị các thông số trong kỹ thuật ñôi khi chúng ta phải giải hệ thống phương trình tuyến tính  a11x1 + a12 x2 a x + a x  21 1 22 2  ⋮ ⋮  ⋮ am1x1 + am 2 x2. + + ⋮ +. ⋯ ⋯ ⋯ ⋯. + a1n xn + a2n xn ⋮ ⋮ + amn xn. = b1 = b2 ⋮ = bm. (I). ðể giải hệ trên người ta ñi theo hai hướng sau: Hướng giải ñúng. Sử dụng các phương pháp giải ñúng ñể tìm ra giá trị chính xác của các nghiệm xj. Một số phương pháp tiêu biểu như : Cramer, Gauss-JordanẦđã ựược khảo sát trong môn toán cao cấp A2, C2. Hướng giải gần ñúng. Sử dụng các phương pháp giải gần ñúng ñể tìm ra giá trị xấp xỉ của các nghiệm xj. Một số phương pháp tiêu biểu như : Lặp ñơn, Seidel,… Nhận xét. Hướng giải ñúng có ưu ñiểm là tìm ra ñược giá trị ñúng của nghiệm trong trường hợp hệ có nghiệm duy nhất và chỉ ra ñược khi nào hệ vô nghiệm hoặc có vô số nghiệm. Tuy nhiên lại rất khó thực hiện trong trường hợp các hệ số aij , bi là các số thập phân. Hướng giải gần ñúng có khuyết ñiểm là chỉ tìm ra ñược giá trị gần ñúng của nghiệm trong trường hợp hệ có nghiệm duy nhất và không chỉ ra ñược khi nào hệ vô nghiệm hoặc có vô số nghiệm. Tuy nhiên lại tỏ ra hiệu quả trong trường hợp các hệ số aij ,bi là các số thập phân. 19.

(20) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Khi mô hình hóa bài toán kỹ thuật bằng hệ phương trình, thường hệ có các hệ số rất lẻ và chỉ có duy nhất nghiệm nên hướng giải gần ñúng chiếm hầu hết khi giải bài toán kỹ thuật. Lưu ý. Các phương pháp giải gần ñúng dưới ñây chỉ giải ñược một số hệ có dạng ñặc biệt. (sẽ ñược chỉ rõ trong từng thuật toán). Nếu không phải là dạng này thì chúng ta phải dùng hướng giải ñúng ñể xử lý. 3.1.2. Các bước giải gần ñúng hệ phương trình tuyến tính. Khi giải gần ñúng nghiệm của hệ phương trình (I) ta cần tuân thủ các bước sau: Bước 1: Kiểm tra (I) có nghiệm ñúng duy nhất Bước 2: Dùng các thuật toán ñể tìm giá trị gần ñúng của nghiệm và ñánh giá sai số 3.1.3. Một số khái niệm toán học cần thiết trong việc thực hiện giải gần ñúng ðể thực hiện bước 1, 2 ta cần nhắc lại và xây dựng một số khái niệm sau. Cho ma trận chữ nhật A cấp m x n:  a11 ⋯ a1n    A= ⋮ ⋯ ⋮  a   m1 ⋯ amn  ðịnh nghĩa1. Ta nói chuẩn của ma trận A là một trong các số sau.  m A 1 = max ∑ aij j =1,.., n   i =1.   . (gọi là chuẩn cột).  n = max  ∑ aij i =1,.., m  j =1 .   . (gọi là chuẩn hàng). A∞. A2=. n m. ∑ ∑ aij. 2. (gọi là chuẩn Euclicd). j =1 i =1.  −1 0 1    Ví dụ. Với A =  4 2 −1  ta có  2 −2 −5   .  A = max{7, 4,7} = 7  1  A ∞ = max{2,7,9} = 9   A 2 = 1 + 0 + 1 + 16 + 4 + 1 + 4 + 4 + 25 = 56. Ghi chú: Người ta thường dùng kí hiệu A chung cho ba chuẩn trên .. Trong không gian véc tơ Rn người ta xây dựng khái niệm chuẩn của véc tơ như sau. 20.

(21) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ðịnh nghĩa2. Trong không gian véc tơ Rn cho vecto x = ( x1, x2 ,..., xn ) . Ta nói chuẩn của vecto x là một trong các số sau.. { } j =1,.., n. x 1 = max. Ghi chú : Khái niệm x. xj. 2. hoặc. x. ∞.  n = ∑ x j  j =1.   hoặc x . 2. =. n. ∑ xj. 2. j =1. của vecto mang ý nghĩa hình học là ñộ dài của vecto ñó. Tính chất của chuẩn (ñọc giáo trình) ðịnh lý4 (Về sự duy nhất nghiệm của hệ (I)).. a11 ⋯ a1n. Xét hệ (I) khi m=n . Nếu. ⋮ ⋯ ⋮ ≠ 0 thì hệ (I) có nghiệm duy nhất an1 ⋯ ann. 3.2. MỘT SỐ THUẬT TOÁN TÌM GIÁ TRỊ GẦN ðÚNG NGHIỆM 3.2.1. Phương pháp lặp ñơn a. Mô tả phương pháp. - Giả sử hệ phương trình (I) có nghiệm duy nhất. - Biến ñổi (I) ñược về dạng  x1 = α11x1 + α12 x2 x = α x + α x  2 21 1 22 2  ⋮ ⋮ ⋮   xn = α n1x1 + α n 2 x2. ðặt. + ⋯ + α1n xn + ⋯ + α 2n xn. + β1 + β2. ⋮ ⋯ ⋮ ⋮ + ⋯ + α nn xn. ⋮ + βn. (II).  x1   α11 ⋯ α1n   β1        x =  ⋮ , α = ⋮ ⋯ ⋮ , β =  ⋮  x  α  β   n  n1 ⋯ α nn   n. Khi ñó (II) ñược viết dưới dạng x = α x + β - Chọn x(0) = β . Tính các xấp xỉ nghiệm x( n +1) theo công thức x( n +1) = α x( n) + β , n = 0,1, 2,... n +1 - đánh giá sai số ∆ n +1 = x( ) − x * với x* là nghiệm ựúng của hệ. ∆ n +1 ≤. α x ( n +1) − x ( n ) 1− α. 21.

(22) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. b. ðiều kiện hội tụ của phương pháp ðịnh lý.. Nếu ma trận α có chuẩn bé hơn 1 thì phương pháp lặp ñơn hội tụ. Ví dụ1. Giải gần ñúng hệ phương trình:. 10x − 2y − 3z = 20  −2x + 20y − 5z = 40  x − 3y + 10z = 8 . (I). Thỏa yêu cầu sau số 10-2 Giải Bước 1: Kiểm tra hệ có nghiệm duy nhất. 10 −2 −3 −2 20 −5 = 1862 ≠ 0 . Vậy hệ có nghiệm duy nhất. Ta có. 1. −3 10. Bước 2: Tính gần ñúng và ñánh giá sai số. Biến ñổi (I) ñược về dạng 0x + 0.2y + 0.3z + 2 x =   y = 0.1x + 0y + 0.25z + 2  z = −0.1x + 0.3y + 0z + 0.8 . (II). ðặt. 0.2 0.3   0  α =  0.1 0 0.25  ,  −0.1 0.3 0  .  2  β =  2  ,  0.8   .  x1    x =  x2  x   3. Khi ñó (II) ñược viết dưới dạng x =αx + β Ta có. ðặt x. α (0). ∞. =max {0.5, 0.35, 0.4} = 0.5 < 1 . Vậy ma trận α thỏa ñiều kiện hội tụ ..  2  = β =  2  . Ta tính nghiệm xấp xỉ x( n +1) , n = 0,1, 2,... theo công thức  0.8   . 22.

(23) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ.  0.2 0.3  2   2   2.64   0  (1)        (0) 0 0.25  2  +  2  =  2.4   x = α x +β =  0.1  −0.1 0.3  0.8   0.8   1.2   0           α∞ 0.5  ∆1 ≤ x (1) − x (0) ≈ × 0.64 = 0.64 > 10−2 ∞ 1− α ∞ 1 − 0.5 .   2.84   (2)   (1)  x = α x +β =  2.564   1.256        α∞  ∆2 ≤ x (2) − x(1) 1− α ∞ . ∞. ≈ 0.2 > 10−2.   2.90516   (4)   (3)  x = α x +β =  2.61026   1.29044     ;    α∞ ∆ 4 ≤ x(4) − x (3) ≈ 0.02 > 10−2 ∞ 1 − α  ∞. V ậy x. (5).   2.8896   (3)   (2)  x = α x +β =  2.598   1.2852     ;    α∞ ∆3 ≤ x (3) − x(2) ≈ 0.05 > 10−2 ∞ 1− α ∞ .   2.909184   (5)   (4)  x = α x +β =  2.613126   1.292562        α∞ ∆5 ≤ x(5) − x(4) ≈ 4 × 10−3 < 10−2 ∞ 1 − α  ∞.  x = 2.909184  hay  y = 2.613126 là xấp xỉ nghiệm thỏa yêu cầu sai số.  z = 1.292562 . Ví dụ2. Giải gần ñúng hệ phương trình:. 19.2x + 2.6y − 1.2z = 20.3  y 4 − 15.3z =  3.7x +  x − 13.5y + z = 8.3 . (I). Thỏa yêu cầu sau số 10-3. 23.

(24) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Ví dụ3. Cho hệ phương trình:. − 21.6y − 1.2z = 2.3  x  y − 15.3z = 7.24  3.7x + 25.7x − 3.5y + z = 8.3 . (I). Tìm nghiệm gần ựúng của hệ sau 4 bước lặp. đánh giá sai số khi nhận xấp xỉ nghiệm này. ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 24.

(25) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 25.

(26) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. 3.2.2. Phương pháp Seidel a. Mô tả phương pháp. - Giả sử hệ phương trình (I) có nghiệm duy nhất. - Biến ñổi (I) ñược về dạng x = α x + β - Chọn x (0) = β . Tính các xấp xỉ nghiệm x( k +1) theo công thức  x( k +1)  1  x2( k +1)   ⋮  ( k +1)  xn. +. α12 x2( k ). + ⋯ +. = α 21x1( k +1). α1n −1xn( k−)1 + α1n xn( k ) + β1. +. α 22 x2( k ). + ⋯ +. α 2n −1xn( k−)1 + α 2n xn( k ) + β 2. ⋮. ⋮. ⋮. ⋮ ⋯ ⋮. ⋮. =. α11x1( k ) ⋮. = α n1x1( k +1). + α n 2 x2( k +1). +. (*). + α nn −1xn( k−+11) + α nn xn( k ) + β n. Nếu phân tích.  α11 α12  α α 22 α =  21  ⋮ ⋮   α n1 α n 2. ⋯ α1n   0 0   ⋯ α 2n   α 21 0 = ⋯ ⋮   ⋮ ⋮   ⋯ α nn   α n1 α n 2 =. L. ⋯ 0   α11 α12   ⋯ 0   0 α 22 + ⋯ ⋮  ⋮ ⋮   ⋯ 0  0 0 +. ⋯ α1n   ⋯ α 2n  ⋯ ⋮   ⋯ α nn  U. Thì (*) sẽ ñươc viết lại là x( k +1) = Lx( k +1) + Ux( k ) + β ⇔ x( k +1) = ( I n − L) −1 Ux( k ) + β   . (**). n +1 - đánh giá sai số ∆ n +1 = x( ) − x * với x* là nghiệm ựúng của hệ. ∆ n +1 ≤. L x ( n +1) − x ( n ) 1− α. Chú ý : Theo biểu diễn (**) ta nhận thấy rằng phương pháp Seidel cũng là một dạng của. phương pháp lặp ñơn . b. ðiều kiện hội tụ của phương pháp ðịnh lý.. Nếu ma trận α có chuẩn bé hơn 1 thì phương pháp lặp Seidel hội tụ.. 26.

(27) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Ví dụ1. Giải gần ñúng hệ phương trình bằng phương pháp Seidel. 10x − 2y − 3z = 20  −2x + 20y − 5z = 40  x − 3y + 10z = 8 . (I). Thỏa yêu cầu sau số 10-2 Giải Bước 1: Kiểm tra hệ có nghiệm duy nhất. 10 −2 −3 −2 20 −5 = 1862 ≠ 0 . Vậy hệ có nghiệm duy nhất. Ta có. 1. −3 10. Bước 2: Tính gần ñúng và ñánh giá sai số (Theo công thức (**)). Biến ñổi (I) ñược về dạng 0x + 0.2y + 0.3z + 2 x =   y = 0.1x + 0y + 0.25z + 2  z = −0.1x + 0.3y + 0z + 0.8 . (II). ðặt.   0 0.2 0.3   0 0 0   0 0.2 0.3         0 0.25  =  0.1 0 0  +  0 0 0.25  α =  0.1   −0.1 0.3 0   −0.1 0.3 0   0 0 0   = L + U   2  x     β =  2  , x =  y       z   0.8  Ta có.  α ∞ = 0.5 < 1, L ∞ = 0.4  0 0  1 0 0  1 0 0  0        0  I 3 − L =  0 1 0  −  0.1 0 0  =  −0.1 1  0 0 1   −0.1 0.3 0   0.1 −0.3 1         . ðặt x. (0).  2  = β =  2  . Ta tính nghiệm xấp xỉ x( n +1) , n = 0,1, 2,... theo công thức  0.8    27.

(28) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ.   2.64   (1)  −1  (0)    x = ( I − L) Ux +β  =  2.464  1.2752     ;    L∞ 0.4 ∆1 ≤ x(1) − x(0) = × 0.64 = 0.512 > 10−2 ∞ 1− α ∞ 1 − 0.5 .   2.87536   (2)  −1  (1)    x = ( I − L) Ux +β  =  2.606336   1.294365        L∞ 0.4 ∆ 2 ≤ x (2) − x(1) ≈ × 0.24 ≈ 0.192 > 10−2 ∞ 1 − 0.5 1− α ∞ . ;.   2.909577    (3) −1  (2)    x = ( I − L) Ux +β  =  2.614549   1.293407        L∞ 0.4 ∆3 ≤ x(3) − x(2) ≈ × 0.034 ≈ 2.72 × 10−2 > 10−2 ∞ α 1 1 0.5 − −  ∞. ;.   2.910932   (4)  −1  (3)    x = ( I − L) Ux +β  =  2.614445   1.29324     ;    L∞ 0.4 ∆ 4 ≤ x(4) − x (3) ≈ × 1.4 × 10−3 ≈ 1.12 × 10−3 < 10−2 ∞ 1 − 0.5 1− α ∞ . V ậy x. (4).  x = 2.910932  hay  y = 2.614445 là xấp xỉ nghiệm thỏa yêu cầu sai số.  z = 1.29324 . Chú ý : Chúng ta cũng có thể tính các xấp xỉ nghiệm theo công thức (*) 28.

(29) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Ví dụ2. Cho hệ phương trình. 10x −2x    x  − x. − 2y − + 20y − − +. 3y y. 3z 5z. − +. t 3t. = 20 = 40. + 10z − t = 8 − 2z + 10t = 10. (I). Bằng phương pháp Seidel (dùng công thức *) tìm nghiệm gần ñúng của hệ sau 3 bước lặp. đánh giá sai số khi nhận xấp xỉ nghiệm này.. ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 29.

(30) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 30.

(31) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. BÀI TẬP Bài1. Giải các hệ phương trình sau bằng phương pháp lặp ñơn , Seidel với sai số 10-3 .  1, 02x − 0,05y − 0,1z = 0,795  1)  −0,11x + 1, 03y − 0, 05z = 0,849  −0,11x − 0,12y + 1,04z = 1,398 . 6,1x + 2, 2y + 1, 2z = 16,55  2) 2, 2x + 5,5y − 1,5z = 10,55 1, 2x − 1,5y + 7, 2z = 16,80 .  1, 02x − 0, 25y − 0,30z = 0,515  3)  −0, 41x + 1,13y − 0,15z = 1,555  −0, 25x − 0,14y + 1, 21z = 2, 780 . 4x − y + z = 8  4) 2x + 5y + 2z = 3  x + 2y + 4z = 11 . 4x + y + 2z = 9  5) 2x + 4y − z = −5  x + y − 3z = −9 . 3x − y + z = 1  6) 3x + 6y + 2z = 0 3x + 3y + 7z = 4 . =9 10x − y  7)  − x + 10y − 2z = 7  − 2y + 10z = 6 . = 24 4x + 3y  8) 3x + 4y − z = 30  − y + 4z = −24 .  0, 42x − 5, 05y − 0,11z = 0, 215  9)  12,5x + 1,02y − 0,05z = 0,743  −0,11x − 0,12y + 2,09z = 1,395 . 2,1x + 2, 2y + 7,5z = 14, 65  10) 5, 2x + 0,5y − 1,5z = 20,15 1, 6x − 4,5y + 1, 2z = −6,18 . Bài2. Giải các hệ phương trình sau bằng phương pháp Seidel sai số 10-5 .  1, 42x1 − 0,5x 2 −  0,5x + 5, 02x −  1 2 1)   0,17x1 + 2,12x 2 +  −0,18x1 − 0,12x 2 +. 0,1x 3 + 0, 2x 4 = 2,525.  0, 42x1 − 5, 05x 2 −  12,5x + 1, 02x − 1,15x 3 − 0,3x 4 = 0, 741  1 2 2)  13,5x 3 + 0, 4x 4 = 5,190  −0,11x1 − 0,12x 2 +  −0,11x1 − 0,12x 2 + 1, 05x 3 − 20,7x 4 = 1,824.  8x1 − x 2 + 2x 3 + x 4 + 2x 5 = 24  2x + 12x − x − 2x + x = 72 2 3 4 5  1 3)  − x1 + 5x 2 + 23x 3 + x 4 + 3x 5 = 46  3x + 2x − 5x − 35x − x = 70 2 3 4 5  1  4x1 − x 2 + x 3 + 2x 4 − 72x 5 = 144. 0,11x 3 + 0,1x 4 = 0, 215 0,05x 3 − 0,5x 4 = 0, 743 2,09x 3 + 0, 4x 4 = 1, 395 1,05x 3 − 5, 2x 4 = 2,092. 25x1 − x 2 + 3x 3 + 2x 4 + x 5 = 75  x + 17x − x − 3x − 4x = 170 1 2 3 4 5  4)  −3x1 + 2x 2 + 35x 3 + x 4 − 5x 5 = 105  4x − 5x − x + 55x + 7x = 330 2 3 4 5  1  − x1 − x 2 + x 3 + 2x 4 + 29x 5 = 580. 31.

(32) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. CHƯƠNG 4 ðA THỨC NỘI SUY VÀ PHƯƠNG PHÁP BÌNH PHƯƠNG BÉ NHẤT MỤC ðÍCH, YÊU CẦU Sau khi học xong chương 4, yêu cầu sinh viên:. 1. Hiểu ñược thế nào là bài toán nội suy và hồi quy. 2. Thực hiện ñược các phương pháp nội suy ña thức. 3. Tìm ñược các hàm xấp xỉ theo phương pháp bình phương bé nhất 4.1. GIỚI THIỆU CHUNG. Khi nghiên cứu các vấn ñề kỹ thuật, kinh tế, xã hội chúng ta thường gặp phải nhu cầu từ các số liệu rời rạc ñã có của các ñại lượng ñang xét, suy ra mối quan hệ toán học giữa chúng, sau ñó sử dụng công cụ toán học nghiên cứu các vấn ñề mà ta quan tâm trên các ñại lượng ñang xét. Ví dụ. Quan sát hai ñại lượng X , Y ta có bảng số liệu:. x. 32. 32.9. 34. 34.5. 35. 36.6. y. 32.4. 33. 33.1. 34.7. 35.2. 33.6. Có rất nhiều câu hỏi liên quan ñến mối quan hệ giữa X,Y mà nếu không sử dụng công cụ toán học thì chúng ta không trả lời ñược ví dụ như: -. Khi X tăng thì Y có tăng hay không ?. -. Khi nào thì Y ñạt cực ñại?. -. Khi X= 36 thì Y là bao nhiêu ?. … Vấn ñề xây dựng mối quan hệ toán học giữa các ñại lượng có thể phát biểu bằng bài toán tổng quát sau Bài toán. Quan sát hai ñại lượng x, y ta ñược bảng số liệu x. x0. x1. …. xn. y. y0. y1. …. yn. Tìm mối liên hệ giữa x,y dưới dạng y = f(x) ?.. 32.

(33) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Khi giải bài toán trên ñiều ñầu tiên chúng ta quan tâm là nên chọn dạng hàm f(x) như thế nào. Các ñịnh lý về xấp xỉ sau ñây của Weierstrass sẽ cho chúng ta gợi ý về dạng hàm f(x). ðịnh lý Weierstrass 1. Cho f (x) là một hàm thực liên tục xác ñịnh trên ñoạn [a,b]. Khi ñó với mọi ε>0 tồn tại một ña thức p m (x) = a 0 + a1x + a 2 x 2 + ... + a m x m với các hệ số thực sao cho với mọi giá trị x thuộc [a,b] ta có |f(x) – pm(x)|<ε. ðịnh lý Weierstrass 2. Cho f (x) là một hàm thực liên tục xác ñịnh trên ñoạn [-π,π] và f(-π) = f(π). Khi ñó với mọi ε>0 tồn tại một ña thức lượng giác a0 m p m (x) = + ∑ a j cos( jx) + b j sin( jx)  2 j=1  với các hệ số thực sao cho với mọi giá trị x thuộc [-π,π] ta có |f(x) – pm(x)|<ε. Như vậy việc chọn ña thức là thích hợp cho dạng hàm f(x). Tiếp theo chúng ta sẽ ñi xác ñịnh các hệ số ai, bj trong ña thức pm(x). Việc xác ñịnh các hệ số thường dựa vào một trong hai dạng yêu cầu: Dạng 1: ða thức pm(x) phải ñi qua các ñiểm (xi ,yi). Tức là pm(xi ) = yi với i=0,1,...,n. Dạng 2: n. ða thức pm(x) ñi gần các ñiểm (xi ,yi) theo nghĩa. ∑ [ pm (xi ) − yi ]. 2. bé nhất.. i =1. Người ta gọi ña thức pm(x) xây dựng theo dạng 1 là ña thức nội suy và ñược dùng khi biết yi = f(xi). ða thức tìm theo dạng 2 gọi là tìm theo phương pháp bình phương bé nhất (hay còn gọi là bài toán hồi quy hoặc hàm hồi quy) nó ñược dùng khi yi ≈ f (x i ) . Chú ý: Khi xây dựng quan hệ giữa y và x theo phương pháp bình phương bé nhất có thể. không phải dạng ña thức. 4.2. ðA THỨC NỘI SUY 4.2.1 ða thức nội suy Lagrange (ñọc giáo trình). 33.

(34) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. 4.2.2. ða thức nội suy Newton Bài toán. Quan sát hai ñại lượng x, y ta ñược bảng số liệu x. x0. x1. …. xn. y. y0. y1. …. yn. Trong ñó : x0 < x1 <…< xn . Tìm mối liên hệ giữa x,y dưới dạng y = p n (x) = a 0 + a1x + a 2 x 2 + ... + a n x n . Thỏa ñiều kiện yi = p n (x i ) . Tìm giá trị y khi x = x* với x* ∉ {x 0 , x1,..., x n } . Giải. a) Các giá trị xi cách ñều : h = xi+1 - xi Bước 1. Tính các hiệu hữu hạn tiến ∆ik.  ∆i = yi +1 − yi  ; k=0,1,...n-1.   k +1 k k ∆i = ∆i +1 − ∆i. Bước 2. Lập ña thức nội suy ða thức nội suy Newton tiến. p Tn (x) = y0 +. ∆0 ∆2 ∆n (x − x 0 ) + 02 (x − x 0 )(x − x1 ) + ... + 0 n (x − x 0 )⋯ (x − x n −1 ) h 2!h n!h. ða thức nội suy Newton lùi. p nL (x) = y n +. ∆n ∆ n −1 ∆2 (x − x n ) + n −22 (x − x n )(x − x n −1 ) + ... + 0n (x − x n )⋯ (x − x1 ) h 2!h n!h. Bước 3 . Tính y khi x = x*. Nếu x* gần x0 thì y ≈ pTn (x* ) . Nếu x* gần xn thì y ≈ p Ln (x* ) .. 34.

(35) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. b) Các giá trị xi không cách ñều Bước 1. Tính các tỉ hiệu hữu hạn tiến f[x i ,...x i + n ]. yi +1 − yi   f[x i , x i +1 ] = x − x i +1 i  ; i = 0,1...; k = 1,2,...   f[x i +1,..., x i + k ] − f[x i ,..., x i + k −1 ] f[x i ,..., x i + k ] = xi + k − xi  Bước 2. Lập ña thức nội suy ða thức nội suy Newton tiến. p Tn (x) = y 0 + f[x 0 , x1 ](x − x 0 ) + ... + f[x 0 ,..., x n ](x − x 0 )⋯ (x − x n −1 ) ða thức nội suy Newton lùi. p Ln (x) = y n + f[x n −1, x n ](x − x n ) + ... + f[x 0 ,..., x n ](x − x n )⋯ (x − x1 ) Bước 3 . Tính y khi x = x*. Nếu x* gần x0 thì y ≈ pTn (x* ) . Nếu x* gần xn thì y ≈ p Ln (x* ) .. Ví dụ1. Cho bảng số liệu. x. 1. 3. 5. 7. 9. 11. y. 3.2. 3.3. 1.7. 2.5. 5.1. 4.3. Xây dựng ña thức nội suy Newton tiến, lùi. Tìm y khi x=1.2345 hoặc x = 9.5437 Giải Bước 1. Tính các hiệu hữu hạn tiến ∆ik. ∆i. ∆i2. ∆ 3i. ∆i4. n. x. y. 0. 1. 3.2. 1. 3. 3.3. 0.1. 2. 5. 1.7. -1.6. -1.7. 3. 7. 2.5. 0.8. 2.4. 4.1. 4. 9. 5.1. 2.6. 1.8. -0.6. -4.7. 5. 11. 4.3. -0.8. -3.4. -5.2. -4.6. ∆ 5i. 0.1 35.

(36) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Bước 2. Lập ña thức nội suy ða thức nội suy Newton tiến. p5T (x). ∆0 ∆ 02 ∆30 = y0 + (x − x 0 ) + (x − x 0 )(x − x1 ) + (x − x 0 )(x − x1 )(x − x 2 ) h 2!h 2 3!h 3 +. ∆ 04. ∆50. 4!h. 5!h 5. (x − x 0 )(x − x1 )(x − x 2 )(x − x 3 ) + 4. (x − x 0 )(x − x1 )(x − x 2 )(x − x 3 )(x − x 4 ). 4.1 0.1 1.7 (x − 1)(x − 3)(x − 5) (x − 1) − (x − 1)(x − 3) + 2 3!23 2!22 4.7 0.1 − (x − 1)(x − 3)(x − 5)(x − 7) + (x − 1)(x − 3)(x − 5)(x − 7)(x − 9) 4 4!2 5!25. p5T (x) = 3.2 +. ða thức nội suy Newton lùi. p5L (x) = y5 + +. ∆14 4!h 4. ∆2 ∆4 ∆3 (x − x 5 ) + 32 (x − x 5 )(x − x 4 ) + 23 (x − x 5 )(x − x 4 )(x − x 3 ) h 2!h 3!h (x − x 5 )(x − x 4 )(x − x 3 )(x − x 2 ) +. ∆50 5!h 5. (x − x 5 )(x − x 4 )(x − x 3 )(x − x 2 )(x − x1 ). 0.8 3.4 5.2 (x − 11) − (x − 11)(x − 9) − (x − 11)(x − 9)(x − 7) 2 2 2!2 3!23 4.6 0.1 − (x − 11)(x − 9)(x − 7)(x − 5) + (x − 11)(x − 9)(x − 7)(x − 5)(x − 3) 4!24 5!25 Bước 3 . Tính y p5L (x)= 4.3 −. Khi x = 1.2345. Vì x gần x0 nên ta có y ≈ p5T (1.2345) = Khi x = 9.5437. Vì x gần x5 nên ta có y ≈ p5L (9.5437) =. Ví dụ2. Cho bảng số liệu. x. 1. 2. 5. 7. 8. 10. y. 3. 2. 1. 4. 5. 1. Xây dựng ña thức nội suy Newton tiến, lùi. Tìm y khi x=2.6375 hoặc x = 8.5722 Giải. 36.

(37) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Bước 1. Tính các tỉ hiệu hữu hạn tiến. n. x. y. TH1. TH2. TH3. TH4. 0. 1. 3. 1. 2. 2. -1. 2. 5. 1. -1/3. 1/6. 3. 7. 4. 3/2. 11/30. 1/30. 4. 8. 5. 1. -1/6. -4/45. -11/630. 5. 10. 1. -2. -1. -1/6. -7/720. TH5. 13/15120. Bước 2. Lập ña thức nội suy ða thức nội suy Newton tiến. 1 1 (x − 1)(x − 2) + (x − 1)(x − 2)(x − 5) 6 30 11 13 − (x − 1)(x − 2)(x − 5)(x − 7) + (x − 1)(x − 2)(x − 5)(x − 7)(x − 8) 630 15120. p5T (x) = 3 − (x − 1) +. ða thức nội suy Newton lùi. 1 p5L (x)= 1 − 2(x − 10) − (x − 10)(x − 8) − (x − 10)(x − 8)(x − 7) 6 7 13 − (x − 10)(x − 8)(x − 7)(x − 5) + (x − 10)(x − 8)(x − 7)(x − 5)(x − 2) 720 15120 Bước 3. Tính y. Khi x = 2.6375. Vì x gần x0 nên ta có y ≈ p5T (2.6375) Khi x = 8.5722. Vì x gần x5 nên ta có y ≈ p5L (8.5722) 4.3. PHƯƠNG PHÁP BÌNH PHƯƠNG BÉ NHẤT Bài toán . Quan sát hai ñại lượng x, y ta ñược bảng số liệu. x. x1. x2. …. xn. y. y1. y2. …. yn. Trong ñó : x1 < x2 <…< xn . Tìm mối liên hệ giữa x, y dưới dạng a) y ≈ p1 (x) = a + bx n. Thỏa ñiều kiện. b) y ≈ p 2 (x) = a + bx + cx 2. ∑ [ p(xi ) − yi ]. 2. c) y ≈ q(x) = a + b cos x + csin x. bé nhất .. i =1. 37.

(38) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Giải n. n. a) ðặt: F(a, b) = ∑ [ p(x i ) − yi ] = ∑ ( a + bx i − yi ) . 2. i =1. 2. i =1. Yêu cầu bài toán tương ñương: Tìm a,b sao cho F(a, b) ñạt cực tiểu. Theo lý thuyết hàm 2 biến ñể F(a, b) ñạt cực tiểu a, b phải thỏa hệ: n  ' Fa (a, b) = 2∑ ( a + bx i − yi ) = 0  i =1  n  ' F (a, b) = 2 ∑ xi ( a + bx i − yi ) = 0  b  i =1 n. n. n n  an + b ∑ x i = ∑ yi  i =1 i =1 ⇔ n n n  2 a x + b x = ∑ i ∑ x i yi  ∑ i  i =1 i =1 i =1. b) ðặt: F(a, b,c) = ∑ [ p(x i ) − yi ] = ∑ ( a + bx i + cx i − yi ) . 2. i =1. 2. i =1. Yêu cầu bài toán tương ñương: Tìm a, b, c sao cho F(a, b,c) ñạt cực tiểu. Theo lý thuyết hàm 3 biến ñể F(a, b,c) ñạt cực tiểu a, b, c phải thỏa hệ n n n n  '  2 2  Fa (a, b,c) = 2∑ a + bx i + cx i − yi = 0  an + b ∑ x i + c∑ x i = ∑ yi   i =1 i =1 i =1 i =1   n n n n  '  n 2 2 3 Fb (a, b,c) = 2∑ x i a + bx i + cx i − yi = 0 ⇔ a ∑ x i + b ∑ x i + c∑ x i = ∑ x i yi   i =1 i =1 i =1 i =1 i =1   n n n n n F' (a, b,c) = 2∑ x 2 a + bx + cx 2 − y = 0 a ∑ x 2 + b ∑ x 3 + c ∑ x 4 = ∑ x 2 y i i i i i i i i  c  i =1 i i =1 i =1 i =1 i =1. (. ). (. ). (. ). n. n. c) ðặt: F(a, b,c) = ∑ [ p(x i ) − yi ] = ∑ ( a + b cos x i + csin x i − yi ) i =1. 2. 2. i =1. Yêu cầu bài toán tương ñương: Tìm a, b, c sao cho F(a, b,c) ñạt cực tiểu. Theo lý thuyết hàm 3 biến ñể F(a, b,c) ñạt cực tiểu a, b, c phải thỏa hệ n n n   an + b ∑ cos x i + c∑ sin x i = ∑ yi  i =1 i =1 i =1  n n n n  2 a cos x + b cos x + c cos x sin x =  ∑ ∑ ∑ i i ∑ yi cos xi i i  i =1 i =1 i =1 i =1  n n n n a ∑ sin x + b ∑ cos x sin x + c∑ sin 2 x = ∑ y sin x i i i i i i  i =1 i =1 i =1 i =1. 38.

(39) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Chú ý. ðể xác ñịnh dạng hàm người thường biểu diễn các cặp ñiểm (xi,yi) lên mặt phẳng. tọa ñộ Oxy, sau ñó nối các ñiểm xem ñường nối có hình dạng nào và chọn dạng hàm theo gợi ý sau: -. Nếu là ñường thẳng (hoặc gần thẳng) chọn dạng a) .. -. Nếu là ñường cong Parabol chọn dạng b).. -. Nếu là ñường tuần hoàn chọn dạng c). Ví dụ1. Quan sát hai ñại lượng x, y ta ñược bảng số liệu. x. 1. 2.3. 3.1. 4. 4.2. 5.5. 6. y. 2. 3.2. 3.1. 3.8. 4. 5. 4.6. Tìm mối liên hệ giữa x, y dưới dạng y = a + bx . Giải.. y = 1.603588 + 0.554593x. Ví dụ2. Quan sát hai ñại lượng x, y ta ñược bảng số liệu. x. 1. 2.3. 3.1. 3.4. 4.1. 5. 5.6. y. 2. 3.2. 3.1. 6.8. 4. 3.5. 1.6. Tìm mối liên hệ giữa x, y dưới dạng y = a + bx + cx 2 . Giải.. y = −1.385703 + 3.55484x − 0.528606x 2. Ví dụ3. Quan sát hai ñại lượng x, y ta ñược bảng số liệu. x. 0. 1.05. 1.57. 3.14. 6.28. y. 2. 3.2. 3. 0. 2. Tìm mối liên hệ giữa x, y dưới dạng y = a + b cos x + csin x . ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 39.

(40) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 40.

(41) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Chú ý.. Khi yêu cầu xây dựng các hàm không phải dạng a), b), c) chúng ta có thể dùng một trong hai cách sau :. Cách 1. Dùng phép ñổi biến ñể ñưa về các dạng a), b), c). n. ∑ [ p(xi ) − yi ]. Cách 2. Dùng ñiều kiện. 2. bé nhất ñể xác ñịnh các hệ số.. i =1. Ví dụ4. Quan sát hai ñại lượng x, y ta ñược bảng số liệu. x. 1. 2.3. 3.1. 4. 4.2. 5.5. 6. y. 2. 3.2. 3.1. 3.8. 4. 5. 4.6. Tìm mối liên hệ giữa x, y dưới dạng. a) y = ax b , a > 0 .. d) y =. 3. e 2x ae2x + be x + c. b) y = aebx , a > 0 .. .. 1. e) y = ax. b + c ln x. c) y =. 1 abln x. − 1 ; a, b > 0 . 3. −1 ; a > 0. b   f) y =  a + + c tan 2 x  cos x  . ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 41.

(42) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 42.

(43) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Ví dụ5. Quan sát hai ñại lượng x, y ta ñược bảng số liệu. x. 1. 2. 3. 4. y. 2. 3.2. 3.1. 3.8. Tìm mối liên hệ giữa x, y dưới dạng a) y = a + x + be x .. b) y = a + b ln x +. c . x. .. ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 43.

(44) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ......................................................................................... ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ......................................................................................... ..…........................................................................................................................................................... ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ............................................................................................... ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................................................................................. ................................................................................................………………………………………… ……………………………………………… 44.

(45) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. BÀI TẬP Bài1. Cho các mốc nội suy sau :. x. 0. 3. 4. 5. 7. f(x). -1. 3. 2. 1. 4. 1) Viết ña thức Newton tiến, tính f(2.5). 2) Viết ña thức Newton lùi, tính f(6.82). Bài 2. Cho các mốc nội suy cách ñều :. x. 1. 2. 3. 4. 5. f(x). 2. -1. 1. 0. 3. 1) Viết ña thức Newton tiến và tính f (1,5). 2) Viết ña thức Newton lùi và tính f (4,5). Bài 3. Từ bảng số liệu ñã cho, bằng pp bình phương bé nhất tìm hàm ñã chỉ ra. x y. 1) y = a ln x + bx + c . Biết. 1 5 x y. 2) y = a + bx + c cos x . Biết. 3) y = ax + b( x 2 + sin x) . Biết. 2 7.7 1 4.6. x y. 4) y = a ( x 2 + 2) + b cos x . Biết. 2 3.7. 1 3.8 x y. 5) y = a (e x − 1) + b ln( x + 1) . Biết. x y. 3 10 3 4. 2 8.9. 3 15. 1 5. 2 4.3. 0 0. 1 5.5. 4 22.2 3 7. 4 15. 3 42.3. Bài 4. Cho bảng dữ liệu. x. 1. 1,5. 2. 2,4. 3. y. 6,62. 3,94. 2,17. 1,35. 0,89. Bằng phép ñổi biến số thích hợp rồi dùng phương pháp bình phương nhỏ nhất tìm các hàm. 1 . ax + b 1 d) y = . 1 + a + bx a) y =. b) y =. ax 2 + bx + c x. 2. e)y = 3 ( b + ax ) . 2. .. c) y = ax 2 + b. f) y=. a cos 2 x. + b tan x + c. 45.

(46) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. CHƯƠNG 5 TÍNH GẦN ðÚNG TÍCH PHÂN XÁC ðỊNH MỤC ðÍCH, YÊU CẦU Sau khi học xong chương 5, yêu cầu sinh viên:. 1. Hiểu ñược thế nào là bài toán tính gần ñúng tích phân xác ñịnh. 2. Thực hiện ñược các phương pháp tính gần ñúng tích phân xác ñịnh, qua ñó biết cách tính giá trị gần ñúng tích phấn xác ñịnh của một hàm bất kỳ. 3. Biết cách áp dụng các phương pháp tính gần ñúng trên vào việc giải các bài toán ngoài thực tế. 4. Biết cách ñánh giá sai số của từng phương pháp. 5.1. GIỚI THIỆU CHUNG. Quá trình tính giá trị các thông số kỹ thuật hoặc giá trị các ñại lượng kinh tế ñôi khi phải tính tích phân b. I = ∫ f (x)dx a. ðể tính ñúng tích phân I ta chỉ có công thức Newton-Leibniz nhưng công thức này sẽ gặp. khó khăn khi hàm f(x) có nguyên hàm phức tạp khó tìm hoặc không có nguyên hàm hoặc cho bởi bảng giá trị rời rạc (ñiều này rất dễ gặp trong kỹ thuật, kinh tế…). Trong thực tế người ta giảm bớt khó khăn khi tính I bằng cách sử dụng các phương pháp tính gần ñúng ñể tìm giá trị xấp xỉ của I rồi dùng nó thay cho giá trị ñúng nhưng với ñiều kiện ñánh giá ñược sai số tuyệt ñối. ða số các phương pháp tính gần ñúng ñều theo các bước sau Bước 1: Chia ñoạn [a,b] thành n ñoạn nhỏ bằng nhau [ x i , x i +1 ]; i=0,...,n-1 . Bước 2: Trên từng ñoạn nhỏ [ x i , x i +1 ] xây dựng ña thức nội suy bậc m pim (x) của f (x) . n −1. x i +1. i =0. xi. Bước 3: Tính I = ∑ Ii ; Ii = *. ∫. pim (x)dx . Kết luận I ≈ I* .. Bước 4: đánh giá sai số ∆ = I − I* .. 46.

(47) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. 5.2. MỘT SỐ CÔNG THỨC TÍNH GẦN ðÚNG 5.2.1. Công thức hình thang. - Giả sử f(x) có ñạo hàm cấp 2 liên tục trên [a,b] - Chia ñoạn [a,b] thành n ñoạn con bằng nhau. - ðặt x 0 = a tính:. - Tính. b−a  h = n  x = x + h; i = 0,1,...,n-1. i  i +1  yi = f (x i ); i = 0,1,...,n.  y + yn  I* = h  0 + y1 + ... + y n −1  . Kết luận I ≈ I* .  2 . - đánh giá sai số : ∆ = I − I* ≤. (b − a)h 2 f ''(c) ; c ∈ [a, b]. . 12. Nhận xét : Giá trị h càng nhỏ thì sai số càng bé. Chú ý : Trong thực hành. - Nếu ñạo hàm cấp 2 của f(x) tính ñược thì sai số ñánh giá theo công thức (b − a)h 2 ∆≤ M với M ≥ max f ''(x) 12 x∈[a,b] - Nếu ñạo hàm cấp 2 của f(x) không tính ñược hoặc quá phức tạp thì sai số ñánh giá theo công thức ∆≤. 1 * * I n − I 2n 3. với I*h , I*2h lần lượt là giá trị I* ñược tính khi chia ñoạn [a,b] thành n và 2n ñoạn có chiều dài bằng nhau.. 47.

(48) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011 1. Th.s ðỗ Hoài Vũ. 2. Ví dụ1. Tính gần ñúng I = ∫ e x dx . 0. đánh giá sai số. Biết chia ựạn [0.1] thành 10 ựoạn có chiều dài bằng nhau Giải. Lập bảng số liệu với h= 0.1. i. 0 1 2 3 4 5 6 7 8 9 10. xi. yi = e 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1. Tính giá trị. -. y +y  I* = h  0 10 + y1 + ... + y9  = 1.467175 2  . x i2. -. y0 = 1 y1 = 1.01005 y 2 = 1.040811 y3 = 1.094174 y 4 = 1.173511 y5 = 1.284025 y 6 = 1.433329 y 7 = 1.632316 y8 = 1.896481 y9 = 2.247908 y10 = 2.718282. đánh giá sai số. Ta có 2. 2. 2. y = e x ⇒ y ' = 2xe x ⇒ y '' = 2e x (1 + 2x 2 ) suy ra Vậ y ∆ ≤. max y ''(x) = 6e .. x∈[0,1]. (1 − 0)0.12 6e = 0.014 12. π 2. Ví dụ2. Tính gần ñúng I = ∫ ln(cos x + sin x)dx . π 6. π π đánh giá sai số. Biết chia ựoạn  ,  thành 12 ựoạn có chiều dài bằng nhau 6 2 2. Ví dụ3. Cho tích phân. ∫. I = ln(1 + e x )dx 1. Nếu dùng phương pháp hình thang thì cần chia ñoạn [1,2] thành ít nhất mấy ñoạn có chiều dài bằng nhau ñể sai số khi tính gần ñúng I không quá 10-3. ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 48.

(49) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 49.

(50) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... .................................................................................................................................................... ..................................................................................................................................................... ...................................................................................................................................................... 50.

(51) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. 5.2.2. Công thức Simpson tổng quát (Simpson 1/3). - Giả sử f(x) có ñạo hàm cấp 4 liên tục trên [a,b] - Chia ñoạn [a,b] thành n =2m ñoạn con bằng nhau. - ðặt x 0 = a tính:. - Tính. I* =. b−a  h = n  x = x + h; i = 0,1,...,n-1. i  i +1  yi = f (x i ); i = 0,1,...,n.. h [ y0 + y2m + 4(y1 + ... + y 2m −1 ) + 2(y 2 + ... + y 2m − 2 )] . 3. - Kết luận I ≈ I* . - đánh giá sai số ∆ = I − I* ≤. (b − a)h 4 (4) f (c) ; c ∈ [a,b]. . 180. Nhận xét : Giá trị h càng nhỏ thì sai số càng bé. Chú ý : Trong thực hành. - Nếu ñạo hàm cấp 4 của f(x) tính ñược thì sai số ñánh giá theo công thức ∆≤. (b − a)h 4 M với M ≥ max f (4) (x) 180 x∈[a,b]. - Nếu ñạo hàm cấp 4 của f(x) không tính ñược hoặc quá phức tạp thì sai số ñánh giá theo công thức ∆≤. 1 * I 2m − I*4m 15. với I*2m , I*4m lần lượt là giá trị I* ñược tính khi chia ñoạn [a,b] thành 2m và 4m ñoạn có chiều dài bằng nhau.. 51.

(52) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Ví dụ1. Tính gần ñúng 1. 2. I = ∫ e x dx . 0. đánh giá sai số. Biết chia ựạn [0.1] thành 10 ựoạn có chiều dài bằng nhau Giải. Lập bảng số liệu với h= 0.1 i 0 1 2 3 4 5 6 7 8 9 10. -. Tính I* =. 2. xi. yi = e x i 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1. y0 = 1 /////////////////////// y1 = 1.01005 y 2 = 1.040811 /////////////////////// /////////////////////// y3 = 1.094174 y 4 = 1.173511 /////////////////////// ////////////////////// y5 = 1.284025 y 6 = 1.433329 /////////////////////// /////////////////////// y 7 = 1.632316 /////////////////////// y8 = 1.896481 /////////////////////// y9 = 2.247908 y10 = 2.718282. h [ y0 + y10 + 4(y1 + ... + y9 ) + 2(y2 + ... + y8 )] = 1.462652 3. - đánh giá sai số ∆ = I − I*. Với mọi x thuộc [0,1] ta có : y = e x. V ậy. 2.  y ' = 2xe x 2 = 2xy ≤ 2e   y '' = 2(y + xy ') ≤ 6e ≤e ⇒ (3)  y = 2(2y '+ xy '') ≤ 20e  (4) (3)  y = 2(3y ''+ xy ) ≤ 76e. max y ''(x) ≤ 76e ≡ M . Kết luận : ∆ ≤. x∈[0,1]. .. (1 − 0) × 0.14 × 76e = 1.148 × 10−4 180. 52.

(53) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Ví dụ2. Tính gần ñúng π 2. I = ∫ cos(e x + 1)dx . 0.  π Không ñánh giá sai số. Biết chia ñoạn 0,  thành 10 ñoạn có chiều dài bằng nhau  2 Ví dụ3.. Cho tích phân 1. I = ∫ ln(1 + x + 1)dx 0. Hỏi phải chia ñoạn [0,1] thành mấy ñoạn con bằng nhau ñể khi tính gần ñúng I bằng công -4. thức simpson bảo ñảm ñược sai số tuyệt ñối < 3.10. ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 53.

(54) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 54.

(55) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. BÀI TẬP Bài1. Tính các tích phân sau bằng công thức hình thang với n = 10. π. I=. 1 ∫0 x + cos x dx ;. 1. xdx ∫0 ln ( 2 + x ) ;. J=. 1. M=. sin x 2dx ∫0 ln ( 2 + x ) ;. 1. K=. 1. N=. e x dx ∫0 sin (1 + x ) ;. x 2dx ∫0 sin (1 + x ) ;. 1. L=. 0. x arcsin xdx ∫0 2x + 1 ;. tgxdx e. 1 1+ x. ;. 1. 1. G=. ∫. H=. e 2x ∫0 1 + cos 3x dx .. Bài2.. Giải lại bài 1 bằng cách sử dụng công thức Simpson 1/3 với n = 10. Bài3. 3,1. Khi tính gần ñúng I =. x3 ∫ x − 1 dx bằng công thức simpson 1/3, cần chia ñoạn [2,1;3,1] 2,1. thành bao nhiêu ñoạn bằng nhau ñể ñạt ñược sai số nhỏ hơn 10-4. Bài4. 1. Khi tính gần ñúng I = ∫ 0. x2 + 1 dx bằng công thức simpson 1/3, cần chia ñoạn [0;1] x+2. thành bao nhiêu ñoạn bằng nhau ñể ñạt ñược sai số nhỏ hơn 0,75.10-4.. 55.

(56) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. CHƯƠNG 6 GIẢI GẦN ðÚNG PHƯƠNG TRÌNH VI PHÂN MỤC ðÍCH, YÊU CẦU Sau khi học xong chương 6, yêu cầu sinh viên:. 1. Hiểu ñược vai trò và tầm quan trọng của bài toán giải gần ñúng phương trình vi phân. 2. Thực hiện ñược các phương pháp tìm nghiệm gần ñúng của phương trình vi phân. 3. Biết cách áp dụng các phương pháp trên vào việc giải quyết các bài toán thực tế. 4. đánh giá ựược sai số của từng phương pháp. 6.1. GIỚI THIỆU CHUNG 6.1.1. ðặt vấn ñề. Nhiều bài toán khoa học kỹ thuật dẫn về việc tìm giá trị của hàm y=y(x) khi x=x*. Biết y là nghiệm riêng của phương trình vi phân thường:  y (a ) = y0  y '( a ) = y 1  (n) . F ( x, y, y ',..., y ) = 0; x ∈ [ a, b] thỏa ñiều kiện ñầu  ⋮   y ( n −1) ( a ) = y n −1 . Giải quyết bài toán trên có hai nhóm phương pháp Phương pháp tìm giá trị thông qua nghiệm chính xác:. Bằng cách dựa vào cách tính tích phân trực tiếp, xác ñịnh ñược dạng tổng quát của nghiệm rồi dựa vào ñiều kiện ban ñầu ñể xác ñịnh nghiệm riêng, sau ñó thay giá trị x* vào nghiệm riêng ñể tìm ra giá trị y. Phương pháp tìm giá trị gần ñúng :. Sử dụng xấp xỉ hàm bằng công thức khai triển Taylor sau ñó dùng phương pháp tính gần ñúng tích phân ñể tính gần ñúng giá trị y. Nhận xét. Hướng giải ñúng có ưu ñiểm là tìm ra ñược giá trị ñúng của y. Tuy nhiên lại rất khó thực hiện bởi không có phương pháp tìm nghiệm riêng tổng quát cho mọi bài toán Hướng giải gần ñúng có khuyết ñiểm là chỉ tìm ra ñược giá trị gần ñúng của nghiệm trong trường hợp hệ có nghiệm duy nhất. Tuy nhiên phương pháp này có thể áp dụng cho một lớp phương trình vi phân rộng hơn rất nhiều so với phương pháp trực tiếp, do ñó ñược dùng nhiều trong thực tế.. 56.

(57) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Trong chương này ta nghiên cứu cách tính gần ñúng và ñánh giá sai số giá trị hàm y thỏa bài toán ñơn giản nhất của phương trình vi phân là bài toán Cauchy ñối với phương trình vi phân cấp 1 Bài toán. Cho hàm y=y(x) thỏa  y ' = f ( x, y ) ;   y ( a ) = y0. x ∈ [ a, b].. (6.1).  x = a, xn = b Tìm giá trị gần ñúng của y khi x = xi thỏa  0 ; i = 0,1,.., n − 1 . đánh giá sai số.  xi +1 = xi + h Giải. Khi giải bài toán trên ta cần tuân thủ các bước sau: Bước 1:. Kiểm tra (6.1) có nghiệm ñúng duy nhất trên [a,b]. Bước 2:. Dùng các thuật toán ñể tìm giá trị xấp xỉ yi của y ( xi ) và ñánh giá sai số. ∆i = yi − y ( xi ) 6.1.2. Một số ñịnh lý cần thiết trong việc thực hiện giải gần ñúng ðể thực hiện bước 1, 2 người ta thường dùng các ñịnh lý dưới ñây ðịnh lý1.. Nếu f ( x, y ) và. ∂f ( x, y ) liên tục trên miền D chứa (x0,y0) thì tồn tại duy nhất hàm y=y(x) ∂y. thỏa (6.1). Ghi chú : Ngoài ðịnh lý 1 còn rất nhiều ðịnh lý nữa nói về sự tồn tại duy nhất nghiệm ðịnh lý2.. Cho hàm số f(x) xác ñịnh và có ñạo hàm ñến cấp n+1 tại x0 và lân cận của x0. Giả sử h là một giá trị sao cho x0 + h cũng thuộc lân cận này. Ta có: f ''( x0 ) 2 f ( n) ( x0 ) n f ( n +1) (c) n +1 f ( x0 + h) = f ( x0 ) + hf '( x0 ) + h + ... + h + h . 2! n! ( n + 1)!. Trong ñó c là hằng số thuộc (x0 , x0+h).. 57.

(58) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. 6.2. MỘT SỐ PHƯƠNG PHÁP TÍNH GẦN ðÚNG 6.2.1. Phương pháp Euler cải tiến. - Kiểm tra bài toán có nghiệm duy nhất. - ðặt x 0 = a, y0 = y(x 0 ) .. - Tính yi +1 :.  y (0) = yi + hf (x i , yi )  i +1 h  (1) (0)   yi +1 = yi + 2 f (x i , yi ) + f (x i +1, yi +1 )   ⋮  h (m −1)   yi(m) +1 = yi +  f (x i , yi ) + f (x i +1, yi +1 )   2. ;. i = 0,1,...,n-1.. Nếu yi(+m1) − yi(+m1−1) ≤ ∆ thì kết luận yi +1 = yi(+m1) .. Ví dụ1. Cho hàm y=y(x) thỏa hệ.  y / = x 2 − y ; x ∈ [0,1]   y (0) = 1 Tìm giá trị xấp xỉ của y(0.1), y(0.2) thỏa yêu cầu sai số 10-3.. Giải Bước 0: Tóm tắt bài toán. i. xi. yi. 0. x0=0. y0=1. 1. x1=0.1. y1=?. 2. x2=0.2. y2=?. ;.  f ( x, y ) = x 2 − y  h = x1 − x0 = x2 − x1 = 0.1. Bước 1: Kiểm tra phương trình có nghiệm duy nhất.  f ( x, y ) = x 2 − y  Ta có  ∂f liên tục trên R2. Vậy phương trình có nghiệm riêng duy nhất.  ∂y ( x, y ) = −1  Bước 2: Tính các giá trị y1 = y(0.1) và y2 = y(0.2).. 58.

(59) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Tính y1.  y (0) = y0 + hf ( x0 , y0 ) = y0 + h( x02 − y0 ) = 1 + 0.1(0 − 1) = 0.9  1  y (1) = y + h  f ( x , y ) + f ( x , y (0) )  = y + h  x 2 − y + x 2 − y (0)  = 0.9055 0 0 0 1 1 0 0 0 1 1   1  2 2   ⇒ y1(1) − y1(0) = 5.5 × 10−3 >10-3    y (2) = y0 + h  f ( x0 , y0 ) + f ( x1, y (1) )  = y0 + h  x02 − y0 + x12 − y (1)  = 0.905225 1  1   1 2 2  ⇒ y1(2) − y1(1) = 2.75 × 10−4 < 10−3  Vậy y1 = 0.905225 .. Tính y2.  y (0) = y1 + hf ( x1, y1 ) = y1 + h( x12 − y1 ) = 0.905225 +0.1(0.12 − 0.905225 )=0.815703  2  y (1) = y + h  f ( x , y ) + f ( x , y (0) )  = 0.821679 1 1 1 2 2   2 2   ⇒ y2(1) − y2(0) = 5.98 × 10−3 >10-3    y (2) = y1 + h  f ( x1, y1 ) + f ( x2 , y (1) )  = 0.82138 2   2 2  ⇒ y2(2) − y2(1) = 2.99 × 10−4 < 10−3  Vậy y2 = 0.82138 .. Ví dụ2. Cho hàm y=y(x) thỏa hệ.  y / = ( y + e x ) cos y ; x ∈ [0.5,1]   y (0.5) = 2 Tìm giá trị xấp xỉ của y(0.55), y(0.58) thỏa yêu cầu sai số 10-3.. ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 59.

(60) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 60.

(61) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. 6.2.2. Phương pháp Runge-Kutta. - Kiểm tra bài toán có nghiệm duy nhất. - ðặt x 0 = a, y0 = y(x 0 ) . Tính yi +1 theo công thức. k1(i) + 2k (i) + 2k 3(i) + k (i) 2 4 yi +1 = yi + 6. k (i) = hf (x i , yi )  1  (i)  k1(i)  h k 2 = hf  x i + , yi +    2 2    . Với    k (i) h  (i) 2  k = hf x + , y +   3  i 2 i 2     (i) (i)  k 4 = hf x i + h, yi + k 3. (. ; i = 0, n − 1. ). Ví dụ1. Cho hàm y=y(x) thỏa hệ.  y / = x 2 y ; x ∈ [1,2]   y (1) = 1 Tìm giá trị xấp xỉ của y(1.1), y(1.12) . Giải Bước 0: Tóm tắt bài toán. i. xi. yi. 0. x0=1. y0=1. 1. x1=1.1. y1=?. 2. x2=1.12 y2=?. ;.  f ( x, y ) = x 2 y   x1 − x0 = 0.1 ≠ x2 − x1 = 0.02. Vì khoảng cách giữa các giá trị x không ñều nên ta phải tính hai lần riêng biệt với hai giá trị h là 0.1 và 0.12 Bước 1: Kiểm tra phương trình có nghiệm duy nhất.  f ( x, y ) = x 2 y  Ta có ∂f liên tục trên R2. Vậy phương trình có nghiệm riêng duy nhất. 2  ∂y ( x, y ) = x  Bước 2: Tính các giá trị y1 = y(1.1) và y2 = y(1.12).. 61.

(62) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. Tính y1 ( với h=0.1).  k (0) = hf ( x0 , y0 ) = h( x02 y0 ) = 0.1  1 2   (0) k1(0)  k1(0)  h h   k = hf x + , y + = h x + y +  0   = 0.115763  2  0   0 0    2 2 2 2         2 (0) (0)  k (0) = hf  x + h , y + k2  = h  x + h   y + k2  = 0.116631  0  0  3  0 2 0 2  2   2      (0) (0) (0) 2 k4 = hf x0 + h, y0 + k3 = h( x0 + h) y0 + k3 = 0.135112. (. Vậy y1 = y0 +. ). (. ). k1(0) + 2k2(0) + 2k3(0) + k4(0) = 1.11665 . 6. Tính y1 (với h=0.12).  k (0) = hf ( x0 , y0 ) = h( x02 y0 ) = 0.12  1 2   (0) k1(0)  k1(0)  h h   k = hf x + , y + = h x + y +     = 0.142923  2  0   0  0 2 0   2 2 2         2 ( 0) (0)  k (0) = hf  x + h , y + k2  = h  x + h   y + k2  = 0.144467  0  0  3  0 2 0 2  2   2      (0) (0) (0) 2 k4 = hf x0 + h, y0 + k3 = h( x0 + h) y0 + k3 = 0.172274 k (0) + 2k2(0) + 2k3(0) + k4(0) Vậy y1 = y0 + 1 = 1.144509 6. (. ). (. ). Ví dụ2. Cho hàm y=y(x) thỏa hệ.  y / = ( y + e x )sin y ; x ∈ [0.5,1]   y (0.5) = 2. Tìm giá trị xấp xỉ của y(0.55), y(0.6) . ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 62.

(63) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 63.

(64) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. BÀI TẬP Bài1:. Giải các phương trình vi phân sau bằng phương pháp Euler cải tiến. Cho ε = 10−4 a ) y / = x 2 − y ; y (0) = 1 trên ñoạn [0;0,5], với bước h = 0,25. b) y / = − xy 2 ; y (0) = 2 trên ñoạn. c) y / = x 2 −. [0;1], với bước h = 0,25.. y ; y (0) = 1 trên ñoạn [0;0,5], với bước h = 0,25. 2. d ) y / = 2 − 3 xy 2 ; y (0) = 1.5 trên ñoạn. e) y / =. [0;1], với bước h = 0,125.. x2 − 2 y ; y (0) = 1 trên ñoạn [0;0,5], với bước h = 0,125. xy + 1. f ) y/ =. 1 + x − 3x 2 y 3 ; y (0) = 2 trên ñoạn 2. [0;1], với bước h = 0,125.. Bài2.. Giải các phương trình vi phân sau bằng phương pháp Runge-kutta . a ) y / = x 2 − y ; y (0) = 1 trên ñoạn [0;0,5], với bước h = 0,25. b) y / = − xy 2 ; y (0) = 2 trên ñoạn [0;1], với bước h = 0,25. c) y / = −. y ; y (0) = 1 trên ñoạn [0;0,5], với bước h = 0,25. 1 + x2. d ) y / = 1 + 3xy 2 ; y (0) = 2 trên ñoạn [0;1], với bước h = 0,25. e) y / =. x2 − y ; y (0) = 1 trên ñoạn [0;0,5], với bước h = 0,25. 1 + 2 xy. 64.

(65) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. MỘT SỐ ðỀ ÔN TẬP THI CUỐI KỲ. ðề số 1 Thời gian : 60 Phút (Không dùng tài liệu) Câu 1 Giải hệ phương trình sau bằng phương pháp Seidel qua 2 bước lặp  2x − 2 y + z + 10t = 10 10x − 2 y − z + t = 10   2 x + 20y − 5 z − 2t = 20  x + 3 y + 20z − 5t = 20. đánh giá sai số khi nhận giá trị xấp xỉ nghiệm ở lẩn ở lần lặp thứ hai. Câu 2 Cho bảng số liệu x y. 0 0. 1 5.5. 3 42.3. Từ bảng số liệu trên, bằng phương pháp bình phương bé nhất tìm hàm có dạng: y = a (e x − 1) + b ln( x + 1) 2. ∫. Cho tích phân I = ln(e x + 2)dx. Câu 3. 1. a) Tính gần ñúng I bằng phương pháp Simpson 1/3. Biết chia ñoạn [1, 2] thành 10 ñoạn có chiều dài bằng nhau . Không ñánh giá sai số . b) Nếu dùng phương pháp hình thang thì cần chia ñoạn [1,2] thành ít nhất mấy ñoạn có chiều dài bằng nhau ñể sai số khi tính gần ñúng I không quá 10-3. Câu 4 Cho hàm y=y(x) thỏa mãn hệ.  y ' = x sin( x + 2 y )   y (0) = 1. a) b). ; x ∈ [0,1]. Dùng phương pháp Euler cải tiến tính giá trị y(0.05) thỏa yêu cầu sai số 10-3. Dùng phương pháp Runge-Kutta tính giá trị y(0.08) không ñánh giá sai số .. .……………………………………………………………………………….. 65.

(66) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ðề số 2 Thời gian : 60 Phút (Không dùng tài liệu) Câu 1 Giải hệ phương trình sau bằng phương pháp Seidel qua 2 bước lặp  10x − y − 2z + t = 10 2 x + 20y − 5 z − 4t = 20    x − 2 y + 2z + 10t = 10  x + 3 y + 20z − 3t = 20. đánh giá sai số khi nhận giá trị xấp xỉ nghiệm ở lẩn ở lần lặp thứ hai. Câu 2 Cho bảng số liệu x y. 0 -1. 1 2.1. 3 12.2. Từ bảng số liệu trên, bằng phương pháp bình phương bé nhất tìm hàm có dạng: y = a( x 2 − 1) + b ln( x + 1) 2. Câu 3. ∫. Cho tích phân I = ln(1 + e x )dx 1. a) Tính gần ñúng I bằng phương pháp Simpson 1/3. Biết chia ñoạn [1, 2] thành 10 ñoạn có chiều dài bằng nhau . Không ñánh giá sai số . b) Nếu dùng phương pháp hình thang thì cần chia ñoạn [1,2] thành ít nhất mấy ñoạn có chiều dài bằng nhau ñể sai số khi tính gần ñúng I không quá 10-3.. Câu 4 Cho hàm y=y(x) thỏa mãn hệ.  y ' = x cos( x + 2 y )   y (0) = 1. ; x ∈ [0,1]. a) Dùng phương pháp Euler cải tiến tính giá trị y(0.05) thỏa yêu cầu sai số 10-3. b) Dùng phương pháp Runge-Kutta tính giá trị y(0.08) không ñánh giá sai số .. .……………………………………………………………………………….. 66.

(67) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ...................................................Phần ghi chép của sinh viên.................................................. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 67.

(68) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 68.

(69) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 69.

(70) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 70.

(71) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 71.

(72) Bài giảng: Toán chuyên ñề 2 –HK3 Năm 2010-2011. Th.s ðỗ Hoài Vũ. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ................................................................................................................................................................. 72.

(73)

phuong phap tinh

Tài liệu liên quan

Tài liệu bạn tìm kiếm đã sẵn sàng tải về