Luận văn Thạc sĩ Toán học: Giá trị đầu của nghiệm bị chặn của phương trình vi phân tuyến tính với hàm ràng buộc tuần hoàn

Chia sẻ: Lavie Lavie | Ngày: | Loại File: PDF | Số trang:47

Thêm vào BST

Báo xấu

74
lượt xem 5
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Luận văn Thạc sĩ Toán học: Giá trị đầu của nghiệm bị chặn của phương trình vi phân tuyến tính với hàm ràng buộc tuần hoàn trình bày tổng quan; biểu diễn nghiệm của phương trình vi phân tuyến tính và sai phân tuyến tính; phân loại đầy đủ của tập những giá trị đầu theo tính bị chặn, tính tuần hoàn và dáng điệu tiệm cận của nghiệm.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Luận văn Thạc sĩ Toán học: Giá trị đầu của nghiệm bị chặn của phương trình vi phân tuyến tính với hàm ràng buộc tuần hoàn

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM THÀNH PHỐ HỒ CHÍ MINH __________________________ Nguyễn Đình Tường Long GIÁ TRỊ ĐẦU CỦA NGHIỆM BỊ CHẶN CỦA PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH VỚI HÀM RÀNG BUỘC TUẦN HOÀN Chuyên ngành: Toán Giải tích Mã số: 60 46 01 LUẬN VĂN THẠC SĨ TOÁN HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS. Lê Hoàn Hoá Thành phố Hồ Chí Minh - 2007
LỜI NÓI ĐẦU Tính chất tiệm cận của nghiệm của phương trình vi phân tuyến tính x  t   Ax  t   f  t  (1), trong đó A  M d   và f :  d là hàm không tầm thường liên tục tuần hoàn chu kì  , đã được nhiều tác giả nghiên cứu. Trong luận văn này, chúng tôi trình bày việc phân loại đầy đủ của tập những giá trị đầu theo dáng điệu tiệm cận của nghiệm phương trình (1). Đặc biệt, chúng tôi chỉ ra việc phân loại đầy đủ của tập những giá trị đầu theo tính chất bị chặn và tính chất tuần hoàn. Nội dung luận văn được trích trong một báo cáo ở Hội nghị Quốc tế về phương trình vi phân tại trường Đại học Khoa học Tự nhiên Tp. Hồ Chí Minh năm 2004. Trong báo cáo này, nhiều chi tiết, bổ đề và định lí chỉ được phát biểu ngắn gọn cho nên chúng tôi cố gắng làm sáng tỏ các chi tiết và chứng minh đầy đủ các định lí và bổ đề. Ngoài ra, luận văn còn chứng minh ý tưởng phát triển mà tác giả đã đề xuất cuối bài báo cáo. Luận văn này ngoài lời nói đầu, phần kết luận, phần tài liệu tham khảo và mục lục sẽ được trình bày trong 5 chương. Chương I là phần tổng quan. Chương II là phần trình bày các kí hiệu và nhắc lại một số kiến thức cần thiết để sử dụng cho các chương sau. Chương III là phần trình bày việc biểu diễn nghiệm của phương trình vi phân tuyến tính và sai phân tuyến tính. Chương IV là phần trình bày việc phân loại đầy đủ của tập những giá trị đầu theo tính bị chặn, tính tuần hoàn và dáng điệu tiệm cận của nghiệm. Chương V là phần phát triển. Cho phép tôi được gửi lời biết ơn sâu sắc đến PGS.TS Lê Hoàn Hoá, người thầy đã tận tình giúp đỡ, hướng dẫn tôi trong quá trình học tập cũng như trong quá trình thực hiện luận văn này. Tôi xin chân thành cảm ơn quý Thầy cô đã tận tình truyền thụ kiến thức cho tôi trong quá trình học tập tại trường và quý Thầy cô trong Hội đồng bảo vệ luận văn đã dành nhiều thời gian để đọc và cho tôi nhiều ý kiến sâu sắc và quí báu. Tôi cũng xin trân trọng cảm ơn đến quý Thầy cô phòng KHCN- SĐH đã nhiệt tình giúp đỡ chúng tôi hoàn tất các thủ tục bảo vệ luận văn. Cuối cùng, cho tôi cảm ơn đến tất cả bạn bè, đồng nghiệp luôn luôn động viên, cổ vũ tôi trong quá trình học tập và nghiên cứu. Tp. Hồ Chí Minh, tháng 4 năm 2007 Tác giả luận văn Nguyễn Đình Tường Long.
Chương 1: TỔNG QUAN Trong luận văn này, chúng tôi sẽ nghiên cứu việc xác định giá trị đầu của nghiệm bị chặn hay nghiệm tuần hoàn của phương trình vi phân tuyến tính sau: x  Ax  f  t  (1) với A  M d   và f :  d là hàm không tầm thường liên tục tuần hoàn chu kì  .  2 , b f   e  s  A Đặt   f  s  ds (2)  0 Nhiều tác giả nghiên cứu tính chất tiệm cận của nghiệm của (1). Chẳng hạn, Massera đã đưa ra kết quả sau: Định lí: Nếu phương trình (1) có nghiệm bị chặn trên 0,+  , thì phưong trình (1) tồn tại nghiệm tuần hoàn chu kì  . Đối với nghiệm tuần hoàn, nó hầu như tầm thường rằng: một nghiệm của phương trình (1) là tuần hoàn chu kì  nếu và chỉ nếu giá trị đầu w  x  0  thoả mãn:  E  e w  b A f . (3) Chẳng hạn, giả sử A   E  N , với   i , i 2  1 và N là ma trận luỹ linh thỏa N m  0  m  1 . Khi đó: m 1 k m 1 k   A  E   N k A  k e e k 0 k ! k 0 k ! và phương trình:  E  e A  w  b f được biểu diễn là: m 1  k    N k w  bf k 1  k !  Từ sự biểu diễn trên, ta có thể tính toán lần lượt mỗi thành phần của nghiệm w của phương trình với b f đã cho. Tuy nhiên, trong trường hợp này nó không rõ ràng để chỉ ra rằng nghiệm w được biểu diễn bởi: 1 m1 k Nw     Bk k 0 k! N kbf (4) với Bk là số Bernoulli. Công thức (4) là một trong những kết quả của chúng tôi sau này. Mục đích của luận văn là đưa ra một sự phân loại đầy đủ của tập những giá trị đầu theo dáng điệu tiệm cận của nghiệm của phương trình (1). Đặc biệt, chúng tôi chỉ ra sự phân loại đầy đủ của tập những giá trị đầu theo tính chất bị chặn và tính chất tuần hoàn.
Phương pháp và ý tưởng cơ bản của luận văn sẽ được giải thích đơn giản trên phương trình vi phân tuyến tính một chiều: dx  ax  t   f  t  , x  0   w dt với a là hằng số và f là hàm vô hướng liên tục tuần hoàn chu kì  . Sau đó, chúng tôi sẽ tổng quát hoá những kết quả đó đến phương trình vi phân tuyến tính (1). Phần còn lại của luận văn bao gồm: chương 2, chương 3, chương 4, chương 5 và chương kết luận. Chương 2: Kiến thức chuẩn bị. Trong chương này, chúng tôi trình bày một số kiến thức cơ bản về nghiệm của phương trình sai phân tuyến tính, nghiệm của phương trình vi phân tuyến tính, số Bernoulli, không gian riêng suy rộng, … để sử dụng cho các chương sau. Chương 3: Biểu diễn nghiệm của phương trình vi phân tuyến tính. Trình bày phương pháp và ý tưởng cơ bản của luận văn trên phương trình vi phân tuyến tính một chiều. Sau đó, dựa trên phương pháp và ý tưởng này để xây dựng cách biểu diễn nghiệm của phương trình vi phân tuyến tính (1). Hơn nữa, chúng tôi cũng chỉ ra rằng cách biểu diễn nghiệm này có sự tương ứng với cách biểu diễn nghiệm của phương trình sai phân tuyến tính. Chương 4: Dáng điệu tiệm cận, tính bị chặn, tính tuần hoàn. Chương này bao gồm hai mục: 4.1 và 4.2. Mục 4.1. Trình bày sự phân loại đầy đủ của tập những giá trị đầu theo tính bị chặn và tính tuần hoàn của nghiệm của phương trình vi phân tuyến tính (1). Sự phân loại này, có sự tương ứng đến sự phân loại của tập những giá trị đầu theo tính bị chặn của nghiệm phương trình sai phân tuyến tính. Mục 4.2. Trình bày sự phân loại đầy đủ của tập những giá trị đầu theo dáng điệu tiệm cận của nghiệm phương trình (1). Chương 5: Phát triển. Chương này trình bày về cách biểu diễn và tính bị chặn của nghiệm của phương trình sai phân tuyến tính tổng quát xn 1  Bxn  b và của phương trình vi phân tuyến tính x  t   A  t  x  t   f  t  có sự tương ứng với nhau. Chương kết luận nêu vắn tắt một số kết quả trình bày trong luận văn và những hướng nghiên cứu dự định trong tương lai. Cuối cùng là tài liệu tham khảo.
Chương 2: KIẾN THỨC CHUẨN BỊ VÀ CÁC KÍ HIỆU Trong chương này, chúng tôi qui ước một số kí hiệu và nêu lại một số kiến thức chuẩn bị cần thiết để sử dụng đến trong các chương sau. 2.1. KÍ HIỆU. 1. Tập các số thực 2. Tập các số phức 3. Tập các số nguyên 4. E Ma trận đơn vị 5. Re- Phần thực, Im- Phần ảo 6. IdV Ánh xạ đồng nhất từ V vào V. 7. N(A) Nhân của A 8. R(A) Miền giá trị của A n 9. f  n   z   f  z  Đạo hàm cấp n của hàm số y  f  z  , f    z   f  z  . 0 z n 10. Bk , k  0,1, 2,... Số Bernoulli 11. M d   Tập các ma trận phức d  d chiều 12. X   xik mn Ma trận cấp m  n mn   X   xik mn / xik  , i  1, m, k  1, n  13. Lloc  I , mn  Tập các ma trận hàm cấp m  n khả tích trên mọi tập compact nằm trong I  . 14.  Kết thúc phần chứng minh. 2.2. SỐ BERNOULLI. Dãy Bk của hệ số của chuỗi luỹ thừa được xác định bởi  z zk 1   Bk với B1   và B2 k 1  0 với k > 0. e  1 k 0 k ! z 2 Số Bernoulli Bk có tính chất sau: 1 n   1 Cnk1Bk  1 , n = 0, 1, 2, … . k n  1 k 0 2.3. NGHIỆM CỦA PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH. Xét hệ phương trình vi phân tuyến tính
x  t   Ax  t   f  t  (2.3.1) x  0  w (2.3.2) với A  M n   , f  t   Lloc  I , n , w n . Bài toán (2.3.1), (2.3.2) có duy nhất một nghiệm. Ta tìm nghiệm của hệ (2.3.1), (2.3.2) theo phương pháp biến thiên hằng số. Nghiệm của hệ (2.3.1), (2.3.2) viết dưới dạng: x  t   etA y  t  (2.3.3) với y  t  cần xác định để (2.3.3) là nghiệm của (2.3.1), (2.3.2). Ta có: w  x  0  y  0 (2.3.4) Lấy đạo hàm hai vế (2.3.3), ta có: x  t   AetA y  t   y  t  etA Thay vào (2.3.1), ta được: y  t   etA f  t  (2.3.5) Vậy y  t  là nghiệm của (2.3.5) thoả điều kiện (2.3.4), suy ra t y  t   w   e sA f  s  ds . 0 Thay vào (2.3.3), ta có:  t  x  t   etA  w   e sA f  s  ds     0  t Hay x  t   etAw   et  s  A f  s  ds . 0 Vậy nghiệm của hệ (2.3.1), (2.3.2) là: t x  t   etAw   e t s A f  s  ds . 0 2.4. NGHIỆM CỦA PHƯƠNG TRÌNH SAI PHÂN TUYẾN TÍNH. Xét phương trình sai phân tuyến tính xn  qxn 1  b , n = 1, 2, … x0  a với b, q  0 là các hằng số.
Nghiệm tổng quát của phương trình sai phân tuyến tính là xn  xn  xn* , trong đó xn là nghiệm tổng quát của phương trình sai phân tuyến tính thuần nhất xn  qxn 1 và có dạng xn  cq n  c  0  , còn nghiệm n 1 riêng xn* của phương trình sai phân tuyến tính không thuần nhất có dạng xn*  c0 q n  b q k . k 0 Do đó nghiệm của phương trình sai phân tuyến tính là: n 1 xn  aq n  b q k k 0 2.5. KHÔNG GIAN RIÊNG SUY RỘNG. 2.5.1. Vec tơ riêng và giá trị riêng của một tự đồng cấu.   Cho  :V  V là một tự đồng cấu của không gian vectơ V trên trương K. Vectơ x  0 của V mà     x   x , với một  nào đó của K, gọi là một vectơ riêng của  ứng với giá trị riêng  của  .  Tập tất cả các vectơ riêng của  ứng với cùng một giá trị riêng   K cùng với vectơ 0 làm thành một không gian vectơ của V, kí hiệu là  và gọi là không gian riêng ứng với giá trị riêng  của  ;  chẳng qua là   N    IV  . Đa thức đặt trưng của tự đồng cấu  kí hiệu là P   và được xác định bởi P     det    IdV  . 2.5.2. Không gian riêng suy rộng của một tự đồng cấu. a. Định nghĩa: Cho tự đồng cấu  của K-không gian vectơ hữu hạn chiều V. Với mỗi   K , xét tập x V / m  0, m   x   0 thì đó là một không gian vectơ con của V. Khi nó khác 0 thì     :    IdV  m nó được gọi là không gian riêng suy rộng của  ứng với   K và kí hiệu M  . b. Tính chất: 1. Rõ ràng nếu  là một giá trị riêng của  thì không gian riêng   M  . Với mọi không gian riêng suy rộng M  ,  phải là một giá trị riêng của  . 2. Với mỗi giá trị riêng  của  , dim  gọi là số chiều hình học của  , dim M  gọi là số chiều đại số của  . Số chiều đại số của  bằng bội của nghiệm  của đa thức đặc trưng . c. Định lí:  là một tự đồng cấu của một không gian vectơ hữu hạn chiều V trên trường K mà đa thức đặc trưng P phân tích được thành tích những nhân tử tuyến tính: P        1  1    2  2 ...    m  m . s s s ( các k (k  1, 2,..., m) khác nhau từng cặp) thì V là tổng trực tiếp các không gian riêng suy rộng:
V  M 1  M 2  ...  M m . 2.6. Định lí 2.6: Nếu f :  nn là một ánh xạ liên tục và tuần hoàn chu kì  thì f bị chặn trên . Chứng minh định lí 2.6. Đặt K  0,   . Khi đó, K là tập compact trong . Vì f là ánh xạ liên tục nên f  K  là tập compact trong nn . Do đó, f bị chặn trên K. Vì f tuần hoàn chu kì  nên suy ra f bị chặn trên .
Chương 3: BIỂU DIỄN NGHIỆM CỦA PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH 3.1. PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH MỘT CHIỀU. Phương pháp và ý tưởng cơ bản của luận văn sẽ được trình bày với phương trình vi phân tuyến tính một chiều: dx  ax  t   f  t  , x  0   w (5) dt với a là hằng số và f là hàm vô hướng liên tục tuần hoàn chu kì  . Nghiệm của phưong trình (5) được biểu diễn: t x  t  : x  t;0,    eatw   eat  s  f  s  ds (6) 0 Tuy nhiên, cách biểu diễn trên không dể thấy dáng điệu tiệm cận của nghiệm. Do đó, chúng ta cần đưa ra một sự biểu diễn khác của nghiệm mà với cách biểu diễn này chúng ta dể dàng nhận biết được dáng điệu tiệm cận của nghiệm. Định lí 3.1: Nghiệm x  t  của phưong trình (5) được biểu diễn như sau: 1) Giả sử e a  1 . Nếu công thức (6) được viết lại:  1  x  t   eat  w   1 e b  u t, bf a f     (7) 1     là nghiệm tuần hoàn chu kì  . t a t s với u t , b f  eat b   e   f  s  ds , thì u t , b f a f 1 e 0 2) Giả sử ea  1 . Nếu công thức (6) được viết lại:  t  x  t   eat  w  b f   v t , b f      (8) t t   với v t , b f  e at  bf   e at  s  f  s  ds , thì v t , b f   là một hàm tuần hoàn chu kì  , 0 tuy nhiên, nó không cần thiết là một nghiệm của (5). Chứng minh định lí 3.1: Ta có: t   t  a  t   s  a  t   s  a  t   s  0 e f  s  ds   e f  s  ds  e f  s  ds 0 
 t  eat  ea  s  f  s  ds   eat  z  f  z  dz ( đặt z = s -  ) 0 0 t at  s   e bf   e at f  s  ds 0 1) Giả sử ea  1 . Khi đó: t  a  t   1 a  t   s  u t     e b  0 e f  s  ds 1  ea f ea t at  s   eat a b f  e b f   e at f  s  ds 1 e 0 t 1 at  s  e at a b f   e f  s  ds  u  t  . 1 e 0 2) Giả sử ea  1 . Khi đó: t  a  t   t  a  t   s  v  t     e bf  0 e f  s  ds  t t  eat b f  eat b f  eat b f   eat  s  f  s  ds  v  t  .  0 Từ sự biểu diễn nghiệm trong định lí 3.1, ta có thể dễ dàng đạt được dáng điệu tiệm cận của nghiệm như sau: Hệ quả 3.1.1: I) Trường hợp ea  1 . 1 1) Giả sử w  b . 1  ea f Khi đó nghiệm x(t;0, w ) là nghiệm tuần hoàn chu kì  , được xác định bởi x  t;0, w   u t , b f .   1 2) Giả sử w  b . 1  ea f   tuần hoàn chu kì  khi t   . a) Nếu Re a  0 , thì nghiệm x(t;0, w ) hội tụ đến nghiệm u t , b f b) Nếu Re a  0 , thì nghiệm x(t;0, w ) hội tụ đến nghiệm u  t , b  tuần hoàn chu kì  khi t   .f c) Nếu Re a  0 , thì nghiệm x(t ;0, w ) là á tuần hoàn. II) Trường hợp ea  1 .
