Tài liệu tham khảo giáo trình Phương pháp tính với C++ - Chương 7 Giải phương trình vi phân | CHƯƠNG 7 GIẢI PHƯƠNG TRÌNH VI PHÂN 1. BÀI TOÁN CAUCHY Một phương trình vi phân cấp 1 có thể viết dưới dạng giải được y f x y mà ta có thể tìm được hàm y từ đạo hàm của nó. Tồn tại vô số nghiệm thoả mãn phương trình trên. Mỗi nghiệm phụ thuộc vào một hằng sô tuỳ ý. Khi cho trước giá trị ban đầu của y là yo tại giá trị đầu Xo ta nhận được một nghiệm riêng của phương trình. Bài toán Cauchy hay bài toán có điều kiện đầu tóm lại như sau cho X sao cho b X a tìm y x thoả mãn điều kiện y x f x y y a a 1 Người ta chứng minh rằng bài toán này có một nghiệm duy nhất nếu f thoả mãn điều kiện Lipschitz f xzy1 - f x y2 I lY1 - y2 với L là một hằng sô dương. Người ta cũng chứng minh rằng nếu fy đạo hàm của f theo y là liên tục và bị chặn thì f thoả mãn điều kiện Lipschitz. Một cách tổng quát hơn người ta định nghĩa hệ phương trình bậc 1 y1 f1 x y1 y2 -. yn y2 f2 x y1 y2 . yn yn fn x y1 y2 . yn Ta phải tìm nghiệm y1 y2 . yn sao cho ỊY x f x X Y a a với Y V y1 y2 . F f f1 f2 . Y v Ầ y1 y2 . . yn . l fn . l y Nếu phương trình vi phân có bậc cao hơn n nghiệm sẽ phụ thuộc vào n hằng sô tuỳ ý. Để nhận được một nghiệm riêng ta phải cho n điều kiện 166 đầu. Bài toán sẽ có giá trị đầu nếu với giá trị Xo đã cho ta cho y xo y xo y xo . Một phương trình vi phân bậc n có thể đưa về thành một hệ phương trình vi phân cấp 1. Ví dụ nếu ta có phương trình vi phân cấp 2 y f x y y y a a y a p Khi đặt u y và v y ta nhận được hệ phương trình vi phân cấp 1 U v v g x u v với điều kiện đầu u a a và v a p Các phương pháp giải phương trình vi phân được trình bày trong chương này là các phương pháp rời rạc đoạn a b được chia thành n đoạn nhỏ bằng nhau được gọi là các bước tích phân h b - a n. 2. PHƯƠNG PHÁP EULER VÀ EULER CẢI TIẾN Giả sử ta có phương trình vi phân y x f x y y a a 1 và cần tìm nghiệm của nó. Ta chia đoạn xo x thành n phần bởi các điểm chia xo xi x2 . xn x Theo công thức khai triển Taylor một hàm lân cận xi ta có x. x 2 . x. x 3 . y xi 1 y xẰ xi 1 - xi y x 1 1 0 xj y xi 1 1 xj y x . 2 6 Nếu .
