Đang tải... (xem toàn văn)
ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA BÁO CÁO BÀI TẬP LỚN PHƯƠNG PHÁP TÍNH Group 4 Giảng viên TS Đậu Thế Phiệt TP HCM 2022 Danh sách thành viên STT Họ Tên MSSV Tiến độ hoàn thành 1 Huỳnh Vă[.]
ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA BÁO CÁO BÀI TẬP LỚN PHƯƠNG PHÁP TÍNH Group Giảng viên: TS Đậu Thế Phiệt TP.HCM - 2022 Danh sách thành viên STT Họ Tên MSSV Tiến độ hoàn thành Huỳnh Văn Quốc Việt 1713947 100% Hoàng Quốc Huy 1810934 100% Đào Văn Tùng 2014996 100% Trần Hữu Hiếu 2013165 100% Ngô Vĩnh Phát 2014080 100% Nguyễn Đăng Hoằng 1913466 100% LỜI NÓI ĐẦU Phương pháp tính mơn học đại cương có tầm quan trọng sinh viên ĐH Bách Khoa nói riêng sinh viên ngành khối khoa học kỹ thuật cơng nghệ nói chung Mục đích mơn học cung cấp đầy đủ nội dung phương pháp tính dùng cho ngành khoa học kỹ thuật Nó giúp sinh viên khối kỹ thuật tiếp thu vấn đề cách nhẹ nhàng trang bị kỹ cho người học tự phát triển khả áp dụng toán học vào tốn thực tế Mơn phương pháp tính bao gồm kiến thức tìm gần nghiệm cho phương trình vi tích phân hàm phức tạp, lý thuyết sai số Cùng với chuẩn đầu ra: Nhắc lại định nghĩa, tính chất, cách tính, vận dụng lý thuyết vào toán áp dụng tốn thực tế , có khả hoạt động nhóm Mục lục I CƠ SỞ LÝ THUYẾT Problem Problem a Cơng thức hình thang b Công thức Simpson 1/3 c Công thức Simpson 8/3 Problem a Phương pháp chụp b Phương pháp sai phân hữu hạn II Code Matlab III Kết 13 Problem 13 Problem 14 Problem 15 I CƠ SỞ LÝ THUYẾT Problem Trường hợp = + −1.5 −0.3 −0.05 + Có giá trị t p(t), lúc ta lập hệ phương trình với ẩn A, B C Trong hệ có cột hàng giá trị ứng với giá trị p(t) −1.5 −0.3 + −1.5 −1.5 −1.5 −0.05 + −0.3 + −0.3 + −0.3 + −0.05 −0.05 () = −0.05 + + + =6 = 4.4 = 3.2 … Để tìm giá trị A, B C ta áp dụng cách chia ma trận, lấy Vector p(t) chia cho ma trận hệ số bên vế trái, ta thu hệ số A, B C Problem a Cơng thức hình thang chia đoạn [a, b] thành n đoạn nhỏ với bước chia ℎ = − Khi a = x0, x1 = x0 + h, …, xk = x0 +kh, …, xn = x0 + nh yk = f(xk), k = 0, 1, …, n Sử dụng cơng thức hình thang cho đoạn [xk, xk+1] ta được: ∫ () =∫ () +∫ () +⋯+∫ −1 ℎ ≈ () ℎ ( + 1)+ ℎ ( + 1)+⋯+ ( −1+ ) b Công thức Simpson 1/3 Tìm giá trị gần tích phân xác định sau: =∫ () Trong ( ) hàm liên tục khả vi đoạn [a;b] Sử dụng công thức Simpson, cụ thể sau: Chia đoạn [a;b] thành đoạn với mốc + 1= ta có giá trị tương ứng: ; = ; 1= ( 1); = ( 2); = ( 3) = Khi ta có được: ℎ − = ( 1+4 ℎ= + 3) (2) Công thức (2) gọi cơng thức Simpson đơn giản Để gia tăng độ xác lấy tích phân khoảng lớn hơn, khoảng (a,b) chia nhỏ thành n bước Việc định trị qua điểm (tương đương với việc có hai khoảng) khoảng địi hỏi phải có số chẳn khoảng Nghĩa khai triển số khoảng dạng n = 2k Ta có: ℎ = − mốc = Giá trị hàm tuong ứng: + ℎ; = [0;2 ] = ( ) = [0;2 ] Ta có phép biến đổi: =∫() = ∫ () +∫ () +⋯+∫ () ℎ 2 −2 ℎ ≈ ( 0+4 ℎ 1+ 2)+ ( 2+4 ( 3+ 4)+⋯+ −2 +4 −1 + ) ℎ ≈ [ + + 4( + + ⋯+ −1) + 2( + + ⋯+ −2)] c Công thức Simpson 8/3 Công thức Simpson 8/3 có dạng sau: (y + 2(y3 + y 6+ +yn−3) + 3(y1 + y + y4 + y 5+ +yn−2 + yn−1) + yn) Problem a Phương pháp chụp Phương pháp Shooting áp dụng để giải tốn phương trình vi phân bậc với điều kiện biên Giả sử, ta có phương trình vi phân bậc với điều kiện biên sau: "( ) = ( ), ( 0) = , ( )= , ≤ ≤ Để giải toán trên, ta đưa phương trình vi phân bậc dạng hệ phương trình vi phân bậc tốn giá trị ban đầu Đặt = → = ” Phương trình vi phân bậc ban đầu trở thành dạng: ′ ′ ′ ′ = = ( ) ( 0)= ′( 0)=? {(0)= Phương pháp chụp yêu cầu ta phải tự đưa dự đoán cho giá trị y’(x0), sau lần lặp giá trị thay đổi theo dự đoán chủ quan người giải Sau đưa dự đốn giá trị, ta áp dụng phương pháp để giải hệ phương vi phân Runge – kutta bậc phương pháp Euler b Phương pháp sai phân hữu hạn Xét toán biên { ( ) "( ) + ( ) ′( ) + ( ) ( ) = ( ), < < ( )= , ( )= Với phương pháp sai phân hữu hạn, chọn số tự nhiên n >0 Chia đoạn [a, b] thành n đoạn điểm chia x = a, xk = x0 +kh, k = 1, 2, …, n-1, xn = b với ℎ = − Tại nút xk, k = 1, 2, …, n-1 bên đoạn [a, b] sử dụng công thức sai phân hướng tâm, ta có: ( +1) −( −1) +1 ′ − 2ℎ ( +1) −2 ( )+( −1 = ( )≈ 2ℎ −1) +1 −2 + −1 = "( ) ≈ 2 ℎ ℎ Thay vào phương trình cho ta +1 −2 + −1 +1 − −1 + + ℎ2 Với 2ℎ ∀ = 1,2,…, −1 = ( Từ điều kiện biên ( ) = , ( ) = , sau biến đổi ta thu hệ phương trình { = 0= ( − ) ℎ2 2ℎ −1 , +( − ℎ2 ) +( + ) ℎ2 2ℎ ), = ( ), =( ), = ( ) = +1 = ,∀ = 1,2,…, − Đây hệ phương trình đại số tuyến tính cấp n-1: AY=B với A ma trận bên 1 1 − + 2 … 2 − 2 − ℎ2 2ℎ + ℎ2 … 2ℎ … … … … ( − [ , ,…, 2 −1 ℎ2 −1 ] −1 −( − … ℎ2 … = 2ℎ ℎ2 ℎ ) ℎ 2ℎ Và = … −2 ( II − −1 −1−( ℎ + 2ℎ ) ) Code Matlab function[]=Project1() disp('>Problem 1Problem 2Problem 3