Phương pháp sai phân giải phương trình vi phân tuyến tính

1 LỜI CẢM ƠN Tôi xin chân thành bày tỏ lòng kính trọng và biết ơn sâu sắc tới TS. Nguyễn Văn Hùng, người thầy đã chỉ ra hướng nghiên cứu, chỉ bảo tận tình, chu đáo, động viên và giúp đỡ tôi trong quá trình thực hiện luận văn. Tôi xin gửi lời cảm ơn chân thành tới Ban giám hiệu, các thầy giáo, cô giáo Phòng Sau đại học, Khoa Toán Trường Đại Học Sư Phạm Hà Nội 2, bạn bè và người thân đã tạo điều kiện, động viên, khuyến khích, giúp đỡ tôi hoàn thành luận văn này. Hà Nội, tháng11 năm 2013 Tác giả Bùi Trung Hiếu 2 LỜI CAM ĐOAN Luận văn này là kết quả của quá trình học tập, nghiên cứu của bản thân dưới sự chỉ bảo, dìu dắt của các thầy giáo, cô giáo, đặc biệt là sự hướng dẫn nhiệt tình và chu đáo của TS. Nguyễn Văn Hùng. Trong khi nghiên cứu, tôi đã kế thừa những thành quả nghiên cứu của các nhà khoa học, nhà nghiên cứu với sự trân trọng và biết ơn. Luận văn với đề tài “Phương pháp sai phân giải phương trình vi phân tuyến tính” không có sự trùng lặp. Hà Nội, tháng 11 năm 2013 Tác giả Bùi Trung Hiếu 3 MỤC LỤC Trang Lời cảm ơn 1 Lời cam đoan 2 Mục lục 3 MỞ ĐẦU 5 NỘI DUNG 7 Chương 1: Một số kiến thức chuẩn bị 7 1.1. Sai số 7 1.2. Số gần đúng 11 1.3. Một số tính chất của phương trình vi phân 11 1.4. Sai phân và tính chất 14 Chương 2: Lược đồ sai phân 19 2.1.Bài toán vi phân 19 2.2. Lưới sai phân 19 2.3. Hàm lưới 20 2.4. Đạo hàm lưới 20 2.5. Qui ước viết vô cùng bé 21 2.6. Công thức Taylor 21 2.7. Liên hệ giữa đạo hàm và hàm lưới 22 2.8. Phương pháp sai phân 23 2.9. Giải bài toán sai phân bằng phương pháp truy đuổi 23 2.10. Sự ổn định của bài toán vi phân 27 2.11. Sự xấp xỉ 27 4 2.12. Sự hội tụ 29 2.13. Trường hợp điều kiện biên loại ba 30 Chương 3: Phương pháp sai phân giải phương trình vi phân tuyến tính 3.1. Bài toán vi phân 35 3.2. Đạo hàm lưới 36 3.3. Phương pháp sai phân 37 3.4. Phương pháp sai phân giải phương trình vi phân tuyến tính 46 3.5. Sự xấp xỉ 58 3.6. Sự ổn định của bài toán sai phân 58 3.7. Bài toán sai phân đối với sai số 75 3.8. Sự hội tụ và sai số 76 KẾT LUẬN 86 TÀI LIỆU THAM KHẢO 87 5 MỞ ĐẦU 1. Lý do chọn đề tài Phương pháp sai phân (hay còn gọi là phương pháp lưới) là phương pháp được áp dụng rộng rãi trong nhiều lĩnh vực khoa học, kỹ thuật. Nội dung của nó là dẫn đối tượng cần xét về việc giải phương trình sai phân (tức là hệ thức hoặc các hệ thức liên hệ các giá trị của các hàm số tại các điểm khác nhau). Một trong những ứng dụng quan trọng của phương pháp sai phân là giải các phương trình vi phân thường và phương trình đạo hàm riêng. Ngày nay, cùng với sự phát triển mạnh mẽ của công nghệ thông tin đã có rất nhiều phần mềm và chương trình có thể giúp chúng ta tìm nghiệm của các phương trình vi phân thường với độ chính xác cao trong một thời gian rất ngắn. Một lớp phương trình vi phân thường rất quan trọng đó là phương trình vi phân tuyến tính. Trong một số ít trường hợp của lớp phương trình vi phân tuyến tính ta có thể tìm được nghiệm tường minh. Chẳng hạn, đó là các phương trình vi phân tuyến tính với hệ số hằng và hàm vế phải có một số dạng đặc biệt. Tuy nhiên, trong những trường hợp tổng quát hơn việc tìm ra nghiệm tường minh cho các phương trình vi phân tuyến tính là một vấn đề rất khó. Chính vì lí do đó người ta phải đưa ra các phương pháp để giải số các phương trình vi phân này. Được sự định hướng của TS Nguyễn Văn Hùng em đã chọn đề tài: “Phương pháp sai phân giải phương trình vi phân tuyến tính” để hoàn thành luận văn đào tạo thạc sĩ chuyên ngành toán gải tích. 2. Mục đích nghiên cứu Đề tài này nhằm nghiên cứu, trình bày về phương trình vi phân và phương pháp sai phân ứng dụng vào giải phương trình vi phân tuyến tính. 6 3. Nhiệm vụ nghiên cứu Với mục đích nghiên cứu ở trên, nhiệm vụ nghiên cứu của luận văn là: - Phương trình vi phân. - Phương pháp sai phân - Phương pháp sai phân giải phuơng trình vi phân tuyến tính 4. Đối tượng và phạm vi nghiên cứu - Đối tượng nghiên cứu: Các kiến thức cơ sở cần thiết, các kết quả về phương pháp sai phân giải phương trình vi phân tuyến tính. - Phạm vi nghiên cứu: Các tài liệu, các bài báo trong và ngoài nước liên quan đến phương pháp sai phân giải phương trình vi phân tuyến tính. 5. Phương pháp nghiên cứu - Thu thập tài liệu và các bài báo viết về phương pháp sai phân và các ứng dụng của nó. - Phân tích, tổng hợp kiến thức. - Tham khảo ý kiến của giáo viên hướng dẫn. 6. Những đóng góp của luận văn. Trình bày một cách có hệ thống về phương pháp sai phân giải phương trình vi phân tuyến tính. Hơn nữa, kết quả thu được có thể mở rộng cho một số một số lĩnh vực khác. 7 CHƯƠNG I: MỘT SỐ KIẾN THỨC CHUẨN BỊ 1.1. Sai số 1.1.1. Làm tròn số và sai số của phép làm tròn số Xét một số thập phân dạng tổng quát   .10 .10 .10 p i p s p i p s a            (1.1) Trong đó , , 0, 0 9 j p j j a a        Nếu   0 p s   thì a là số nguyên Nếu   0 p s k k     thì a có phần lẻ gồm k chữ số Nếu s   thì a là số thập phân vô hạn Làm tròn số a là bỏ đi 1 số các chữ số bên phải của số a để được số a gọn hơn và gần đúng với số a . Quy tắc làm tròn : Xét số a ở dạng (1.1) và ta sẽ giữ lại đến bậc thứ i , phần bỏ đi là  thì :   1 1 .10 .10 .10 p i i p i i a            Trong đó : 1 1 1 ; 0 .10 .10 khi 2 ; 2 2 1 1 ; .10 .10 khi 2 1; 2 2 i i i i i i i i i l l l l                                 8 Ta kí hiệu sai số của phép làm tròn là a  , như vậy a a a    , rõ ràng 1 .10 2 i a   . Vì * * a a a a a a a a          , do đó khi làm tròn thì sai số tuyệt đối tăng thêm a  . 1.1.2. Chữ số có nghĩa, chữ số chắc Xét số a ở dạng (1.1) nghĩa là được viết dưới dạng thập phân, khi đó chữ số có nghĩa là mọi chữ số khác 0 và những chữ số 0 bị kẹp giữa hai số khác 0 hoặc nó là những chữ số 0 ở hàng được giữ lại. Xét số a ở dạng (1.1)   .10 .10 .10 p i p s p i p s a            Chữ số j  ở (1.1) của số a là số chắc nếu : .10 ; j a     là tham số cho trước Tham số  sẽ được chọn để sao cho một chữ số vốn là chắc thì sau khi làm tròn vẫn là chắc, rõ ràng i a là chữ số chắc thì 1 i a  cũng là chữ số chắc. 1.1.3. Sai số tính toán Giả sử phải tìm đại lượng y theo công thức   1 2 , , , n y f x x x  Gọi     * * * * * * 1 2 , , , ; n x x x x y f x   là giá trị đúng còn   1 2 , , , ;  n x x x x 9    y f x là giá trị gần đúng * y , * i i i x x x    . Giả sử   1 2 , , , n f x x x là hàm số khả vi liên tục thì :     * * * * 1 2 1 2 , , , , , , n n y y y f x x x f x x x      * 1 n i i i i f x x x      Với i f x  là đạo hàm theo i x tính tại điểm trung gian. Vì f là khả vi liên tục , i x  khá bé nên:   1 2 1 , , , n i n i i y f x x x x x       (1.2) Vậy 1 ln n y i i i y f x y x          . (1.3) a. Sai số của phép toán cộng trừ Nếu 1 n i i y x    thì 1 i x y   , vì vậy ta có : 1 n i i y x      Chú ý rằng : Nếu tổng đại số 1 n i i y x    bé về giá trị tuyệt đối thì y y  lớn, phép tính sẽ kém chính xác. Ta khắc phục bằng cách tránh công thức đưa đến hiệu quả của hai số gần nhau. b. Sai số của phép toán nhân, chia 10 Giả sử 1 1 n i i n p i i x y x       . Áp dụng (1.2) và (1.3) ta có 1 y q x x       và y y y    c. Sai số của phép tính lũy thừa Xét   , 0 y x x       , khi đó y x     Như vậy, Nếu 1   thì độ chính xác là giảm đi Nếu 1   thì độ chính xác tăng lên Nếu 1    ( phép nghịch đảo ) thì độ chính xác là không đổi Nếu * 1 , k k     ( phép khai căn ) thì độ chính xác tăng lên. d. Sai số của phép tính logarit Xét ln y x  , ta có y x    . 1.1.4. Bài toán ngược của sai số Giả sử đại lượng y được tính theo công thức :   1 2 , , , n y f x x x  . Yêu cầu đặt ra là cần tính i x  như thế nào để y    , với  là cho trước. [...]... (1.9) là phương trình vi phân tuyến tính thuần nhất, ngược lại thì gọi là phương trình vi phân tuyến tính không thuần nhất Phương trình thuần nhất có vế trái trùng với vế trái của phương trình không thuần nhất (1.9) gọi là phương trình thuần nhất tương ứng với phương trình không thuần nhất (1.9) 1.3.2.2 Phương trình vi phân tuyến tính thuần nhất hệ số hằng Để giải phương trình vi phân tuyến tính thuần... là sai số tuyệt đối của a a   là sai số tương đối của a a Rõ ràng a,  a càng nhỏ càng tốt 1.3 Một số khái niệm về phương trình vi phân 1.3.1 Phương trình vi phân tuyến tính cấp một Phương trình vi phân tuyến tính cấp một có dạng 12 y  a  x  y  b  x  (1.4) Nếu b  x  là hàm hằng 0, thì gọi (1.4) là phương trình vi phân tuyến tính thuần nhất, ngược lại thì gọi là phương trình vi phân tuyến. .. gọi tắt sai phân hữu hạn cấp k là sai phân cấp k , còn sai phân cấp 1 gọi tắt là sai phân Định nghĩa 2: Ta gọi sai phân cấp 2 của hàm xn là sai phân của sai phân cấp 1 của xn và nói chung, sai phân cấp k của hàm xn là sai phân của sai phân cấp k  1 của hàm số đó Như vậy, sai phân cấp 2 của hàm xn là  2 xn    xn   xn 1  xn  xn 2  xn 1   xn 1  xn   xn 2  2 xn 1  xn ; Sai phân cấp... x  C  x  e  , a x dx thay vào phương trình đã cho được C  x   b  x  e  , tích phân hai vế, cuối cùng nhận được nghiệm tổng quát của phương trình không thuần nhất (1.4) là : a x dx  a x dx y    b  x e  dx  C  e      1.3.2 Phương trình vi phân tuyến tính cấp cao 1.3.2.1 Phương trình vi phân tuyến tính cấp n Phương trình vi phân tuyến tính cấp n có dạng (1.8) 13 n n 1 n... lẽ đó ta nói toán tử sai phân Lh xấp xỉ toán tử vi phân L tới cấp O h 2 Hơn nữa, vì  0  y0      0 và  N  y N      0 nên ta cũng nói: Bài toán sai phân (2.9) – (2.10) xấp xỉ bài toán vi phân (2.1) –(2.2) 29 2.12 Sự hội tụ Định nghĩa: Gọi y  x  là nghiệm của bài toán vi phân (2.1) –(2.2) và  i là nghiệm của bài toán sai phân (2.9) – (2.10) Nói phương pháp sai phân (2.9) – (2.10) hội... tại các nút lưới, ta có bảng kết quả sau: Nghiệm gần i xi đúng Nghiệm đúng y ( xi ) vi 0 -0.5 0 0 1 -0.25 -0.19943 -0.21031 2 0 -0.26393 -0.27415 3 0.25 -0.20295 -0.21031 4 0.5 0 0 Sai số đạt: max  vi  y ( xi )   0.01 i  0,4 35 CHƯƠNG III: PHƯƠNG PHÁP SAI PHÂN GIẢI PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH 3.1 Bài toán vi phân Cho 2 số a và b với a  b Tìm hàm y  y  x  xác định tại a  x  b thỏa mãn:... (2.8) 2.8 Phương pháp sai phân Ta tìm cách tính gần đúng giá trị của nghiệm đúng y  xi  tại các nút xi  h Gọi các giá trị gần đúng đó là  i Muốn có  i ta thay bài toán vi phân (2.1) – (2.2) bởi bài toán sai phân:   xi  qi i  fi Lh   a x 0 ,  N   h  trong đó: ai  p  xi   , qi  q  xi  , f i  f  xi  2  2.9 Giải bài toán sai phân bằng phương pháp truy đuổi Vi t cụ thể... tuyến tính thuần nhất, ngược lại thì gọi là phương trình vi phân tuyến tính không thuần nhất Phương trình vi phân tuyến tính cấp 1 thuần nhất y  a  x  y  0 (1.5) là phương trình phân ly biến số được và có nghiệm tổng quát là y  C.e   a  x dx (1.6) Để tìm nghiệm tổng quát của phương trình không thuần nhất ta dùng phương pháp biến thiên hằng số Lagrange : Thay C bởi hàm số C  x  để y  C... cos x e x sin  x, xe x sin  x, , x k 1e x sin  x đều là nghiệm của (1.10) Điều đó có nghĩa là bằng cách giải phương trình đặc trưng tìm được đủ n nghiệm thì ta có được ngay hệ thống n nghiệm độc lập tuyến tính của phương trình vi phân (1.10) và do đó có nghiệm tổng quát 1.4 Sai phân và tính chất 1.4.1 Các khái niệm cơ bản Xét dãy số  xn  , dạng khai triển của nó là:  x0 , x1 , x2 , , xn , ... vào phương trình ta được   e x  n  a1 n1  a2 n2   an1  an  0 Giải phương trình đặc trưng  n  a1 n1  a2 n2   an1  an  0 (1.11) - Nếu  là nghiệm thực bội k của phương trình (1.11), thì các hàm e x , x.e x , , x k 1e x là k nghiệm riêng độc lập tuyến tính của (1.10) - Nếu     i là nghiệm phức bội k của phương trình (1.11) thì   i cũng là nghiệm bội k của phương . phương pháp sai phân giải phương trình vi phân tuyến tính. - Phạm vi nghiên cứu: Các tài liệu, các bài báo trong và ngoài nước liên quan đến phương pháp sai phân giải phương trình vi phân tuyến. 3.3. Phương pháp sai phân 37 3.4. Phương pháp sai phân giải phương trình vi phân tuyến tính 46 3.5. Sự xấp xỉ 58 3.6. Sự ổn định của bài toán sai phân 58 3.7. Bài toán sai phân đối với sai. nghiên cứu của luận văn là: - Phương trình vi phân. - Phương pháp sai phân - Phương pháp sai phân giải phuơng trình vi phân tuyến tính 4. Đối tượng và phạm vi nghiên cứu - Đối tượng nghiên