Luận văn: Sử dụng Matlab mô hình hóa Vật lý THPT và giải bài toán ngược trọng lực

Luận văn ứng dụng phần mềm Matlab mô hình hóa bài tập Vật lý THPT và giải bài toán ngược trọng lực, xác định ranh giới phân chia mật độ.

• Trường học: Đại Học Quốc Gia Hà Nội, Trường Đại Học Giáo Dục
• Chuyên ngành: Sư Phạm Vật Lý
• Đề tài: Sử Dụng Phần Mềm Matlab Mô Hình Hóa Một Số Bài Tập Vật Lý THPT Và Giải Bài Toán Ngược Trọng Lực Xác Định Độ Sâu Ranh Giới Phân Chia Mật Độ
• Loại tài liệu: Khóa Luận Tốt Nghiệp
• Năm xuất bản: 2018
• Địa điểm: Hà Nội

Tóm tắt

I. Giới thiệu về Matlab trong Mô hình hóa Vật lý THPT

Matlab (Matrix Laboratory) là một phần mềm mạnh mẽ được phát triển bởi công ty MathWorks, chuyên biệt trong mô hình hóa vật lý và tính toán khoa học kỹ thuật. Phần mềm này được sử dụng rộng rãi tại các trường đại học, viện nghiên cứu trên toàn thế giới. Với thư viện hàm phong phú, Matlab cho phép giải quyết các bài toán phức tạp với đối tượng khác nhau như đại lượng vô hướng, vector, hệ phương trình tuyến tính và phi tuyến. Khả năng xử lý dữ liệu và biểu diễn đồ họa trong không gian 2D và 3D giúp người dùng quan sát kết quả một cách trực quan. Matlab có thể chạy trên nhiều hệ điều hành và loại máy tính khác nhau, từ máy vi tính đến siêu máy tính. Tính năng tích hợp với các ngôn ngữ lập trình như C, C++, Fortran, Java làm cho Matlab trở thành công cụ lý tưởng cho giáo dục Vật lý hiện đại.

1.1. Các ứng dụng chính của Matlab trong giảng dạy

Matlab được ứng dụng rộng rãi trong mô phỏng các hiện tượng vật lý từ chuyển động ném xiên, dao động điều hòa đến thương bắn dẫn. Phần mềm hỗ trợ giáo viên trực quan hóa các bài tập phức tạp thông qua đồ thị sinh động. Các chương trình tính toán trong Matlab cho phép thay đổi dễ dàng các tham số đầu vào, giúp học sinh hiểu rõ hơn mối quan hệ giữa các đại lượng vật lý và kết quả của chúng.

1.2. Ưu điểm của Matlab so với phần mềm khác

Khác với các phần mềm thương mại cứng nhắc, Matlab cho phép người dùng chủ động xây dựng các chương trình tính toán theo nhu cầu riêng. Tính linh hoạt này giúp tối ưu hóa quá trình dạy học, vượt qua những giới hạn của phần mềm có sẵn. Matlab cung cấp môi trường làm việc thân thiện với không gian làm việc (workspace), cửa sổ lệnh (command window) và trình soạn thảo (editor) tích hợp.

II. Mô hình hóa các bài tập Vật lý THPT bằng Matlab

Mô hình hóa vật lý là quy trình chuyển đổi các hiện tượng tự nhiên thành các mô hình toán học có thể mô phỏng trên máy tính. Matlab cung cấp các công cụ mạnh mẽ để mô phỏng chuyển động của vật ném xiên, một trong những bài toán kinh điển của Vật lý THPT. Thông qua biểu diễn đồ họa, học sinh có thể quan sát quỹ đạo chuyển động, tốc độ và gia tốc tại các thời điểm khác nhau. Dao động điều hòa là một chủ đề khác được mô hình hóa hiệu quả, cho phép học sinh thấy được mối liên hệ giữa vị trí, vận tốc và gia tốc. Ngoài ra, Matlab còn được dùng để mô phỏng thương bắn dẫn và phương trình mắt súng, giúp học sinh hiểu sâu hơn về cơ học cổ điển. Các chương trình mô hình hóa có thể được sửa đổi dễ dàng để khám phá tác động của các thông số khác nhau đến kết quả.

2.1. Chuyển động ném xiên và phân tích quỹ đạo

Chuyển động ném xiên được mô hình hóa bằng cách giải hệ phương trình vi phân mô tả quỹ đạo. Matlab vẽ đồ thị quỹ đạo chi tiết, hiển thị vị trí x, y theo thời gian. Người dùng có thể thay đổi góc ném, vận tốc ban đầu để thấy ảnh hưởng đến tầm xa và độ cao cực đại. Điều này giúp học sinh nắm vững các khái niệm về chuyển động projectile.

