›
›

Bài giảng Phương trình vi phân - Nguyễn Thị Phương Lan

pdf 32 trang ngocly 2460
Download
Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Phương trình vi phân - Nguyễn Thị Phương Lan", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pdfbai_giang_phuong_trinh_vi_phan_nguyen_thi_phuong_lan.pdf

Nội dung text: Bài giảng Phương trình vi phân - Nguyễn Thị Phương Lan

  1. Trường ĐHQN Khoa Toán BÀI GIẢNG PHƯƠNG TRÌNH VI PHÂN (Số đvht: 3) Dành cho sinh viên : Khoa Hóa Hệ : Sư phạm Khóa : 32 Năm học : 2011-2012 Giảng viên : Nguyễn Thị Phương Lan - 1 -
  2. Chương I: PHƯƠNG TRÌNH VI PHÂN CẤP MỘT §1 PHƯƠNG TRÌNH VI PHÂN CẤP MỘT Trong nhiều lĩnh vực kỹ thuật, vật lý, khoa học xã hội ta thường gặp các bài toán dẫn đến việc xác định một hàm thỏa mãn phương trình có chứa một hay nhiều đạo hàm của hàm đó. Các phương trình như vậy gọi là phương trình vi phân (PTVP). PTVP là phương trình liên hệ giữa biến độc lập, hàm cần tìm và các đạo hàm của nó. - Nếu hàm cần tìm chỉ phụ thuộc vào một biến độc lập thì ta có PTVP thường. - Nếu hàm cần tìm phụ thuộc vào hai hoặc nhiều biến độc lập thì ta có phương trình đạo hàm riêng. - Cấp của PTVP là cấp cao nhất của đạo hàm có mặt trong phương trình đó. - Nghiệm của PTVP là mọi hàm thỏa mãn phương trình ấy. Trong học phần này ta chỉ xét đến PTVP thường (còn gọi là PTVP). 1 Ví dụ: y' y 0 là PTVP cấp một, y'' cos x là PTVP cấp hai. x u  u x y 0 là phương trình đạo hàm riêng cấp một. x  y 2u  2 u 0 là phương trình đạo hàm riêng cấp hai. x2  y 2 1.1 Định nghĩa: PTVP cấp một có dạng: F x, y , y ' 0 (1) Nếu giải được đối với y' thì PTVP cấp một có dạng dy y', f x y hay f x, y (2) (dạng chuẩn) hoặc dx P x, y dx Q x , y dy 0 (3) ( dạng vi phân) dy Ví dụ: y' 2y , ex cos x , xydx x2 y 2 dy 0 là các PTVP cấp một. x dx 1.2 Nghiệm của PTVP cấp một: là hàm thỏa mãn phương trình ấy. - Nghiệm tổng quát của PTVP cấp một là nghiệm có chứa một hằng số tùy ý. y x,, C C const . Ví dụ hàm y Cx2 , C const là nghiệm tổng quát của PT y' 2 y . x Về mặt hình học nghiệm tổng quát xác định một họ đường (cong) gọi là họ đường tích phân. - Nghiệm riêng của PTVP cấp một là nghiệm nhận được từ nghiệm tổng quát bằng cách chọn hằng số phù hợp. Chú ý: - Đôi khi giải PTVP ta không tìm được nghiệm tổng quát dưới dạng tường minh - 2 -
  3. y x,, C C const mà được một hệ thức dạng  x, y , C 0, C const nó xác định nghiệm tổng quát dưới dạng ẩn. Hệ thức ấy được gọi là tích phân tổng quát. Hệ thức  x, y , C0 0 được gọi là tích phân riêng. - PTVP có thể có một số nghiệm không nằm trong họ nghiệm tổng quát đó là những nghiệm kỳ dị. 1.3 Bài toán Cauchy (bài toán đầu): PTVP dạng y', f x y cùng với điều kiện y x0 y 0 lập nên bài toán Cauchy (bài toán đầu) của PTVP cấp một. Điều kiện y x0 y 0 với x0, y 0 là các hằng số cho trước được gọi là điều kiện đầu. Ví dụ: Tìm nghiệm riêng thỏa mãn điều kiện đầu y 1 2 của phương trình y' 2 y . x 1.4 Sự tồn tại và duy nhất nghiệm của bài toán Cauchy: Xét phương trình y', f x y f Định lý: Nếu các hàm f x, y và liên tục trong hình chữ nhật D có chứa điểm y x0, y 0 thì tồn tại một lân cận của điểm x0 sao cho PTVP y', f x y có một nghiệm duy nhất thỏa mãn điều kiện y x0 y 0 , nghĩa là bài toán Cauchy y x0 y 0 của PTVP y', f x y có một nghiệm duy nhất. §2 CÁCH GIẢI MỘT SỐ PHƯƠNG TRÌNH CẤP MỘT 2.1 Phương trình phân ly biến số (tách biến) 1. Phương trình dạng: A x dx B y dy (1) trong đó A x là hàm số liên tục của biến x, B y là hàm số liên tục của biến y được gọi là phương trình tách biến. Để giải (1) ta chỉ cần tích phân hai vế. dy Ví dụ 1: a) Giải phương trình vi phân: 3x2 y 2 * dx b) Tìm nghiệm bài toán Cauchy y 0 0 của (*) . Ví dụ 2: Thực nghiệm chỉ ra các chất phóng xạ như uranium có tốc độ phóng xạ tỉ lệ với khối lượng M t tại thời điểm đang xét. Ta có thể viết công thức để tính khối lượng tại bất kỳ thời điểm nào bằng cách giải phương trình dM kM . dt 2. Phương trình dạng: y' f ax by c (2) được đưa về (1) bằng cách z ax by c, z z x . - 3 -
  4. 2.2 Phương trình đẳng cấp và gần đẳng cấp: 2.2.1 Phương trình đẳng cấp: y Phương trình dạng: y' f , x 0 (3) x có thể đưa (3) về phương trình tách biến bằng cách đặt y du u , x 0, u u x y ux,' y x u . x dx du Thay vào (3) ta được x f u u (4). dx dx du du - Nếu f u u 0 thì 4 lnx  u ln C , C 0 x f u u f u u x Ce u , C 0. y y  Thay u , x 0 ta được tích phân tổng quát của (1) là x Ce x , C 0 . x - Nếu f u u tại u u0 thì có thể kiểm tra hàm y u0 x cũng là nghiệm của (3). Đó là một nghiệm riêng. dy y - Nếu f u  u thì (3) có dạng là phương trình tách biến. Nghiệm tổng quát của dx x nó là y Cx, C const . Ví dụ: Giải các phương trình: x y y2 2 xy a) y' b) y' . x y x2 2.2.2 Phương trình gần đẳng cấp: ax by c Phương trình dạng: y' f (5) a1 x b 1 y c 1 trong đó a,,,,, b c a1 b 1 c 1 là các hằng số. - Nếu c c1 0 thì (5) là phương trình đẳng cấp. - Nếu ít nhất một trong các hằng số c hoặc c1 khác 0 thì a b a) Nếu ab1 a 1 b 0 thì (5) có thể đưa về phương trình đẳng cấp bằng a1 b 1 cách đặt: x X ax by c 0 , trong đó ,  là nghiệm của hệ . y Y  a1 x b 1 y c 1 0 dY aX bY Khi đó (5) f là phương trình đẳng cấp. dX a1 X bY 1 - 4 -
  5. a b a1 b 1 b) Nếu 0  a1  a , b 1  b thì a1 b 1 a b dy ax by c 5 f . dx ax by  c1 đặt z ax by, z z x thì ta được phương trình tách biến. Ví dụ: Giải phương trình: 2x y 1 dx x 2 dy 0. 2.3 Phương trình vi phân toàn phần, thừa số tích phân: 2.3.1 Phương trình vi phân toàn phần: Phương trình dạng: P x, y dx Q x , y dy 0 (6) trong đó P P x,,, y Q Q x y là các hàm số liên tục và có các đạo hàm riêng liên tục PQ  trong miền mở, đơn liên D  R2 và thỏa mãn điều kiện ,,  x y D thì (6) y  x được gọi là PTVP toàn phần. UU  Khi đó sẽ tồn tại hàm U x, y D sao cho: PQ , và Pdx Qdy dU . x  y và (1), dU 0 x y . Vậy U x,, y C C const là tích phân tổng quát của (6) . Cách giải: Giả sử xy y x x0, y 0 D Uxy ,,,,, Pxydx 0 Qxydy Qxydy 0 Pxydx . x0 y 0 y 0 x 0 Ví dụ: Giải phương trình: 3x2 6 xy 2 dx 6 x 2 y 4 y 3 dy 0. 2.3.3 Thừa số tích phân: PQ  Xét PTVP P x, y dx Q x , y dy 0 (6). Nếu thì (6) không y  x phải là PTVP toàn phần. Tuy nhiên có thể tìm được hàm   x, y 0 sao cho phương trình Pdx  Qdy 0 (7) là PTVP toàn phần. Hàm   x, y được gọi là thừa số tích phân của (6). Khi đó nếu U x,, y C C const là tích phân tổng quát của (7) cũng đồng thời là tích phân tổng quát của (6). - 5 -
  6. Cách tìm thừa số tích phân: Vì (7) là PTVP toàn phần nên:   PQ     PQPQ    y  x  y  y  x  x 1 PQ  a) Nếu   x và F1 x chỉ phụ thuộc vào x thì có thể tìm được Q  y  x F x dx thừa số tích phân  e 1 . 1 QP  b) Nếu   y và F2 y chỉ phụ thuộc vào y thì có thể tìm được P  x  y F y dy thừa số tích phân  e 2 . Ví dụ: Giải các phương trình: 3 2y 2 2 a) 2xy x y dx x y dy 0 b) y 1 xy dx xdy 0 . 3 2.4 Phương trình vi phân tuyến tính cấp một: Định nghĩa: Phương trình dạng: y' p x y q x (8) trong đó p x , q x là các hàm số liên tục. - Nếu q x  0 thì (8) được gọi là PTVP tuyến tính thuần nhất. - Nếu q x 0 thì (8) được gọi là PTVP tuyến tính không thuần nhất. Cách giải 1: Phương pháp biến thiên hằng số 1) Xét phương trình thuần nhất tương ứng: y' p x y 0 (9) - y 0 là nghiệm của (9). dy p x dx - y 0 thì (2) p x dx y Ce , C 0. y Ngoài ra nghiệm y 0cũng được ghép vào nghiệm tổng quát ứng với C 0 . Vậy p x dx nghiệm tổng quát của (9) là y Ce ,  C const (10) 2) Để tìm nghiệm tổng quát của (8) ta dùng phương pháp biến thiên hằng số. Xem C C x ta tìm C x để (10) là nghiệm tổng quát của (8). Ta có pxdx pxdx dC pxdx yCxe'' Cxpxe e pxy dx thay vào (8) p x dx p x dx dC q x e dx C q x e dx K, K const . - 6 -
  7. Vậy nghiệm tổng quát của (8) là pxdx pxdx pxdx pxdx pxdx y qxe dxKe Ke e qxe dxKconst, Chú ý: Công thức nghiệm: Nghiệm tổng quát của phương trình không thuần nhất = nghiệm tổng quát của phương trình thuần nhất tương ứng + Nghiệm riêng của phương trình không thuần nhất. Cách giải 2: Ta tìm nghiệm tổng quát của (8) dưới dạng y u. v trong đó u u x , v v x mà một trong hai hàm đó có thể chọn tùy ý. Thay vào (8) ta được vu''. v p x v u q x (11). Tìm v x từ điều kiện v' p x v 0 . Thay vào (11) có thể tìm u x từ phương trình vu' q x . Vậy có thể tìm được nghiệm tổng quát của (8). y Ví dụ: 1) Giải bài toán Cauchy: y' 2 4 x , y 1 2 . x 2) Giải phương trình: ey dx xe y 1 dy 0 . 2.5 Phương trình Bernoulli: Định nghĩa: Phương trình dạng: y' p x y y q x (12) trong đó R, p x , q x là các hàm số liên tục. Nếu 0 hoặc 1 thì (12) là PTVP tuyến tính cấp một. Nếu 0 và 1thì (12) được gọi là phương trình Bernoulli. Cách giải: - y 0 là nghiệm của (12). - y 0 thì (12) y y' p x y1 q x (13) 1 Đặt zyzzxy 1 , ' yzz ' ' 1 yy '. Thay vào (13) ta được 1 z' 1 p x z 1 q x (14) là PTVP tuyến tính cấp một. Ví dụ: Giải phương trình: xy' 4 y x2 y . 2.6 Phương trình Clairaut: Định nghĩa: Phương trình dạng: y xy'' f y (15) trong đó f là hàm số khả vi. - 7 -
  8. Cách giải: Đặt y' t , ta có y xt f t . Lấy đạo hàm hai vế đối với x, ta được dt dt dt y'' t x f t t hay x f' t 0 . dx dx dx dt - Nếu 0 thì t là hằng số, ta được họ đường thẳng D phụ thuộc tham số t có dx t phương trình y tx f t . - Nếu x f' t thì y tf' t f t , đó là phương trình tham số của đường tích phân kỳ dị E. Dễ thấy đường E tiếp xúc với mọi đường tích phân Dt . 1 Ví dụ: Giải phương trình: y xy'' y 2 . 4 2.7 Phương trình Lagrange: Định nghĩa: Phương trình dạng y x g y'' f y (16) trong đó f và g là các hàm số khả vi. Cách giải: Đặt y' t , ta có y x g t f t . Lấy đạo hàm hai vế đối với x, ta được dt dt dx y''' g t xg t f t t hay g t t g' t x f ' t 0. dx d x dt Đó là phương trình tuyến tính đối với x t . Nếu nghiệm tổng quát của nó là x C t  t , trong đó C là hằng số tùy ý thì y C t  t g t f t . Ta được phương trình tham số của các đường tích phân. Ví dụ: Giải phương trình: y xy''2 y 2 . §3 PHƯƠNG PHÁP XẤP XỈ PICARD (PHƯƠNG PHÁP GẦN ĐÚNG LIÊN TIẾP) Phương pháp cho ta nghiệm gần đúng của bài toán đầu y',; f x y y x0 y 0 khi giả thiết bài toán có nghiệm duy nhất trong khoảng nào đấy có chứa x0 . Sau khi tích phân bài toán trở thành x y x y0 f t, y t dt x0 xác định một dãy các hàm như sau: x x x yx1 y 0 ftydtyx , 0 , 2 y 0 ftydt , 1 , , yxn y 0 fty , n 1 dt . x0 x 0 x 0 - 8 -
  9. Từ đó ta xây dựng được y1 x từ y0 và f x, y ; y2 x từ y1 x và f x, y xác định mỗi hàm từ hàm ngay trước nó và f x, y . Ta đưa ra sơ đồ xấp xỉ Picard. x yn x y0 f t, y n 1 dt . x0 Với các điều kiện đặt lên hàm f x, y mà ta sẽ xét đến trong định lý tồn tại và duy nhất nghiệm. Có thể chứng minh dãy y0, y 1 , , yn , hội tụ về nghiệm thực y x . Do đó sơ đồ Picard là một công cụ lý thuyết hữu hiệu để chứng minh định lý tồn tại và duy nhất nghiệm của PTVP. Ví dụ: Tìm nghiệm gần đúng của bài toán đầu sau: y' 1 y2 , y 0 0 . Giải: Áp dụng sơ đồ Picard ta tính được x x x3 2 y00 ; y 1 0 dt x ; y 2 0 1 t dt x ; 0 0 3 x t3 x 32 x 5 x 7 ; y3 0 1 t dt x 0 3 3 15 63 Để so sánh kết quả ta tìm nghiệm tổng quat của phương trình y' 1 y2 là arctan y x C . Với điều kiện đầu đã cho thì C 0 và y tan x là nghiệm của bài toán đầu đã xét. Khai triển Maclaurin của tan x ở lân cận x 0 có dạng x32 x 5 17 x 7 tanx x 3 15 315 §4 GIỚI THIỆU CÁC PHƯƠNG PHÁP SỐ Phương pháp số để giải bài toán đầu là một cách xác định nghiệm gần đúng tại các điểm riêng biệt nào đấy mà chỉ cần dùng đến các phép toán cộng, trừ, nhân, chia và tính giá trị hàm. Mọi phương pháp số đều dẫn đến tìm nghiệm gần đúng tại x0, x 1 , , trong đó hiệu giữa hai giá trị x bằng hằng số, tức là xn 1 x n h . Ta mô tả ba phương pháp để tìm nghiệm gần đúng của bài toán đầu y',, f x y y x0 y 0 . 3.1 Phương pháp Euler. Giả sử h nhỏ, ta dùng gần đúng yxh yx hyx', yx hfxy . Đặt xi x0 ih và tính y0 y x 0 , y1 yhfxyy 0 0, 0 , 2 yhfxy 1 1 , 1 , , yn 1 y n hfxy n , n . - 9 -
  10. Vậy bước thứ n của phương pháp Euler có dạng yn 1 y n h f x n, y n . Về mặt hình học nghiệm gần đúng nhận được như một đường gấp khúc mà đoạn đầu tiên là tiếp tuyến với đường cong nghiệm tại x0 . Ví dụ: Xác định 5 xấp xỉ của phương pháp Euler giải bài toán đầu sau đây với h 0,2 y' x y ; y 0 0 . Nghiệm gần đúng yn 1 y n 0,2 x n y n . Nghiệm chính xác y ex x 1. n xn yn Giá trị y đúng 0 0,0 0,0 0,0 1 0,2 0,0 0,021 2 0,4 0,04 0,091 3 0,6 0,128 0,222 4 0,8 0,274 0,425 5 1,0 0,489 0,718 3.2 Phương pháp Euler cải tiến. Đây là phương pháp biến thể của phương pháp Euler. Tại mỗi bước tính giá trị phụ * yn 1 y n h f x n, y n rồi tính giá trị mới h * yn 1 y n f x n,, y n f x n 1 y n 1 . 2 Kết hợp hai biểu thức ta viết bước thứ n của phương pháp Euler cải tiến h yn 1 y n fxy n,,, n fxhyhfxy n n n n . 2 h Về mặt hình học, trong khoảng x, x ta gần đúng y theo đường thẳng qua n n 2 xn, y n với hệ số góc f xn, y n rồi tiếp tục dọc theo đường thẳng với hệ số góc * f xn 1, y n 1 . Ví dụ: Xác định 5 xấp xỉ của phương pháp Euler cải tiến của bài toán đầu nêu trên. * Nghiệm gần đúng yn 1 y n 0,2 x n y n . yn 1 y n0,1 x n y n x n 0,2 y n 0,2 x n y n  - 10 -
  11. Do vậy yn 1 y n 0,22 x n y n 0,02 n xn yn Giá trị y đúng 0 0,0 0,0 0,0 1 0,2 0,0200 0,0214 2 0,4 0,0884 0,0918 3 0,6 0,2158 0,2221 4 0,8 0,4153 0,4255 5 1,0 0,7027 0,7183 3.3 Phương pháp Runge-Kutta. Phương pháp được thiết lập bằng cách lấy trung bình có trọng số của f x, y tại các điểm xác định trong khoảng xn, x n 1 . h y y A 2 B 2 C D n 1 n6 n n n n trong đó h h An f x n,,,, y n B n f x n y n A n 2 2 h h Cfxn n ,,, y n B n DfxhyhC n n n n 2 2 Ví dụ: Xác định 5 xấp xỉ của phương pháp Runge-Kutta của bài toán đầu nêu trên. An x n y n, B n 1,1 x n y n 0,1 Cn 1,11 x n y n 0,11 ; D n 1,222 x n y n 0,222 yn 1 y n 0,2214 x n y n 0,0214 n xn yn Giá trị y đúng 0 0,0 0,0 0,0 1 0,2 0,02140 0,02140 2 0,4 0,09181 0,09182 3 0,6 0,22210 0,22211 4 0,8 0,42552 0,42557 5 1,0 0,71825 0,71828 Cấp của phương pháp số. Phương pháp số có cấp n với n nguyên dương, nếu phương pháp chính xác đến đa thức cấp n của h. Phương pháp Euler là cấp một, phương pháp Euler cải tiến là cấp hai, phương pháp Runge-Kutta là cấp bốn. - 11 -
  12. Chương II: PHƯƠNG TRÌNH VI PHÂN CẤP CAO VÀ HỆ PHƯƠNG TRÌNH VI PHÂN §1 PHƯƠNG TRÌNH VI PHÂN CẤP CAO 1.1 Định nghĩa: PTVP cấp n là phương trình có dạng: F x, y , y ', , y n 0 Nếu giải được đối với đạo hàm cấp n thì PTVP cấp n có dạng: y n f x, y , y ', , y n 1 (1) 1.2 Bài toán Cauchy (bài toán đầu): Bài toán Cauchy đối với phương trình (1) là bài toán tìm nghiệm y y x của k k phương trình (1) thỏa mãn các điều kiện đầu y x0 y 0 , k 0, n 1, trong đó ' n 1 x0, y 0 , y 0 , , y 0 là các hằng số cho trước và được gọi là các giá trị đầu. Định lý (Sự tồn tại và duy nhất nghiệm của bài toán Cauchy): Xét phương trình y n f x, y , y ', , y n 1 n 1 f  n 1 Nếu các hàm f x, y , y ', , y và k ,k0, n 1 liên tục trong miền D R có y ' n 1 chứa điểm x0, y 0 , y 0 , , y 0 thì tồn tại duy nhất một nghiệm y y x của bài toán Cauchy của phương trình (1) trong một lân cận nào đó của điểm x0 . 1.3 Nghiệm tổng quát: Định nghĩa: Hàm số y y x, C1 , , Cn , C i const , i 1, n phụ thuộc vào n hằng số tùy ý được gọi là nghiệm tổng quát của phương trình (1) nếu thỏa mãn - Hàm y x, C1 , , Cn thỏa mãn (1) với mọi giá trị CC1, , n . ' n 1 - Với mọi giá trị x0, y 0 , y 0 , , y 0 cho trước, bài toán Cauchy bao giờ cũng giải được. Các khái niệm khác như nghiệm riêng, nghiệm kỳ dị, tích phân tổng quát, tích phân riêng được định nghĩa tương tự như đối với PTVP cấp một. - 12 -
  13. §2 HẠ THẤP CẤP PTVP CẤP CAO 2.1 Phương trình dạng: F x, y n 0. Nếu giải được đối với y n thì ta có phương trình y n f x (1) n 1 Đặt z y, z z x 2 z ' f x v à z x f x dx C1 . phương trình (2) có dạng như phương trình (1) nhưng cấp thấp hơn một đơn vị. Tích phân n lần ta được kết quả. Ví dụ: Giải phương trình: y''' x2 . 2.2 Phương trình dạng: F y n 1 , y n 0 (3) Đặt z y n 1 , z z x ta có 3 F z , z ' 0 là PTVP cấp một. n 1 Giả sử phương trình này có nghiệm z f x,, C1 y f x C 1 là phương trình có dạng (1) nhưng cấp thấp hơn một đơn vị. Ví dụ: Giải phương trình: y''' y '' 1. 2.3 Phương trình dạng: F x, y n 1 , y n 0 (4) Đặt z y n 1 , z z x ta có 4 F x , z , z ' 0 là PTVP cấp một. y'' Ví dụ: Giải phương trình: y''' x . x Đặt z y n 1 , z z x 4 F x , z , z ' 0 là PTVP cấp một. 2.4 Phương trình dạng: F y, y ', y '' 0 hoặc y'' f y , y ' (5) dy dz dz dy dz Đặt z y', z z x , xem z là hàm của y . Ta có y'' . z . dx dx dy dx dy dz Thay vào (5) ta được F y, z , z 0 là PTVP cấp một, trong đó y được xem là biến dy độc lập, z là hàm của y. Ví dụ: Giải phương trình: 2yy '' y '2 0 . - 13 -
  14. §3 PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH CẤP HAI 3.1 Định nghĩa: Phương trình dạng: y'' p x y ' q x y f x 1 trong đó p x , q x v à f x là các hàm số liên tục. - Nếu f x  0 thì (1) được gọi là PTVP tuyến tính cấp hai thuần nhất. - Nếu f x 0 thì (1) được gọi là PTVP tuyến tính cấp hai không thuần nhất. - Nếu p x , q x là các hằng số thì (1) được gọi là PTVP tuyến tính cấp hai với hệ số hằng. Ví dụ: y'' x2 y ' ex y 0 PTVP tuyến tính cấp hai thuần nhất. 1 y'' 2 xy ' y ex sin x PTVP tuyến tính cấp hai không thuần nhất. x y'' 2 y ' y x2 PTVP tuyến tính cấp hai với hệ số hằng. 3.2 PTVP tuyến tính cấp hai thuần nhất: Phương trình dạng: y'' p x y ' q x y 0 2 trong đó p x , q x là các hàm số liên tục. 3.2.1 Định nghĩa: Hai hàm y1 x , y 2 x được gọi là phụ thuộc tuyến tính nếu tồn tại các hằng số CC1, 2 không đồng thời bằng 0 sao cho C1 y 1 C 2 y 2 0. Trường hợp ngược lại thì chúng được gọi là độc lập tuyến tính. y1 x Nhận xét: Hai hàm y1 x , y 2 x là độc lập tuyến tính nếu C const . y2 x 3.2.2 Định lý: Nếu các hàm y1 y 1 x , y 2 y 2 x là hai nghiệm riêng độc lập tuyến tính của (2) thì hàm y C1 y 1 C 2 y 2 , trong đó CC1, 2 là các hằng số tùy ý là nghiệm tổng quát của (2). Chú ý: Nếu các hàm y1 y 1 x , y 2 y 2 x là hai nghiệm riêng phụ thuộc tuyến tính của phương trình (2) thì y1 k y 2, k const y C 1 y 1 C 2 y 2 kC 1 C 2 y 2 thực chất chỉ phụ thuộc vào một hằng số nên y không phải là nghiệm tổng quát của (2). Nhận xét: - Từ định lý muốn tìm nghiệm tổng quát của (2) chỉ cần tìm hai nghiệm riêng độc lập tuyến của nó (các nghiệm đó được gọi là hệ nghiệm cơ bản). 3.2.3 Công thức Liouville: Nếu đã biết một nghiệm riêng y1 y 1 x của (2) thì có thể tìm nghiệm riêng độc lập tuyến tính y2 y 2 x của nó theo công thức Liouville: 1 p x dx y2 y 1. u , trong đó u u x 2 e dx . y1 x - 14 -
  15. Ví dụ: Tìm nghiệm tổng quát của phương trình: 2x 2 y'' y ' y 0 1 x2 1 x 2 biết một nghiệm riêng y1 x x . 3.3 PTVP tuyến tính cấp hai không thuần nhất: Phương trình dạng: y'' p x y ' q x y f x 3 trong đó p x , q x v à f x là các hàm liên tục, f x 0. Phương trình y'' p x y ' q x y 0 2 được gọi là phương trình thuần nhất tương ứng của (3). 3.3.1 Công thức nghiệm: Nếu gọi y là nghiệm tổng quát của phương trình thuần nhất (2), Y là nghiệm riêng của phương trình không thuần nhất (3) thì y y Y là nghiệm tổng quát của (3). 3.3.2 Nguyên lý cộng nghiệm ( chồng chất nghiệm ): Xét phương trình không thuần nhất y'' pxyqxy ' fx1 fx 2 4 Nếu Y1 là nghiệm riêng của phương trình y'' p x y ' q x y f1 x , Y2 là nghiệm riêng của phương trình y'' p x y ' q x y f2 x thì YYY 1 2 là nghiệm riêng của (4). Kết quả này còn được mở rộng đối với vế phải của (4) là tổng của hữu hạn hàm. 3.3.3 Phương pháp biến thiên hằng số: Giả sử y C1 y 1 C 2 y 2 , trong đó CC1, 2 là các hằng số là nghiệm tổng quát của phương trình thuần nhất tương ứng (2). Khi đó nếu C1 C 1 x , C 2 C 2 x là những hàm số thỏa mãn hệ phương trình: '' C1 y 1 C 2 y 2 0 '''' C1 y 1 C 2 y 2 f x thì hàm y C1 x y 1 C 2 x y 2 là nghiệm tổng quát của phương trình không thuần nhất (3). §4 PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH CẤP HAI VỚI HỆ SỐ HẰNG 4.1 Phương trình tuyến tính thuần nhất: Phương trình dạng: y'' py ' qy 0 1 , trong đó p,q là các hằng số. - 15 -
  16. Cách giải: Ta tìm nghiệm riêng của (1) dưới dạng y ekx , trong đó k là hằng số. Thay y, y ', y '' vào (1) ta được phương trình đại số bậc hai k2 pk q 0 2 và gọi là phương trình đặc trưng, nó có đúng hai nghiệm trong trường số phức £ . Ta có các khả năng xảy ra như sau: k1 x k 2 x - Nếu phương trình đặc trưng (2) có hai nghiệm thực k1 k 2 thì y1 e, y 2 e là hai nghiệm riêng độc lập tuyến tính của (1). Do đó nghiệm tổng quát của (1) k1 x k 2 x là y C1 e C 2 e , trong đó CC1, 2 là các hằng số. kx - Nếu phương trình đặc trưng (2) có nghiệm thực kép k1 k 2 k thì y1 e là một kx nghiệm riêng của (1). Nghiệm riêng y2 xe tìm được theo công thức Liouville. kx Do đó nghiệm tổng quát của (1) là y C1 xC 2 e , trong đó CC1, 2 là các hằng số. - Nếu phương trình đặc trưng (2) có hai nghiệm phức liên hợp k1 i , k2 i  thì i  x x i  x x i  x x i  x x ye1 eee cos xixye sin  , 2 ee e cos  xix sin  là hai nghiệm riêng phức của (1). Ta có y y y y y 1 2 excos x , y 1 2 e x sin  x 12 2 2i là hai nghiệm riêng thực độc lập tuyến tính của (1). Do đó nghiệm tổng quát của (1) là x y e C1cos x C 2 sin  x , trong đó CC1, 2 là các hằng số. Ví dụ : 1) Giải các phương trình: a) y'' 3 y ' 2 y 0 , b) y'' 4 y ' 4 y 0. 2) Giải bài toán Cauchy: y''2'4 y y 0; y 0 1,'0 y 1. 4.2 Phương trình tuyến tính không thuần nhất: 4.2.1 Phương trình dạng: y'' py ' qy f x 3 trong đó p,q là các hằng số, f x là hàm số liên tục. Có thể tìm nghiệm của (3) bằng phương pháp biến thiên hằng số. 1 Ví dụ : Giải phương trình: y'' y . sin x 4.2.2 Phương pháp hệ số bất định (Phương pháp Lagrange): Nếu vế phải f x của (3) có dạng đặc biệt thì có thể tìm nghiệm riêng của (3) theo phương pháp hệ số bất định. - 16 -
  17.  x a) Trường hợp: Vế phải f x e Pn x , trong đó  R , Pn x là đa thức bậc n. - Nếu  không trùng với nghiệm của phương trình đặc trưng (2) thì nghiệm riêng  x của (3) có dạng Y e Qn x . - Nếu  trùng với nghiệm đơn của phương trình đặc trưng (2) thì nghiệm riêng  x của (3) có dạng Y xe Qn x . - Nếu  trùng với nghiệm kép của phương trình đặc trưng (2) thì nghiệm riêng 2  x của (3) có dạng Y x e Qn x . trong đó Qn x là đa thức bậc n mà hệ số của nó được xác định theo phương pháp hệ số bất định.  x b) Trường hợp: Vế phải f x e Pn x cos x P m x sin  x , trong đó ,  R ; Pn x , P m x là các đa thức bậc n,m ; l max m , n - Nếu  i  không trùng với cặp nghiệm phức của phương trình đặc trưng (2) thì  x nghiệm riêng của (3) có dạng Y e Ql x cos x R l x sin  x . - Nếu  i  trùng với cặp nghiệm phức của phương trình đặc trưng (2) thì  x nghiệm riêng của (3) có dạng Y xe Ql x cos x R l x sin  x . trong đó Ql x , R l x là các đa thức bậc l mà hệ số của chúng được xác định theo phương pháp hệ số bất định. Đặc biệt nếu f x  Acos x B sin  x, trong đó A,B là các hằng số. - Nếu i không trùng với cặp nghiệm phức của phương trình đặc trưng (2) thì nghiệm riêng của (3) có dạng: Y  Mcos x N sin  x. - Nếu i trùng với cặp nghiệm phức của phương trình đặc trưng (2) thì nghiệm riêng của (3) có dạng: Y x Mcos x N sin  x. trong đó M,N là các hằng số được xác định theo phương pháp hệ số bất định. Ví dụ : Giải các phương trình: a) y'' 3 y ' 2 y ex 3 4 x b) y'' y 4 x sin 2 x c) y'' y ' 5 ex cos x . §5 PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH CẤP CAO 5.1 Định nghĩa: PTVP tuyến tính cấp n là phương trình có dạng: n n 1 y axyn 1 axyaxyfx 1 ' 0 (1) trong đó f x , a0 x , a 1 x , , an 1 x là các hàm số liên tục . - 17 -
  18. - Nếu f x  0 thì (1) được gọi là PTVP tuyến tính cấp n thuần nhất. - Nếu f x 0 thì (1) được gọi là PTVP tuyến tính cấp n không thuần nhất. - Nếu a0 x , a 1 x , , an 1 x là các hằng số thì (1) được gọi là PTVP tuyến tính cấp n với hệ số hằng. 5.2 Biểu diễn dưới dạng toán tử: Nếu ký hiệu vế trái của (1) là L y và gọi là toán tử vi phân tuyến tính cấp n thì (1) có dạng L y f x . Phương trình thuần nhất tương ứng là L y 0 (2). Toán tử L y có các tính chất: - L Cy CL y , C const . - L y1 y 2 L y 1 L y 2 . m m - L  Ck y k  C k L y k , C k const . k 1 k 1 Từ các tính chất của toán tử L y ta thấy nếu các hàm y1, y 2 , , yn là các nghiệm của phương trình thuần nhất (2) thì tổ hợp tuyến tính y C1 y 1 C 2 y 2 Cn y n , trong đó CCC1, 2 , , n là các hằng số cũng là nghiệm của (2). 5.3 Sự độc lập tuyến tính, phụ thuộc tuyến tính - Định thức Wronski - Hệ nghiệm cơ bản, nghiệm tổng quát của phương trình thuần nhất: 5.3.1 Sự độc lập tuyến tính, phụ thuộc tuyến tính: Định nghĩa: Hệ n hàm yyxy1 1 , 2 yx 2 , , yn yxxab n , , được gọi là phụ thuộc tuyến tính nếu tồn tại các hằng số CCC1, 2 , , n không đồng thời bằng 0 sao cho C1 y 1 C 2 y 2 Cn y n 0. - Các hàm trên được gọi là độc lập tuyến tính nếu chúng không phụ thuộc tuyến tính. 5.3.2 Định thức Wronski: Giả sử các hàm yyxy1 1 , 2 yx 2 , , yn yxxab n , , là các nghiệm của phương trình thuần nhất (2). Định thức Wronski của các nghiệm này được xác định bởi: y1 y 2 yn ''' y1 y 2 yn W = n 1 n 1 n 1 y1 y 2 yn - 18 -
  19. Định lý: Tập hợp n nghiệm yyxy1 1 , 2 yx 2 , , yn yxxab n , , của phương trình thuần nhất (2) là độc lập tuyến tính trong khoảng a, b x0 a, b :W 0. 5.3.3 Hệ nghiệm cơ bản, nghiệm tổng quát của phương trình thuần nhất: Định nghĩa: Hệ gồm n nghiệm độc lập tuyến tính của phương trình thuần nhất (2) được gọi là hệ nghiệm cơ bản của phương trình ấy. Định nghĩa: Nếu hệ n hàm y1 y 1 x , y 2 y 2 x , , yn y n x là hệ nghiệm cơ bản của phương trình (2) thì hàm y C1 y 1 C 2 y 2 Cn y n , trong đó CCC1, 2 , , n là các hằng số là nghiệm tổng quát của (2). Công thức nghiệm: Nghiệm tổng quát của phương trình không thuần nhất (1) là y y Y , trong đó y là nghiệm tổng quát của phương trình thuần nhất tương ứng, Y là nghiệm riêng của (1). 5.4 Phương pháp biến thiên hằng số: Giả sử y C1 y 1 C 2 y 2 Cn y n , trong đó CCC1, 2 , , n là các hằng số là nghiệm tổng quát của phương trình thuần nhất. Nếu C1 C 1 x , C 2 C 2 x , , Cn C n x là những hàm số thỏa mãn hệ phương trình: ''' C1 y 1 C 2 y 2 Cn y n 0 '''''' C1 y 1 C 2 y 2 Cn y n 0 ''' n 1 n 1 n 1 C1 y 1 C 2 y 2 Cn y n f x thì hàm y C1 x y 1 C 2 x y 2 Cn x y n là nghiệm tổng quát của phương trình không thuần nhất. 5.5 PTVP tuyến tính cấp n với hệ số hằng: PTVP tuyến tính cấp n với hệ số hằng là phương trình dạng: n n 1 y an 1 y a 1 y ' a 0 y f x trong đó a0, a 1 , , an 1 là các hằng số, f x là hàm số liên tục. 5.5.1 Phương trình tuyến tính thuần nhất cấp n: Phương trình có dạng: n n 1 y an 1 y a 1 y ' a 0 y 0 (3) trong đóa0, a 1 , , an 1 là các hằng số. - 19 -
  20. Phương trình đặc trưng của (3) là: n n 1 k an 1 k a 1 k a 0 0 (4) (4) là phương trình đại số bậc n, nó có đúng n nghiệm trong trường số phức £ . Ta có: - Nếu (4) có n nghiệm thực phân biệt thì nghiệm tổng quát của (3) là: k1 x k 2 x kn x y C1 e C 2 e Cn e ; CCC1, 2 , , n là các hằng số - Nghiệm thực k bội m m n của (4) cho nghiệm của (3) là: kx m 1 e C1 xC 2 x Cm ; CCC1, 2 , , m là các hằng số. - Cặp nghiệm phức liên hợp i  của (4) cho nghiệm của (3) là: x e C1cos x C 2 sin  x ; CC1, 2 là các hằng số. - Cặp nghiệm phức liên hợp i  bội m m n của (4) cho nghiệm của (3) là: x m 1 ecos x C1 xC 2 x Cm ; CCC1, 2 , , m là các hằng số. x m 1 esin x D1 xD 2 x Dm ; DDD1, 2 , , m là các hằng số. Chú ý: Nghiệm tổng quát của (3) phải chứa đúng n hằng số. Ví dụ : Giải các phương trình: a) y 4 3 y ''' 3 y '' y ' 0 b) y 4 y 0 . 5.5.2 Phương trình tuyến tính không thuần nhất cấp n: Giải tương tự như phương trình không thuần nhất cấp hai. Ví dụ : Giải phương trình: y 4 2 y '' y cos2 x . $6 PHƯƠNG TRÌNH EULER 6.1 Phương trình Euler cấp hai thuần nhất: Phương trình dạng: x2 y'' Axy ' By 0 1 trong đó A, B là các hằng số. Cách giải: dz 1 Cách 1: Đặt x ez , z z x , x 0 z ln x . Ta có: dx x dy dy dz1 dy y'. . dx dz dx x dz d d 1 dy 1 dy 1 d dy 1 dy 1 d2 y dz 1 dy 1 d 2 y y'' y ' . . dx dx x dz x2 dz x dx dz x 2 dz x dz 2 dx x 2 dz x 2 dz 2 . - 20 -
  21. Thay vào (1) ta được 2 2 2 1 d y dy 1 dy d y dy x 2 2 Ax. By 0 2 A 1 By 0 x dz dz x dz dz dz là PTVP tuyến tính cấp hai với hệ số hằng, trong đó z được xem là biến độc lập, y y z là hàm cần tìm. Cách 2: Ta tìm nghiệm riêng của (1) dưới dạng: y xk, x 0, k const y ' kx k 1 , y '' k k 1 x k 2 . Thay vào (1) ta được: k xkk 1 AkB 0 kk 1 AkB 0 k2 A 1 k B 0 (2). và gọi là đặc trưng của (1). k1 k 2 - Nếu (2) có hai nghiệm thực k1 k 2 thì y1 x, y 2 x là hai nghiệm riêng độc lập k1 k 2 tuyến tính của (1). Do đó nghiệm tổng quát của (1) là y C1 x C 2 x , trong đó CC1, 2 là các hằng số. k - Nếu (2) có nghiệm thực kép k1 k 2 k thì y1 x là một nghiệm riêng của (1). k Nghiệm riêng y2 xln x tìm được theo công thức Liouville. Do đó nghiệm tổng k quát của (1) là y C1 C 2 ln x x , trong đó CC1, 2 là các hằng số. - Nếu (2) có hai nghiệm phức liên hợp k1 i , k 2 i  thì nghiệm tổng quát của (1) là y x C1cos  ln x C 2 sin  ln x , trong đó CC1, 2 là các hằng số. Ví dụ : Giải phương trình: x2 y'' 4 xy ' 6 y 0. 6.2 Phương trình Euler cấp hai không thuần nhất: Phương trình dạng: x2 y'' Ax y ' B y f x trong đó A, B là các hằng số. Để giải phương trình Euler-Cauchy cấp hai không thuần nhất có thể dùng phương pháp biến thiên hằng số. Một số trường hợp đặc biệt có thể dùng phương pháp hệ số bất định. Ví dụ : Giải phương trình: x2 y'' 5 xy ' 12 y x ln x . 6.3 Phương trình Euler cấp cao: Có thể giải tương tự như phương trình Euler-Cauchy cấp hai. Phương trình thuần nhất có dạng: n n n 1 n 1 x y An 1 x y A 1 xy ' A 0 y 0 (3) trong đó AAA0, 1 , , n 1 là các hằng số. - 21 -
  22. Phương trình đặc trưng của (3) là một phương trình đại số bậc n, nó có đúng n nghiệm trong trường số phức £ kk 1 kn 1 Akkn 1 1 kn 2 AkA 1 0 0 Ví dụ : Giải phương trình: x3 y''' 5 x 2 y '' 18 xy ' 26 y 0 . 6.4 Phương trình dạng: n n n 1 n 1 axby An 1 axb y AaxbyAy 1 ' 0 0 (4) trong đó AAA0, 1 , , n 1 là các hằng số. Có thể đưa (4) về (3) bằng cách đặt t ax b hoặc giải (4) bằng cách đặt ax b ez , z z x . §7 HỆ PHƯƠNG TRÌNH VI PHÂN CẤP MỘT 7.1 Hệ chuẩn tắc cấp một là hệ n PTVP cấp một dạng dy 1 f x, y , y , , y dx 1 1 2 n dy 2 f x, y , y , , y dx 2 1 2 n (1) dy n f x, y , y , , y dx n1 2 n n 1 trong đó f1, f 2 , , fn là các hàm số liên tục trong miền mở G R , x là biến độc lập, dy dy dy y, y , , y là các hàm cần tìm, 1, 2 , , n là các đạo hàm cấp một của chúng. 1 2 n dx dx dx dy z cos x dx Ví dụ: Hệ phương trình: là hệ chuẩn tắc cấp một. dz 3y 4 z 4cos x sin x dx 7.1.1 Định nghĩa: Tập hợp n hàm y1 x , y 2 x , , yn x khả vi, liên tục trong khoảng n 1 a, b  R sao cho điểm x, y1 x , y 2 x , , yn x G  R và dy i f x, y x , y x , , y x , i 1, n , x a, b là nghiệm của hệ chuẩn tắc (1). dx i1 2 n - 22 -
  23. 7.1.2 Bài toán Cauchy: là bài toán tìm nghiệm y1 x , y 2 x , , yn x của hệ (1) thỏa 0 0 0 mãn các điều kiện yi x0 y i , i 1, n trong đó x0, y 1 , , yn là những số cho trước. 7.1.3 Định lý: ( về sự tồn tại và duy nhất nghiệm của bài toán Cauchy) Giả sử các hàm số f1, f 2 , , fn ở vế phải của các phương trình của hệ (1) liên tục f  f  f và có các đạo hàm riêng i, i , , i ,i 1, n liên tục trong miền DG . Khi đó tồn y1  y 2  yn tại duy nhất n hàm y1 x , y 2 x , , yn x là nghiệm của (1) trong một lân cận U nào đó 0 0 0 của điểm x0 thỏa mãn các điều kiện yi x0 y i , i 1, n trong đó x0, y 1 , , yn D. 7.1.4 Nghiệm tổng quát: Giả sử DG là miền thỏa mãn các điều kiện của định lý. Tập hợp n hàm yi y i x, C1 , C 2 , , C n , i 1, n (2) phụ thuộc vào n tham số CCC1, 2 , , n được gọi là nghiệm tổng quát của hệ (1) nếu 1) Tập hợp các hàm (2) là nghiệm của hệ (1) với mọi hằng số CCC1, 2 , , n . 0 0 2) Với mọi giá trị x0, y 1 , , yn Dcho trước, bài toán Cauchy bao giờ cũng giải được. Các khái niệm khác như nghiệm riêng, nghiệm kỳ dị, tích phân tổng quát, tích phân riêng được định nghĩa tương tự như đối với PTVP cấp một. Ví dụ: 1) Chứng minh rằng hệ hàm: x 3 x y1 x C 1 e C 2 e x 3 x y2 x C 1 e 3 C 2 e cos x là nghiệm tổng quát của hệ phương trình: ' y1 cos x y 2 ' (*) y2 4cos x sin x 3 y 1 4 y 2 2) Giải bài toán Cauchy đối với hệ (*) với điều kiện đầu y1 0 1, y 2 0 2. 7.2 Cách giải hệ chuẩn tắc cấp một: 7.2.1 Phương pháp đưa về phương trình vi phân cấp cao (phương pháp khử): Là phương pháp đưa về một phương trình vi phân cấp cao đối với một hàm số chưa biết bằng cách khử các hàm số chưa biết còn lại từ những phương trình của hệ. - 23 -
  24. Ví dụ: Giải các hệ phương trình: y2 y' y' z y' y z z a) b) c) z' y x z' y z x 1 z' y 2 7.2.2 Phương pháp tổ hợp tích phân: là phương pháp tổ hợp một số phương trình vi phân của hệ, sau đó qua một số phép biến đổi và lấy tích phân ta được nghiệm của hệ. Ví dụ: a) Giải hệ phương trình: dy y dx2 y 3 z dz z dx2 y 3 z b) Tìm nghiệm của bài toán Cauchy y 0 1, z 0 2. y' y2 yz Ví dụ: Giải hệ phương trình: 2 . z' z yz §8 HỆ PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH CẤP MỘT VỚI HỆ SỐ HẰNG 8.1 Định nghĩa: Hệ PTVP tuyến tính cấp một với hệ số hằng là hệ phương trình dạng: dy 1 a y a y a y f x dx 11 1 12 2 1n n 1 dy 2 a y a y a y f x dx 21 1 22 2 2n n 2 (1) dy n a y a y a y f x dx n1 1 n 2 2 nn n n trong đó aij , i 1, n , j 1, n là các hằng số, fi x , i 1, n là các hàm số liên tục, yi y i x , i 1, n là các hàm cần tìm. Nếu các hàm số fi x 0, i 1, n thì hệ (1) được gọi là hệ thuần nhất. Nếu ít nhất một trong các hàm số fi x 0, i 1, n thì hệ (1) được gọi là hệ không thuần nhất. - 24 -
  25. 8.2 Hệ thuần nhất (Phương pháp Euler): Có thể giải hệ thuần nhất mà không cần đưa về PTVP cấp cao. Để thuận tiện ta xét hệ gồm hai phương trình: dy a y a z dx 11 12 (2) dz a y a z dx 21 22 trong đó aij ; i , j 1,2 là các hằng số, y x , z x là các hàm cần tìm. Cách giải: y ekx  Ta tìm nghiệm riêng của (2) dưới dạng kx , trong đó ,,k R . z  e a k a  0 Thay vào (2) ta được 11 12 (3) a21 a 22 k  0 (3) là hệ phương trình đại số tuyến tính thuần nhất. Nó có nghiệm khác 0 khi a11 k a 12 0 a11 k a 22 k a 12 a 21 0 (4) a21 a 22 k (4) là phương trình đại số bậc hai đối với k, nó có đúng hai nghiệm trong trường số phức £ và được gọi là phương trình đặc trưng của (2). - Nếu (4) có hai nghiệm thực k1 k 2 thì thay k1 vào (3) ta được các nghiệm 1,  1 . k1 x k 1 x Khi đó y1 1 e, z 1  1 e là hai nghiệm riêng độc lập tuyến tính tương ứng với k1 . Tương tự thay k2 vào (3) ta được các nghiệm 2,  2 . Khi đó k2 x k 2 x y2 2 e, z 2  2 e cũng là hai nghiệm riêng độc lập tuyến tính tương ứng với k2 . Các nghiệm riêng độc lập tuyến tính y1,,, y 2 z 1 z 2 được gọi là hệ nghiệm cơ bản. Do đó nghiệm tổng quát của (2) là: y C y C y y C ek1 x C e k 2 x 1 1 2 2 1 1 2 2 , k1 x k 2 x z C1 z 1 C 2 z 2 z 1 C 1 e  2 C 2 e trong đó CC1, 2 là các hằng số. - Nếu (4) có nghiệm thực kép k1 k 2 k . Khi đó ta tìm nghiệm riêng của (2) dưới dạng y Ax B ekx, z Cx D e kx , trong đó A, B, C, D là các hằng số được xác định theo phương pháp hệ số bất định. - 25 -
  26. - Nếu (4) có hai nghiệm phức liên hợp k1,2 p iq thì ứng với k1 p iq thay vào (3) ta được các nghiệm 1,  1 và các nghiệm riêng phức ứng vớik1 là p iq x y e 1 . p iq x z 1 e Dùng công thức Euler tách phần thực, phần ảo ta được hệ nghiệm cơ bản. Ví dụ: Giải các hệ phương trình: dy dy dy y 2 z 2y z 3y z dx dx dx a) b) c) dz dz dz 3y 4 z y 2 z 4y z dx dx dx 8.3 Hệ không thuần nhất (phương pháp biến thiên hằng số): Xét hệ dy a y a z f x dx 11 12 1 (5) dz a y a z f x dx 21 22 2 trong đó f1 x , f 2 x là các hàm số liên tục, được gọi là hệ không thuần nhất. Có thể giải (5) bằng phương pháp biến thiên hằng số. y C y C y Giả sử hệ thuần nhất tương ứng (2) có nghiệm 1 1 2 2 , z C1 z 1 C 2 z 2 trong đó CC1, 2 là các hằng số. y C x y C x y Ta tìm nghiệm tổng quát của (5) dưới dạng 1 1 2 2 , z C1 x z 1 C 2 x z 2 C1'' y 1 C 2 y 2 f 1 x trong đó C1 x , C 2 x là các hàm số thỏa mãn hệ . C1'' z 1 C 2 z 2 f 2 x y1 y 2 Hệ này có nghiệm duy nhất vì 0 (do y1,,, y 2 z 1 z 2 là hệ nghiệm cơ bản). z1 z 2 dy 2y 4 z 1 4 x dx Ví dụ: Giải hệ phương trình sau: . dz 3 y z x2 dx 2 - 26 -
  27. Chương III KHÁI NIỆM PHƯƠNG TRÌNH ĐẠO HÀM RIÊNG - PHÂN LOẠI PHƯƠNG PHÁP TÁCH BIẾN - NGUYÊN LÝ CỘNG NGHIỆM 1. Khái niệm phương trình đạo hàm riêng (PTĐHR) Phương trình đạo hàm riêng (PTĐHR) là phương trình chứa hàm cần tìm của hai hoặc nhiều biến với các đạo hàm riêng theo các biến này. Cấp của PTĐHR là cấp cao nhất của đạo hàm riêng có mặt trong phương trình. Nghiệm của PTĐHR là hàm thỏa mãn phương trình. Nghiệm tổng quát là nghiệm có chứa số hàm bằng số cấp của phương trình. Nghiệm riêng là nghiệm có thể nhận được từ nghiệm tổng quát bằng cách chọn thích hợp hàm tùy ý. Nghiệm đặc biệt là nghiệm không thể nhận được từ nghiệm tổng quát bằng cách chọn thích hợp hàm tùy ý. 2u Ví dụ: 2x y (1) là PTĐHR cấp hai. x  y Bằng cách thế vào phương trình ta thấy 1 uxy , xy2 xy 2 Fx Gy 2 là nghiệm của (1). Nó chứa hai hàm độc lập tùy ý F x và G y . Vậy nó là nghiệm tổng quát. Trường hợp riêng F x sin 2 x , G y 2 y5 3 ta được một nghiệm riêng 1 u x, y x2 y xy 2 sin 2 x 2 x 5 3. 2 Ví dụ: Phương trình truyền sóng. Xét một sợi dây đàn hồi bị kích động tại đầu x 0. Có thể viết phương trình mô tả dao động của sợi dây đàn hồi 2u  2 u a2 t2  x 2 trong đó u x, y là chuyển vị của dây tại vị trí x và thời điểm t, a2 const là hệ số đàn hồi. Ví dụ: Phương trình truyền nhiệt. Phương trình biểu thị sự truyền nhiệt hay tỏa nhiệt trong thanh có dạng TT 2 a2 t  x2 trong đó T x, y là nhiệt độ trí x và thời điểm t, a2 const là hệ số truyền nhiệt. Bài toán biên của PTĐHR là bài toán tìm kiếm các nghiệm của PTĐHR trong miền xác định nào đấy thỏa mãn các điều kiện trên biên của miền gọi là điều kiện biên. - 27 -
  28. Định lý liên quan đến tồn tại và duy nhất của nghiệm gọi là định lý tồn tại và duy nhất. Bài toán Cauchy của PTĐHR. Cho PTĐHR với các đạo hàm riêng của hàm cần tìm theo biến t, x, y, z nu  u  u 2 u n F t,,,,, x y z u , ,2 , . t  t  x  x Với giá trị t t0 cho trước các giá trị của hàm cần tìm u và các đạo hàm theo t đến cấp n 1 (điều kiện đầu). ku k k x, y , z , k 0, n 1 . t t t0 Bài toán Cauchy bao gồm tìm nghiệm của PTĐHR đã cho thỏa mãn các điều kiện ban đầu. Trong chương này ta xét các PTĐHR tuyến tính cấp hai. 2. Phân loại PTĐHR tuyến tính cấp hai của hàm hai biến u x, y có dạng 2u  2 u  2 u  u  u A B C D E Fu G (1) x2  x  y  y 2  x  y trong đó A, B, , G có thể là hàm của x, y nhưng không phụ thuộc u. PTĐHR cấp hai của hàm hai biến không có dạng trên gọi là PTĐHR phi tuyến. Nếu G 0 thì (1) được gọi là phương trình thuần nhất. Nếu G 0 thì (1) được gọi là phương trình không thuần nhất. Tùy thuộc vào dấu của BC2 4A ta phân loại PTĐHR như sau: B2 4 AC 0 - phương trình loại elliptic. B2 4 AC 0 - phương trình loại hypebolic. B2 4 AC 0 - phương trình loại parabolic. 3. Phương pháp tách biến (Phương pháp Fourier). Có nhiều phương pháp giải bài toán biên của PTĐHR tuyến tính. Phương pháp tách biến là một trong những phương pháp quan trọng nhất. Ta tìm nghiệm tổng quát sau đó cho thỏa mãn điều kiện biên. Các định lý sau là cơ sở cho phương pháp. Định lý (nguyên lý cộng nghiệm) Nếu u1, u 2 , , un là các nghiệm của PTĐHR tuyến tính thuần nhất thì C1 u 1 C 2 u 2 Cn u n trong đó CCC1, 2 , , n là các hằng số, cũng là nghiệm của phương trình. Định lý: Nghiệm tổng quát của PTĐHR tuyến tính không thuần nhất nhận được bằng cách cộng nghiệm riêng của phương trình tuyến tính không thuần nhất vào nghiệm tổng quát của phương trình tuyến tính thuần nhất tương ứng. - 28 -
  29. Nếu ABF, , , là những hằng số, thì nghiệm tổng quát của phương trình tuyến tính thuần nhất có thể tìm được bằng cách đặt u eax by trong đó a, b là các hằng số cần xác định. Phương pháp tách biến: Giả sử nghiệm được biểu diễn dưới dạng tích của các hàm chưa biết mà mỗi hàm chỉ phụ thuộc vào một biến độc lập. Kết quả của phương pháp là có thể viết phương trình ở hai vế mà mỗi vế chỉ phụ thuộc vào một biến, vì vậy mỗi vế phải bằng hằng số. Ta lần lượt giải cho từng hàm. Hợp các nghiệm này cho ta nghiệm cần tìm. Ví dụ: Bằng phương pháp tách biến giải bài toán biên u  u 4 ,u 0, y 8 e 3 y . x  y Giải: Đặt u x,. y X x Y y vào phương trình đã cho. Ta có XY'' X' Y 4 XY ' c , 4XY trong đó c là hằng số (vì X X x , Y Y y ). Do đó X Ae4cx,;, Y Be cy A B là các hằng số. Nghiệm của phương trình đã cho uxy ,., XYABe c 4 x y Ke c 4 x y K là hằng số. Từ điều kiện biên u 0, y Kecy 8 e 3 y điều này xảy ra khi K 8 và c 3. Vậy nghiệm cần tìm là u x, y 8 e 3 4 x y 8 e 12y 3 y . - 29 -
  30. BÀI TẬP PHƯƠNG TRÌNH VI PHÂN 1.Giải các phương vi phân tách biến sau: a) xydx x 1 dy 0 b) y'cot g x y 2 c) x 1 y2 y 1 x 2 y ' 0 d) x2 1 y ' 2 xy 2 0 e)2x2 yy ' y 2 2 f) y' y 2 x 1 2. Giải các phương vi phân đẳng cấp sau: y a)2x3 y ' y 2 x 2 y 2 b) xy' y x tg x c) 3x2 y 2 y y 2 x 2 xy ' 0 d) x2 3 y 2 dx 2 xydy 0 3. Giải các phương vi phân sau: 2 y 2 a) y' 2 b) 2x 4 y 6 dx x y 3 dy 0 x y 1 y 2 y 2 x c) y' tg d) 1 x2 y ' y 2 1 0 x 1 x 1 4) Giải các phương trình vi phân toàn phần: a) x 2 x2 y 2 dx y x 2 2 y 2 dy 0 b) x3 3 xy 2 2 dx 3 x 2 y y 3 dy 0 3 y 2 x c) dx yln xdy 0 d) 3x 1 ln y dx 2 y dy 0 x y 5. Tìm các thừa số tích phân và giải các phương trình vi phân sau: a) 1 x2 y dx x 2 y x dy 0 b) ydx 2 x yey dy 0 c) x2 y dx xdy 0 d) 2x2 y 2 y 5 dx 2 x 3 2 x dy 0 6. Giải các phương trình vi phân tuyến tính cấp một: 2y a) y' x2 ex b) xy' x 1 y 3 x2 e x x xy c) y' x arcsin x d) xy' y x2 arctan x 1 x2 dy 1 e) y2 dx 2 xy 3 dy 0 f) dx xcos y sin 2 y - 30 -
  31. 7) Giải các phương trình Bernoulli: a) xy' y y2 ln x b) xy2 y' x 2 y 3 c)3y ' y 1 2 x y4 d) y' y xy 1 2xy y dy 3 e) y' 2 4 arctan x f) xsin y x 2 y 0 1 x 1 x2 dx 8) Dùng xấp xỉ Picard giải bài toán đầu. Tìm nghiệm chính xác rồi so sánh: a) y' xy ; y 0 1 b) y' yex ; y 1 1 c) y' xy 2 x x3 ; y 0 0 9) Sử dụng phương pháp Euler, Euler cải tiến, Runge-Kutta để giải gần đúng bài toán đầu: d a) 0,1  20 ;  0 100 h 1 dt dy b) x2 y; y 0 3 h 0,1 dx 10) Giải các phương trình vi phân cấp cao khuyết: y' a) y'' x x 1 0 b) y'' 2 y ' 1 2 y 0 c) xy'' y ' x2 ln x x 1 d) yy'' y '2 y ' 3 0 e) yy'' y '2 y 3 0 11) Giải các PTVP tuyến tính cấp hai bằng phương pháp biến thiên hằng số: x2 2 x 2 3 1 a) y'' 2 y ' y b) y'' 9 y c) y'' y x3 cos3x x d) y'' 2 y ' y 3 e x x 1 e) y'' y tan2 x 12) Giải các PTVP tuyến tính cấp cao với hệ số hằng: a) y'' y ' 2 y x2 e 2x sin x b) y'' 4 y cos2 x c) y'' 3 y ' 4 y e 4x xe x d) y'' y sin x cos3 x e) y''' 4 y '' 13 y ' x2 cos3 x f) y''' y '' y ' y 4 xex cos x 13) Giải các PTVP tuyến tính cấp hai thuần nhất nếu cho trước một nghiệm riêng: 2 2 1 1 a) x y'' xy ' x y 0, y1 x cos x 4 x 2x b) 2x 1 y '' 4 x 2 y ' 8 y 0, y1 x e 2 c) 2x 1 y '' 4 xy ' 4 y 0 , y1 x x . - 31 -
  32. 14) Giải các phương trình Euler: a) x2 y'' xy ' 4 y x 2 ln x b) x2 y'' 4 xy ' 4 y x 2 x 4 , x 0 1 2 c) x2 y'' 2 xy ' 2 y d) 2x 1 y ''42 x 1'8 y y 8 x 4 x2 15) Giải các hệ phương trình vi phân bằng cách đưa về phương trình vi phân cấp cao: dy dz dy dy 4 3y sin x 3y 4 z 0 y 8 z ex dx dx dx dx a) b) c) dz dz dz z cos x 2y 5 z 0 2y z e 3x dx dx dx 16) Giải các hệ phương trình vi phân bằng phương pháp tổ hợp tích phân: x dy 2 2 2 y' y z z y y' z yz dx a) b) c) 2 x dz z y z' y z' 2yz y2 dx 17) Giải các hệ phương trình vi phân tuyến tính thuần nhất bằng phương pháp Euler: dy dy 2y 3 z y z dx dx a) b) dz dz y y 3 y dx dx dx x y z dy dt 4y 3 z dx dy c) d) x y z dz dt 3y 4 z dx dz x y z dt 18) Giải các hệ PTVP tuyến tính không thuần nhất bằng phương pháp biến thiên hằng số: dy dy dy y z 1 3y 2 z 3 e2x y z cos x dx dx dx a) b) c) dz dz dz 4y z x y 2 z ex 2y z cos x sin x dx dx dx - 32 -