Hàm mũ ma trận và ứng dụng với hệ phương trình vi phân tuyến tính cấp một

TRƯỜNG ĐẠI HỌC SƯ PHẠM HÀ NỘI 2KHOA TOÁN Nguyễn Thị Lan Anh HÀM MŨ MA TRẬN VÀ ỨNG DỤNG ĐỐI VỚI HỆ PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH CẤP MỘT KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC Hà Nội – Năm 2016..

TRƯỜNG ĐẠI HỌC SƯ PHẠM HÀ NỘI 2

KHOA TOÁN

Nguyễn Thị Lan Anh

HÀM MŨ MA TRẬN VÀ ỨNG DỤNG ĐỐI VỚI

HỆ PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH CẤP MỘT

KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC

Hà Nội – Năm 2016

TRƯỜNG ĐẠI HỌC SƯ PHẠM HÀ NỘI 2

KHOA TOÁN

Nguyễn Thị Lan Anh

HÀM MŨ MA TRẬN VÀ ỨNG DỤNG ĐỐI VỚI

HỆ PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH CẤP MỘT

Chuyên ngành: Toán giải tích

KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC

NGƯỜI HƯỚNG DẪN KHOA HỌC:

TS Trần Văn Bằng

Hà Nội – Năm 2016

Em xin bày tỏ lòng biết ơn sâu sắc TS Trần Văn Bằng - Người đãtận tình hướng dẫn và giúp đỡ em để em hoàn thành bài khóa luận củamình Đồng thời em cũng xin trân thành cảm ơn các thầy cô trong tổGiải tích và các thầy cô trong khoa Toán - Trường Đại học Sư phạm HàNội 2 Ban chủ nhiệm khoa Toán đã tạo điều kiện cho em hoàn thành tốtbài khóa luận này.

Trong khuôn khổ của một khóa luận, do điều kiện thời gian, do trình

độ có hạn và cũng là lần đầu tiên nghiên cứu khoa học cho nên khôngtránh khỏi những hạn chế những thiếu sót nhất định Vì vậy, em kínhmong nhận được sự đóng góp của các thầy cô và các bạn

Em xin chân thành cảm ơn!

Hà Nội, ngày 2 tháng 05 năm 2016

Sinh viênNguyễn Thị Lan Anh

Khóa luận này là kết quả của bản thân em đạt được trong quá trìnhhọc tập và nghiên cứu, dưới sự chỉ dẫn của TS Trần Văn Bằng và sựgiúp đỡ của các Thầy, Cô trong khoa Toán Trường Đại học Sư phạm HàNội 2 và của các Thầy, Cô trực tiếp giảng dạy em.

Trong khi nghiên cứu, hoàn thành bản khóa luận này em đã tham khảomột số tài liệu đã ghi trong phần tài liệu tham khảo

Em xin khẳng định kết quả của đề tài "Hàm mũ ma trận và ứngdụng đối với hệ phương trình vi phân tuyến tính cấp một"

là kết quả của việc nghiên cứu, học tập và nỗ lực của bản thân, không có

sự trùng lặp với kết quả của các đề tài khác

Hà Nội, ngày 2 tháng 05 năm 2016

Sinh viênNguyễn Thị Lan Anh

1.4 Hệ phương trình vi phân tuyến tính cấp một 11

2 Hàm mũ ma trận và ứng dụng đối với hệ phương trình

Lời mở đầu

Phương trình vi phân là một trong những chuyên ngành của toán học

nó có nhiều ứng dụng trong các lĩnh vực khoa học và công nghệ, nó được

coi là cầu nối giữa lí thuyết và ứng dụng Trong đó lí thuyết hệ phương

trình vi phân tuyến tính cấp một là một lí thuyết quan trọng trong lí

thuyết phương trình vi phân Việc giải một hệ phương trình vi phân

dù là tuyến tính thì nói chung cũng không đơn giản Trong khóa luận

này tôi muốn tìm hiểu phương pháp hàm mũ ma trận để giải hệ phương

trình vi phân tuyến tính Đây là một phương pháp cho ta những công

thức biểu diễn nghiệm của hệ rất gọn và đẹp

Xuất phát từ nhận thức trên và lòng ham mê môn học với sự hướng

dẫn tận tình của thầy giáo T.S Trần Văn Bằng, em mạnh dạn chọn

đề tài: "Hàm mũ ma trận và ứng dụng đối với hệ phương trình

vi phân tuyến tính cấp một" để thực hiện khóa luân tốt nghiệp của

mình Nội dung của khóa luận bao gồm:

Chương 1 Kiến thức chuẩn bị

Chương 2 Hàm mũ ma trận và ứng dụng đối với hệ phương

trình vi phân tuyến tính cấp một

Do lần đầu tiên thực tập nghiên cứu, thời gian có hạn và năng lực

bản thân còn hạn chế nên chắc chắn bài nghiên cứu này khó tránh khỏi

