(Luận văn) chéo hóa phương trình vi phân tuyến tính

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC SƯ PHẠM NGUYỄN BÌNH DƯƠNG lu an CHÉO HĨA PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH n va p ie gh tn to oa nl w d LUẬN VĂN THẠC SĨ TOÁN HỌC nf va an lu z at nh oi lm ul z m co l gm @ an Lu THÁI NGUYÊN - 2020 n va ac th si ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC SƯ PHẠM Nguyễn Bình Dương CHÉO HĨA PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH lu an va Chun ngành: Tốn giải tích n p ie gh tn to Mã số: 8460102 d oa nl w LUẬN VĂN THẠC SĨ TOÁN HỌC nf va an lu Người hướng dẫn khoa học lm ul PGS TSKH Đoàn Thái Sơn z at nh oi z m co l gm @ an Lu Thái Nguyên - 2020 n va ac th si Lời cam đoan Tôi xin cam đoan nội dung trình bày luận văn trung thực không trùng lặp với đề tài khác Nguồn tài liệu sử dụng cho việc hoàn thành luận văn nguồn tài liệu mở Các thông tin, tài liệu luận văn lu an ghi rõ nguồn gốc va n Thái Nguyên, ngày tháng năm 2020 tn to Người viết luận văn p ie gh nl w d oa Nguyễn Bình Dương nf va an lu Xác nhận lm ul khoa chuyên môn Xác nhận người hướng dẫn z at nh oi z PGS TSKH Đoàn Thái Sơn m co l gm @ an Lu n va i ac th si Lời cảm ơn Luận văn hoàn thành khoa Toán, Trường Đại học Sư phạm – Đại học Thái Nguyên hướng dẫn PGS TSKH Đồn Thái Sơn Tơi xin cảm ơn thầy hướng dẫn tận tình hiệu kinh nghiệm q trình nghiên cứu hồn thiện luận văn lu an Tôi xin chân thành cảm ơn Ban chủ nhiệm khoa Tốn, thầy giáo va n tổ môn giảng dạy tạo điều kiện thuận lợi cho gh tn to trình học tập thực luận văn thạc sĩ p ie Cuối cùng, xin gửi lời cảm ơn gia đình, bạn bè quan tâm giúp đỡ, w động viên tơi suốt q trình làm luận văn oa nl Tuy nhiên, luận văn tránh khỏi thiếu sót, tơi mong d nhận góp ý từ thầy giáo để luận văn hoàn thiện lu nf va an Tôi xin chân thành cảm ơn! Thái Nguyên, ngày tháng năm 2020 lm ul Người viết luận văn z at nh oi z m co l gm @ Nguyễn Bình Dương an Lu n va ii ac th si Mục lục lu an Lời cam đoan i Lời cảm ơn ii Mục lục iii Danh mục ký hiệu, chữ viết tắt iv Mở đầu va n Chương Khái niệm phổ nhị phân mũ cho phương trình vi phân gh tn to tuyến tính p ie 1.1 Tính nhị phân mũ nl w 1.2 Phổ nhị phân mũ 15 d oa 1.3 Ví dụ an lu Chương Chéo hóa phương trình vi phân tuyến tính khơng nf va gian hữu hạn chiều 17 lm ul 2.1 Một số kết chuẩn bị tích phân Lebesgue hàm liên tục 17 2.2 Phép biến đổi tương đương 23 2.3 Chéo hóa 25 z at nh oi tuyệt đối z @ l gm Kết luận 33 m co Tài liệu tham khảo 32 an Lu n va iii ac th si Danh mục ký hiệu, chữ viết tắt lu an n va tập số thực Rn×n tập ma trận vng cấp n A : R → Rn×n Hàm ma trận khả tích địa phương ∅ tập rỗng ⊕ tổng trực tiếp ie gh tn to R p GLN (R) Nhân phép chiếu P d an lu oa ker P Ảnh phép chiếu P nl im P w Nhóm ma trận tuyến tính khả nghịch cấp N kết thúc chứng minh nf va z at nh oi lm ul z m co l gm @ an Lu n va iv ac th si Mở đầu Xét hệ phương trình vi phân tuyến tính x˙ = A(t)x (1) với hàm ma trận A : R → RN ×N liên tục bị chặn Nếu (1) phương lu trình vi phân với hệ số không phụ thuộc thời gian, tức an va A ∈ RN ×N n x˙ = Ax, gh tn to việc thay đổi biến x 7→ T −1 x biến phương trình thành phương trình ie có dạng p x˙ = T −1 ATx w oa nl Do đó, cách lựa chọn ma trận khả nghịch T cách phù hợp ta d đơn giản hóa phương trình tuyến tính với hệ số khơng phụ thuộc thời nf va chuẩn Jordan an lu gian Ví dụ ta chọn ma trận khả nghịch T cho T −1 AT dạng lm ul Trong khuân khổ nội dung luận văn này, tìm hiểu kết tương tự cho hệ phương trình tuyến tính phụ thuộc thời gian z at nh oi (1) Ở sử dụng khái niệm khả quy định nghĩa Coppel [3] Cụ thể, hệ (1) gọi khả quy tồn phép đổi biến z khả nghịch phụ thuộc thời gian biến đổi hệ (1) thành hệ đường chéo @ gm khối với số chiều đường chéo khối nhỏ hẳn N co l Để mở rộng kết dạng chuẩn Jordan cho hệ tuyến tính phụ thuộc vào thời gian cần nhắc tới lý thuyết phổ phù hợp hệ (1) m an Lu Khái niệm phổ coi khái niệm khái quát giá trị riêng theo n va ac th si cách thích hợp Nhắc lại rằng, lịch sử có nhiều khái niệm phổ cho phương trình (1), ví dụ khái niệm phổ Lyapunov cho hệ quy, khái niệm phổ Morse cho hệ động lực tích lệch (xem Colonius Kliemann [4]), khái niệm phổ Bohl với mục đích mơ tả tất tốc độ tăng trưởng hệ với thời gian dương (xem Daleckii Krein [5]) Tuy nhiên, khái niệm phổ Sacker-Sell, hay gọi phổ nhị phân mũ, dường khái niệm phổ phù hợp để giải câu hỏi mở rộng Nhằm trình bày cách có hệ thống kết mở rộng định lý dạng chuẩn Jordan cho hệ tuyến tính phụ thuộc thời gian Siegmund [9], luận văn bao gồm nội dung sau: Chương 1: Khái niệm phổ nhị phân mũ cho phương trình vi phân tuyến lu an tính n va Chương 2: Chéo hóa phương trình vi phân tuyến tính khơng gian Cuối phần kết luận tóm tắt kết đạt danh mục gh tn to hữu hạn chiều p ie tài liệu tham khảo d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu n va ac th si Chương Khái niệm phổ nhị phân mũ cho phương trình vi phân tuyến tính lu an Xét hệ phương trình vi phân tuyến tính khơng ơtơnơm có dạng n va (1.1) tn to x˙ = A(t)x ie gh với A : R → RN ×N hàm đo khả tích địa phương, tức với p đoạn [a, b] ⊂ R ta có oa nl w Z b kA(t)k dt < ∞ a d an lu Mục đích chương trình bày khái niệm phổ nhị phân mũ nf va cho phương trình (1.1) định lý phổ nhị phân mũ cho phương trình z at nh oi lm ul Kết chương chứng minh [9] 1.1 Tính nhị phân mũ z Ta kí hiệu Φ : R × R → RN ×N , (t, τ ) 7→ Φ(t, τ ) toán tử tiến hóa gm @ phương trình (1.1), tức Φ(., τ )ξ giải toán giá trị ban đầu (1.1) với l χ(τ ) = ξ Một phép chiếu bất biến (1.1) định nghĩa hàm m co P : R → RN ×N phép chiếu P (t), t ∈ R, cho n va (1.2) an Lu P (t)Φ(t, s) = Φ(t, s)P (s), t, s ∈ R ac th si Lưu ý P liên tục đồng thức P (·) ≡ Φ(., s)P (s)Φ(s, ) Chúng ta nói (1) có tính nhị phân mũ có phép chiếu bất biến P số K ≥ 1, α > cho kΦ(t, s)p(s)k ≤ Ke−α(t−s) với t ≥ s kΦ(t, s)[I − P (s)]k ≤ Keα(t−s) với t ≤ s Tiếp theo nghiên cứu lớp phương trình dịch chuyển sau x˙ = (A(t) − γI)x (1.3) lu γ ∈ R Dễ dàng thấy Φγ (t, s) := e−γ(t−s) Φ(t, s) toán tử an tiến hóa Nếu với giá trị γ mà phương trình dịch chuyển va n (1.3) có tính nhị phân mũ với họ phép chiếu bất biến P P tn to phép chiếu bất biến x˙ = A(t)x, tức (1.2) thỏa mãn Hơn p ie gh nữa, ta có bất đẳng thức sau oa nl w kΦ(t, s)P (s)k ≤ Ke(γ−α)(t−s) với t ≥ s d kΦ(t, s)[I − P (s)]k ≤ Ke(γ+α)(t−s) với t ≤ s an lu nf va Nhận xét 1.1 Nếu x˙ = [A(t) − γI]x có tính nhị phân mũ với phép chiếu lm ul bất biến P ≡ I x˙ = [A(t) − ζI]x có tính nhị phân mũ với phép chiếu bất biến cho ζ > γ Khẳng định với z at nh oi ζ < γ P ≡ Tiếp theo ta trình bày khái niệm mô tả tốc độ tăng trưởng mũ z gm @ hàm số an Lu n va m (b) γ − - bị chặn supt≤0 kg(t)k e−γt < ∞ co (a) γ + - bị chặn supt≥0 kg(t)k e−γt < ∞ l Định nghĩa 1.2 Giả sử γ ∈ R Một hàm liên tục g : R → RN ac th si [a,b] liên tục tuyệt đối Bổ đề sau với hàm liên tục tuyệt đối khả vi hầu hết điểm bổ đề mối liên hệ đạo hàm hàm với phép lấy tích phân Bổ đề 2.2 Với a, b ∈ R, với a ≤ b ta có khẳng định sau: (A) Nếu f : [a, b] → RN khả tích Lebesgue Z t f (σ)dσ, t ∈ [a, b] F (t) := a liên tục tuyệt đối Dt F (t) = f (t) tồn với hầu hết giá trị lu t ∈ [a, b] an n va (B) Nếu F : [a, b] → RN liên tục tuyệt đối Dt F (t) tồn với hầu hết ie gh tn to giá trị t ∈ [a, b] , Dt F có tính chất khả tích Lebesgue [a, b] Z t F(t) − F(a) = Dt F (σ)dσ a p Chứng minh Xem Craven [1] nl w d oa Tích f.g hai hàm liên tục tuyệt đối f, g : I ⊂ R → R liên an lu tục tuyệt đối hàm hợp f ◦ g hàm liên tục tuyệt đối chưa liên tục tuyệt đối Mệnh đề sau giả thiết thêm để nf va đảm bảo cho tính liên tục tuyệt đối phép lấy hợp thành lm ul Bổ đề 2.3 Giả sử gi : R → R, i = 1, , N hàm liên tục tuyệt đối z at nh oi f : RN → R khả vi liên tục Khi hàm hợp f ◦ (g1 , , gN ) : R → R, t 7→ f (g1 (t), , gN (t)) z gm @ liên tục tuyệt đối tục f tồn số C > với m N co l Chứng minh Giả sử [a, b] ⊂ R khoảng compact Khi tính khả vi liên i=1 n va 18 an Lu |f (x1 , , xn ) − f (¯ x1 , , x¯N )| ≤ C max |xi − x¯i | , xi , x¯i ∈ gi ([a, b]) ac th si Bây giờ, cho trước ε > Do tính liên tục tuyệt đối gi tồn δ > cho N max i=1 n X |gi (βk ) − gi (αk )| < ε/C k=1 với a ≤ α1 < β1 ≤ α2 < β2 ≤ · · · ≤ αn < βn ≤ b với Pn k=1 (βk − αk ) < δ Điều suy n X |f (g1 (βk ) − gm (βk )) − f (g1 (αk ) − gm (αk ))| < ε k=1 có điều phải chứng minh Đặt Sym := {A ∈ RN ×N : AT = A} biểu thị tập ma trận đối xứng N × N Đặt L(Sym) tập ánh xạ tuyến tính Sym đặt lu an P d := {A ∈ Sym : A xác định dương} Giờ xét hàm sau va n σ : Sym → Sym, A 7→ A2 gh tn to Dễ dàng nhận thấy σ |P d : P d → P d có tính song ánh Giờ p ie tính đạo hàm w Bổ đề 2.4 Ánh xạ σ : P d → Sym, A 7→ A2 khả vi liên tục với đạo hàm an lu d oa nl Dσ : P d → L(Sym), A 7→ Dσ(A) nf va Dσ(A) : Sym → Sym, X 7→ XA + AX lm ul phép đẳng cấu tuyến tính với A ∈ P d z at nh oi Chứng minh Lấy A ∈ P d X ∈ Sym Khi 2 (A + X) − A − [XA + AX] ≤ kXk z điều σ khả vi P d Để chứng minh tính liên tục @ gm Dσ đặt A, B ∈ P d X ∈ Sym Khi bất đẳng thức co l kDσ(A)X − Dσ(B)Xk = kXA + AX − XB − BXk m = kX(A − B) + (A − B)Xk n va 19 an Lu ≤ kA − Bk · kXk ac th si kDσ(A) − Dσ(B)kL(Sym) ≤ kA − Bk tính liên tục Dσ chứng minh Để kết thúc chứng minh, ta cần Dσ ∈ L(Sym) Tính tồn ánh: Với A ∈ P d Q ∈ Sym phương trình Lyapunov (hoặc Sylvester) XA + AX = Q có nghiệm X ∈ Sym từ ta có tính tồn ánh Tính đơn ánh : Với A ∈ P d Q ∈ Sym với Dσ(A)X = Chúng ta X = 0, có nghĩa Ker Dσ(A) = {0} Đặt v vectơ riêng A với giá trị riêng a > 0, có nghĩa Av = av Do AX = −XA nên ta có lu an AXv = −XAv = −aXv va Xv vectơ riêng A với giá trị riêng −a Ma trận A xác n tn to định dương nên khơng có giá trị riêng giá trị âm Do X = p ie gh Bổ đề 2.5 Ánh xạ σ: Pd → Pd, A 7→ A2 vi phôi C nghịch đảo √ σ −1 : Pd → Pd, A 7→ A liên tục oa nl w Chứng minh Chúng ta σ: Pd → Pd, A 7→ A2 , vi phôi C d (2.1) lu nf va an Dựa vào Bổ đề 2.4, σ khả vi liên tục với A ∈ P d đạo hàm Dσ(A) phép đẳng cấu tuyến tính Định lý ánh xạ nghịch lm ul đảo cho kết σ vi phôi C lân cận A z at nh oi σ(A) Do σ : P d → P d có tính song ánh, (2.1) chứng minh Giờ z chứng minh phát biểu bổ đề Với tính khả vi liên tục √ σ −1 : Pd → Pd, A 7→ A , ta có với A ∈ P d X ∈ Sym bất đẳng thức √ √ −1 A + X − A − Dσ (A)X ≤ R(X) kXk gm @ co l Với limX→0 R(X) = Từ ta thu chứng minh bổ đề m Bổ đề 2.6 Đặt P0 ∈ RN ×N phép chiếu đối xứng X : R → n va 20 an Lu GLN (R), t 7→ X(t), liên tục tuyệt đối Khi ta có: ac th si (A) Ánh xạ ˜ : R → RN ×N , t 7→ P0 X(t)T X(t)P0 +[I−P0 ]X(t)T X(t)[I−P0 ] (2.2) R ˜ liên tục tuyệt đối R(t) xác định dương, đối xứng với t ∈ R Hơn nữa, có hàm liên tục tuyệt đối R : R → RN ×N ma trận đối xứng xác định dương R(t), t ∈ R thỏa mãn ˜ R(t)2 = R(t), P0 R(t) = R(t)P0 lu (B) Ánh xạ an S : R → RN ×N , t 7→ X(t)R(t)−1 va n liên tục tuyệt đối ma trận S(t) khả nghịch thỏa mãn tn to ie gh S(t)P0 S(t)−1 = X(t)P0 X(t)−1 p thỏa mãn √ (2.3) d oa nl w kS(t)k ≤ lu an r 2 −1 −1 −1 S(t) ≤ X(t)P0 X(t) + X(t)[I − P0 ]X(t) nf va z at nh oi lm ul với t ∈ R (2.4) z Chứng minh (A) Chứng minh chia làm ba bước ˜ : R → RN ×N xác định (2.2) liên tục tuyệt đối Bước 1: Ánh xạ R ˜ xác định dương đối xứng R(t) gm @ ˜ Do P0 có tính chất đối xứng, nên rõ ràng R(t) có tính chất đối l xứng với t ∈ R Kết luận P0 X(t)T X(t)P0 P0 = m co xác định dương P0 6= 0, với v ∈ RN ta có D E T v, P0 X(t) X(t)P0 v = hX(t)P0 v, X(t)P0 vi > an Lu n va 21 ac th si ˜ xác định dương tương tự vậy, ta có kết luận thứ hai R(t) ˜ có tính tuyệt Do X có tính tuyệt đối liên tục , ta suy R đối liên tục Bước 2: Tồn hàm liên tục tuyệt đối R : R → RN ×N ma trận xác định dương, đối xứng R(t) , t ∈ R ˜ R(t)2 = R(t) ˜ ˜ Do R(t) ∈ P d, R(t) liên hợp với ma trận đường chéo với giá trị ˜ dương Do mà tồn bậc hai R(t) R(t) lu an ma trận có tính đối xứng xác định dương Mọi giá trị tính n va tn to R tổng hàm vi phân liên tục với giá trị liên tục tuyệt ˜ Một tổng liên tục tuyệt đối hàm ma trận R đối R liên tục tuyệt đối gh p ie Bước 3: Với t ∈ R ma trận R(t) R(t)−1 giao hoán với phép w chiếu P0 oa nl Xem Coppel [2] (B) Việc chứng minh chia làm hai bước d an lu Bước 1: Ánh xạ S(t) = X(t)R−1 (t) liên tục tuyệt t ∈ R nf va ma trận S(t) khả nghịch thỏa mãn z at nh oi lm ul S(t)P0 S(t)−1 = X(t)P0 X(t)−1 X(t) R−1 (t) khả nghịch, S(t) khả nghịch Mỗi phần tử X R−1 hàm liên tục tuyệt đối t ∈ R Do tổng tích z hàm liên tục tuyệt đối có tính liên tục tuyệt đối nên ta @ gm suy giá trị tính S liên tục tuyệt đối t ∈ R m co với R(t) với t ∈ R nên đồng thức l Chúng ta có S(t)P0 S(t)−1 = X(t)R(t)−1 P0 R(t)X(t)−1 Do P0 giao hoán n va 22 an Lu S(t)P0 S(t)−1 = X(t)P0 X(t)−1 ac th si thỏa mãn Bước 2: Với t ∈ R ước lượng (2.3) (2.4) Xem Coppel [2] 2.2 Phép biến đổi tương đương Xét hệ phương trình vi phân tuyến tính x˙ = A(t)x (2.5) với hàm ma trận khả tích địa phương A : R → RN ×N , N ∈ N Trước hết giới thiệu khái niệm tương đương hai phương trình vi phân lu tuyến tính có hệ số phụ thuộc vào thời gian an n va Định nghĩa 2.7 Xét hệ (2.5) hệ to (2.6) gh tn z˙ = B(t)z p ie với B : R → RN ×N khả tích địa phương Khi (2.5) (2.6) gọi oa nl cho w tương đương tồn hàm liên tục tuyệt đối S : R → GLN (R) −1 −1 |S| := sup kS(t)k < ∞ S := sup S(t) < ∞ (2.7) d t∈R z(t) = S(t)−1 x(t) nf va an lu t∈R lm ul Bổ đề 2.8 Với hàm liên tục tuyệt đối S : R → GLN (R) thỏa mãn z at nh oi điều kiện tính bị chặn (2.7) phát biểu sau tương đương (A) Hệ (2.5) (2.6) tương đương theo phép biến đổi.S z (B) Với tốn tử tiến hóa ΦA ΦB (2.5) (2.6) thỏa mãn gm @ ΦA (t, τ )S(τ ) = S(t)ΦB (t, τ ) l m co không đổi với t, τ ∈ R an Lu (C) Hàm S giải phương trình vi phân S˙ = A(t)S − SB(t) n va 23 ac th si Chứng minh Xem Bổ để 2.2 Daleckii Krein [5, Chứng minh bổ đề 2.1, trang 158] Bổ đề 2.9 Giả sử hệ (2.5) (2.6) tương đương theo phép biến đổi S Nếu với γ ∈ R hệ x˙ = [A(t) − γI]x có tính nhị phân mũ với hàng số K ≥ 1, α > phép chiếu bất biên P (t), t ∈ R, hệ z˙ = [B(t) − γI]z có tính nhị phân mũ với số |S| · S −1 · K ≥ 1, α > phép chiếu lu bất biến S(t)−1 P (t)S(t), t ∈ R an n va Chứng minh Đặt γ ∈ ρ(A) x˙ = [A(t) − γI]x nhận phép nhị phân tn to theo số mũ với số K ≥ 1, α > phép chiếu bất biến P Chúng gh ta z˙ = [B(t) − γI]z nhận phép nhị phân theo số mũ p ie Chúng ta xác định phép chiếu Q(t) := S(t)−1 P (t)S(t) Sử dụng Bổ đề nl w 2.8 ta dễ dàng thấy d oa Q(t)ΦB (t, s) = ΦB (t, s)Q(s) với t, s ∈ R lu an Q tốn tử bất biến với (2.6) Các bất đẳng thức nhị phân nf va với x˙ = [A(t) − γI]x lm ul