2.2. Dao động điều hòa và dao động tắt dần

Dao động điều hòa được mô phỏng thông qua phương trình vi phân mô tả hệ thống lò xo-khối lượng. Matlab hiển thị đồ thị li độ, vận tốc, gia tốc theo thời gian. Dao động tắt dần cũng được mô hình hóa để cho thấy tác động của lực cản. Các thí nghiệm ảo này giúp học sinh hiểu rõ các định luật dao động.

III. Bài toán ngược trọng lực và Phép biến đổi Fourier

Bài toán ngược trọng lực là một ứng dụng nâng cao của Matlab trong việc xác định ranh giới phân chia mật độ dựa trên dữ liệu trường trọng lực. Phương pháp này sử dụng phép biến đổi Fourier để chuyển đổi dữ liệu từ miền không gian sang miền tần số. Phép biến đổi Fourier là một công cụ toán học quan trọng trong xử lý tín hiệu, cho phép phân tích cấu trúc tần số của dữ liệu trọng lực. Tính chất đối xứng và tính chất tuyến tính của phép biến đổi Fourier được ứng dụng để xử lý dữ liệu hiệu quả. Sự chuyển dịch trong miền tần số cũng được khai thác để xác định vị trí của các ranh giới mật độ. Phương pháp Cranser và phương pháp Gauss-Boussinesq được sử dụng để giải bài toán ngược, tìm ra chiều sâu của các ranh giới dựa trên dị thường trọng lực quan sát được.

3.1. Phép biến đổi Fourier và các tính chất cơ bản

Phép biến đổi Fourier chuyển đổi dữ liệu từ miền không gian sang miền tần số, giúp xác định các thành phần tần số trong dữ liệu. Tính chất đối xứng đảm bảo các phần tử thực có đối xứng liên hợp. Tính chất tuyến tính cho phép xử lý tổng hợp các tín hiệu. Trong Matlab, hàm fft() thực hiện biến đổi Fourier nhanh, cho phép xử lý dữ liệu lớn hiệu quả.

3.2. Phương pháp Cranser và Gauss Boussinesq

Phương pháp Cranser sử dụng mối quan hệ giữa dị thường trọng lực và mô hình mật độ để xác định ranh giới. Phương pháp Gauss-Boussinesq cải tiến phương pháp này bằng cách tính đến ảnh hưởng của các lớp sâu. Matlab giải các phương trình tích phân này để tìm độ sâu ranh giới tối ưu.

IV. Ứng dụng thực tế và Kết quả mô hình hóa

Các mô hình được xây dựng trong Matlab không chỉ phục vụ mục đích giáo dục mà còn có ứng dụng thực tế trong địa vật lý và kỹ thuật địa chất. Khu vực San Jacinto, Mỹ đã được sử dụng làm trường hợp nghiên cứu để xác định ranh giới phân chia mật độ ở độ sâu khoảng 30-40 km. Dữ liệu trọng lực đo đạc được xử lý bằng phép biến đổi Fourier và các phương pháp Cranser-Gauss-Boussinesq để tìm ra cấu trúc địa chất sâu. Khu vực An Giang, Việt Nam là một ví dụ khác, nơi mô hình hóa được áp dụng để nghiên cứu cấu trúc địa chất của thềm lục địa Sông Cửu Long. Các kết quả tính toán cho thấy sự thay đổi mật độ địa chất tại các độ sâu khác nhau, hỗ trợ công tác thăm dò tài nguyên và nghiên cứu địa học. Thành công của các ứng dụng này chứng minh tính hữu ích và hiệu quả của Matlab trong giải quyết các bài toán ngược.

4.1. Thông số mô hình và cấu hình dữ liệu

Thông số mô hình bao gồm mật độ các lớp địa chất, độ sâu ranh giới, kích thước lưới tính toán. Dữ liệu trọng lực quan sát được chuẩn hóa và loại bỏ nhiễu trước khi xử lý. Matlab cho phép nhập dữ liệu từ file và lưu kết quả dễ dàng. Các thông số được tối ưu hóa để phù hợp với dữ liệu thực tế.

4.2. Áp dụng thực tế tại San Jacinto và An Giang

