Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 29 trang
THÔNG TIN TÀI LIỆU
Cấu trúc
HỆ PHƯƠNG TRÌNH VI PHÂN CẤP 1
ĐỊNH NGHĨA
BÀI TOÁN CAUCHY
PHƯƠNG PHÁP KHỬ
Slide 5
Cách khử cho hệ 2 pt (tuyến tính)
HỆ PTVP TUYẾN TÍNH CẤP 1 HỆ SỐ HẰNG
Ví dụ
Slide 9
PP TRỊ RIÊNG GIẢI HỆ KHÔNG THUẦN NHẤT
Slide 11
Slide 12
Slide 13
Slide 14
PPTRỊ RIÊNG TÌM NGHIỆM HỆ THUẦN NHẤT
Slide 16
Slide 17
Slide 18
Slide 19
Slide 20
Cấu trúc nghiệm hệ tt không thuần nhất
PP biến thiên hằng số tìm Xr
Slide 23
Slide 24
Slide 25
Slide 26
Slide 27
Slide 28
Slide 29
Nội dung
HỆPHƯƠNGTRÌNHVIPHÂNCẤP ĐỊNH NGHĨA F1(t,x1,x2,…, xn, x1’,x2’,…,xn’) = Hệ tổng qt … Fn(t,x1,x2,…, xn, x1’,x2’,…,xn’) = x1’ = f1(t,x1,x2,…, xn) Hệ tắc … xn’ = fn(t,x1,x2,…, xn) t : biến x1, x2 , …, xn : ẩn hàm BÀI TỐN CAUCHY Tìm nghiệm hệ x1’ = f1(t,x1,x2,…, xn) ……………………… xn’ = fn(t,x1,x2,…, xn) x1(t0) = α1 Thỏa điều kiện ………… xn(t0) = αn Hệ n ptvp cấp tương đương ptvp cấp n nên hệ nghiệm có n số tự PHƯƠNG PHÁP KHỬ B1: xây dựng ptvp cấp n theo hàm chọn trước B2: giải ptvp cấp n vừa tìm rút hệ với (n – 1) hàm Vd: x ' = x '(t ) = y + e t t y ' = y '(t ) = − x + 3y − e (1) (2) y′′ = − x '+ 3y '− e t y′′ = −2 y − e t + 3y '− e t (3) ⇒ ⇒ t t x ' = y + e x ' = 2y + e (3) ⇔ y "− 3y '+ y = −2e t t Tt cấphệ số 2t ⇔ y = C1e + C2e + 2te t (2) ⇒ x = − y '+ 3y − e t t = − C1e − 2C2e 2t t t 2t t −2(t + 1)e + 3(C1e + C2e + 2te ) − e t 2t = 2C1e + C2e + (4t − 3)e x = 2C1et + C2e t + (4t − 3)e t t 2t t y = C1e + C2e + 2te t t Cách khử cho hệ pt (tuyến tính) x ′ = a1x + b1y + f1 (t ) y ′ = a2 x + b2 y + f2 (t ) (1) (2) Lấy đạo hàm pt (1) theo t (3) Thay y’ từ pt (2) vào (3) (4) Rút y từ (1) thay vào (4) Pt kết pt cấp theo ẩn hàm x biến t Nếu xuất phát từ pt (2), ta có pt cấp theo y HỆ PTVP TUYẾN TÍNH CẤPHỆ SỐ HẰNG X’(t) = AX(t) + F(t) x1′ (t ) M ÷ ÷ x′ ( t ) ÷ n x1 (t ) M ÷ ÷ x (t) ÷ n f1 (t ) M÷ ÷ f (t) ÷ n (Hệ ẩn hàm ) a11 L A =L L a n1 L a1n L ÷: ma trậ n vuô ng cấ pn ÷ ann ÷ Ví dụ x ' = x '(t ) = y + e t 1/ t y ' = y '( t ) = − x + y − e x(t ) X( t ) = ÷ y ( t ) et F( t ) = −e t ÷ ÷ 2 A = ÷ − x ' = x + y + 2z + t + sin t 2 / y ' = 2x + 4y + z + t t z ' = y − z + e − ln t t + sin t 1 ÷ ÷ ⇔ X ( t ) = X( t ) + t , ÷ ÷ −2 ÷ et − ln t ÷ x(t ) ÷ X( t ) = y ( t ) ÷ z( t ) ÷ PP TRỊ RIÊNG GIẢI HỆ KHƠNG THUẦN NHẤT A chéo hóa X’ = AX + F(t) ⇔ X’ = PDP-1X + F(t) ( ⇔ ∃ P: P-1AP = D (chéo) ) Đặt Y = P-1X: ⇔ P-1X’ = DP-1X + P-1F(t) ⇔ Y’ = DY + G(t) y1′ λ1 K y1 g1 (t ) y′ λ K y g ( t ) 2= 2 + ′ yn 0 K λn yn gn (t ) PPTRỊ RIÊNG TÌM NGHIỆM HỆ THUẦN NHẤT X’(t) = AX(t) ⇔ Y’ = DY y1 ' λ1 K y1 y ' λ K y = 2 y n ' 0 K λn y n y1 '(t ) = λ1y1 (t ) y '(t ) = λ y (t ) 2 ⇔ yn '(t ) = λn yn (t ) y1 '(t ) = λ1y1 (t ) y '(t ) = λ y (t ) 2 ⇔ yn '(t ) = λn yn (t ) n y1 ( t ) = c1eλ1t λ2 t y ( t ) = c 2e ⇔ λn t y n ( t ) = c ne ⇒ X = PY = ∑ c k e Pk { Xk = e k =1 λk t λk t (Pk cột thứ k P) } Pk , k = 1, , n : hệnghiệ mđltt củ a hệthuầ n nhấ t Định lý: Hệ X’ = AX(t), ma trận A có n giá trị riêng thực λ1, λ2 … λn (kể trị riêng bội), n vector riêng P1, P2 , … , Pn độc lập tuyến tính ⇒ Nghiệm tổng qt pt nhất: n X ( t ) = [ x1 ( t ) , x ( t ) , K , xn ( t ) ] = ∑ c k e Pk T k =1 λk t Vd: x1′ = x1 + x + 2x 1 2 x′2 = x1 + x + 2x ⇔ X′ = 1 ÷X ÷ x′ = x + x + x 2 4÷ 3 1− λ A − λI = 1− λ λ1 = ⇔ λ2 = 2 4−λ A = λ (6 − λ ) = 1 p1 ÷ ÷ ⇔ 1 p =0 ( A − λ1I)P = ÷ ÷ 2 ÷ p ÷ 1 2 Chọn vector riêng: P1 = −1÷, P2 = ÷ ÷ ÷ 0÷ −1÷ p1 −5 ÷ ÷ ( A − λ2I)P = ⇔ −5 p2 = ÷ ÷ 2 −2 ÷ p ÷ 1 Chọn VTR: P3 = ÷ ÷ 2÷ λ1t λ1t λ2 t X1 = e P1, X2 = e P2 , X3 = e P3 = e6 tP2 ⇒ X = ∑ Ck Xk k =1 1 2 1 0t 0t 6t ÷ ÷ ÷ = C1e −1 + C2e + C3e ÷ ÷ ÷ 0÷ −1÷ 2÷ 6t x1 C1 + 2C2 + C3e ÷ 6t ÷ ⇔ x = −C1 + C3e ÷ ÷ x ÷ 6t ÷ ÷ −C2 + 2C3e Cấu trúc nghiệm hệ tt khơng X = X0 + Xr X0 : nghiệm tổng qt hệ pt X’(t) = AX(t) (1) Xr : nghiệm riêng hệ pt khơng Cấu trúc nghiệm tổng qt hệ X0 = C1X1 + C2X2 + …+ CnXn { Xk , k = 1, ,n }: hệ nghiệm độc lập tuyến tính (1) PP biến thiên số tìm Xr Xr = C1(t)X1 + …+ Cn(t)Xn Ci tìm từ hệ pt: C’1(t)X1 + …+ C’n(t)Xn = F(t) Ví dụ x1′ = x + et (1) t x′2 = − x1 + 3x − e Hệ nhất: x1′ = x (2) x′2 = − x1 + 3x 2 A = , ÷ −1 et F( t ) = −e t ÷ ÷ Trị riêng VTR A: 2 1 λ1 = 1, P1 = ÷, λ1 = 2, P2 = ÷, 1 1 Các nghiệm đltt hệ t 2 t 1 X1 = e ÷, X = e ÷ 1 1 Nghiệm tổng qt hệ t 2 t 1 X0 = C1X1 + C2 X2 = C1e ÷+ C2e ÷ 1 1 Tìm Xr pp biến thiên số: Trong X0 xem C1 C2 hàm cố theo t Tìm C1 C2 từ hệ: C’1(t)X1 + …+ C’n(t)Xn = F(t) t e t 2 t 1 ⇔ C1′ e ÷+ C′2e ÷ = ÷ 1 1 −et ÷ C1′ 2et + C′2e t = e t ⇔ t 2t t ′ ′ C1e + C2e = −e C1 (t ) = 2t Chọn: −t C2 (t ) = 3e C1′ = ⇔ −t C′2 = −3e t 2 t 1 X0 = C1X1 + C2 X2 = C1e ÷+ C2e ÷ 1 1 C1 (t ) = 2t −t C2 (t ) = 3e ⇒ Xr = C1 (t ) X1 + C2 (t ) X Nghiệm tổng qt: t 2 − t t 1 = 2te ÷+ 3e e ÷ 1 1 4te t + 3et = ÷ 2te t + 3et ÷ X = X0 + Xr Ví dụ x1′ (t ) = 3x1 + x + e t x′2 (t ) = x1 + x + t F = 2t + A = 3t ÷ ÷ 4 Chéo hóa A 1 −1 / −1 / 0 P= P = D= ÷ ÷ ÷ −1 1/ 1/ 5 Đặt : y1 / −1 / x1 Y = P X ⇔ ÷= ÷ ÷ y2 / / x2 −1 et − t t / − / e 3÷ −1 P F= ÷ ÷= ÷ / / t et + t ÷ 3 3 Hệ viết lại theo y1, y2 Y′ = DY + P −1F et − t y1′ 2y1 3÷ ⇔ ÷= + ÷ ÷ y′2 5y et + t ÷ 3 3 y′ = y + e t − t 3 ⇔ t t y′2 = 5y + e + 3 y = − et + t + + C e2t 12 ⇒ y = − e t − t + + C e5 t 12 15 75 ⇒ X = PY x1 1 y1 hay ÷ = ÷ ÷ x −1 y ... 1 1 = 1, P1 = ÷, 1 = 2, P2 = ÷, 1 1 Các nghiệm đltt hệ t 2 t 1 X1 = e ÷, X = e ÷ 1 1 Nghiệm tổng qt hệ t 2 t 1 X0 = C1X1 + C2 X2 = C1e ÷+ C2e ÷ 1 1 ... nghiệm hệ x1’ = f1(t,x1,x2,…, xn) ……………………… xn’ = fn(t,x1,x2,…, xn) x1(t0) = 1 Thỏa điều kiện ………… xn(t0) = αn Hệ n ptvp cấp tương đương ptvp cấp n nên hệ nghiệm có n số tự PHƯƠNG PHÁP KHỬ B1:... 1 1 y1 1 x1 Y = P X ⇔ ÷= P ÷ y2 x2 1 t t 2e 1 e 1 P F( t ) = = t ÷ ÷ ÷ −3et ÷ 1 −e ÷ 1 (1) ⇔ Y′ = DY + P F(t ) t ′ y y e 1