1) Nếu b f  0 , thì mỗi nghiệm x(t ;0, w ) là một nghiệm tuần hoàn chu kì  . 2) Nếu b f  0 , thì tất cả nghiệm x(t;0, w ) là không bị chặn. Ý tưởng cơ bản đã được phát biểu một cách đơn giản trong định lí 3.1. Chúng tôi sẽ tổng quát hoá 1 những kết quả ở trên đến phương trình (1). Khi đó vấn đề là tìm những hàm tưong ứng với eat b , 1  ea f t eat b f trong (7) và (8).  Để làm được như vậy, ta viết lại nghiệm của (6): x  t;0, w   eatw  z  t   g  t  t với g  t   z  t    eat  s  f  s  ds, và sẽ tìm một hàm z(t) sao cho g(t) trở nên một hàm tuần hoàn chu kì  . 0 t  t a  t   s  a t  s  g(t) tuần hoàn chu kì   z (t   )  0 e f  s  ds  z  t    e f  s  ds 0 t  t Vì 0 e at   s  f  s  ds  eat b f   eat  s  f  s  ds nên ta có: 0 z  t     z  t   eat b f Đặt  z  t  : z  t     z  t   e at b f (9) Nếu z(t) là một nghiệm của phương trình (9), khi đó, với 0  s   , n  1, 2,..., ta có: z  s  n   z  s   n  1   e  a s   n 1  bf  z  s   n  2    e  a s   n  2   bf  e  a s   n 1  bf n 1 = . . . = z  s    ea s  j b f . j 0 n 1 Vậy z  s  n   z  s   e as  eaj b f . (10) j 0  Giả sử ea  1 , khi đó từ (10) cho ta: eas  ea s  n  n 1 1  ean z  s  n   z  s   1  ea b f ( vì  j 0 e aj  1  ea ) eas ea s  n   z s  b  b . 1  e a f 1  ea f
eas eat Do đó, nếu chọn z  s   b với 0  s   , thì z  t   b (t  s  n ) . 1  ea f 1  e a f Nếu ea  1 , khi đó từ (10) cho ta:  s  n s z  s  n   z  s   neasb f  z  s      easb f     s s  n as  n  z  s   easb f  e b f (vì ea  1 ).   s t Do đó, nếu chọn z  s    easb f với 0  s   , thì z  t    et b f ( t  s  n ).   n 1 Lưu ý rằng, trong mối quan hệ (10), tổng  j 0 eaj b f là nghiệm tổng quát của phương trình sai phân: xn  e a xn 1  b f , x0  0 . 3.2. TỔNG VÔ HẠN. Cho   0 và c :  d là một hàm liên tục. Ta tìm nghiệm liên tục của phương trình:  z  t   z  t     z  t   c  t  , t  (11) đó là, tổng vô hạn z  t   1c  t  . Nói chung, z(t) bao gồm một hằng số tuần hoàn, nghĩa là, nếu z(t) là một tổng vô hạn, khi đó z(t) + c(t) cũng là một tổng vô hạn đối với bất kì hàm c(t) tuần hoàn chu kì  nào. Bổ đề 3.1.1: Cho c :  d là một hàm liên tục, và hàm  : t0   , t0   d liên tục thỏa mãn:   t0     t0     c  t0    . Khi đó phương trình (11) với điều kiện đầu z  t     t  , t  t0   , t0  , có duy nhất trên một nghiệm liên tục z(t). Mặt khác, nghiệm z(t) sẽ được xác định như sau: Cho s  t0   , t0  . Khi đó n 1 z  s  n     s    c  s  k , (n = 0, 1, 2,… ). k 0 n z  s  n     s    c  s  k , ( n = 1, 2,… ). k 1 Chứng minh bổ đề 3.1.1: Vì z  t     t  , t  t0   , t0  , nên   t0     z  t0    và   t0   z  t0  Do đó, từ   t0     t0     c  t0    ta được: z  t0   z  t0     c  t0   
Vậy phương trình (11) có nghiệm z(t) trên . Khi đó, với s  t0   , t0  , n = 1, 2, …, ta có: z  s  n   z  s   n  1   c  s   n  1   z  s   n  2    c  s   n  2    c  s   n  1  n 1 n 1  ...  z  s    c  s  k     s    c  s  k  k 0 k 0 z  s  n   z  s   n  1   c  s  n   z  s   n  2    c  s   n  1   c  s  n  n n  ...  z  s    c  s  k     s    c  s  k  . k 1 k 1 Vậy nghiệm z(t) được xác định là: với s  t0   , t0  , n 1 z ( s  n )    s    c  s  k ,(n  0,1, 2...) k 0 n z  s  n     s    c  s  k ,  n  1, 2,... k 1 Vì  liên tục trên t0   , t0  và c liên tục trên nên nghiệm z(t) liên tục trên . * Chứng minh tính duy nhất nghiệm: Giả sử phương trình (11) còn có nghiệm z  t  . Đặt z  t   z  t   z  t  , t  . Chứng minh tương tự trên, nghiệm z  t  được xác định như sau: với s  t0   , t0  , n 1 z( s  n )    s    c  s  k ,(n  0,1, 2...) k 0 n z  s  n     s    c  s  k ,  n  1, 2,... k 1 Suy ra z  t   0 , t  hay z  t   z  t  , t  . Chúng ta xét trường hợp c(t) được xác định bởi: z t     z t   B t  b , t  (12) với B(t), t  , là một ma trận hàm liên tục và không suy biến thỏa
B  s  k   B  s  B k   , k  . Định lí 3.2: Cho xm  0  , m  , là nghiệm của phương trình sai phân: xm  B   xm1  b, x0  0 . (13) Khi đó đối với một hàm ban đầu liên tục  : t0   , t0   d thỏa mãn:   t0     t0     B  t0    b Phương trình (12) với điều kiện đầu z  t     t  , t  t0   , t0  , có một nghiệm duy nhất z(t) trên , được xác định bởi: Với s  t0   , t0  , z  s  n     s   B  s  xn  0  , (n = 0, 1, 2,…) (14) z  s  n     s   B  s  x n  0  , (n = 1, 2,… ) (15) Chứng minh định lí 3.2: Phương trình sai phân (13) có nghiệm tổng quát là: n 1 xn  0    B k   b với n = 1, 2, … k 0 n x n  0    B  k   b với n = 1, 2, … k 1 Áp dụng bổ đề 3.1.1, phương trình (12) có nghiệm duy nhất z(t) trên và n 1 z  s  n     s    B  s  k b (n = 0, 1, 2, …) k 0 n 1    s   B  s   B k   b k 0    s   B  s  xn  0  n z  s  n     s    B  s  k b (n = 1, 2, …) k 1 n    s   B  s   B  k  b k 1    s   B  s  x n  0  .  Chú ý: B  t   etA , A  M d   , thỏa mãn điều kiện B  s  k   B  s  B k   , k  , trong định lí 3.2.
3.3. PHƯƠNG TRÌNH VI PHÂN TUYẾN TÍNH. Ta xem xét trở lại việc biểu diễn nghiệm của phương trình x  Ax  t   f  t  (1) Để áp dụng ý tưởng ở trên đối với phương trình (1) với x  0   w là điều kiện đầu của nghiệm, ta sẽ tịnh tiến nghiệm x  t;0, w  của phương trình (1) như sau: x  t;0, w   e Atw  z  t   g  t  (16) t với g  t   z  t    e At  s  f  s  ds (17) 0 Ta sẽ tìm một hàm z(t) sao cho g(t) là một hàm liên tục và tuần hoàn chu kì  . Rõ ràng, g(t) là một hàm liên tục và tuần hoàn chu kì  khi và chỉ khi z(t) là một nghiệm liên tục của phương trình: z  t     z  t   e At b f , t  , đó là  z  t   1 e At b f , t   . Vì g(t) là một hàm tuần hoàn chu kì  và số hạng thứ hai của vế phải (17) được xác định trên , cho nên ta chỉ cần xét z(t) trên 0,   .   * Ta tính 1 e At b f , t  o,   . Giả sử   A là phổ của A, bao gồm r phần tử i với chỉ số mi , i = 1,…,r. Đặt M i  N  A   E  i mi là không gian riêng suy rộng của A ứng với giá trị riêng i . Khi đó ta có sự phân tích tổng trực tiếp d  M 1  M 2  ...  M r d Đặt Pi là phép chiếu từ lên M i cảm sinh từ sự phân tích này. Để tính toán đơn giản, ta gọi  là một trong các giá trị riêng của A và m là số chỉ số của  . Xét phương trình P z  t     P z  t    P etAb f (18) Đặt z  s     s  , s    ,0 . Nếu t  s  n , s    ,0 , n = 1, 2, …, khi đó từ định lí 3.2, phương trình (18) có nghiệm duy nhất P z  t  và được biểu diễn bởi P z  t   P z  s  n   P  s   P e sA xn  0  (19)
với xn  0  là nghiệm của phương trình sai phân xn  e A xn 1  b f , x0  0 (20) và điều kiện đối với hàm ban đầu trở thành P  0   P     e A P b f . 3.3.1. TRƯỜNG HỢP RỜI RẠC. Xét phương trình sai phân xn  e A xn 1  b , x0  w (21) Nghiệm tổng quát của phương trình (21) là xn  en A x0  Sn  e A  b n 1 với Sn  e A    ek A , n  1 . k 0  Ta tính mỗi thành phần P xn của nghiệm xn . P xn thỏa mãn phương trình P xn  en A P x0  Sn  e A  P b  ti etA P  et e A  E  t P  et   A   E  P i Ta có: i 0 i ! m 1 i t Vì  A   E  P  0 với i  m nên etA P  et   A   E  P . i i i 0 i !  k  i n 1 n 1 m 1 Sn  e  P   e P   e   A  E A k A k i Do đó P k 0 k 0 i 0 i! m 1 n 1   i     k i ek   A  E i P . i 0  k 0  i! n 1 Đặt E in  z  :  k i ekz . k 0   i kz * Nếu e  1 , ta có: k e  i e . i kz z n 1  i kz  i n 1  i enz  1 Suy ra E  z  :  i e  i i n  e  i z kz . k  0 z z k 0 z e  1 1 Đặt   z  : , e 1 z
Khi đó Eni  z  : i z i    z   e nz  1  i  zi   z  e  nz  i z i   z i i Vì z i   z  e nz   j 0 Ci j n j enz i  j   z  nên i Eni  z    Ci j n j e nz i  j   z    i   z  . j 0 Suy ra m 1  i  i Sn  e A  P     Ci j n j en   i j      i     A  E i P i 0  j 0  i! m 1 i i m 1 i   en  Ci j n j i  j     A  E P    i    A  E i i P i 0 j 0 i! i 0 i! Mặt khác m 1 i i  Ci j n j      A  E  i j i P i  0 j 0 i! m 1 i 1  m2 i      n  A  E   1!  A   E   n  A  E i i i P   i P i 0 i! i 0 i!  2  m 3 i +       A   E   ni  A   E  P 2 2 i 2! i 0 i!  m1 + ... +       A   E  P m 1 m 1  m  1! Vì  A   E  P  0 với i  m nên i m 1 i i  Ci j n j      A  E  i j i P i  0 j 0 i! m 1 i  m 1 i =     ni  A   E  P        A   E   ni  A  E i 1 i P i 0 i! 1! i 0 i!  m1 m 1 i + … +  m 1    A   E   ni i !  A   E  m 1 i P  m  1 ! i 0 m 1 i m 1 i        A   E   ni  A   E  i i i = P . i 0 i! i! i 0 Do đó, ta có m 1 i m 1 i m 1 i Sn  e A  P =        A  E  ni  A   E  P   i     A  E i i i i P i 0 i! i 0 i! i 0 i!
m 1 i = en X   A  ni  A  E P  X   A P . i i 0 i! m 1 i với X   A    i     A   E  ,    A . i i 0 i! Vậy m 1 i P xn  en A P x0  en X   A  ni  A  E P b  X   A P b i i 0 i!  en A P x0 + en A X   A P b  X   A P b  en A  P x0  X   A P b   X   A P b . * Nếu e z  1 , với số Bernoulli Bi , i  0,1, 2,..., ta có E 0n  z   n  B0 n n  n  1 1 1 1 1 En1  z     n  n 2  C21 B1n  C22 B0 n 2 2 2 2 2 2 1 1 1 1 1 1 En2  z   n  n 2  n3  C31B2 n  C32 B1n 2  C33 B0 n3 6 2 3 3 3 3 ……………………………….. i 1 Bi 1 j j Đến bước thứ i = m – 1 ta được: Eni  z    Ci j 1 n . j 1 i 1 Do đó, ta có: m 1  i 1 Bi 1 j j   i Sn  e A  P     Ci j 1 n   A   E  P i i 0  j 1 i  1  i !  B1 B  Hay Sn  e A  P = C11B0 nP   C21 n  C22 0 n 2   A   E  P  2 2  3 B0 3   2  1 B2 2 B1 2 n   A   E  P 2 +  C3 n  C3 n  C3  3 3 3  2!  Bm1 B B   m1  A   E  P m 1 +…+  Cm1 n  Cm2 m2 n 2  ...  Cmm 0 n m   m m m   m  1! m 1 ni 1  i  n  m  2 i 1 i = B0   A   E  P  B1  A   E    A   E  P i i i 0 i  1 i ! 1! i 0 i  1 i !  m1  A  E m 1 + … + Bm 1 nP .  m  1!
Từ  A   E  P  0 với i  m , ta có: i ni 1  i m 1  n  m 1 i 1 i Sn  e  P = B0   A   E  P  B1  A   E    A   E  P A i i i 0 i  1 i ! 1! i 0 i  1 i !  m1 m 1 ni 1  i  A  E   A   E  P m 1 i + … + Bm1  m  1! i 0 i  1 i ! m 1 i m 1 ni 1  i  Bi  A  E   A   E  P . i i = i 0 i! i 0 i  1 i ! m 1 ni 1  i = Y  A   A   E  P . i i 0 i  1 i ! m 1 i Y  A   Bi  A   E  ,    A . i với i 0 i! Vậy m 1 ni 1  i 1 n  m 1 i 1 i P xn    A   E  P x0  Y  A   A   E  P b  P x0 i 1 i i 0 i  1 i ! i 0 i  1 i ! m 1 ni 1  i   A   E    A   E  P x0  Y  A P b   P x0 . i i 0 i  1 i ! Tóm lại các kết quả trên, ta có: Định lí 3.3: 1) Nếu e  1 , khi đó P xn  en A  P x0  X   A P b   X   A P b m 1 i X   A    i    i !  A   E  ,    A . i với i 0 2) Nếu e  1 , khi đó m 1 ni 1  i P xn    A   E    A   E  P x0  Y  A P b   P x0 i i 0 i  1 i ! m 1 i Y  A   Bi  A   E  ,    A . i với i 0 i! 3.3.2. BIẾN ĐỔI ĐẾN TRƯỜNG HỢP LIÊN TỤC. Để tính toán đơn giản, đặt b f  b . Gọi xn  0  là nghiệm của xn  e A xn 1  b với điều kiện đầu x0  0 .
Áp dụng định lí 3.3 để nhận được hàm thành phần P z  s  n   P  s   P e sA xn  0  1) Giả sử e  1 . Ta có: P xn  0   en A X   A P b  X   A P b Suy ra P e sA xn  0   e sA  en A X   A P b  X   A P b   e sA X   A P b  e s  n  A X   A P b. Đặt X  X   A P b. Khi đó P z  t   P  s   P e sA xn  0  trở thành P z  s  n    P  s   e sA X   e s  n  A X , Hay P z  s  n   e s  n  A X  P  s   e sA X . Do đó c  t  : P z  t   etA X là tuần hoàn chu kì  , và P z  t   etA X  c  t  . * Bây giờ ta kiểm tra hàm ban đầu P  s  : e sA X  c  s  , s    ,0  thỏa mãn điều kiện: P  0   P     P e Ab của hàm ban đầu trong định lí 3.2. Từ P  0    X  c  0  và c  0   c    , điều kiện P  0   P     P e Ab tương đương với điều kiện X  e A X  e A P b (22) Hay e A X  X  P b Ta chứng minh (22) luôn đúng. Thật vậy, Áp dụng định lí 3.3, trường hợp x0  0 , ta có: P x1  e A X  X Mặt khác, từ phương trình xn  e A xn 1  b , ta có: x1  b . Do đó P b  e A X  X hay e A X  X  P b . Vậy (22) luôn đúng.  Cho nên ta có P z  t   1   P etAb   etA X   A P b  c  t  với c(t) là một hằng số tuần hoàn. Suy ra P x  t;0, w   e At Pw  P z  t   P g  t 