Tại khu vực San Jacinto, mô hình cho thấy ranh giới Mohorovičić ở độ sâu khoảng 32 km. Ở An Giang, Việt Nam, mô hình hóa xác định các lớp mẫu cột nước và mật độ trầm tích, hỗ trợ khoan thăm dò dầu khí. Các ứng dụng thực tế này chứng minh giá trị của Matlab trong công tác địa vật lý ứng dụng.

Trích đoạn nội dung tài liệu

Chương 1: Matlab và ứng dụng trong việc mô hình hóa một số bởi tập Vật dý THPT, ~_ Chương 2: Xác định ranh giới phân chia mật độ -_ Chương 3: Mỏ hành hóa và áp đụng. DANH MUC HINA ANH Tình 1.1: Ném xiên trơng trường hợp góc nêm œ 30° Hình 1.2: Ném xiên trơng trường hợp góc tiêm œ 45° Tình 1.3: Ném xiên trong trường hợp góc tiêm ø 60° Hình 1.4: Dao déng điều hỏa Hình 1.5: Dao dông tắt dan Hình 1.6 Trường trọng lục gây bởi quả cầu Tình 1.7: Bản đỏ c đường đồng mức trường trọng lực gầy bởi quả cầu.8: Hình ảnh mặt sóng Tình 31: Mô hình ranh giới phâu chứa mật độ Tình 3.1: Độ sâu mồ hình Hình 3.2: Kết quả xáo định đi thường trọng lực gây bởi ranh giới phân chia mật độ trên Hình 3.3: Kết quả xác định độ sâu tính toàn theo mô hình Tình 3.4: Sự thay đổi của mật độ theo độ sâu Tình 3.5: Dị thường trọng lực Tình 3.6: Kết quả xác định độ sâu Hình 3.7: Kết quả dị thường trọng lực Hình 3.8: Sự thay đối của mật độ theo độ sâu Hình 3.9: Dị thường trọng lực khu vực (.uong Phuoe Toạn và Dang Van Tel [19]) Tình 3.10: Kết quả xác định độ sâu Tình 3.11: Kết quả dị thường trọng lực DANH MUC HINA ANH Tình 1.1: Ném xiên trơng trường hợp góc nêm œ 30° Hình 1.2: Ném xiên trơng trường hợp góc tiêm œ 45° Tình 1.3: Ném xiên trong trường hợp góc tiêm ø 60° Hình 1.4: Dao déng điều hỏa Hình 1.5: Dao dông tắt dan Hình 1.6 Trường trọng lục gây bởi quả cầu Tình 1.7: Bản đỏ c đường đồng mức trường trọng lực gầy bởi quả cầu.8: Hình ảnh mặt sóng Tình 31: Mô hình ranh giới phâu chứa mật độ Tình 3.1: Độ sâu mồ hình Hình 3.2: Kết quả xáo định đi thường trọng lực gây bởi ranh giới phân chia mật độ trên Hình 3.3: Kết quả xác định độ sâu tính toàn theo mô hình Tình 3.4: Sự thay đổi của mật độ theo độ sâu Tình 3.5: Dị thường trọng lực Tình 3.6: Kết quả xác định độ sâu Hình 3.7: Kết quả dị thường trọng lực Hình 3.8: Sự thay đối của mật độ theo độ sâu Hình 3.9: Dị thường trọng lực khu vực (.uong Phuoe Toạn và Dang Van Tel [19]) Tình 3.10: Kết quả xác định độ sâu Tình 3.11: Kết quả dị thường trọng lực (Figure) dé hién thi cac két qua dang dé thi, hinh ảnh, một cửa số soạn tháo (Editor) để soạn thảo chương trình, các thu viện ham ngoai Đề đóng Matlab: kích chuột trải vào [x] ở cửa số lệnh Matlab. Làm việc ữ chế độ hội thoai: Trong chế độ hội thoại MalLab cho phép ta nhập trực tiếp các lệnh từ bàn phím và nhận kết quả xử lý các lệnh này. Chế dộ hội thoại chỉ nên sử dụng khi giải các bài toán nhỏ, cỏ cấu trúc đơn giấn, chí sử dụng 1 lần.

Mỗi lần Matlab chỉ xứ lý từng động lệnh nhập mội + Lại dẫu mời >> của MatLab trong cửa số lệnh ta nhập dòng lệnh cần xứ lý, rồi bam Enter dé thuc hiện dong lệnh này. Trên một déng lệnh có thể nhập nhiều lênh, các lệnh cách nhau bồi đấu phẩy <., > hoặc dân chấm phảy < >. Với một lệnh. đài hoặc dòng lệnh dải, để thông bảo cho MatLab biét déng hiện thời cén tiếp tục xuống dong dudi, ta đùng 3 đâu châm liên nhau < bam Enter để xuống dang, rối lại gõ tiếp +Khi gặp lệnh thực hiện được, nêu kết quả có dạng số, dạng ký tự hoặc biểu thức toán va sau lệnh này không có đâu châm phầy Ahi MatLab sé dura két qua va @ offa 36 lénh: niéu két qua 14 hinh ãnh, nó sẽ được đua kết quả ra cửa số để hoạ.

