Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 37 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
37
Dung lượng
594,07 KB
Nội dung
Khóa luận tốt nghiệp Trường ĐHSP Hà Nội LỜI CẢM ƠN Trong q trình thực khố luận tốt nghiệp em nhận nhiều quan tâm, giúp đỡ thầy cô giáo bạn sinh viên Em xin chân thành cám ơn thầy khoa tốn trường ĐHSPHà Nội 2, thầy tận tình dạy dỗ em năm học vừa qua tạo điều kiện giúp đỡ em hồn thành thành khóa luận tốt nghiệp Em xin bày tở lòng biết ơn sâu sắc tới TS Bùi Kiên Cường người trực tiếp hướng dẫn, bảo tận tình cho em suốt q trình thực khóa luận tốt nghiệp Do cịn hạn chế thời gian nên đề tài em khơng tránh khỏi thiếu sót Em mong nhận giúp đỡ góp ý quý thầy cô bạn để đề tài em hoàn thiện Em xin chân thành cám ơn! Hà nội, tháng 05 năm 2013 Sinh viên Vũ Thị Cúc Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội LỜI CAM ĐOAN Khóa luận tốt nghiệp cơng trình nghiên cứu em hướng dẫn TS Bùi Kiên Cường Trong trình nghiên cứu thực đề tài em có tham khảo số tài liệu (đã nêu phần tài liệu tham khảo) Hà nội, tháng 05 năm 2013 Sinh viên Vũ Thị Cúc Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội MỤC LỤC LỜI MỞ ĐẦU 1.Lý chọn đề tài 2.Mục đích nghiên cứu Chương 1: Phương Trình Truyền Sóng 1.1 Một vài khái niệm mở đầu 1.1.1 Chuyển động hàm điều hòa đơn giản 1.1.2 Sóng đứng sóng dịch chuyển 1.1.3 Điều hòa chồng chất âm 10 1.2 Nguồn gốc phương trình truyền sóng 11 1.2.1 Bài toán 11 1.2.2 Giải toán 11 1.3 Nghiệm phương trình truyền sóng 14 1.3.1 Nghiệm phương trình truyền sóng 14 1.3.2 Sự chồng chất sóng đứng 18 1.4 Ví dụ 22 Chương 2: Phương trình truyền nhiệt 24 2.1 Nguồn gốc phương trình truyền nhiệt 24 2.2 Trạng thái ổn định phương trình truyền nhiệt đĩa 26 2.3 Ví dụ 30 Kết luận 35 Tài liệu tham khảo 36 Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội LỜI MỞ ĐẦU 1.Lý chọn đề tài Trong giải tích cổ điển phương trình đạo hàm riêng, sinh viên tiếp cận giải tích Fourier ứng dụng Qua tìm hiểu ta thấy rằng, chất giải tích Fourier bắt nguồn từ tượng vật lý bên thực tiễn Để làm sáng tỏ chất làm tài liệu cho khóa sau, em chọn đề tài “ Căn nguyên giải tích Fourier” để làm tài liệu cho khóa luận tốt nghiệp 2.Mục đích nghiên cứu Qua đề tài nghiên cứu ta làm rõ được: Cơ sở vật lý giải tích Fourier Bản chất vật lý phương pháp tách biến Bản chất vật lý phương pháp D’Alembert 3.Nhiệm vụ nghiên cứu Nghiên cứu nhằm đạt mục đích đề 4.Đối tượng, phạm vi nghiên cứu Đối tượng nghiên cứu: giải tích Fourier Phạm vi nghiên cứu: phương trình truyền sóng, phương trình truyền nhiệt 5.Phương pháp nghiên cứu Nghiên cứu lý luận, phân tích tổng hợp đánh giá Vũ Thị Cúc K35E – Sp Toán Khóa luận tốt nghiệp Trường ĐHSP Hà Nội 6.Cấu trúc khóa luận Ngồi phần mở đầu, kết luận, tài liệu tham khảo, khóa luận gồm chương: Chương 1: Phương trình truyền sóng Chương 2: Phương trình truyền nhiệt Trong suốt trình nghiên cứu ẹm nhận giúp đỡ tận tình thầy tổ giải tích khoa tốn trường ĐHSP Hà Nội đặc biệt thầy giáo Bùi Kiên Cường Em xin bày tỏ lòng biết ơn sâu sắc tới thầy Em mong nhận đóng góp ý kiến quý báu quý thầy cô bạn sinh viên để đề tài hoàn thiện Em xin chân thành cám ơn! Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Chương Phương Trình Truyền Sóng 1.1 Một vài khái niệm mở đầu 1.1.1 Chuyển động hàm điều hòa đơn giản Chuyển động hàm điều hịa đơn giản mơ tả có trạng thái hệ thống dao động đơn giản (gọi dao động hàm điều hịa đơn giản), sở tự nhiên để bắt đầu nghiên cứu dao động Xét vật khối lượng m gắn vào lò xo nằm ngang Một đầu gắn vào tường, giả thiết hệ nằm bề mặt khơng có ma sát Chọn trục tọa độ có gốc trùng với trọng tâm vật vật trạng thái nghỉ (khi lị xo khơng bị kéo hay nén), hình Kéo vật rời khỏi vị trí cân ban đầu m y m y y(t) Hình Chuyển động hàm điều hịa đơn giản thả ra, chuyển động hàm điều hòa đơn giản Chuyển động mơ tả tốn học ta sử dụng phương trình vi phân để thể thay đổi chuyển động vật Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Gọi y t vị trí vật thời điểm t Ta giả sử rằng, lò xo trạng thái lý tưởng để định luật Hooke’s thỏa mãn: Lực phản hồi F tác dụng lò xo lên vật cho F ky t Ở k đại lượng vật lý gọi số lò xo Áp dụng định luật Newton’s (lực = khối lượng x gia tốc), ta có: ky t my t , Với y đạo hàm bậc hai hàm y , biến số t Với c k m , phương trình vi phân trở về: y t c y t (1) Nghiệm tổng quát phương trình (1 ) cho công thức: y t a cos ct b sin ct , với a, b số Tất hàm số có cơng thức thỏa mãn phương trình (1), rõ ràng, hàm nghiệm phương trình cho Trong biểu y t trên, với đại lượng c cho, a b đồng thời số thực Để xác định nghiệm riêng phương trình, ta phải đặt điều kiện ban đầu cho hai số chưa biết a b Chẳng hạn, cho y y vị trí ban đầu vận tốc ban đầu vật nghiệm tốn vật lý cho cơng thức: y t y cos ct Vũ Thị Cúc y sin ct c K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Đơn giản ta kiểm tra rằng: tồn số A R cho: a cos ct b sin ct A cos ct Theo cách lý giải vật lý trên, gọi A a b biên độ dao động, c tần số tự nhiên, pha (xác định từ bội số nguyên 2 ) 2 c chu kỳ chuyển động Đồ thị tiêu biểu hàm số A cos ct minh họa hình 2, hình đại diện đường lượn sóng, thu tịnh tiến dãn (hoặc co lại) đồ thị hàm số cos t Ta quan sát hai ý phần ví dụ chuyển động hàm điều hịa đơn giản Chú ý diễn tả toán học phần lớn hệ dao động bản, gọi dao động điều hòa đơn giản, bao gồm hàm lượng giác sở sin t , cos t Hình Đồ thị hàm số A cos ct Điều quan trọng sau ta nhắc lại mối quan hệ hàm với số phức, công thức Euler’s eit cos t i sin t Nhận xét thứ hai hàm điều hòa đơn giản coi hàm theo biến thời gian với hai điều kiện ban đầu, để xác định vị trí, để xác định Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội vận tốc (ví dụ, thời điểm t ) Nhiều hệ dao động tổng qt có chung tính chất này, thấy bên 1.1.2 Sóng đứng sóng dịch chuyển Như ta thấy, dao động sợi dây xem chuyển động sóng chiều Bây giờ, ta mô tả hai loại chuyển động có biểu diễn đồ thị đơn giản Loại thứ nhất, ta khảo sát sóng đứng Đó chuyển động mô tả đồ thị y u x, t theo biến thời gian t , cho y u x, t hình Nói cách khác, tồn hình dạng ban đầu y x biểu diễn sóng thời điểm t , phần tử khuyêch đại t , phụ thuộc vào thời gian t , cho y u x, t , với: u x, t x t Sự tự nhiên sóng đứng gợi ý cho ta ý tưởng phương pháp tách biến phương trình đạo hàm riêng Loại chuyển động thứ hai sóng thường thấy tự nhiên sóng dịch chuyển Biểu diễn đặc biệt đơn giản: Hình Sóng đứng hai thời điểm khác t t t0 Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội tồn hình dạng ban đầu F x cho u x, t F x t Khi t tăng, biên dịch chuyển sang phải đoạn ct , với c số dương, hay: u x, t F x ct Về đồ thị, vị trí mơ tả hình Hình Sóng dịch chuyển hai thời điểm khác t t t0 Vì sóng dịch chuyển theo thời gian t với tốc độ c , nên vận tốc sóng số Hàm F x ct sóng dịch chuyển chiều chuyển động sang bên phải Tương tự, u x, t F x ct sóng dịch chuyển chiều chuyển động sang bên trái 1.1.3 Điều hòa chồng chất âm Hiện tượng vật lý cuối mà ta muốn đề cập đến (không vào chi tiết) tượng mà nhạc công phát từ lâu Đó hịa âm Âm nghe hợp âm cách kết hợp âm bội sinh tương ứng từ âm sắc nhạc cụ Ý tưởng kết hợp hay chồng chất âm thực toán học dựa khái niệm tuyến tính, ta thấy bên Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Do đó, dự đoán Am hệ số Fourier sine thứ n f , là: (6) A n 2 f x sin nxdx Chúng ta quay trở lại công thức công thức tương tự sau Ta chuyển câu hỏi chuỗi hàm sin Fourier đoạn 0, sang câu hỏi tổng quát đoạn , Nếu ta biểu diễn hàm f dạng chuỗi hàm sin đoạn 0, , khai triển đoạn , ta thác triển f đoạn cách biến thành hàm lẻ Tương tự, đặt câu hỏi hàm chẵn g x đoạn , biểu diễn dạng chuỗi hàm cos hay không, cụ thể: g x Am cos mx m 0 Tổng quát hơn, hàm F đoạn , biểu diễn thành f g , f hàm lẻ g hàm chẵn, nên đặt câu hỏi F viết lại dạng m 1 m 0 F x A m sin mx Am cos mx, Hay cách áp dụng đồng thức Euler’s eix cos x i sin x , ta mong muốn F có dạng: F x ameimx m 1 Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Tương tự với công thức (6) , ta sử dụng kết : 2 e imx 0, e inx dx 1, nm n m để nhận thấy an 2 F x e inx dx Đại lượng a n gọi hệ số Fourier thứ n F Ta viết lại tốn đặt sau: Câu hỏi: Cho hàm F xác định , , với số Fourier xác định trên, liệu điều sau F x (7) ae m m imx hay khơng? Cơng thức tốn này, có dạng số mũ phức, ta áp dụng hầu hết trường hợp Quay trở lại với phương trình truyền sóng Để xây dựng cơng thức hóa tốn cách xác, phải đặt hai điều kiện ban đầu gặp chuyển động điều hịa đơn giản sóng đứng Địi hỏi, u thỏa mãn phương trình vi phân hai điều kiện: u x,0 f x Vũ Thị Cúc u x,0 g x , t K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Ở f g hàm cho trước Lưu ý điều phù hợp với (4), f với g biểu diễn dạng: f x A m 11 m sin mx g x mB m 11 m sin mx 1.4.Ví dụ: Sợi dây kéo Bây ta áp dụng lập luận vào toán kéo sợi dây Để đơn giản, ta chọn đơn vị cho sợi dây nằm đoạn 0, , thỏa mãn phương trình truyền sóng với c Giả sử sợi dây kéo đến độ cao h điểm p với p ; vị trí ban đầu Nghĩa là, ta kéo sợi dây cho vị trí sợi dây có dạng hình tam giác cho bởi: xh p , f x h x , h 0 x p p x , Được minh họa đồ thị 8: Hình Vị trí ban đầu sợi dây bị kéo Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Ta chọn vận tốc ban đầu g(x) Khi đó, ta tính hệ số Fourier f , giả sử kết câu hỏi cho đưa trước cơng thức (5) dương, ta có: f x Am sin mx với Am m 1 2h sin mp m p p Do u x, t Am cos mt sin mx (8) m 1 lưu ý chuỗi hội tụ tuyệt đối Nghiệm biểu diễn dạng sóng dịch chuyển Với: u x, t (9) f x t f x t Ở đây, f x xác định x sau: ta thác triển f đoạn , cách làm cho f lẻ sau thác triển f toàn chu kỳ 2 nghĩa f x 2k f x với tất số nguyên k Ta thấy (8) kéo theo (9) theo cách đồng thức lượng giác cos v sin u Vũ Thị Cúc sin u v sin u v 2 K35E – Sp Toán Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Chương Phương Trình Truyền Nhiệt 2.1 Nguồn gốc phương trình truyền nhiệt Xét kim loại vơ hạn mà ta xem như mặt phẳng R , giả sử có phân bố nhiệt lên thời điểm t Kí hiệu nhiệt độ điểm x, y thời điểm t hàm u x, y , t Xét hình vng nhỏ có tâm điểm x0 , y0 có cạnh song song với trục tọa độ có độ dài cạnh h ( hình 9) Năng lượng nhiệt S thời điểm t cho công thức: H t u x, y, t dxdy, S với số, gọi hệ số dẫn nhiệt chất liệu Do đó,lượng nhiệt truyền tới S là: H u dxdy, S t t xấp xỉ bằng: h2 u ( x0 , y0 , t ), t diện tích S h Giờ ta áp dụng định luật Newton’s hạ nhiệt, trạng thái mà nhiệt chuyển từ nơi có nhiệt độ cao vào nơi có nhiệt độ thấp theo tỉ lệ,tức gradient Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Hình Nhiệt truyền qua hình vng nhỏ Do đó, lượng nhiệt truyền qua cạnh thẳng đứng bên phải h u x0 h 2, y0 , t , x gọi độ dẫn vật liệu Lập luận tương tự với cạnh khác hình vng, ta tổng nhiệt tỏa hình vng S là: u u u x0 h 2, y0 , t x0 , y0 h 2, t x0 , y0 h 2, t y y x h Áp dụng định luật giá trị trung bình cho h tiến dần , ta được: u 2u 2u ; t x y Phương trình gọi phương trình truyền nhiệt phụ thuộc vào thời gian, thường gọi tắt phương trình truyền nhiệt Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội 2.2 Phương trình trạng thái cân nhiệt đĩa Sau thời gian, khơng có thêm biến đổi nhiệt, hệ đạt tới cân nhiệt u t Trong trường hợp này, phương trình nhiệt phụ thuộc vào thời gian rút gọn thành phương trình trạng thái cân nhiệt: 2u u x y (10) 2 Tốn tử x y có vai trị quan trọng tốn học vật lý thường viết tắt gọi tốn tử Laplace Do đó, phương trình cân nhiệt viết dạng: u 0, nghiệm phương trình gọi hàm điều hịa Xét đĩa đơn vị mặt phẳng D x, y R : x y2 , biên đường tròn đơn vị C Trong hệ tọa độ cực r , , với r 2 , ta có: D r, : r 1 C r , : r 1 Bài toán gọi toán Dirichlet (đối với toán tử Laplace hình trịn đơn vị), giải phương trình cân nhiệt đĩa đơn với điều kiện biên u f C Điều tương ứng với cố định hàm phân phối nhiệt độ xác định trước hình trịn, đợi thời gian dài, nhìn vào tỏa nhiệt bên đĩa Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Hình 10 Bài tốn Dirichlet đĩa Khi phương pháp tách biến áp dụng vào giải (10) có vấn đề khó giải điều kiện biên khơng dễ biểu diễn theo dạng tọa độ thình chữ nhật Vì điều kiện biên biểu diễn tốt tọa độ cực r , , cụ thể u 1, f ta viết lại toán tử Laplace dạng tọa độ cực Áp dụng quy tắc dây chuyền ta có : 2u u 2u u r r r r Chúng ta nhân hai vế với ,vì u r2 2u u 2u r r r Dùng phương pháp tách biến, tìm nghiệm dạng u , F G , ta được: r F r rF r F r G G Điều xảy hai vế phương trình số, số, gọi Vì Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội G G 0, r F r rF r F r Do G tuần hoàn với chu kỳ 2 , ta biết trước m2 m số nguyên, đó: sin m G A cos m B Áp dụng đồng thức Euler’s eix cox i sin x , ta viết lại G dạng số mũ phức: G Aeim Be im Với m2 m hai nghiệm đơn giản phương trình F F r r m F r r m Nếu m F r F r log r hai nghiệm Nếu m ý r m dần tới trạng thái không bị chặn r dần tới Do F G khơng bị chặn gốc Tương tự, m F r log r Chúng ta bỏ qua nghiệm mâu thuẫn với trực giác Do đó, cịn lại hàm đặc biệt sau: um r , r m e im , với m Z Ta thấy nhận xét quan trọng công thức (10) tuyến tính, tương tự sợi dây dao động, ta chồng chất nghiệm đặc biệt vào để thu công thức nghiệm tổng quát : u r , Vũ Thị Cúc a m m m r eim K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Nếu biểu diễn tất nghiệm phương trình trạng thái cân nhiệt với hàm f phù hợp ta có: u 1, a m m eim f Trong điều kiện ta đặt câu hỏi : Cho hàm f xác định đoạn 0,2 với f f 2 có tìm hệ số am thỏa mãn: f ae m m im hay khơng? Chú thích lịch sử: D’Alembert (năm 1747) lần giải phương trình chuỗi dao động cách sử dụng phương pháp sóng dịch chuyển Nghiệm Euler chỉnh sửa sau Năm 1753, D.Bernoulli giả định phương pháp dựa vào chuỗi Fourier đưa (4), Euler khơng hồn tồn đồng ý q chung chung hợp lý hàm tùy ý mở rộng chuỗi Fourier D’Alembert nhà tốn học khác hồi nghi phương pháp Mối hoài nghi Fourier dập tắt (năm 1807) ơng nghiên cứu phương trình nhiệt Nghiên cứu ông chứng thuyết phục hàm biểu diễn chuỗi Fourier Vũ Thị Cúc K35E – Sp Toán Khóa luận tốt nghiệp Trường ĐHSP Hà Nội 2.3 Ví dụ: Giải tốn Drichlet sau: u xx u yy 0, x2 y2 u x y 1 x x xy Bài giải Đặt: r 1 x r cos 2 y r sin u u u u r r r r Khi nghiệm tốn có dạng: a0 n r an cos n bn sin n n1 a u 1, 1n an cos n bn sin n n 1 g cos cos3 2cos sin u r , cos cos 2cos 1 cos cos cos 2cos 1 cos cos3 cos cos cos3 4 4 Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Sử dụng phương pháp đồng hệ số hai vế ta được: a0 a 1 a3 a an n 4, bn n 1 1 u r , r cos r cos3 4 4 Theo cách đăt ta có: x r cos y r sin r x y arctan y x Nghiệm tốn có dạng: 1 u x, y x y cos arctan y x Vũ Thị Cúc 1 x y cos3 arctan y x 4 K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Kết luận Những vấn đề mà ta đặt giải chương điểm khởi đầu cho giả tích Fourier mà Joseph Fourier (1786-1830) người tin hàm số F biểu diễn dạng cơng thức (7) Hay nói cách khác, ý tưởng ơng hàm tổ hợp tuyến tính hàm lượng giác sin mx cos mx , với m số nguyên.Mặc dù ý tưởng ngầm xuất cơng trình trước đó, Fourier cho người tiền nhiệm chưa có nhìn tổng qt Ông sử dụng ý tưởng nghiên cứu khuếch xạ nhiệt; điều khơi mào cho chủ đề “Giải tích Fourier” Một lần em xin cám ơn thầy Bùi Kiên Cường tồn thể thầy khoa tốn Đặc biệt thầy tổ giải tích tạo điều kiện cho em hồn thành khóa luận Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Trường ĐHSP Hà Nội Tài liệu tham khảo Nguyễn Xuân Liêm, Giải tích –tập 2, Nxb Giáo Dục, 2006 Nguyễn Đình Trí-Nguyễn Hồ Quỳnh, Bài tập toán học cao cấp-tập 2, Nxb Giáo Dục 2001 Trần Đức Long-Nguyễn Đình Sang-Hồng Quốc Tồn, Giáo trình giải tích 2, Nxb ĐH Quốc gia Hà Nội, 2001 Nguyễn Thừa Hợp, Phương trình đạo hàm riêng, Nxb ĐH Quốc Gia, 2006 Elias M.Stein, Rami Shakarchi, Fourier Analysis, An Introduction, Psinceton leeture in Analysis, 2002 Vũ Thị Cúc K35E – Sp Tốn Khóa luận tốt nghiệp Vũ Thị Cúc Trường ĐHSP Hà Nội K35E – Sp Toán ... cận giải tích Fourier ứng dụng Qua tìm hiểu ta thấy rằng, chất giải tích Fourier bắt nguồn từ tượng vật lý bên thực tiễn Để làm sáng tỏ chất làm tài liệu cho khóa sau, em chọn đề tài “ Căn nguyên. .. chọn đề tài “ Căn nguyên giải tích Fourier? ?? để làm tài liệu cho khóa luận tốt nghiệp 2.Mục đích nghiên cứu Qua đề tài nghiên cứu ta làm rõ được: Cơ sở vật lý giải tích Fourier Bản chất vật lý... tượng, phạm vi nghiên cứu Đối tượng nghiên cứu: giải tích Fourier Phạm vi nghiên cứu: phương trình truyền sóng, phương trình truyền nhiệt 5.Phương pháp nghiên cứu Nghiên cứu lý luận, phân tích