Bài tập lớn Phương pháp tính 3

Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 v Bài tập lớn Phương pháp tính 3 v v v v v vBài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3 Bài tập lớn Phương pháp tính 3

TỔNG LIÊN ĐOÀN LAO ĐỘNG VIỆT NAM TRƯỜNG ĐẠI HỌC TÔN ĐỨC THẮNG KHOA CÔNG NGHỆ THÔNG TIN

BÀI TẬP LỚN PHƯƠNG PHÁP TÍNH

BÁO CÁO CÂU 3 VÀ 4

Người hướng dẫn: ThS VÕ ĐỨC VĨNH Người thực hiện: THÁI TRUNG TÍN - 51503315

Lớp : 15050301

Khóa : 19

THÀNH PHỐ HỒ CHÍ MINH, NĂM 2017

TỔNG LIÊN ĐOÀN LAO ĐỘNG VIỆT NAM TRƯỜNG ĐẠI HỌC TÔN ĐỨC THẮNG KHOA CÔNG NGHỆ THÔNG TIN

BÀI TẬP LỚN PHƯƠNG PHÁP TÍNH

BÁO CÁO CÂU 3 VÀ 4

Người hướng dẫn: ThS VÕ ĐỨC VĨNH Người thực hiện: THÁI TRUNG TÍN - 51503315

Lớp : 15050301

Khóa : 19

THÀNH PHỐ HỒ CHÍ MINH, NĂM 2017

LỜI CẢM ƠN

Trong suốt quá trình làm bài tập lớn, em đã gặp rất nhiều khó khăn từ cách tiếp cận và trình bày ý tưởng nhưng nhờ có ThS Võ Đức Vĩnh - Khoa Công nghệ thông tin - Trường đại học Tôn Đức Thắng - đã tận tình hướng dẫn đã giúp em nhìn nhận vấn đề cụ thể, tiếp cận đề tài dễ dàng.

Em xin chân thành cảm ơn thầy vì những lời chỉ bảo vô cùng quý báu của thầy đã giúp em có những thu hoạch quý giá để hoàn thành quá trình làm bài tập lớn này.

Bài thu hoạch này được thực hiện trong khoảng thời gian gần

2 tuần Do vậy, không tránh khỏi những thiếu sót là điều chắc chắn, em rất mong nhận được những ý kiến đóng góp quý báu của quý Thầy Cô và các bạn học cùng lớp để kiến thức của em trong lĩnh vực này được hoàn thiện hơn.

Một lần nữa em xin chân thành cảm ơn.

TP.Hồ Chí Minh,ngày 15 tháng 04 năm 2017

Tác giả Thái Trung Tín

ĐỒ ÁN ĐƯỢC HOÀN THÀNH

TẠI TRƯỜNG ĐẠI HỌC TÔN ĐỨC THẮNG

Tôi xin cam đoan đây là sản phẩm đồ án của riêng tôi và được sự hướng dẫn của ThS Võ Đức Vĩnh Các nội dung nghiên cứu, kết quả trong đề tài này là trung thực và chưa công

bố dưới bất kỳ hình thức nào trước đây Những số liệu trong các bảng biểu phục vụ cho việc phân tích, nhận xét, đánh giá được chính tác giả thu thập từ các nguồn khác nhau có ghi rõ trong phần tài liệu tham khảo.

Ngoài ra, trong đồ án còn sử dụng một số nhận xét, đánh giá cũng như số liệu của các tác giả khác, cơ quan tổ chức khác đều có trích dẫn và chú thích nguồn gốc.

Nếu phát hiện có bất kỳ sự gian lận nào tôi xin hoàn toàn chịu trách nhiệm về nội dung đồ án của mình Trường đại học Tôn Đức Thắng không liên quan đến những vi phạm tác quyền, bản quyền do tôi gây ra trong quá trình thực hiện (nếu có).

TP.Hồ Chí Minh,ngày 15 tháng 04 năm 2017

Tác giả Thái Trung Tín

PHẦN XÁC NHẬN, ĐÁNH GIÁ CỦAGIẢNG VIÊN

Phần xác nhận của giáo viên hướng dẫn

TP.Hồ Chí Minh,ngày 15 tháng 04 năm 2017

Thái Trung Tín Phần đánh giá của giáo viên chấm bài

TP.Hồ Chí Minh,ngày 15 tháng 04 năm 2017

Thái Trung Tín

TÓM TẮTTìm hiểu về đạo hàm và tích phân, mô phỏngchuyển động của 1 con lắc đôi bằng phương thứcRunge-Kutta Tìm hiểu về phương trình vi phân thường

và một số cách giải, phương thức Euler lùi, sự ổn địnhcủa phương thức này, có thêm kiến thức về tính hộicủa phương thức Newton

Mục lục

1.1 Phần a 9

1.2 Phần b 10

2 Câu 4 15 2.1 Phần a 15

2.2 Phần b 16

2.3 Phần c 16

2.3.1 i 16

2.3.2 ii 17

2.3.3 iii 18

2.3.4 iv 18

7

8 MỤC LỤC

Hình 1.2: Question3b.m

1.2 PHẦN B 11

Hình 1.3: Question3b.m

Hình 1.4: Đồ thị biểu diễn câu 3b phần 1

1.2 PHẦN B 13

Hình 1.6: question3b2.m

∆h θ 2 (t = 100) 0.05 -1.052068113668714 0.025 -1.052071573182922 0.0125 -1.052071785738407 0.00625 -1.052071798892810 0.0001 -1.052071799764450 Bảng 1.2: Kết quả câu 3b phần 2

hệ số góc = 4.0152 hoặc −0.4853

14 CHƯƠNG 1 CÂU 3

Hình 1.7: Đồ thị câu 3b phần 2

ωi+1(k) = ω(k−1)i+1 −F (ω(k−1)i+1 )

F0(ωi+1(k−1)) = ωi+1(k−1)−ω(k−1)i+1 −ωi−hf (ti+1,ωi+1(k−1))

1−hf y (t i+1 ,ωi+1(k−1))

15

2.3 PHẦN C 17

2.3.2 ii

Hình 2.2: question4cii.m

Hình 2.3: Xanh lá : 790 ; Xanh dương: 646

Đường màu xanh lá cho giá trị đúng với y tại 2000xấp xỉ bằng 1

1 − hλ

22 CHƯƠNG 3 CÁC HÀM BỔ TRỢ

• rk4.m: phương thức Runge-Kutta

Hình 3.2: rk4.m

Chương 4

Tài liệu tham khảo

• Lab 6,7,8 Numerical Analysis , sakai.it.tdt.edu.vn

• https : //vi.wikipedia.org/wiki/EulerM ethod

• https : //en.wikipedia.org/wiki/N umericalStability

23

• rk4.m: phương thức Runge-Kutta

Hình 3. 2: rk4.m