những thiếu sót Em mong nhận được sự đóng góp, ý kiến của thầy cô

và các bạn để đề tài này được hoàn chỉnh và đạt được kết quả cao hơn

Em xin chân thành cảm ơn!

Kiến thức chuẩn bị

Chương này trình bày một số kiến thức cơ bản về ma trận, sự hội tụ

trong không gian định chuẩn Mat(n ∗ n, K), hệ phương trình vi phân cấpmột và hệ phương trình vi phân tuyến tính cấp một nhằm thuận tiện

Nếu x = (x1, , xn) là một vectơ n chiều thì ta có thể xem x như làmột ma trận n hàng, một cột và do đó

Dễ dàng thấy chuẩn của ma trận có các tính chất sau:

(i) kA + Bk ≤ kAk + kBk

(ii) kABk ≤ kAk kBk

(iii) kAxk ≤ kAk kxk

Ma trận đơn vị n chiều được kí hiệu là In (hay đơn giản là I khikhông sợ nhầm lẫn) Đa thức det(Iλ − A) bậc n của λ gọi là đa thức đặc

trưng của ma trận A, nghiệm của nó gọi là giá trị riêng của ma trận A

Hai ma trận vuông cấp n là A và B được gọi là đồng dạng nếu tồn

tại một ma trận vuông cấp n không suy biến P sao cho B = P AP−1.Nếu A và B đồng dạng thì chúng có chung một đa thức đặc trưng vì

det(λI−B) = det(P (λI−A)P−1) = detP det(λI−A)detP−1 = det(λI−A)

Đặc biệt các hệ số của lũy thừa λ của đa thức det(λI − A) là bất

biến đối với phép biến đổi đồng dạng Hai bất biết quan trọng nhất đối

với phép đồng dạng là detA và trace(A) (tức là định thức và vết của ma

trận A).

Giả sử cho f là một phép biến đổi tuyến tính của không gian n chiều

Rn trên trường C vào chính nó

f : Rn → Rn.Giả sử h = {h1, , hn} là một cơ sở của không gian Rn Khi đó

tính f trong Rn tồn tại một cơ sở sao cho ma trận A của phép biến đổi

trong đó Ki(i = 1, 2, , p) là ma trận có dạng

(tức là các phần tử thuộc đường chéo chính đều bằng λi và các phần tử

kề trên đều bằng 1, các phần tử khác đều bằng 0 Cấp của ma trận này

ta kí hiệu là ni)

Dạng trên đây được gọi là dạng chính tắc Jordan của ma trận A

Do đó

(i) Với bất kì ma trận vuông B nào cũng tồn tại ma trận vuông không

suy biến P sao cho P BP−1 = A trong đó A là ma trận dạng chínhtắc Jordan

(ii) Giả sử A là ma trận phức và λ1, λ2, , λk là các giá trị riêng (phức)phân biệt của A với bội tương ứng là m1, , mk Khi đó A đồngdạng với ma trận J :

trong đó

Nếu λi là các giá trị riêng đơn với ∀i = 1, k thì ma trận A đồngdạng với ma trận chéo

ma trận Ji có dạng Ji = λiImi + Zi trong đó Imi là ma trận đơn vịcấp mi và

tức là so với Zi đường chéo đơn vị trong Zi2 bị dịch về bên phải mộtđơn vị, còn các phần tử còn lại đều bằng không Suy ra Zmi

Kí hiệu: lim

m→∞Am = A hay Am → A khi m → ∞

Nhận xét 1.1 Dãy {Am} hội tụ khi và chỉ khi mỗi dãy các phần tử

của nó hội tụ Hơn nữa nếu Am = (a(m)ij ), A = (aij) thì

lim

m→∞Am = A ⇔ lim a(m)ij = aij(∀i, j =1, n) Thật vậy

(i) Nếu mỗi dãy các phần tử của {Am} hội tụ thì hiển nhiên dãy matrận đó hội tụ

(ii) Ngược lại, giả sử ta có lim

m→∞Am = A = (aij)m×n thế thì(∀ε > 0)(∃Nε) sao cho ∀m ≥ Nε thì kAm − Ak < ε

a(m)ij − aij

< ε ⇔

a(m)ij − aij

... data-page="21">

Hàm mũ ma trận ứng dụng hệ phương trình vi phân tuyến tính cấp một< /h2>

Chương trình bày khái niệm hàm mũ ma

trận tính chất mệnh đề từ đưa cấu trúc

nghiệm để giải hệ phương. .. độc lập tuyến tính (1.6) Khi đó, ma trận vng X(t) có

các cột x1, x2, , xn gọi ma trận hệ (1.6)

Dễ thấy X ma trận hệ (1.6) C ma trận. .. bởi

x(t, t0, x0) = Xt0(t)x0.Tổng quát, ta gọi nghiệm phương trình ma trận X(t) =A(t)X ma trận nghiệm hệ (1.6)