“khi gặp lệnh không thực hiện được, matlab sẽ thông báo lỗi ở cửa sẽ lệnh và đừng lại ở lệnh này. Làm việc ö chế độ tập trình: Các lệnh được soạn thảo trong các tệp văn bản có phần mở rộng là <. Các tập này ứm-ile) dược chia làm 2 loại: Các tếp lệnh (m-seripL file) và các tộp ham (m-function file), Cé thể chạy các tệp nảy từ cứa số lệnh hoặc từ cửa số soạn thio. Các ham nếu có biển vào thì chỉ có thể chạy được từ cửa số lệnh.

Cúc toắn từ cho ma trận --Ö : cộng ma trận A và l3 (2 ma trận cùng kich thuớc) * A-B: trừ ma trận A và B (2 ma trên cứng kích thuớe) * A*B: nhân rna trận (số cột của Á bằng số hàng của B) + A.*D: nhân từng phântử cia A va B (A, I3 củng kích thước) SƠ ĐÒ KHÓI CAU TRÚC LỆNH for <Biến for>=<Ma trận for> < Nhóm lệnh cẳn lặp > Lệnh lặp While.end SƠ ĐÒ KHÓI CAU TRUC LENH while <Biéu thie logic> LL7. Dohoa Vẽ đỗ thị hảm số y=f(x) Plot(x,y, ‘mau sic’) Vẽ đỗ thị hảm 2 bien U=f(x,y) dưới dang các đường đồng mức vả mặt cong, Contour(x,y,U,n) (n là số đường đồng mức) 'Mesh(x,y,U) hoặc Surf(x,y,U) 1. Mỡ và lưu file số liệu Nd file: Tniportdaita Vi dur: sl = importdata('data.dat') DANH MUC HINA ANH Tình 1.1: Ném xiên trơng trường hợp góc nêm œ 30° Hình 1.2: Ném xiên trơng trường hợp góc tiêm œ 45° Tình 1.3: Ném xiên trong trường hợp góc tiêm ø 60° Hình 1.4: Dao déng điều hỏa Hình 1.5: Dao dông tắt dan Hình 1.6 Trường trọng lục gây bởi quả cầu Tình 1.7: Bản đỏ c đường đồng mức trường trọng lực gầy bởi quả cầu.8: Hình ảnh mặt sóng Tình 31: Mô hình ranh giới phâu chứa mật độ Tình 3.1: Độ sâu mồ hình Hình 3.2: Kết quả xáo định đi thường trọng lực gây bởi ranh giới phân chia mật độ trên Hình 3.3: Kết quả xác định độ sâu tính toàn theo mô hình Tình 3.4: Sự thay đổi của mật độ theo độ sâu Tình 3.5: Dị thường trọng lực Tình 3.6: Kết quả xác định độ sâu Hình 3.7: Kết quả dị thường trọng lực Hình 3.8: Sự thay đối của mật độ theo độ sâu Hình 3.9: Dị thường trọng lực khu vực (.uong Phuoe Toạn và Dang Van Tel [19]) Tình 3.10: Kết quả xác định độ sâu Tình 3.11: Kết quả dị thường trọng lực CHUONG 1: MATLAB VA UNG DUNG TRONG VIEC M6 HINA HOA MỘT SỐ HÀI FẬP VẬT LÝ THPT 1. Gisi thiéu phan mém Matlab.

Matlab (matrrix laboralory) là mội. bộ chương trình phần tiểm lớn của công, ty MathWork, lẻ ngôn ngữ có tính năng cao trong tính toán khoa học, kỹ thuật và được sử đụng rộng rãi tại nhiều trường, đại học, viện nghiền cứu trong nước và trên thể giới. Matlab vo thé chay dược dưới nhiều hệ điều hành, trên nhiều loại máy tính, từ các máy vị tính đến cáo siêu máy tín. Với thư viện các hàm phong phủ, Miatlab cho phép giải quyết các bài toán với các đối tượng khác nhau là đại lượng vô tưởng, veolor, các hệ phương trình tuyển lính, phí tuyến, hay các xâu kỉ Lự,.

