Luận văn phân tích độ nhạy liên hợp đối với phương trình vi phân đại số (tt)

BỘ GIÁO DỤC VÀ ĐÀO TẠO UBND TỈNH THANH HÓA TRƯỜNG ĐẠI HỌC HỒNG ĐỨC LÊ THỊ OANH PHÂN TÍCH ĐỘ NHẠY LIÊN HỢP ĐỐI VỚI PHƯƠNG TRÌNH VI PHÂN ĐẠI SỐ Chun ngành: Tốn Giải tích Mã số: 60 46 01 02 TÓM TẮT LUẬN VĂN THẠC SĨ THANH HĨA – 2016 Luận văn hồn thành Trường Đại Học Hồng Đức Người hướng dẫn: TS.Hoàng Nam Phản biện 1: ………………………… Phản biện 2: …………………………… Luận văn bảo vệ Hội đồng chấm luận văn thạc sĩ khoa học tại: Trường Đại Học Hồng Đức Vào hồi: ….giờ… ngày …tháng… năm2016 Có thể tìm hiểu luận văn -Thư viện trường Đại học Hồng Đức, -Bộ mơn: Giải Tích,Trường Đại học Hồng Đức MỞ ĐẦU Lý chọn đề tài Trong thực tiễn nay, kể khoa học ứng dụng có nhiều vấn đề, nhiều tốn, chẳng hạn mô tả động lực, mô tả hệ thống mạng điện, lý thuyết điều khiển, đòi hỏi phải quan tâm giải hệ phương trình vi phân đại số Từ cuối năm 70 đầu năm 90 kỉ XX có nhiều nhà tốn học giới nghiên cứu phương trình vi phân đại số, số nhà tốn học thuộc Đại học Humbodt Berlin, nhóm nhà toán học Nga, Mỹ, Ba Lan số nước khác Ở nước ta, vào năm 90 kỉ XX có số nhà tốn học thuộc Đại học sư phạm Hà Nội, Đại học khoa học Tự nhiên - Đại học Quốc gia Hà Nội nghiên cứu phương trình vi phân đại số Trong thời gian qua có nhiều kết thu phương trình vi phân đại số, chẳng hạn kết nghiệm, tính ổn định, ổn định tiệm cận, ổn định mũ phương trình vi phân đại số, tính ổn định hệ có nhiễu nhỏ, tính nhị phân, phương trình liên hợp phương trình vi phân đại số,… góp phần bước hoàn thiện phát triển lý thuyết phương trình vi phân đại số đẩy mạnh việc ứng dụng chúng thực tiễn Các mơ hình tốn học sử dụng để khảo sát, điều tra tượng vật lý ngày trở nên thực tế Các đặc tính mơ hình họ thường sử dụng thông số mà giá trị khơng biết xác, cách gọi để phân tích độ nhạy tham số Lĩnh vực áp dụng, bao gồm tối ưu hóa, tham số ước lượng, điều khiển tối ưu, mơ hình đơn giản hóa, độ nhạy q trình, phân tích khơng chắn thiết kế thí nghiệm cho loạt vấn đề khoa học kỹ thuật,… Nghiên cứu gần phương pháp phần mềm để phân tích độ nhạy phương trình vi phân đại số (Cao et al., 2003; Feehery et al., 1997; Li and Petzold, 2000; Li and Petzold, 1999; Li at al., 2000 and Maly and Petzold, 1997) chứng tỏ dãy giá trị nhạy cảm tính tốn đáng tin cậy hiệu thông qua vi phân kết hợp với kỹ thuật nghiệm phương trình vi phân đại số Tuy nhiên, có số vấn đề khó khăn ta phân tích độ nhạy hệ phương trình vi phân đại số với số lượng lớn tham số so với số lượng biến Theo phương pháp này, ta tính độ nhạy tham số cách sử dụng hệ liên hợp tuyến tính, khơng cần thiết phải mở rộng hệ ban đầu, điều nâng cao hiệu tính tốn vấn đề xem xét Phương pháp sử dụng loạt thí nghiệm, nơi kết so sánh với ước lượng độ nhạy phương pháp hữu hạn khác nhau, sử dụng rộng rãi cho đơn giản việc thực Bởi vậy, Luận văn nghiên cứu đề tài: “Phân tích độ nhạy liên hợp phương trình vi phân đại số” Mục đích nghiên cứu - Mục đích nghiên cứu luận văn nhằm giới thiệu phương pháp dựa hệ liên hợp để phân tích độ nhạy tham số phương trình vi phân đại số Phương pháp sử dụng loạt thí nghiệm kết so sánh với việc ước lượng độ nhạy phương pháp khác hữu hạn, sử dụng rộng rãi tính đơn giản việc thực chúng xét tính ổn định hệ liên hợp Nhiệm vụ nghiên cứu - Trình bày phương pháp độ nhạy liên hợp hệ phương trình vi phân đại số phụ thuộc vào tham số Một phương pháp dựa hệ liên hợp để phân tích độ nhạy tham số phương trình vi phân đại số - Nghiên cứu đề xuất mở rộng phương pháp dựa việc sử dụng hàm so sánh cho phép chuyển hệ liên hợp ban đầu hệ đại số làm cho phù hợp với nhiều phương pháp số hiệu - Trình bày tính ổn định, tính ổn định hệ liên hợp (đối với phương trình vi phân đại số bán tường minh) hệ liên hợp bổ sung (đối với phương trình vi phân đại số hoàn toàn ẩn) Đối tượng phạm vi nghiên cứu - Phương trình vi phân đại số phương trình liên hợp; - Độ nhạy liên hợp phương trình vi phân đại số; - Tính ổn định, ổn định số phương trình vi phân Phương pháp nghiên cứu - Tìm kiếm ,tổng hợp tài liệu từ giáo trình ,sách phường vi phân đại số độ nhạy liên hợp phương trình vi phân đại số Sau phân tích ,tổng hợp để trình bày rõ ràng ,hợp logic vần đề Cấu trúc luận văn Luận văn cấu trúc gồm phần Mở đầu, hai chương nội dung sau: Chương 1: Phương trình vi phân đại số phương trình liên hợp Chương trình bày số khái niệm, kết phương trình vi phân đại số phương trình vi phân đại số liên hợp tính ổn định nghiệm phương trình vi phân Chương 2: Độ nhạy liên hợp phương trình vi phân đại số Chương giới thiệu phương pháp dựa hệ liên hợp để phân tích độ nhạy tham số phương trình vi phân đại số; mở rộng phương pháp dựa việc sử dụng hàm so sánh, cho phép chuyển hệ liên hợp ban đầu hệ đại số làm cho phù hợp với nhiều phương pháp số hiệu Cuối phần Kết luận phần Danh mục tài liệu tham khảo 3 CHƯƠNG PHƯƠNG TRÌNH VI PHÂN ĐẠI SỐ VÀ PHƯƠNG TRÌNH LIÊN HỢP 1.1 Một số khái niệm phương trình vi phân đại số Định nghĩa 1.3 Số tự nhiên k gọi số ma trận A số tự nhiên bé thỏa mãn kerAk  kerAk 1 Kí hiệu số ma trận A ind  A  , ind  A   {k : kerAk  kerAk 1} 1 Định nghĩa 1.5 Nếu cặp  A, B quy det (cA  B) 0 với c  C ind {  cA  B  A} 1 gọi số cặp ma trận  A, B Như ind {A, B} ind {  cA  B  A} với c  C Định nghĩa 1.6 Phương trình A( Px) ' (B  AP ') x  q , (1.2) A, B : I  L(C m , C m ), f : I  C m ma trận hàm thỏa mãn giả thiết sau  T1 dim  imA  t   r  m, t  I ;  T 2 Cặp ma trận  A  t  , B  t   quy số với t  I ;  T 3 Tồn phép chiếu Q  C1 ( I , L(C m , C m ) lên kerA , gọi phương trình vi phân đại số tuyến tính số chuyển Định lý 1.1 [10] Giả sử (1.1) phương trình vi phân đại số quy số  Khi x(t) nghiệm  thỏa mãn điều kiện đầu x(0)  x  ker A(0)  N (0) ; (1.4) x(t ) PS (t )u (t )  Q (t )G  (t )  B (t )u (t )  f (t ) , t    , u(t)là nghiệm u '(t ) P '(t )u (t )  P(t )( I  P '(t ))G  (t ) B(t )u (t )  f (t )) toán giá trị đầu  Nếu ta sử dụng phép u (0) P(0) x chiếu tắc PS I  QG  B lên S(t) dọc N(t) cơng thức (1.5) (1.6) viết lại sau  x(t ) PS (t )u (t )  Q(t )G  f (t )  1 1 u '(t ) ( P '(t ) PS (t )  P (t )G (t ) B(t )u (t )  P(t )( I  P '(t ))G (t ) f (t )  x(t ) PS (t )u (t ) f(t)=0  1 u '(t ) ( P '(t ) PS (t )  P (t )G (t ) B (t )u (t ) 1.2 Phương trình liên hợp phương trình vi phân đại số  1.7   1.8 Định nghĩa 1.8 Phương trình (1.9) ( A )  B s với A, B : I  L(C m , C m ), s : I  C m thỏa mãn giả thiết  T1 dim  imA  t   r  m, t  I  T 2 Cặp ma trận  A  t  , B  t   quy số với t  I  T 3 Tồn phép chiếu Q  C1 ( I , L(C m , C m ) lên kerA , gọi phương trình liên hợp phương trình vi phân đại số (1.2) số   m Định lý 1.2[3] Với giả thiết  T1  ,  T2  ,  T3  , với   C m s  C I , C bất kì, tốn giá trị ban đầu  L s   A t   t    0          (1.10)  m  có nghiệm   C A I , C cho A  t0    t0    0   m Định lý 1.3[3] Với giả thiết  T1  ,  T2  ,  T3  , với   C m s  C I , C nghiệm toán giá trị ban đầu (1.10)  A    B   p Xét phương trình   (1.11)   m m A, B : I  L C , C ma trận hàm liên tục thỏa mãn giả thiết  T1    T3  Định nghĩa 1.9 Cho phương trình A  t  x t   B  t  x  t  q, t  I , (1.13)     , detA t  0, rankA  t  r, A  t   C1 ( I , L(C m , C m )), B  t  C I , L C m t  I , x C gọi phương trình vi phân đại số số (index- tractable) I (i) dimN1  t   const , (ii) N1  t   S  t   C m , t  I , Q  t  phép chiếu lên N  t  ker A  t  , t  I , P  t  I  Q  t  ; A1  t   A  t   B0  t  Q(t ), B0  t  B  t   A  t  P(t ), N  t  ker A  t  ,     S  t   z  C m : B  t  z  imA  t  , N1  t  ker A1  t  , S1  t   z  C m : B  t  P  t  z  imA1  t  Xét phương trình vi phân dạng Hessenberg, nghĩa hệ dạng  x1 ' B11 x1  B12 x2 q1  B12 x1 =q  đó, x ( x1T , xT2 )T , x1  R m1 , x2  R m2 , m m1  m2 Định lý 1.5[14] Giả sử (1.13) phương trình vi phân đại số giải có số với Q hàm khả vi liên tục Khi đó, kết sau (i) Bài toán giá trị ban đầu (1.8), (1.2) có nghiệm   C1N ( I , R m ) : x  C ( I , R m ) : Px  C ( I , R m ) với q  C ( I , R m ), Q1G2 1q  C1 ( I , R m ) (ii) Nếu x(.) nghiệm phương trình x(t )  M (t ) :im (t )  S (t ), t  I (iii) Với x*  M (t* ) qua có nghiệm phương trình thời điểm t*  I Không gian nghiệm M(t) không gian riêng S(t) dim M (t ) m  dim N  dim(( N  S (t )) Bổ đề 1.6[14] Giả sử phương trình (1.13) chuyển với số 2, Q1 Q1 A2 B1 khả vi liên tục Khi phương trình (1.13) tương đương với hệ  1  z   PP1 A2 B1   PP1  z 0,  y 0    v    QQ1    QP1 A2 B1  z ,      (1.15) z PP1 x, y PQ1 x v Qx Hơn nữa, z  R  PP1  t  với số t   t0 ,   , 1   nghiệm z toán giá trị đầu z    PP1 A2 B1   PP1    z 0 , z  t  z thoã mãn   z  t   im  PP1  t   , t   t0 ,   1.3 Sự ổn định (Lyapunov) hệ phương trình vi phân đại số Xét hệ phương trình vi phân đại số tuyến tính sau A(t ) x t   B  t  x  t  0 (1.16)   n n n I  R , A, B  L , det A 0, q  C I , R Giả sử phương trình (1.16) có số KerA  t  trơn Định nghĩa 1.10 Nghiệm tầm thường x  t  0 hệ (1.16) gọi ổn định (theo nghĩa Lyapunov) với số   cho trước với t  I tồn   (t0 ,  )  cho x0  R n thoả mãn P  t0  x0   x  t ; t0 , x0    với t t0 Định nghĩa 1.11 Nghiệm tầm thường x  t  0 hệ (1.16) gọi ổn định tiệm cận ổn định tồn số  (t0 )  cho P  t0  x0    t0  x  t ; t0 , x0   t   Định nghĩa 1.12 Nghiệm tầm thường x  t  0 hệ (1.16) gọi ổn định tiệm cận mũ tồn số dương  với số   cho trước tồn số   (t0 ,  )  cho   t t x  t ; t0 , x0    e   với t t0 x0  R n thoả mãn P  t0  x0   Định nghĩa 1.13 Nghiệm tầm thường x  t  0 phương trình vi phân đại số có số (1.16) gọi ổn định (theo nghĩa Lyapunov) với số   cho trước với t  I tồn   (t0 ,  )  cho x0  R n thoả mãn P  t0  P1 (t0 ) x0   x  t ; t0 , x0    với t t0 Định nghĩa 1.14 Nghiệm tầm thường x  t  0 phương trình vi phân đại số có số (1.16) gọi ổn định tiệm cận ổn định tồn số  (t0 )  cho P  t0  P1 (t0 ) x0   x  t ; t0 , x0   t   Định nghĩa 1.15 Nghiệm tầm thường x  t  0 phương trình vi phân đại số có số (1.16) gọi ổn định tiệm cận mũ tồn số dương  với số   cho trước tồn số   (t0 ,  )  cho x0  R n thoả mãn P  t0  P1 (t0 ) x0     t t x  t ; t0 , x0    e   với t t0 Định lý 1.7 Giả sử Qs QA1 B0 bị chặn  t0 ,   Khi đó, nghiệm tầm thường (1.11) ổn định tiệm cận mũ nghiệm tầm thường phương trình   u   PA1 B0  P u 0, t   t0 ,   ổn định tiệm cận mũ imP  t  Định lý 1.8 [14] Giả sử phương trình (1.13) chuyển với số 2, khả vi liên tục đồng thời 1  QQ1   QP1 A2 B1 Q1 Q1 A2 1B1 bị chặn Khi đó, nghiệm tầm thường phương trình( 1.13) ổn định tiệm cận mũ nghiệm tầm thường nghiệm phương trình   z  PP1 A2 1B1   PP1   z 0, t   t0 ,   ổn định tiệm cận mũ im  PP1  t   CHƯƠNG ĐỘ NHẠY LIÊN HỢP CỦA PHƯƠNG TRÌNH VI PHÂN ĐẠI SỐ Trong chương này, ta giới thiệu phương pháp dựa hệ liên hợp để phân tích độ nhạy tham số phương trình vi phân đại số Ngoài ta mở rộng phương pháp dựa việc sử dụng hàm so sánh, cho phép chuyển hệ liên hợp ban đầu hệ đại số làm cho phù hợp với nhiều phương pháp số hiệu Phương pháp độ nhạy liên hợp trình bày phương trình vi phân đại số phụ thuộc tham số Hệ liên hợp suy với điều kiện ban đầu quán phương trình vi phân Hessenberg số 2.1 Đạo hàm hệ liên hợp độ nhạy Xét phương trình cân có dạng tổng quát      F  x, x, t , p  0,     x   x  p  ,  (2.1) tạo nên tổ hợp phương trình vi phân phương trình đại số, x  R nx vec tơ biến trạng thái p  R n p vec tơ tham số Vấn đề tính tốn độ nhạy tham số hệ phương trình vi phân đại số dạng (2.1) có dạng sau: Đối với phụ thuộc tham số hệ phương trình vi phân đại số dạng (2.1), tìm dx thời gian T, với j 1, 2, , n p dp j Các nghiệm chúng địi hỏi nghiệm tương thích hệ phương trình vi phân đại số ban đầu với với n p - hệ nhạy cảm, thu cách lấy đạo hàm phương trình vi phân đại số ban đầu tương ứng theo tham số Ta qua tâm đến việc tính tốn độ nhạy dG hàm mục tiêu dp nghiệm nhiễu tham số T  g ( x, t , p)dt , G  x, p   loại trù độ nhạy dg hàm g ( x, t , p) xác định thời điểm T Bây dp ta xét hệ phương trình vi phân đại số tuyến tính cho cơng thức sau Ax  Bx 0, A, B ma trận hàm phải thỏa mãn đủ trơn, phương trình vi phân đại số liên hợp phương trình ( A )  B 0, 2.1.1 Độ nhạy G  x, p  Trước hết, để giải toán độ nhạy hàm G  x, p  , cho công thức (2.2) Để giới thiệu hệ số Lagrange  , ta hình thành hàm mục tiêu bổ sung I ( x, p) sau I  x, p  G  x, p   T   F ( x, x, t , p)dt ,  đó,  hệ số Lagrange F ( x, x , p, t ) 0 xác định (2.1) Lấy đạo hàm hai vế theo p , độ nhạy G  x, p  p dG  dp T  g Ta giả sử P    * Fp dt  T    g x    * Fx  ( * Fx ) x p dt   * Fx x p  T ( * Fx )   * Fx  g x , (2.7) (2.8) phương trình vi phân đại số có số có số 1, để đơn giản ta lấy  (t ) Fx 0 t T dG  dp T  g Phương trình (2.7) trở thành P    * Fp dt  ( * Fx x p ) T0 (2.10) Để tìm điều kiện ban đầu (tại t T ) hệ liên hợp, ta phải thực việc xem xét cấu trúc hệ phương trình vi phân đại số Đối với trưịng hợp hệ phương trình vi phân đại số có số số 1, ta chọn   Fx t T 0 , dẫn tới phương trình độ nhạy dG  dp T  g P (2.11) dG dp      Fp dt    Fx x p  t 0 Cách chọn khơng đủ cho trường hợp hệ phương trình vi phân đại số Hessenberg có số Xét hệ d d a  x  f ( x , x , p ),  d 0  f ( x , p ), đây, x d x a thành phần nghiệm vi phân thành phần nghiệm đại số tương ứng, A f f f , , ma trận tích CB khả nghịch Ở đây, hệ liên hợp cho B  C  x d x a x d công thức  d*   d* A   a* C  g d , x  d  B  g xa ,   d   biến vi phân biến đại số tương ứng ta giả sử  d (T )  *C t T ,   chưa xác định Khi đó, phương trình độ nhạy (2.10) trở thành dG  dp T  ( g P   d f p1    f p2 )dt  ( d x dp ) t 0  g xa (CB )  f p2 Ta sử dụng ký hiệu hệ liên hợp bổ sung tương ứng với (2.6)  (t )   * Fx  g x  *  (t )   Fx 0  (T ) 0  (219) đây, ta giới thiệu hàm  (t ) Từ (2.19) trở đi, dễ thấy phương trình vi phân đại số liên hợp T phải giải ngược lại thời gian Để tính số hạng tích phân ( g p   * Fp )dt phương trình (2.8), ta sử dụng biến cầu phương β, mở rộng lần với hệ liên hợp  (t )   * F  g xx  *  (t )   Fx 0  *    ( g   Fp ) 0   (T ) 0   (T ) 0  T * Trong trường hợp này, với t = 0,  (0)  ( g p   Fp )dt thỏa mãn phương trình độ nhạy dG dp dG  (0)   * Fx x p dp t 0 2.1.2 Độ nhạy g  x, T , p  Bây ta xem xét việc tính toán dg dg d dG d dG   Từ (2.8), ta có dp dp dp dT dT dp 10 dg ( g p    Fp )(T )  dp T ký hiệu T   T Fp dt   T Fx x p ) t 0  (  d (  Fx x p ) dT  T Đối với phương trình vi phân đại số có số có số d (  Fx x p ) t T dT 0 Đối với phương trình vi phân đại số Hessenberg có số dạng (2.11) d (  Fx x p ) t T dT Phương trình liên hợp tương ứng  d ( g x (CB )  f p2 ) dt  t T  (T Fx )  (T FX ) 0 Đối với phương trình vi phân đại số có số số (như trình bày trên, trường hợp số lại khác), để tìm điều kiện biên cho phương trình này, ta viết   (t , T ) , phụ thuộc vào t T Khi  * (T , T ) Fx  (T Fx ) t T t T 0 Ta có điều kiện biên   dFx        (T , T )   Fx   dt   (2.22) t T Đối với trưịng hợp phương trình vi phân đại số có số 1, từ (2.8) (2.22) dẫn đến  (T Fx ) t T  g x    Fx  t T Đối với trường hợp phương trình vi phân thường qui ẩn, Fx khả nghịch, đó, ta có    T , T  0 điều dẫn đến   T     T  T Ta thảo luận điều kiện biên thích hợp phương phương trình vi phân đại số Hessenberg số phần sau Ta có dg ( g p   (t , T ) Fp ) dp t T  T  T (t , T ) Fp dt  T (t , T ) Fx x p biến liên hợp  (t , T ) phụ thuộc vào t T T (t , T ) biểu thị t 0 , (*) d  (t , T ) dT Phương trình liên hợp tương ứng  (t )  T (t , T ) Fx 0   (t )  T (t , T ) Fx 0   (T ) ( g x   Fx ) t T giải tương tự hệ (2.10) Một lần nữa, để tính số hạng tích phân phương trình (*), ta sử dụng biến  ta xác định hệ liên hợp sau 11  (t )  T (t , T ) Fx 0   (t )  T (t , T ) Fx 0   (T ) ( g x   Fx ) t T     (t , T ) F 0 T p  Cuối cùng, biểu thị phương trình độ nhạy dg ( g p   (t , T ) Fp ) dp t T dg cho dp   (0)  T (t , T ) Fx x p t 0 2.2 Một số ví dụ 2.2.1 Phương trình vi phân thường dạng chuẩn tắc Cho phương trình vi phân thường với điều kiện đầu  x  f ( x, t )   x(0) x0 ( p) hàm g  x  t T , hệ liên hợp tương ứng T  f xT   T (T ) g x độ nhạy cho dx dg T (0) x0 p , x0 p  Phương trình tương tự với hệ liên dp dp hợp nói chung cách xác định phương trình vi phân thường 2.2.2 Phương trình vi phân thường dạng ẩn Cho hệ  F  x , x, t  0,   x(0) xo ( p), với A  F F ma trận không suy biến, B  hàm g  x  t T , hệ liên hợp tương  x x ứng    ( A T )  B T 0  AT (T ) g x  độ nhạy cho dg T (0) A(0) x0 p dp 2.2.3 Phương trình vi phân đại số có số bán tường minh Cho hệ 12  x  f ( x d , x a , t ),  d a   f (x , x ) ,  d d  x (0) x0 ( p), với A  f f f f d , , , ma trận không suy biến hàm g x B  C  D  d a d a x x x x   t= T, hệ liên hợp tương ứng Td  ATd  C Ta ,   d  a  B T  D T ,  Td g  xd  độ nhạy cho  dx p dg Td (0) Nếu hàm g phụ thuộc hai biến x d x a , dp dp  g g x d  T  , x  T  , phương trình liên hợp tương tự trên, điều kiện biên bây   cho công thức Td (T ) g xd  C * D* 1 g xa 2.2.4 Phương trình vi phân đại số Hessenberg có số Cho hệ phương trình sau  x d  f ( x d , x a , t ),   f ( x d ),   d d  x (0) x0 ( p), với A  f f f , , CB khả nghịch, giả sử hàm g  x, T , p  phụ thuộc B  C  x d x a x d vào thành phần vi phân x d x , hệ liên hợp tương ứng  Td  ATd  C Ta ,   BTd ,   d    1  * T (T ) ( I  C ( B C ) B ) g xa  độ nhạy cho dx0p dg d T (0) Nếu hàm g  x, T , p  phụ thuộc vào x d x a , trước hết ta dp dp giải điều kiện biên từ  d  d  a    A   C   g xd ,  B d  g *d ,  x  d    1    (T )  C ( B C ) g xa  13  d   ATd  C Ta , Hệ liên hợp  T điều kiện biên cho Td  d 0 B T ,  dC       Td  T , T  P*  g xd  A d (T , T )  ( B C ) g xa  dt   2.3 Tính ổn định Trong phần này, ta nghiên cứu tính ổn định hệ liên hợp Giả sử rằng, hệ ban đầu ổn định, liệu hệ liên hợp có ổn định khơng ? Nếu ta xét phương trình liên hợp (2.8), liệu khơng ổn định Xét ví dụ et x  et x 0 Hệ tương đương với x  (2.24) x 0 Rõ ràng, hệ ổn định Nhưng hệ liên hợp (2.8) (2.24) et   et   et  0 , tương đương với    0 (2.26) Chú ý rằng, ta giải hệ liên hợp ngược Do đó, hệ liên hợp (2.26) khơng ổn định Ký hiệu  Fx* , ta có dạng hệ liên hợp bổ sung (2.8)    Fx  g x    Fx   0  Nếu giải hệ liên hợp bổ sung (2.27) thay cho (2.26),  thỏa mãn (2.27)    0 ổn định ngược Ta biểu diễn rằng, trường hợp tổng quát hệ liên hợp bổ sung (2.27)  ổn định hệ ban đầu ổn định Vì [6] ta giải phương trình (2.27) thay (2.8) 2.3.1 Phương trình vi phân thường dạng tường minh Xét hệ phương trình vi phân thường tuyến tính x  A  t  x Hệ phương trình liên hợp tương ứng   A* Do, hệ liên hợp giải theo chiều ngược, ta thực phép đổi biến đổi sau  T  t Hệ liên hợp giải ngược chuyển thành   A* Khi đó, ta có kết sau 14 Định lý 2.1.[7] Nếu hệ phương trình vi phân x  A  t  x ổn định, hệ liên hợp ổn định 2.3.2 Phương trình vi phân đại số có số bán tường minh Xét hệ phương trình vi phân đại số có số dạng bán tường minh d d a  x  A(t ) x  B(t ) x ,  d a  C (t ) x  D(t ) x ,  d  A d  C   , đó, D(t) khả nghịch Hệ liên hợp tương ứng   d    B   D  Hệ ban đầu ổn định phương trình vi phân thường sở (Essential Underlying ODE - EUODE) x d  Ax d  B ( D)  Cx d ổn định Phương trình vi phân thường sở tương ứng   với hệ liên hợp cho  d  A d  C  D 1   B d ,  d  ( A  B D 1 C )  d 2.3.3 Phương trình vi phân đại số Hessenberg có số Xét hệ phương trình vi phân đại số Hessenberg có số tuyến tính  x d  A(t ) x d  B(t ) x a  q  t  ,  C (t ) x d  r (t ),  đây, ma trận C  t  B  t  không suy biến Nếu B đủ trơn, tồn ma trận hàm trơn, bị chặn R  t   R  nd  na  nd , độc lập tuyến tính, chuẩn hóa theo hàng thành sở không gian không B* Ở nd số chiều biến vi phân x d na số chiều biến đại số x a  R Do R  t  B  t  0   khả nghịch Đặt biến u Rx d , t T , x d cho C 1  R  u  x     Su  Fr , C  r d S  t   R nd  nd  n  thỏa mãn điều kiện RS I , CS 0 Phương trình vi phân thường sở  )u  ( RAF  RF  )r (t )  Rq  t  tương ứng u ( RSA  RS (2.28) Hệ phương trình ban đầu ổn định biến vi phân phương trình vi phân thường sở (2.28) ổn định Bây giờ, ta xét hệ liên hợp  d  A (t ) d  C  (t )   qˆ  t  ,    B (t ) d  r (t ) Khi đó, phương trình vi phân thường sở phương trình vi phân đại số liên hợp  v S  R v  S  C  B C   1  r  S  A Rv  S  AC  BC   1 r 15 RS=I dẫn đến S  R  S  R  , thành phần phương trình phương trình   RAS * v vi phân thường sở cho v  S * R *  S * A R v  RS     2.3.4 Phương trình vi phân đại số hồn tồn ẩn Bổ đề 2.2 [7] Cho hệ phưong trình vi phân đại số tuyến tính phụ thuộc A  t  x  B  t  x  f  t  ma trận khả vi, phụ thuộc thời gian, khơng suy biến P(t) nhân phương trình phương trình vi phân đại số Q(t) đổi biến, hệ liên hợp phương trình vi phân đại số biến đổi hệ biến đổi phưong trình vi phân đại số liên hợp Định lý 2.3.[7] Đối với hệ phương trình vi phân đại số tổng quát có số số 1, hệ phương trình vi phân đại số ban đầu ổn định hệ phương trình vi phân đại số liên hợp bổ sung (2.27)  ổn định 2.4 Tính ổn định số Xét phương trình vi phân đại số tuyến tính tổng quát A  t  x  t   B  t  x  t  0 (2.29) Hệ liên hợp tương ứng mà chúng giải theo hướng ngược ( A )  B 0  xn 1  xn  Rời rạc hóa hệ liên hợp (2.29), ta thu An 1    Bn 1 xn 1 0 Suy h   n  An  hBn  1 Do đó, ta có 0* N AN  AN   hBN   Giả sử Cn  AN   hBN   1 An*1n 1 1 AN   A1  hB1  AN   A1  hB1  1 1 A1  A0  hB0  1 (2.33) A1 Tính ổn định số hệ ban đầu kéo theo C N bị chặn Từ (2.33), hệ liên hợp ta có 1 0 N AN C N  A0  hB0  Do đó, tất phương trình vi phân thường tuyến tính phương trình vi phân đại số, tính ổn định tiệm cận số phương pháp Euler theo hướng ngược lại hệ ban đầu kéo theo tính ổn định tiện cận số hệ liên hợp 16 KẾT LUẬN Trên sở báo [6,7] sau thời gian nghiên cứu, tìm tịi với giúp đỡ tận tình TS Hồng Nam, tơi hồn thành Luận văn thạc sỹ với tiêu đề: “Phân tích độ nhạy liên hợp phương trình vi phân đại số” theo yêu cầu nội dung kế hoạch đề Luận văn thu số kết sau: Trình bày phương pháp độ nhạy liên hợp hệ phương trình vi phân đại số phụ thuộc vào tham số Phương pháp sử dụng loạt thí nghiệm kết so sánh với việc ước lượng độ nhạy phương pháp hữu hạn khác, sử dụng rộng rãi Nghiên cứu đề xuất mở rộng phương pháp dựa việc sử dụng hàm so sánh cho phép chuyển hệ liên hợp ban đầu hệ đại số làm cho phù hợp với nhiều phương pháp số hiệu Trình bày tính ổn định, tính ổn định hệ liên hợp (đối với phương trình vi phân đại số bán tường minh) hệ liên hợp bổ sung (đối với phương trình vi phân đại số hồn tồn ẩn) Tuy nhiên, thời gian có hạn nên luận văn khơng tránh khỏi khiếm khuyết, mong thầy cô bạn góp ý để luận văn hồn thiện 17