Luận văn Phương trình sóng phi tuyến với một đầu biên chứa số hạng memory

Trong luận văn này, chúng tôi xem xét phương trình sóng phi tuyến với điều kiện biên và điều kiện đầu dưới đây u t u K u u u u f x t x t T tt xx t t             p2 q2  , , 0 1 , 0  (1.1)           0 0, 0, , 1, 0 t x  t u t k t s u s ds u t     (1.2) u x u x u x u x     ,0 , ,0   0 1 t     (1.3) trong đó p q K , 2, , 0    là các hằng số cho trước và , , , , f k u u 0 1 là các hàm cho trước thỏa các điều kiện mà ta sẽ chỉ ra sau. Bài toán (1.1) – (1.3) có nhiều ý nghĩa trong Cơ học mà nhiều nhà Toán học quan tâm nghiên cứu trong thời gian gần đây [3] – [10]. Về mặt toán học, ta có thể nói rằng bài toán (1.1) – (1.3) là sự kết hợp những đặc điểm quan trọng của hai bài báo đã được công bố trước đây [3] và [4]

pdf44 trang | Chia sẻ: duongneo | Lượt xem: 951 | Lượt tải: 0download
Bạn đang xem trước 20 trang tài liệu Luận văn Phương trình sóng phi tuyến với một đầu biên chứa số hạng memory, để xem tài liệu hoàn chỉnh bạn click vào nút DOWNLOAD ở trên
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM THÀNH PHỐ HỒ CHÍ MINH ___________________ Hoàng Quốc Công PHƯƠNG TRÌNH SÓNG PHI TUYẾN VỚI MỘT ĐẦU BIÊN CHỨA SỐ HẠNG MEMORY 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: TS. NGUYỄN THÀNH LONG Thành phố Hồ Chí Minh, năm 2010 LỜI CẢM ƠN Lúc đầu khi nhận đề tài này, với vốn kiến thức còn hạn hẹp tôi đã gặp phải rất nhiều khó khăn. Tuy nhiên, với sự hướng dẫn tận tình và mang tính khoa học của Thầy hướng dẫn, TS. Nguyễn Thành Long, tôi đã dần khắc phục được các khó khăn trên để hoàn thành đề tài này. Trước tiên, tôi xin chân thành cảm ơn Thầy Nguyễn Thành Long, người đã tận tình dìu dắt tôi vượt qua nhiều trở ngại trong suốt thời gian thực hiện đề tài. Tôi xin chân thành cảm ơn các Thầy Cô trong và ngoài trường Đại học Sư Phạm Tp. Hồ Chí Minh, những người đã tạo mọi điều kiện thuận lợi cũng như tận tình truyền đạt kiến thức và kinh nghiệm quý báu của mình cho tôi trong suốt thời gian tôi học tập và làm việc. Tôi xin trân trọng cảm ơn các Thầy Cô trong Hội đồng đã dành thời gian để đọc và cho những nhận xét rất có giá trị khoa học đối với luận văn của tôi. Cuối cùng, tôi xin gửi lòng biết ơn của mình đến gia đình và các bạn tôi, những người đã tạo điều kiện tốt nhất cho tôi hoàn thành luận văn này và cho tôi những lời khuyên, lời động viên vô cùng hữu ích. Hoàng Quốc Công Chương 1 PHẦN TỔNG QUAN Trong luận văn này, chúng tôi xem xét phương trình sóng phi tuyến với điều kiện biên và điều kiện đầu dưới đây    22 , , 0 1 , 0qptt xx t tu t u K u u u u f x t x t T          (1.1)           0 0, 0, , 1, 0 t xt u t k t s u s ds u t    (1.2)        0 1,0 , ,0tu x u x u x u x  (1.3) trong đó , 2, , 0p q K   là các hằng số cho trước và 0 1, , , ,f k u u là các hàm cho trước thỏa các điều kiện mà ta sẽ chỉ ra sau. Bài toán (1.1) – (1.3) có nhiều ý nghĩa trong Cơ học mà nhiều nhà Toán học quan tâm nghiên cứu trong thời gian gần đây [3] – [10]. Về mặt toán học, ta có thể nói rằng bài toán (1.1) – (1.3) là sự kết hợp những đặc điểm quan trọng của hai bài báo đã được công bố trước đây [3] và [4]. Trong [3], Nguyễn Thành Long, Lê Văn Út, và Nguyễn Thị Thảo Trúc đã nghiên cứu sự tồn tại và duy nhất nghiệm cũng như tính chính quy, tính ổn định và khai triển tiệm cận nghiệm của bài toán  ( ) , , 0 1, 0 ,tt xx tu t u Ku u f x t x t T         (1.4) 1 1 0 ( ) (1, ) ( ) (1, ) ( ) (1, ) ( ) ( ) (1, ) , (1, ) 0, t x tt u t K t u t t u t g t k t s u s ds u t          (1.5) 0 1( ,0) ( ) , ( ,0) ( ),tu x u x u x u x  (1.6) trong đó ,K  là các hằng số, và 1 1 0 1, , , , , , ,f K g k u u  là các hàm cho trước. Như vậy, số hạng 22 qp t tK u u u u   trong (1.1) là sự tổng quát hóa từ số hạng tKu u trong (1.4). Các điều kiện biên (1.2) – (1.3) cũng chính là (1.5) – (1.6) sau khi đã hoán đổi 2 đầu biên 0x  và 1x  , đồng thời làm triệt tiêu các hàm K1, λ1 và g. Sự đặc biệt hóa này tưởng chừng sẽ mang lại thuận lợi cho chúng ta khi nghiên cứu (1.1) – (1.3), nhưng thật ra nó lại khiến chúng ta gặp đôi chút khó khăn hơn trong các ước lượng khi mà điều kiện cực tiểu cho hàm λ1 lúc này không còn là một giá trị dương nữa. Trong [5], Nguyễn Thành Long, Alain Phạm Ngọc Định, và Trần Ngọc Diễm đã khảo sát sự tồn tại và duy nhất nghiệm, tính chính quy, và khai triển tiệm cận của nghiệm của bài toán sau đây  , , 0 1, 0 ,tt xx t tu u K u u u u f x t x t T         (1.7)    0, ,xu t P t (1.8)      1 11, 1, 1, 0,x tu t K u t u t   (1.9)        0 1,0 , ,0 ,tu x u x u x u x  (1.10) trong đó 1 1, , , , ,K K    là các hằng số cho trước; 0 1, ,f u u là các hàm cho trước; hàm  ,u x t cần tìm và giá trị biên chưa biết  P t thỏa mãn bài toán Cauchy cho phương trình vi phân thường dưới đây           2 0 1 0, , 0 0 , 0 ttP t P t hu t t T P P P P        (1.11) với 0 10, 0, ,h P P   là các hằng số cho trước. Từ (1.11), nếu ta giải phương trình vi phân tuyến tính cấp 2 hệ số hằng, ta sẽ có           0 0, 0, t P t g t hu t k t s u s ds    trong đó                0 0 1 110 cos 0 sin , sin .g t P hu t P hu t k t h t        Ta nhận xét thấy phương trình (1.1) mà ta đặt ra là một trường hợp tổng quát hơn của phương trình (1.7) trong bài toán (1.7) – (1.11), có được bằng cách thêm hàm hệ số  t vào trước số hạng ttu . Hai điều kiện biên (1.3) và (1.10) là như nhau, trong khi điều kiện biên (1.2) lại là trường hợp đặc biệt của (1.8) – (1.9), (1.11). Nội dung chính của luận văn gồm các chương mục được trình bày theo thứ tự sau: Chương 1 là phần mở đầu tổng quan về bài toán mà ta sẽ khảo sát trong luận văn, chỉ ra vài kết quả quan trọng đã có trước đó, đồng thời nêu bố cục của luận văn. Chương 2 trình bày một số kết quả chuẩn bị, bao gồm việc nhắc lại một số không gian Sobolev, một số kết quả về các phép nhúng compact giữa các không gian. Chương 3 nghiên cứu sự tồn tại và duy nhất nghiệm yếu của bài toán (1.1) – (1.3). Chương 4 nghiên cứu tính trơn của nghiệm đối với dữ kiện ban đầu. Chương 5 nghiên cứu tính ổn định của nghiệm đối với dữ kiện ban đầu, tức là hàm    0 1 0 1, , , , , , , ,f k u u u f k u u  , nghiệm yếu của bài toán (1.1) – (1.3), là liên tục theo nghĩa mà ta sẽ chỉ ra khi xem xét vấn đề này. Chương 6 nghiên cứu bài toán nhiễu theo hai tham số bé  ,K  :                   22 , 0 0 1 ( ) ( , ), 0 1, 0 , ( ) 0, 0, , 1, 0, ,0 , ,0 , qp tt xx t t t K x t u t u K u u u u f x t x t T P t u t k t s u s ds u t u x u x u x u x                    trong đó các tham số 0 1, , , ,f k u u cho trước. Cụ thể như sau a/ Nghiên cứu dáng điệu tiệm cận của nghiệm yếu ,Ku u  của bài toán  ,KP  khi 0K  và 0  . b/ Nghiên cứu khai triển tiệm cận của nghiệm yếu ,Ku u  của bài toán  ,KP  theo hai tham số bé K, λ, có nghĩa là có thể xấp xỉ nghiệm ,Ku  bởi một đa thức theo hai biến K, λ:     , , , i jijK i j N u x t U x t K      theo nghĩa cần phải chỉ ra các hàm   ,ijU x t và thiết lập đánh giá      12 2, * , N i j ijK N i j N u U x t K C K         theo một chuẩn thích hợp *|| || , với các tham số dương K, λ đủ bé và hằng số NC độc lập đối với các tham số bé này. Chương 7 trình bày một bài toán cụ thể để minh họa phương pháp tìm nghiệm tiệm cận của bài toán nhiễu ở trên. Sau cùng là phần kết luận và danh mục các tài liệu tham khảo. Chương 2 MỘT SỐ KẾT QUẢ CHUẨN BỊ Ta đặt các kí hiệu:    0,1 , 0, , 0TQ T T     . Chúng ta bỏ qua các định nghĩa của các không gian hàm      ,, , m p m pC L W   . Ta có thể xem trong [1]. Để tạo thuận lợi khi trình bày và làm cho luận văn gọn gàng, ta sẽ quy ước một vài kí hiệu vắn tắt như sau:        , , 0, ,2, , , 1 , 0,1,...m p m p p p p m mW W L L W H W p m            Chuẩn trong L2 kí hiệu là || || . Kí hiệu ,  chỉ tích vô hướng trong L2 hay tích đối ngẫu của một phiếm hàm tuyến tính liên tục với một phần tử của không gian hàm. Kí hiệu || ||X dùng để chỉ chuẩn trong một không gian Banach X, và X / là không gian đối ngẫu của X. Ta kí hiệu  0, ;pL T X là không gian Banach các hàm thực đo được  : 0,u T X sao cho     1 0, ; 0 p T pp L T X X u u t dt       với 1 p   và     0, ; 0 supL T X Xt Tu ess u t     khi p   Ta viết              , , , , t tt x xxu t u t u t u t u t u t u t   để chỉ  , ,u x t        2 22 2, , , , , , ,u u u ux t x t x t x tt t x x         , theo thứ tự đó. Ta lại đặt     1 0,1 : 1 0V v H v   và xét tích vô hướng      1 0 , u va u v x x dx x x     Khi đó, V là một không gian con của H1, và trên V, chuẩn  , xVv a v v v  tương đương với chuẩn 1Hv cảm sinh trên V. Từ đây ta có bổ đề sau Bổ đề 2.1. Phép nhúng V ↪  0C  là compact và   0 0,1 , C Vv v v V   Bổ đề 2.2. (Bổ đề Gronwall) Giả sử  và u là các hàm liên tục trên [ , ],a b trong đó  xác định không âm,  là một hằng số dương (không phụ thuộc t). Khi đó, nếu        , ,t a u t s u s ds t a b     thì      exp , ,t a u t s ds t a b       Bổ đề 2.3. [2] (Bổ đề về tính compact của Lions) Cho ba không gian 0 1, ,X X X thỏa 0 1,X X phản xạ, 0X ↪ X compact, X ↪ 1X liên tục Với 0 , 1 ,T p q      , ta đặt     0 10, ; : 0, ;p qW u L T X u L T X   và trang bị cho W một chuẩn    0 10, , 0, ,p qW L T X L T Xv v v  Khi đó, phép nhúng W ↪  0, ;pL T X là compact. Chương 3 SỰ TỒN TẠI VÀ DUY NHẤT NGHIỆM Sau đây chúng ta sẽ nghiên cứu sự tồn tại và duy nhất nghiệm của bài toán (1.1) – (1.3) với những giả thiết về các dữ kiện đầu vào như sau: (GT1) , 0 ; , 2K p q   (GT2) 20 1, u H V u V   (GT3)  2,1 0,k W T (GT4)    2 0, ; 0oH T t     (GT5)  2, t Tf f L Q Định lý 3.1. Nếu các giả thiết (GT1) – (GT5) được thỏa thì bài toán (1.1) – (1.3) có duy nhất một nghiệm yếu u với các tính chất      2 1 20, ; ; 0, ; ; 0, ;t ttu L T V H u L T H u L T L      và    1,0, 0,u W T  Chứng minh. Gọi  jw là một cơ sở đếm được của không gian 2V H . Ta sẽ tìm nghiệm xấp xỉ theo phương pháp Galerkin của bài toán (1.1) – (1.3) dưới dạng 1 ( ) ( ) , m m mj j j u t c t w    trong đó, các hàm hệ số  mjc t thỏa mãn hệ phương trình vi phân thường             0 , , 0 0, t m j mx jx j mu t w t u t w w k t s u s ds             2 2 , , ; 1p qm m m m j jK u t u t u t u t w f t w j m         0 0 1 ,0 m m m mj j j u u w u       mạnh trong 2H   1 1 1 ,0 m m m mj j j u u w u        mạnh trong 1H (3.1) Với các giả thiết (GT1) – (GT5) được thỏa, hệ phương trình vi phân (3.1) có nghiệm um trên  0, mT tương ứng. Các ước lượng sau cho phép ta chọn mT T với mọi m. Ước lượng tiên nghiệm 1. Nhân lần lượt các phương trình thứ j của (3.1) 1 với  mjc t , sau đó cộng các vế tương ứng với nhau rồi lấy tích phân theo biến thời gian từ 0 đến t, ta thu được đẳng thức sau              2 0 0 0 2 0, 0, 2 0 0, t t m m m m mS t S u t k t u d k u s ds                2 0 0 0 2 0, 0, t t s mx m ms u s ds u s k s u d ds             0 2 , t mf s u s ds  (3.2) trong đó            2 2 0 2 2p q t p q m m mx m mL L KS t u t t u t u t u s ds p        (3.3) Từ (3.3) và (GT4), ta thấy rằng      2 20m m mxS t u t u t  (3.4) Từ (3.4) và bổ đề 1, ta có đánh giá          0 0 0,1 , mm m mxC S t u x t u t u t      (3.5) Dưới đây là các đánh giá vắn tắt cho từng số hạng trong vế phải của (3.2).             2 2 0 00 0 12 0, 0, 0, 0, t t m m m mu t k t u d u t k t u d                    22 0 0 0 1 t t m mS t k d S s ds      (3.6)        2 00 0 2 0 2 0 0, t t m m k k u s ds S s ds  (3.7)        2 00 0 t t mx m s s u s ds S s ds      (3.8)       0 0 2 0, 0, t s m mu s k s u d ds          22 00 0 0 0 1t t s s m mS s ds k d S d ds                   22 00 0 0 t t t m m tS s ds k d S d          (3.9)            2 0 0 0 0 2 , 2 t t t t m m mf s u s ds f s u s ds f s ds S s ds        (3.10) Từ (3.2) và (3.6)–(3.10), đồng thời chọn tham số 12  ta có đánh giá          2 00 0 2 0 2 2 t t m m m s S t S f s ds S s ds               2 22 0 00 0 4 04 3 t t m k k t k d S s ds              (3.11) Do H1 ↪  0C  , nên từ (3.1)2,3 và (3.3) ta nhận được   12 0mS M với mọi m (3.12) trong đó hằng số M1 chỉ phụ thuộc vào μ, K, u0 và u1. Theo (3.12) và (GT5) ta có   21 0 2 t TM f s ds A  với mọi  0,t T (3.13) Mặt khác, theo (GT3) ta lại có       2 22 0 00 4 04 3 t T k k t k d B       với mọi  0,t T (3.14) Với các đánh giá (3.11), (3.13) và (3.14) vừa thu được, ta suy ra rằng       00 2t m mT T s S t A B S s ds          với mọi  0,t T (3.15) Áp dụng bổ đề Gronwall cho (3.15), rồi dựa vào (GT4) ta có kết quả quan trọng sau     00 2 exp t m TT T s S t A B ds C               với mọi  0,t T (3.16) Ước lượng tiên nghiệm 2. Đạo hàm 2 vế của (3.1)1 theo biến thời gian, ta có          , , ,m j mx jx mx jxu t w t u t w t u t w                 0 0 0, 0 0 0, t j m j mw k t s u s ds w k u t   (3.17)              2 21 1 , , ; 1p qm m m m j jK p u t u t q u t u t w f t w j m           Nhân lần lượt các phương trình thứ j với  mjc t , sau đó cộng lại rồi lấy tích phân 2 vế theo biến thời gian từ 0 đến t, ta có             0 10 2 , 2 0 ,m m mx mx mx mxX t X t u t u t u u                2 0 0 2 , 3 t t mx mx mxs u s u s ds s u s ds                   0 0 0 2 0, 0, 2 0 0, 0, t s t m m m mu s k s u d ds k u s u s ds                 1 2 0 0 2 1 , , , t p m m mK p u x s u x s u x s dxds           0 2 , t mf s u s ds   (3.18) với               12 2 2 2 0 0 2 1 , , t q m m mx m mX t u t t u t q u x s u x s dxds           Từ đây và bổ đề 1, ta có các đánh giá      2 20m m mxX t u t u t   (3.19)            0 0 0,1 , mm m mx mC V X t u x t u t u t u t          (3.20) Lấy tích phân từng phần cho 2 số hạng thứ 6 và thứ 7 trong vế phải của (3.18), ta có                    20 10 2 0 0 0 0 0, 2 0 0, 0,m m m m m m mX t X k u u k u t k u t u t            0 12 0 , 2 ,mx mx mx mxu u t u t u t            2 2 0 0 2 0 0, 3 t t m mxk u s ds s u s ds                 0 0 2 , 2 0, 0, t t mx mx m ms u s u s ds u t k t s u s ds                 2 0 0 0 0, 2 0, 0, t s m m mk u t u s k s u d ds                1 2 0 0 2 1 , , , t p m m mK p u x s u x s u x s dxds           0 2 , t mf s u s ds   (3.21) Dựa vào (3.4), (3.5), (3.16), (3.19), (3.20) và bổ đề 2.1, ta có những ước lượng về các số hạng trong vế phải của (3.21) như sau          22 0 2 0 0, 0, 0m m mT Ck u t u t k X t   (3.22)    0 1 0 12 0 , 2 0 . .mx mx m mV Vu u u u   (3.23)           22 0 2 , Tmx mx m Ct u t u t X t t      (3.24)        2 00 0 2 0 2 0 0, t t m m k k u s ds X s ds   (3.25)        2 00 0 33 t t mx ms u s ds s X s ds     (3.26)          22 00 0 0 2 , t t t mx mx m TCs u s u s ds s ds X s ds       (3.27)           2 2 00 0 2 0, 0, t t m m m TCu t k t s u s ds X t k d             (3.28)      2 2 0 0 0 0,m Tk Ck u t    (3.29)           2 2 00 0 0 0 2 0, 0, t s t t m m m TtCu s k s u d ds k d X s ds                (3.30)         1 2 0 0 2 1 , , , t p m m mK p u x s u x s u x s dxds               1 12 2 0 0 0 2 1 , , , p t p m m m TCK p u x s u x s u x s dxds                     2 1 1 0 00 2 1 , p t m mp TC K p X s u x s dxds                 2 1 00 2 1 p t m mp TC K p X s u s ds              2 1 00 2 1 p t mp TC K p X s ds      (3.31)        2 0 0 0 2 , t t t m mf s u s ds f s ds X s ds      (3.32) Từ (3.21) – (3.32), đồng thời chọn tham số 14  , ta có                20 1 0 14 0 2 0 0 0 0 0, 8 0 . .m m m m m m mV VX t X k u u k u t u u             2 22 02 0 0 4+ 4 0 0 4 4 t TC k k t k d                 22 2 20 0 0 0 t t t T s ds f s ds t k d C                       2 1 0 00 0 8 0 12 8 1 12 p t mp TCk s K p X s ds                  (3.33) Từ (GT2) – (GT4) và (3.19), ta có kết quả                   2 0 1 0 1 2 4 0 2 0 0 0 0 0, 8 0 . . 4 0 (3.34) m m m m m m m m V VX t X k u u k u t u u u M        trong đó hằng số M2 chỉ phụ thuộc vào μ, k, K, u0 và u1. Hơn nữa, từ (3.1) và (1.2), bằng phương pháp tích phân từng phần cho số hạng chưa μ của vế trái, ta có            2 2 20 0 0 1 10 0 , 0 , 0 0 , 0p qm mxx m m m m m m mu u u K u u u u u f u         Do đó      2 20 0 0 1 10 0 0p qm mxx m m m mu u K u u u u f       Kết quả này cùng với (GT5) và (3.1)2,3 cho ta   30mu M  (3.35) trong đó hằng số M3 chỉ phụ thuộc vào μ, f, K, λ, u0 và u1. Mặt khác, do H1 ↪  0C  , nên ta có   24 Tt D  a.e., trên  0,T (3.36) Tổng hợp các kết quả (3.33) – (3.36), ta có đánh giá     0 ( ) , t m T mX t E q s X s ds   (3.37) trong đó          2 1 0 0 0 8 0 12 ( ) 8 1 12 p p TCk s q s K p            và ET là một hằng số chỉ phụ thuộc vào T. Từ (GT3), (GT4), (3.37), theo bổ đề Gronwall, ta thu được kết quả sau 0 ( ) exp ( ) , t m T TX t E q s ds F  với mọi  0,t T (3.38) Qua giới hạn. Từ (3.4), (3.5), (3.16), (3.19), (3.20) và (3.38) ta suy ra các dãy hàm dưới đây bị chặn trên các không gian hàm L tương ứng. Từ đó, ta có thể trích ra một dãy con của dãy  mu vẫn kí hiệu là  mu sao cho   trong 0, ;mu u L T V yếu *  1 trong 0, ;mu u L T H  yếu *  2 trong 0, ;mu u L T L  yếu *      1,0, 0, trong 0,mu u W T   yếu * (3.39) Nhờ bổ đề về tính compact của Lions [2], từ (3.38), ta lại rút ra được một dãy con của dãy  mu vẫn kí hiệu là  mu sao cho mu u mạnh  2trong TL Q mu u  mạnh  2trong TL Q    0, 0,mu u   mạnh   0trong 0,C T (3.40) Mặt khác, nhờ vào bất đẳng thức    1 , , , 0, 0x x y y R x y x y R R R            (3.41) nên từ (3.5), (3.16), (3.20), (3.38) và (3.40)1–2, ta suy ra rằng m mu u u u   mạnh  2trong TL Q m mu u u u      mạnh  2trong TL Q (3.42) Trong (3.1), ta cho m  . Khi đó, dựa vào các kết quả (3.39), (3.40) và (3.42) ta chỉ ra rằng có hàm u thỏa mãn bài toán                       0 2 2 1 , , 0 0, , , , t x x p q u t v t u t v v k t s u s ds K u t u t u t u t v f t v v H               (3.43)   00u u và   10u u  (3.44) Từ (3.43) và các giả thiết (GT4) – (GT5) ta có kết quả  2 2 21 0, ;p qxxu u K u u u u f L T L            (3.45) Do đó  20, ;u L T V H  . Việc chứng minh sự tồn tại của nghiệm yếu hoàn tất. Tính duy nhất của nghiệm. Giả sử u1, u2 là hai nghiệm yếu của bài toán (1.1) – (1.3). Khi đó, ta có           2 1 2 1, 0, ; ; 0, ; ; 0, ; 0, 0, 1;2 i i i i u L T V H u L T H u L T L u W T i              Ta đặt 1 2u u u  . Khi đó, u thỏa mãn bài toán sau                             0 2 2 1 1 2 2 2 2 1 1 1 2 2 , , 0 0, , , 0, t x x p p q q u t v t u t v v k t s u s ds K u t u t u t u t v u t u t u t u t v v H  
Luận văn liên quan