với kết quả nhanh chóng vả chính xác. Matlab cho phép xử lý đữ liệu, biêu điễn đỗ họa một cách linh hoạt, đơn giân và chính xác trang không gian hai chiêu cũng như ba chiêu, giúp người sử dụng có thể quan sát kết quả một cách trực quan và dưa ra giải pháp tốt nhất. Do được tích hợp cùng với một số ngôn ngữ lập trình thông dụng khác rằn: Ơ, ƠI 1, Fortram, Java. rên những ứng dụng của Matlab co thé dé dang chuyển dỏi sang các ngôn ngữ dỏ.

Mãi trường làm việc trong Matlab. Matlab là một ngôn ngữ thông địch khi làm việc trong môi trường của nó. Matlab thực hiện việc quãn ly tép én va tép ham (cae m-file) ở thự mục hiện hành. và các thư mục trên đường dan path Cac thư mục nảy thuộc miễn quản lý của Matlab.

Chúng ta có thể bổ sung thôi áo Thụ trục mới vảo miễn quên lý của Matlab bằng addpath [1] Khởi động tà đồng Matlab. Để khối dong; ta kich đúp chuột trái vào biểu tượng Matlab. Môi trường Matlab gồm các phân chính: Không gian làm việc (work space) để lưu trữ các biển trong qué. tinh Liuh loàn;, một cửa số lệnh (conunand wirrdow) để nhập lệnh vả hiển thị kết quả dạng số, xâu ký tự và biểu thức, các cửa sỏ dỗ hoạ Luu file: Save Vi du: save(‘data.dal’,’sI','-asen') 1.

Sử dụng phan mém Matlab mã hình hỏa một số bài tập Vật lý THPT. Phần mềm Matlab dược các nhà Vật lý sử dụng như mội. công cụ tiên ích cho khảo sát, và xử lý số liệu Vật lý. Trong chương trinh Vật lý phổ thông, hầu hết các tài toán đều có thể mô hình hóa hoặc mô phông nhưng trong phạu vị của khóa hiận tôi chọn những bài toán diễn hình nhất dễ mô tả thông qua việc xây dựng các chương trinh cho máy tính.

Công việc này cũng đồng thời tạo tiền để, cơ sở cho những nghiên cứu ở chương sau 1. Chuyễn động của vật ném xiên. Bài toán Từ một đỉnh tháp cao H người ta ném một hèn đá xuống đất với vận tốc v théo phương hợp với mặt phẳng ngàng một góc ơ. Bồ qua sức cân của không khí Với H, v, ø không bị giới hạn.

Trong bải toán dưới đây, chủng tôi chọn vận tốc ban. đầu là 25 mưa, vị trí ban đầu theo phuống ngang là Ô mụ vị trí bạn đầu theo phương thẳng đứng là 40 m, góc ném là ơ thay dối. Chương trình Chúng tôi xây dụng chương trình mô phỏng hiện tượng ném xiên của vật nặng với các thông số như trên. Tuy nhiên, chúng tôi lưuý rằng tất cá cáo thông số đầu vào của vật năng có thể thay đổi thông qua việc chỉnh sửa code chương trình.

Qua trình nêm xiên được mô phỏng trực quan và sinh động như kết quã rữnư hình vẽ tbên dưới. 6 MỦ ĐẦU Hiện nay, giáo dục dang đổi mới mạnh mẽ (heo đà phát triển của xã hội, đồi hỏi các nhà giáo dục phải đổi mới không những nội dung mà cén cả vẻ phương pháp giảng dạy. Với sự phát triển không ngừng của khoa học công nghệ, máy tỉnh đã trở thành một công cụ dắc lực giúp giải các bài tập cũng như mồ phóng, mô bình hóa các hiện tượng Vật lý, đem lại hiệu quả cao trong giảng đạy ỏ các trường Đại học, Cao đẳng, các trường Trung học phố thông (THPT). Nhiều phần mềm thương mại dä dược thiết kế, xây dựng phục vụ cong tac ging dạy.

Tuy nhiên, các phần mềm này thường cứng nhắc, khó thay đổi được nội đưng, các thông số đầu vào hoặc tiểu có thấy đối dược thì cũng bị giới hạn trong một phạm vì nhất định. Hạn chế đỏ có thể được khắc phục bằng việc chủ động xây dung các chương trình tính toán, mỗ phỏng, mê hình hóa trên cáo phần mềm, ngôn ngữ lập trình. Một trong những công cụ mạnh mẽ dược tôi lựa chọn là phần mềm Matlab.