SỬ DỤNG PHƯƠNG PHÁP TÁCH BIẾN FOURIER ĐỂ TÌM DAO ĐỘNG CỦA CÁC MÀNG CÓ HÌNH DẠNG ĐẶC BIỆT

Chia sẻ: Nguyen Nhi | Ngày: | Loại File: PDF | Số trang:11

Thêm vào BST

Báo xấu

339
lượt xem 44
download

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

Phương pháp tách biến Fourier đã được sử dụng để nghiên cứu dao động của các màng có hình dạng đặc biệt. Hai bài toán dao động tổng quát của màng chữ nhật và màng tròn đã được đưa ra như một ví dụ minh hoạ về tính hiệu quả và đơn giản của phương pháp tính toán này khi giải các bài toán vật lý có liên quan đến phương trình vi phân đạo hàm riêng.

Chủ đề:

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

Lưu

Nội dung Text: SỬ DỤNG PHƯƠNG PHÁP TÁCH BIẾN FOURIER ĐỂ TÌM DAO ĐỘNG CỦA CÁC MÀNG CÓ HÌNH DẠNG ĐẶC BIỆT

SỬ DỤNG PHƯƠNG PHÁP TÁCH BIẾN FOURIER ĐỂ TÌM DAO ĐỘNG CỦA CÁC MÀNG CÓ HÌNH DẠNG ĐẶC BIỆT Nguyễn Thị Minh, Phạm Hữu Kiên Khoa Vật l ý trường Đại học Sư pham Thái Nguyên Tóm tắt: Phương pháp tách biến Fourier đã được sử dụng để nghiên cứu dao động của các màng có hình dạng đặc biệt. Hai bài toán dao động tổng quát của màng chữ nhật và màng tròn đã được đưa ra như một ví dụ minh hoạ về tính hiệu quả và đơn giản của phương pháp tính toán này khi giải các bài toán vật lý có liên quan đến phương trình vi phân đạo hàm riêng. Trong nghiên cứu này chúng tôi đã trình bày lời giải chi tiết về dao động của màng chữ nhật và màng tròn. Kết quả tìm được cũng được chúng tôi thảo luận và nhận xét tương đối đầy đủ trong công trình này. 1. Giới thiệu Dao động tự do của các vật có hình dạng đặc biệt (vuông, chữ nhật, tròn và quạt…) đã được nghiên cứu, tính toán bằng nhiều phương pháp khác nhau. Nhưng đơn giản và hiệu quả hơn cả là đưa bài toán về dạng các phương trình vi phân đạo hàm riêng, sau đó sử dụng các công cụ toán học để giải các phương trình vi phân này như: biến đổi ảnh Laplace, tách biến Fourier,... Mặc dù, đã có rất nhiều giáo trình viết về những phương pháp và kết quả tính toán đối với các dạng dao động này. Xong trên thực tế những giáo trình này vẫn chưa được giảng dạy phổ biến và có tính chất hệ thống tại các trường Đại học ở Việt Nam. Tuy nhiên, những vấn đề này lại thu hút được sự quan tâm rất lớn của những nhà nghiên cứu về Vật lý lý thuyết cũng như các bạn Sinh viên khoa Vật lý ở các trường Đại học Sư phạm. Vì vậy, với mong muốn tìm hiểu, trang bị những kiến thức bổ ích, chúng tôi đã lựa chọn hướng nghiên cứu “Sử dụng phương pháp tách biến Fourier để tìm dao động của các màng có hình dạng đặc biệt”. Với mong muốn tìm lời giải cho các bài toán dao động tự do của các loại màng có hình dạng đặc biệt, chúng tôi đã sử dụng phương pháp tách biến Fourier để tìm nghiệm của các phương trình vi phân đạo hàm riêng mô tả dao động của màng trên. Kết quả tính toán được, chúng tôi đã thảo luận và nhận xét rất rõ ràng. Chúng tôi hy vọng rằng nghiên cứu này có thể làm tài liệu tham khảo tốt cho Sinh viên đã và đang theo học Vật lý ở các trường Đại học Sư phạm, Đại hoc Khoa học Tự nhiên và Học viên ôn thi Cao học…. 2. Phương pháp tính toán Chúng tôi đã sử dụng phương pháp tách biến Fourer để tính toán trong quá trình nghiên cứu. Giả sử tìm nghiệm phương trình [1-3]: 1
U tt − a 2 Δu = G ( x, y, z...t ) , ,, (1) nghiệm thoả mãn phương trình vi phân (1), được tìm bằng cách phân tích hàm U ( x, y, z ,...t ) thành tích các hàm chứa các biến độc lập với nhau, cụ thể là chúng ta đặt: U ( x, y, z ,...t ) = T(t ) X ( x )Y( y ) Z ( z ) ⋅⋅⋅, (2) sau đó thay phương trình (2) vào phương trình (1), kết hợp với điều kiện biên và điều kiện ban đầu sẽ tìm được nghiệm của bài toán. Cụ thể chúng tôi xét chuyển động dao động của một màng mỏng được căng ra trên mặt phẳng (x,y), các dao động chuyển động ngang vuông góc với mặt phẳng của màng dưới sự tác dụng của ngoại lực (hình 1) [1-3, 13]. Đưa vào các đại lượng sau: U = U ( x, y, t ) là đại lượng mô tả dao động ngang của màng tại ( x, y, t ) ρ là mật độ khối lượng trên một đơn vị diện tích của màng. T là sức căng tính trên một đơn vị diện tích. w = w ( x, y, t ) là ngoại lực tác dụng tính trên một đơn vị diện tích. Fe = w ( x, y, t ) ΔxΔy là ngoại lực tác dụng lên yếu tố diện tích ΔxΔy . ∂u Fd = β ΔxΔy là lực hãm hay lực tắt dần tỷ lệ với vận tốc dao động trên yếu tố diện ∂t tích Δx ⋅ Δy , β là hệ số tắt dần. U(x,y,z) y u x X Hình 1. Dao động của màng 2
T ⋅ Δx T ⋅ Δy β2 α2 β1 α1 T ⋅ Δx T ⋅ Δy y − Δy / 2 yx Δy / 2 + x − Δx / 2 x + Δx / 2 Hình 2. Hình chiếu của các lực tác dụng lên yếu tố màng Xét độ dịch chuyển nhỏ, trong đó các góc α1 , α 2 , β1 , β 2 là nhỏ và tổng các lực tác dụng lên yếu tố của màng theo hướng x và y là hình chiếu của sức căng lên mặt phẳng, ta có: T T Tcosα = ≅ ≈ T ≡ To 1+tg 2α 2 ⎛ ∂u ⎞ 1+ ⎜ ⎟ ⎝ ∂x ⎠ Tcosα 2 = Tcosα1 = To , Tcosβ 2 = Tcosβ1 = To . Áp dụng định luật II-Newton, lấy tổng các lực theo hướng vuông góc với yếu tố của màng: ∂ 2u = T Δy ( sin α 2 − sin α1 ) + T Δx ( sin β 2 − sin β1 ) ρΔxΔy ∂ 2t 2 (3) ∂u + w ( x, y, t ) ΔxΔy − βΔxΔy . ∂t Chú ý rằng: ∂u 2 tgα ∂u ⎛ ∂u ⎞ ∂x sin α = ≈ ≈ ;⎜ ⎟ ≈ 0 , ∂x ⎝ ∂x ⎠ 1 + tg 2α 2 ⎛ ∂u ⎞ 1+ ⎜ ⎟ ⎝ ∂x ⎠ do đó, ta có: ⎡ ∂u ⎛ Δx ⎞ ∂u ⎛ Δx ⎞⎤ T Δy ( sin α 2 − sin α1 ) = To Δy ⎢ ⎜ x + , y, t ⎟ − ⎜ x − , y, t ⎟ ⎥ ; ⎣ ∂x ⎝ ⎠ ∂x ⎝ 2 2 ⎠⎦ ⎡ ∂u ⎛ Δy ⎞ ⎤ Δy ⎞ ∂u ⎛ T Δx ( sin β 2 − sin β1 ) = To Δx ⎢ ⎜ x, y + , t ⎟ − ⎜ x, y − ,t ⎟ . 2 ⎠⎥ ⎣ ∂y ⎝ 2 ⎠ ∂y ⎝ ⎦ 3
Chia hai vế phương trình (3) cho To ΔxΔy , rồi lấy giới hạn khi Δx → 0 , ta được phương trình: ρ ⎛ ∂ 2u β ∂u ⎞ ∂ 2u ∂ 2u 1 + w ( x, y , t ) = ⎜ 2 + + ⎟, ρ ∂t ⎠ ∂ 2 x ∂ 2 y To To ⎝ ∂t 1 ⎛ ∂ 2u ∂u ⎞ ∂ 2u ∂ 2u 1 + 2 + w ( x, y , t ) = 2 ⎜ 2 + k ⎟ , hay (4) ∂ x ∂ y To a ⎝ ∂t ∂t ⎠ 2 β To trong đó: a 2 = ;k = ... ρ ρ Phương trình (4) gọi là phương trình dao động của màng. Xét các trường hợp đặc biệt: ●Trường hợp 1: Cho β = 0, w = 0 , thu được phương trình sóng hai chiều ∂ 2u ∂ 2 u 1 ∂ 2 u + = . (5) ∂ 2 x ∂y 2 a 2 ∂t 2 ●Trường hợp 2: Dao động của màng trong trạng thái dừng, suy ra phương trình Poison: ∂ 2u ∂ 2 u 1 + 2 = − w ( x, y ) . (6) ∂x ∂y 2 To r Gọi C là đường cong kín, nó là biên của màng nằm trên mặt phẳng ( x, y ) , gọi n là vectơ pháp tuyến đơn vị hướng ra ngoài đường cong trên biên. Điều kiện biên của dao động của màng có các dạng sau: Điều kiện biên Dirichlet được gọi là điều kiện biên gắn chặt U ( x, y, t ) x , y∈c = f ( x, y ) . (7) Điều kiện biên Neumann được gọi là điều kiện biên tự do ∂u r = gradu n x , y∈c = f ( x, y ) . (8) ∂n Điều kiện hỗn hợp bao gồm cả điều kiện biên gắn chặt và điều kiện biên tự do, có dạng ⎧u ( x, y, t ) x , y∈c = f ( x, y ) ⎪ 1 (9) ⎨ ∂u r ⎪ = gradu.n x , y∈c2 = g ( x, y ) ⎩ ∂n Điều kiện ban đầu dao động của màng là hình dạng ban đầu và vận tốc ban đầu của mọi điêm trên bề mặt của màng khi t = 0 , có dạng: ⎧u ( x, y, o ) = f ( x, y ) ⎪ ⎨ ∂u ( x, y, 0 ) (10) = g ( x, y ) ⎪ ∂t ⎩ 4
3. Kết quả và thảo luận Vận dụng các phương pháp tính toán được trình bày trong mục 2, chúng tôi đã đưa ra lời giải tương đối chi tiết hai bài toán dao động tổng quát có hình dạng đặc biệt chữ nhật và màng tròn. Từ các kết quả thu được, chúng tôi đã đưa ra những nhận xét về mặt vật lý, trên cơ sở đó giải thích một vài hiện tượng quan sát trong thực tế. Bài toán thứ nhất: Tìm dao động ngang của màng chữ nhật 0 ≤ x ≤ L1 , 0 ≤ y ≤ L2 có các mép gắn chặt, với điều kiện biên ban đầu: u ( x, y, 0 ) = Axy ( L1 − x ) ( L2 − y ) ; ut ( x, y, 0 ) = 0 [1-3, 7]. y L2 x O L1 Hình 3. Dao động của màng chữ nhật Bài toán trở thành tìm nghiệm của phương trình vi phân hai biến: utt = a 2 (u xx + u yy ) , (11) =u =u =u =0, thoả mãn các điều kiện biên: u (12) x =0 y =0 y = L2 x=L1 và điều kiện ban đầu: u ( x, y, 0 ) = Axy ( L1 − x ) ( L2 − y ) , ut ( x, y, 0 ) = 0. (13) Hàm dao động xác định trong khoảng 0 ≤ x ≤ L1 , 0 ≤ y ≤ L2 . Dao động của màng chữ nhật này có dạng như sau: ( 2n + 1) π x sin ( 2m + 1) π y , ∞ ∞ u ( x, y, t ) = ∑∑ ( Amn cos ( aλnmt ) + Bmn sin ( aλmnt ) ) sin L1 L2 m =0 n =0 từ điều kiên ban đầu (12) xác định được các hệ số Amn, Bmn , ta có : (( f , X ) ,Y ) = L1 L2 nπξ nπη 4 ∫ ∫ f (ξ ,η ) sin n m d ξ dη . Amn = sin 2 2 L1 L2 L1 L2 Xn Ym 00 ( 2n + 1) π x sin ( 2m + 1) π y dxdy L1 L2 4 ∫ ∫ Axy ( L − x )( L 2 − x ) sin Amn = 1 L1 L2 L1 L2 00 ( 2n + 1) π x ( 2m + 1) π y L1 L2 4A 4 ∫ x ( L1 − x ) sin dx ∫ y ( L2 − y ) sin dy = Ι1Ι 2 , = L1 L2 L1 L2 L1 L2 0 0 5
tính Ι1 , áp dụng tích phân từng phần L1 ( 2n + 1) π x + L1 ( 2n + 1) π x dx L1 L1 Ι1 = − x ( L1 − x ) ∫ ( L − 2 x) cos cos ( 2n + 1) π ( 2n + 1) π 1 L1 L1 0 0 ( 2n + 1) π x dx L1 L1 ∫ ( L1 − 2 x) cos = ( 2n + 1) π L1 0 Áp dụng tích phân từng phần, ta có: L1 ( 2n + 1) π x ( 2n + 1) π x dx L1 L12 2 L12 ( L1 − 2 x ) ∫ sin Ι1 = + sin ( 2n + 1) π ( 2n + 1) π 2 2 L1 L1 2 2 0 0 4 L13 = ( 2n + 1) π 3 3 tương tự, áp dụng tích phân từng phần, ta có: ( 2m + 1) π y L2 4 L23 y ( L2 − y ) sin ∫ = Ι2 = ( 2m + 1) π 3 3 L2 0 64 AL13 L23 Như vậy: Amn = Bmn = 0 . và π 3 ( 2n + 1) ( 2m + 1) 3 3 Như vậy, dạng dao động cụ thể của màng chữ nhất là: ( 2n + 1) π x sin ( 2m + 1) π y sin ∞ 22 64 AL L L1 L2 ∑ u ( x, y , t ) = × × 12 ( 2n + 1) ( 2m + 1) π 3 3 6 m ,n =0 (14) ⎧ ( 2n+1) + ( 2m+1) ⎫ . 2 2 ⎪ ⎪ × cos ⎨π at ⎬ 2 2 L1 L2 ⎪ ⎪ ⎩ ⎭ Từ kết quả tính được trong biểu thức (14), chúng tôi có một số nhận xét về mặt Vật lý như sau: Biểu thức (14) mô tả dao động tự do của màng chữ nhật, vì hàm u ( x, y, t ) khả vi hai lần theo t , hai lần theo x , hai lần theo y và thoả mãn phương trình vi phân (11) cùng với các điều kiện biên (12) và điều kiện ban đầu (13). Mọi điểm ( x, y ) đều dao động điều hoà với biên độ ( 2n + 1) π x sin ( 2m + 1) π y sin ∞ 22 64 AL L L1 L2 ∑ × 12 ( 2n + 1) ( 2m + 1) π 3 3 6 m,n =0 ( 2n + 1) ( 2m + 1) 2 2 với tần số dao động riêng: ω = π a + . 2 L2 L1 2 6
πx πy 64 AL1 L2 2 Nếu m = n = 0 , biên độ dao động là 2 sin sin , các điểm trên biên đều bất π 6 L1 L2 động, khi đó ta nói biên của màng là đường nút. Do 0 ≤ x ≤ L1 , 0 ≤ y ≤ L2 nên các hàm πx πy ≥ 0 , các điểm của màng đều nằm về phía này hay phía kia của mặt phẳng oxy . sin sin L1 L2 L1 L Độ lệch lớn nhất đạt được tại x = và y = 2 . Đó là tâm của màng. Đó là điểm 2 2 bụng của sóng đứng. 3π x πy 64 AL1 L2 Nếu n = 1, m = 0 biên độ dao động là sin sin , các điểm ( x, y ) của 27π 6 L1 L2 L1 màng mà 0 < x < và 0 < y < L2 có biên độ dao động âm, các điểm ( x, y ) của màng mà 2 L1 < x < L1 và 0 < y < L2 có biên độ dao động dương. Do đó nửa trái và nửa phải của màng 2 trong quá trình dao động luôn uốn về hai phía trái nhau. Độ lệch lớn nhất đạt được tại hai ⎛ L L ⎞ ⎛ 3L L ⎞ điểm ⎜ 1 , 2 ⎟ ; ⎜ 1 , 2 ⎟ . ⎝4 2⎠⎝ 4 2⎠ L2 Nếu n = 0, m = 1 , các đường nút là x = 0; x = L1 , y = 0, y = , y = L2 và hai điểm bụng 2 ⎛ L L ⎞ ⎛ L 3L ⎞ là ⎜ 1 , 2 ⎟ ; ⎜ 1 , 2 ⎟ . ⎝2 4⎠⎝2 4 ⎠ Bài toán thứ hai: Khí lý tưởng choán đầy giữa hai mặt cầu đồng tâm S r1 và S r2 (Hình 4). Bán kính của mặt cầu bên trong S r1 thay đổi thêo quy luật: r ( t ) = r1 + ε sin ωt , −∞ < t < +∞, 0 < ε < r2 − r1 , còn bán kính mặt cầu bên ngoài Sr2 không đổi. Tìm dao động của chất khí giữa các mặt cầu [2, 3, 14]. Gọi u ( r , t ) là thế vận tốc của các hạt khí lý tưởng.Phương trình dao động của khí lý ⎛ ∂ 2u 2 ∂u ⎞ ∂ 2u = a2 ⎜ 2 + tưởng có dạng: ⎟, (15) ∂t 2 ⎝ ∂r r ∂r ⎠ với các điều kiện ban đầu: r ( t ) = r1 + ε sin ωt , −∞ < t < +∞, 0 < ε < r2 − r1 , (16) ∂u ∂u = εω cos ωt , = 0. và điều kiện biên: (17) ∂r ∂r r = r1 r = r2 7
r1 r2 Hình 4. Khí lý tưởng choán đầy giữa hai mặt cầu đồng tâm Thế vận tốc u ( r , t ) của các hạt khí lý tưởng mô tả dao động điều hoà được thiết lập với tần số ω . Bằng phương pháp tách biến, đặt: u ( r , t ) = R ( r ) cos ωt . (18) Thay (18) vào (15) ta có phương trình: dR 2 ( r ) 2 dR ( r ) ω 2 + 2 R (r ) = 0 . + (19) dr 2 r dr a Phương trình (19) có dạng là phương trình Betsen, có các tính chất sau: 2m ⎛x⎞ ⎜⎟ ∞ ⎝2⎠ , J 0 ( x ) = ∑ ( −1) m ( m !) 2 m =0 2 m +1 ⎛ x⎞ ⎜⎟ ∞ ⎝2⎠ J1 ( x ) = ∑ ( −1) m , ( m !)( m + 1)! m =0 x ∫ ξ J (ξ ) dξ = xJ ( x ) , 0 1 0 J1 ( x ) = − J 0 ( x ) , ′ x J 0 ( x ) + ( x 3 − 4 x ) J1 ( x ) . ∫ ξ J (ξ ) d ξ = 2 x 2 0 0 8
Sử dụng các tính chất của hàm Betsen và các điều kiện (16), (17) chúng ta xác định được: ⎛ω ωr ⎞ ⎛ω ωr ⎞ ωr 1 ωr 1 2aε ⎜ cos 2 − sin 2 ⎟ cos ω r 2aε ⎜ cos 2 − sin 2 ⎟ sin ω r ⎝a a r2 a⎠ ⎝a a r2 a⎠ a+ a R (r ) = ω ( r1 − r2 ) ω ( r1 − r2 ) ⎛1 1⎞ ⎛1 1⎞ r r ⎜ − ⎟ cos ⎜ − ⎟ cos r1 r2 ⎠ a ⎝ r1 r2 ⎠ a ⎝ Như vậy, dao động của chất khí giữa các mặt cầu: ⎧ ⎫ ⎛ω ωr ⎞ ω r2 1 − sin 2 ⎟ cos ω r ⎪ ⎪ 2aε ⎜ cos ⎝a a r2 a⎠ ⎪ a +⎪ ω ( r1 − r2 ) ⎪ ⎛1 1⎞ ⎪ r ⎪ ⎜ − ⎟ cos ⎪ ⎪ ⎝ r1 r2 ⎠ a ⎪ u ( r, t ) = ⎨ ⎬ cos ωt. (20) ⎪ 2aε ⎛ ω cos ω r2 − 1 sin ω r2 ⎞ ωr ⎪ ⎜ ⎟ sin ⎪ ⎪ ⎝a a r2 a⎠ a⎪ ⎪+ ω ( r1 − r2 ) ⎛1 1⎞ r⎪ ⎪ ⎜ − ⎟ cos ⎪ ⎪ ⎝ r1 r2 ⎠ a ⎩ ⎭ Từ kết quả tính được, theo trong biểu thức (20) chúng tôi có một số nhận xét về mặt Vật lý như sau: Dao động của chất khí giữa các mặt cầu là một dao động điều hoà với biên độ: ωr 1 ωr ω r 2aε ⎛ ω cos ω r2 − 1 sin ω r2 ⎛ω ⎞ ⎞ ωr 2aε ⎜ cos 2 − sin 2 ⎟ cos ⎜ ⎟ sin ⎝a a r2 a ⎝a a r2 a ⎠ ⎠ a+ a R (r ) = ω ( r1 − r2 ) ω ( r1 − r2 ) ⎛1 1⎞ ⎛1 1⎞ r r ⎜ − ⎟ cos ⎜ − ⎟ cos ⎝ r1 r2 ⎠ a ⎝ r1 r2 ⎠ a và tần số góc dao động là ω . Khi r càng lớn thì R ( r ) càng giảm, điều này chứng tỏ các điểm càng xa tâm màng thì dao động cáng nhỏ và đến biên màng thì đứng yên. 4. Kết luận Nghiên cứu này đã thu được một số kết quả sau: i) Đã trình bày một cách có hệ thống và lôgic phương pháp tách biến Fourier và sử dụng phương pháp này để tìm các phương trình dao động tổng quát và dao động cụ thể (Biết điều kiện và điều kiện ban đầu) của màng chữ nhật và tròn. 9
ii) Từ những kết quả tính toán được, chúng tôi có những thảo luận và nhận xét quan trọng về mặt Vật lý của bài toán. Vấn đề tiếp theo là nghiên cứu và tính toán sự truyền nhiệt trong các màng chữ nhật và tròn. Lời cảm ơn Để hoàn thành bản đề tài này, chúng tôi xin được gửi lời cảm ơn đến các thầy cô cùng toàn thể các bạn Sinh viên trong Khoa lý đã trao đổi và động viên chúng tôi. 10
Tài liệu tham khảo Đỗ Đình Thanh, Phương pháp toán lý, NXBGD, 2002. [1] Phan Huy Thiện, Phương trình toán lý, NXBGD, 2007. [2] Phan Huy Thiện, Tuyển tập bài tập phương trình toán lý, NXBGD, 2008. [3] Nguyễn Đình Trí, Tạ Văn Đĩnh, Nguyễn Hồ Quỳnh, Toán cao cấp tập 1. NXBGD, [4] 1997. Nguyễn Đình Trí, Tạ Văn Đĩnh, Nguyễn Hồ Quỳnh, Toán cao cấp tập 2. NXBGD, [5] 1997. Nguyễn Đình Trí- Tạ Văn Đĩnh, Nguyễn Hồ Quỳnh, Toán cao cấp tập 3. NXBGD, [6] 1997. Vũ Thị Kim Liên, Phạm Hữu Kiên, Bài tập toán Cho Vật lý, Thái Nguyên, 2009. [7] Phạm Hữu Kiên, Khoá luận tốt nghiệp Đại học, Thái Nguyên, 2003. [8] Phạm Hữu Kiên, Luận văn tốt nghiệp Cao học, Hà nội, 2006 [9] Đào Văn Phúc, Điện động lực, NXBGD,1984 [10] Đỗ Thị Liên, Khoá luận tốt nghiệp Đại học, Thái Nguyên, 2008 [11] Phạm Đức Long, Đề tài nghiên cứu khoa học, Thái Nguyên, 2009. [12] Nguyễn Thị Thu Hằng, Đề tài nghiên cứu khoa học, Thái Nguyên,2009. [13] Phạm Hữu Kiên, Nguyễn Thị Thu Hằng, Tạp chí khoa học và công nghệ ,số 2(50), năm [14] 2009 trang 43-47. Nguyễn Thi Bảo Ngọc, Dao động sóng, NXBĐHQGHN, 1996. [15] 11