(LUẬN văn THẠC sĩ) giải số phương trình vi phân đại số bằng phương pháp đa bước

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN NGUYỄN THỊ HẢI DUNG GIẢI SỐ PHƯƠNG TRÌNH VI PHÂN ĐẠI SỐ BẰNG PHƯƠNG PHÁP ĐA BƯỚC LUẬN VĂN THẠC SĨ KHOA HỌC Chuyên ngành: Toán ứng dụng Mã số: 60 46 01 12 NGƯỜI HƯỚNG DẪN KHOA HỌC PGS TSKH VŨ HOÀNG LINH Hà Nội - 2014 TIEU LUAN MOI download : skknchat@gmail.com Mục lục Mục lục i Lời nói đầu Lời cảm ơn iii v Giới thiệu 1.1 Phương trình vi phân đại số 1.1.1 Phương trình vi phân đại số số 1.1.2 Phương trình vi phân đại số số cao 1.2 Phương pháp đa bước giải phương trình vi phân 1.2.1 Điều kiện ổn định phương pháp đa bước 1.2.2 Một số phương pháp đa bước cụ thể Phương pháp đa bước giải phương trình vi phân đại số số 2.1 Sự hội tụ phương pháp đa bước 2.2 Một số phương pháp đa bước giải phương trình vi phân đại số 2.2.1 Phương pháp Euler ẩn 2.2.2 Phương pháp BDF 2.3 Thử nghiệm số 1 5 19 19 25 25 27 27 Phương pháp đa bước giải phương trình vi phân đại số số 30 3.1 Sự tồn lời giải số 31 3.2 3.3 Ảnh hưởng nhiễu Sai số địa phương Sự hội tụ phương pháp BDF i TIEU LUAN MOI download : skknchat@gmail.com 33 35 MỤC LỤC 3.3.1 Sai số địa phương 3.3.2 Sự hội tụ BDF 3.4 Phương pháp đa bước tổng quát 3.5 Giải hệ phương trình phi tuyến 3.6 Thử nghiệm số Kết luận Tài liệu tham khảo Phụ lục phép lặp Newton ii TIEU LUAN MOI download : skknchat@gmail.com 35 35 39 42 43 46 47 48 Lời nói đầu Phương trình vi phân thường nghiên cứu từ lâu, lý thuyết phương trình vi phân ẩn, có phương trình vi phân đại số, quan tâm mạnh mẽ vòng 30 năm trở lại Phương trình vi phân đại số tốn đặt khơng chỉnh, có nhiều điểm đặc biệt mà ta khơng thể tìm thấy phương trình vi phân thường Ví dụ ma trận hệ số ma trận suy biến, tồn nghiệm phụ thuộc vào vế phải, , khiến việc nghiên cứu vấn đề định tính giải số phương trình vi phân đại số trở nên phức tạp nhiều so với phương trình vi phân thường Phương trình vi phân đại số có nhiều ứng dụng rộng rãi, chúng mô hệ động lực có ràng buộc, chẳng hạn hệ học, hệ mạch điện, hệ kỹ thuật hóa học, lý thuyết điều khiển, động lực học chất lỏng nhiều lĩnh vực khác Động thái chuyển động đối tượng vật lý thường mơ hình hóa qua hệ phương trình vi phân Nhưng trạng thái hệ thống vật lý chịu số ràng buộc (về vị trí, lượng, ) hạn chế mơ tả phương trình (ràng buộc) đại số Những hệ bao gồm phương trình vi phân phương trình đại số, gọi hệ phương trình vi phân đại số Khái niệm số sử dụng lý thuyết phương trình vi phân đại số để đo độ phức tạp phương trình vi phân đại số phương trình vi phân thường Chỉ số số nguyên không âm, cung cấp thơng tin hữu ích cấu trúc tốn học phức tạp việc phân tích hệ phương trình vi phân đại số Trong luận văn này, chúng tơi trình bày số phương pháp đa bước để giải phương trình vi phân số phương trình vi phân số 2, cụ thể phương pháp Euler ẩn phương pháp BDF k bước Ngoài phần mở đầu, kết luận tài liệu tham khảo, luận văn chia thành ba chương: iii TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc MỤC LỤC Chương 1: Giới thiệu Trình bày phương trình vi phân đại số số phương trình vi phân đại số số cao Trình bày số phương pháp đa bước cụ thể giải phương trình vi phân điều kiện ổn định phương pháp đa bước Chương 2: Phương pháp đa bước giải phương trình vi phân đại số số Trình bày dạng tổng quát phương trình vi phân đại số số 1, phương pháp đa bước áp dụng cho toán, hội tụ tốn nhiễu suy biến Lấy ví dụ minh họa thử nghiệm số Chương 3: Phương pháp đa bước giải phương trình vi phân đại số số Trình bày dạng tổng quát phương trình vi phân đại số số 2, phương pháp đa bước áp dụng cho toán, ảnh hưởng nhiễu, hội tụ phương pháp BDF phương pháp đa bước nói chung Lấy ví dụ minh họa thử nghiệm số iv (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Lời cảm ơn Trước trình bày nội dung luận văn, tơi xin cảm ơn Ban chủ nhiệm khoa Toán - Cơ - Tin học tồn thể thầy giáo, giáo khoa Tốn - Cơ - Tin học, phòng Sau đại học, trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, giảng dạy tận tình tạo điều kiện thuận lợi để tơi hồn thành tốt luận văn Đặc biệt, xin gửi lời cảm ơn chân thành tới thầy giáo PGS.TS Vũ Hoàng Linh, người trực tiếp hướng dẫn tận tình bảo tơi suốt q trình tơi học tập thực luận văn Nhân dịp này, xin cảm ơn gia đình ln ủng hộ động viên suốt thời gian học tập Cuối cùng, xin cảm ơn tất bạn, anh, chị, em lớp cao học Tốn khóa 2010 - 2012 khóa 2011 - 2013 tận tình giúp đỡ động viên tơi q trình học tập Xin chân thành cảm ơn! Hà Nội, ngày 10 tháng 10 năm 2014 Học viên Nguyễn Thị Hải Dung (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Giới thiệu 1.1 1.1.1 Phương trình vi phân đại số Phương trình vi phân đại số số Xét phương trình εz 00 + (z − 1)z + z = (1.1) Ta thay đồng thức   d z3 εz + (z − 1)z = εz + ( − z) dx {z } | 00 :=y vào (1.1), ta có y = −z =: f (y, z), z3 εz = y − ( − z) =: g(y, z) (1.2) Đặt ε = (1.2) ta toán đơn giản y = −z =: f (y, z), z3 = y − ( − z) =: g(y, z) (1.3) Trong việc giải (1.2) không đơn giản (1.3) dễ dàng giải y = −z = (z − 1)z , (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Giới thiệu từ suy ln |z| − z2 = x + C Phương trình (1.3) gọi phương trình vi phân đại số Ta thấy, phương trình vi phân đại số kết hợp phương trình vi phân phương trình đại số 1.1.2 Phương trình vi phân đại số số cao Dạng tổng quát phương trình vi phân đại số phương trình vi phân ẩn F (x, u, u0 ) = 0, (1.4) u : R → Rm lời giải, F : R × Rm × Rm → Rm hàm số, suy biến ∂F ∂u0 Định nghĩa 1.1 Phương trình vi phân đại số (1.4) có số vi phân d = m m số nhỏ vi phân F (u0 , u) = 0, dm F (u0 , u) dF (u0 , u) = 0, , = 0, dx dxm (1.5) cho từ phương trình (1.5) rút hệ phương trình vi phân thường u0 = ϕ(u) Hệ số Xét phương trình vi phân đại số y = f (y, z), (1.6) = g(y, z) (1.7) (khơng có z ) Ta lấy đạo hàm (1.7), thu gy (y, z)f (y, z) + gz (y, z)z = 0, suy z = −gz−1 (y, z).gy (y, z)f (y, z), (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Giới thiệu gz khả nghịch lân cận lời giải Vì tốn (1.6), (1.7) có số vi phân gz khả nghịch Hệ số Xét phương trình vi phân đại số y = f (y, z), (1.8) = g(y) (1.9) Trong z khơng có mặt ràng buộc đại số Lấy đạo hàm (1.9) ta thu "ràng buộc ẩn" = gy (y)f (y, z) (1.10) Nếu gy (y)fz (y, z) khả nghịch lân cận lời giải phương trình (1.8), (1.10) phương trình số Lấy vi phân phương trình (1.10) cho ta phương trình vi phân z, phương trình (1.8), (1.9) phương trình vi phân số Nếu giá trị ban đầu thỏa mãn = g(y0 ) = gy (y0 )f (y0 , z0 ) ta gọi chúng "tương thích" Chỉ trường hợp này, phương trình (1.8) (1.9) có lời giải địa phương Chỉ số nhiễu Quan niệm thứ hai số, giải thích số tiêu chuẩn (đơn vị đo) độ nhạy cảm lời giải nhiễu toán cho trước Định nghĩa 1.2 Phương trình (1.4) có số nhiễu p = m dọc theo lời giải u(x) [0, x), m số nguyên nhỏ cho lời giải u b(x) phương trình có nhiễu F (x, u b, u b0 ) = δ(x), (1.11) tồn [0, x) và có đánh giá  (m−1) kb u(x) − u(x)k ≤ C kb u(0) − u(0)k + max kδ(ξ)k + + max δ (ξ) , 0≤ξ≤x 0≤ξ≤x (1.12) với biểu thức vế phải đủ nhỏ, C số (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Giới thiệu Hệ số Để tính tốn số nhiễu phương trình (1.6), (1.7), ta xét hệ bị nhiễu yb0 = f (b y , zb) + δ1 (x), (1.13) = g(b y , zb) + δ2 (x) (1.14) Ta thấy hiệu zb−z đánh giá nhờ định lý hàm ẩn mà không cần đạo hàm lời giải Vì gz khả nghịch, từ phương trình (1.14), (1.7), ta có kb z (x) − z(x)k ≤ C1 (kb y (x) − y(x)k + kδ2 (x)k) , (1.15) với vế phải (1.15) đủ nhỏ Trừ (1.13) cho (1.6), lấy tích phân từ → x, sử dụng điều kiện Lipschitz cho f ước lượng cho zb(x)−z(x) cho ta e(x) = kb y (x) − y(x)k thỏa mãn Zx Zx e(x) ≤ e(0) + C2 e(t)dt + C3 0 x Z kδ2 (t)k dt + δ (t)dt Trong ước lượng này, lấy chuẩn phần tích phân cho δ2 , phần ngồi tích phân cho δ1 Điều trường hợp nhiễu phương trình đại số (1.7) quan trọng nhiễu phương trình vi phân (1.6) Cuối ta áp dụng Bổ đề Gronwall Rξ Rx kb y (x) − y(x)k ≤ C4 (kb y (0) − y(0)k + kδ2 (t)k dt + max δ1 (t)dt ) 0≤ξ≤x 0 ≤ C5 (kb y (0) − y(0)k + max kδ2 (ξ)k + max kδ1 (ξ)k) 0≤ξ≤x 0≤ξ≤x Bất đẳng thức với bất đẳng thức (1.15) số nhiễu toán Hệ số Xét nhiễu phương trình (1.8), (1.9) sau: yb0 = f (b y , zb) + δ(x), (1.16) (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Phương pháp đa bước giải phương trình vi phân đại số số ước lượng muốn có thành phần y b) Vì zn phụ thuộc vào yn−k , , yn−1 không phụ thuộc vào giá trị z trước nên ta áp dụng Định lý 3.2 với ybj = y(xj ), zbj = z(xj ), δ = O(hp ), θ = Ta suy k CX kzn − z(xn )k ≤ kgy (y (xn )) (yn−j − y (xn−j )) + O (hp )k (3.46) h j=1 Sử dụng (3.44) ynl = y (xn ) + O (hp ) suy n−k+1   P l+1 l − y g (y (x )) y kgy (y (xn )) (yn−j − y (xn−j ))k = y n n−j n−j l=0   n−k+1
 l P  2p+1 l+1 l ) kgy yn yn−j − yn−j + O(h ≤ l=0 = O(hp+1 ) suy kzn − z(xn )k ≤ C6 hp (3.47) c) Tổng quát, số C3 , C5 , C6 phụ thuộc vào C0 , C1 , C2 với giả thiết (3.32) Với p ≥ ta thu gọn cỡ bước h C5 hp−1 ≤ C0 , C3 hp−1 ≤ C1 , C6 hp−1 ≤ C2 lời giải số không vi phạm điều kiện (3.32) khoảng xét Với p = (phương pháp Euler ẩn), thay đánh giá (3.34) k∆zn+k k ≤ C k∆Yn k (3.48) ∆yn+1 − ∆yn = O(h k∆yn k) (3.49) Thay (3.35) ta có số ẩn số hạng O( ) phụ thuộc vào C0 , C1 , C2 h đủ nhỏ Phương pháp tiêu chuẩn cho ta kết hội tụ 3.4 Phương pháp đa bước tổng quát Cho phương pháp đa bước tổng quát (3.3) với đa thức ban đầu ρ(ζ) = k X αi ζ i i=0 39 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Phương pháp đa bước giải phương trình vi phân đại số số σ(ζ) = k X βi ζ i i=0 Ta có kết hội tụ Định lý 3.4 Xét phương trình vi phân đại số số (3.1) thỏa mãn điều kiện (3.2) Giả thiết phương pháp đa bước ổn định ổn định mạnh ∞ (phần tử σ(ζ) thuộc đường tròn đơn vị |ζ|< 1) Nếu có cấp xác p≥ sai số tồn cục thỏa mãn yn − y(xn ) = 0(hp ), zn − z(xn ) = 0(hp ), xn = nh ≤ const, với giá trị ban đầu thỏa mãn (j = 0, , k − 1) yj − y(xj ) = 0(hp+1 ), zj − z(xj ) = 0(hp ) (3.50) Chứng minh Trong (3.31) ta giả thiết kzen − zbn k ≤ C3 h Thay vào (3.32), (3.33) ta có k∆yn+k k ≤ C(k∆yn k + h k∆Zn k , k−1 c X k∆Zn+k k ≤ ( kgy (ym+k )∆yn+j k + h k∆Yn k + h k∆Zn k) h j=0 (3.34), (3.35) trở thành k X i=0 αi ∆yn+1 = h k X βi fz (b y n+k , zbn+k )∆zn+i + 0(h k∆Yn k + h2 k∆Zn k , i=0 = gy (n+k )∆yn+k + 0(h k∆Yn k + h2 k∆Zn k) (3.51) 40 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Phương pháp đa bước giải phương trình vi phân đại số số Sử dụng phép truy hồi cho ∆Zn ta thu kết sau: nhân phương trình (3.51) với ((gy fz )−1 gy )(b yn+k , zbn+k) lấy k k P P h βi ∆Zn+i = αi ((gy fz )−1 gy )(b yn+k , zbn+k )∆yn+i +0(h k∆Yn k+h2 k∆zn k) i=0 i=0 (3.52) với phép chiếu Pn , Qn (3.36) vectơ Un , Vn ta thu (3.39), (3.40) cộng thêm số hạng 0(h2 k 4Zn k) Từ (3.52) ta có h∆Zn+1 = (B ⊗ I)h∆Zn + 0(h kUn k + kVn k + h2 k∆Zn k),  0 −Bk−1 −B1 −B0  0    B=     với Bj = Bj/Bk Trong phần ta sử dụng chuẩn kB ⊗ Ik ≤ K < (do phương pháp ổn định chặt ∞) Do ta có bất đẳng thức      kUn k + 0(h) 0(h) 0(h) kUk   kVn+1 k  ≤  0(h) f + 0(h) 0(h)   kVn k h k∆Zn k 0(h) 0(1) K + 0(h) h k∆Zn+1 k từ k∆yn k ≤ C3 hp+1 , k∆Zn k ≤ C7 (ρn + k n + h)hp , kgy (n+k ∆yn+j k ≤ C4 (ρn + k n + h)hp+1 với j = 0, , k − Cộng sai số lan truyền (3.45) lại ta thu ước lượng cho thành phần y z 41 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Phương pháp đa bước giải phương trình vi phân đại số số 3.5 Giải hệ phương trình phi tuyến phép lặp Newton Hệ phi tuyến (3.3), (3.4) thường giải phép lặp Newtơn có ích cho việc nghiên cứu hội tụ Sử dụng η (3.10) thay đổi h Khi hệ phi tuyến trở thành (bỏ qua số) y − η − hf (y, z) = 0, (3.53) g(y) = Ma trận Jacobian hệ phi tuyến (3.53)   I − hfy −hfz J= gy (3.54) nghịch đảo có dạng   P + 0(h) fz (gy fz )−1 + 0(h) −1 J = −h−1 (gy fz )−1 g(y) + 0(1) −h−1 (gy fz )−1 + 0(1) P = I − fz (gy fz )−1 gy hình chiếu (3.36) Ta xét phương pháp Newtơn phép lặp với hàm   y  y − η − hf (y, z) −1 φ(y, z) = − J0 (3.55) g(y) z Chỉ số J0 cho biết đối số đạo hàm (3.54) b cho lời giải (3.53) Ta ước lượng vài xấp xỉ cố định (b η , ζ) sử dụng kí hiệu {fy }0 thay cho fy (b η , ςb) Tính tốn trực tiếp Φ0 (y, z) ta o  n 

−1 fz (gy fz ) {gy }0 − gy + O (h) h{P }0 fz + O h  n o0 n o 
−1 −1 −1 h (gy fz ) {gy }0 − gy + O (1) (gy fz ) gy ({fz }0 − fz ) + O (h) 0 Nếu ta giả thiết (b η , ςb) xấp xỉ với điểm cố định (3.55) với sai số O (h), ta có điểm cố định 42 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Phương pháp đa bước giải phương trình vi phân đại số số
 O (h) O h2 Φ (y, z) = O (1) O (h) Với ma trận D = diag(I, hI) ta có Dφ (y, z)D−1 = 0(h) Do chuẩn kyk + h kzk ta đạt hệ số h phép lặp Newtơn đơn giản  3.6 Thử nghiệm số Ví dụ 3.1 Giải tốn sau phương pháp Euler ẩn y10 = λy1 − y4 y20 + y30 = (2λ − sin2 t)(y2 + y3 ) + (y2 − y3 )2 = y2 − y3 − 2(sin t)(y1 − 1) = y2 + y3 − 2(y1 − 1)2 , λ tham số, y1 (0) = 2, y2 (0) = Giải Ta đưa tốn phương trình vi phân số cách đặt x1 = y1 , x2 = 12 (y2 + y3 ), z1 = 12 (y2 − y3 ), z2 = y4 Ta có x01 = λx1 − z2

x02 = 2λ − sin2 t x2 + sin2 t (x1 − 1)2 = x2 − (x1 − 1)2 Đây phương trình vi phân đại số số 43 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Phương pháp đa bước giải phương trình vi phân đại số số Dùng phương pháp Euler ẩn ta có bảng kết sau:(lấy h = 0.1) Ví dụ 3.2 Xét tốn Ví dụ 3.1, giải phương pháp BDF bước, bước, bước Kết thử nghiệm số phương pháp BDF bước cho phương trình vi phân đại số số (lấy h = 0.1) 44 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Chương Phương pháp đa bước giải phương trình vi phân đại số số Kết thử nghiệm số phương pháp BDF bước cho phương trình vi phân đại số số (lấy h = 0.1) Kết thử nghiệm số phương pháp BDF bước cho phương trình vi phân đại số số (lấy h = 0.1) 45 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Kết luận Luận văn trình bày giải số phương trình vi phân đại số phương pháp đa bước Việc giải số phương trình vi phân đại số phương pháp đa bước trình bày từ giới thiệu phương trình vi phân đến số phương pháp đa bước cụ thể như: phương pháp Adams, phương pháp Nystrom, phương pháp BDF k bước Tiếp theo, luận văn trình bày khái niệm phương trình vi phân đại số số 1, phương trình vi phân đại số số vấn đề liên quan: hội tụ, sai số, ảnh hưởng nhiễu, Bên cạnh đó, luận văn trình bày số chứng minh chi tiết liên quan đến sai số, ảnh hưởng nhiễu, hội tụ Đặc biệt, luận văn có trình bày số ví dụ, toán minh họa với việc thử nghiệm số tốn (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Tài liệu tham khảo Phạm Kỳ Anh, 2008,Giải tích số, Nhà xuất Đại học quốc gia Hà Nội E Hairer, S P Norsett, G Wanner (1993), Solving Ordinary Differential Equations I: Nonstiff Problems, Springer-Verlag, second edition E Hairer, S P Norsett, G Wanner (1991), Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, Springer-Verlag, second revised edition J.C.Butcher (2003), Numerical Methods for Ordinary Differential Equations, The Uninversity of Auckland, New Zealand, Wiley Publishers Linda R.Petzold, UriM.Ascher (1997), Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations L Shampine, M W Reichelt (1997), The MATLAB ODE Suite W Wasow (1965),Asymptotic expansions for ordinary differential equations, Interscience, John Wiley and Sons, New York 64c (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Phụ lục Chương trình chạy matlab cho phần thử nghiệm số Ví dụ 2.1 Phương pháp Euler ẩn cho phương trình vi phân đại số số 1 c l e ar a l l syms y t z y01 = [ ] ; z01 = [ ] h=0.1; f=y−y∗ z+cos ( t )−sin ( t )+exp ( t ) ∗( sin ( t ) ) ^2+( sin ( t ) ) ^ ; nghiem = [ ; ] ; for i =1:10 T=h∗ i ; 10 % phuong t r i n h e u l e r an 11 g=[y−y01 ( i )−h∗ subs ( f , t ,T) ;− z+y^2−2∗exp (T) ∗ sin (T) −exp (2∗T) ] ; 12 dcx =1; 13 % Phuong phap newton 14 while dcx >0.00001 15 diemdau=nghiem ; 16 b=subs ( g , { y z } , { nghiem ( , ) nghiem ( , ) } ) ; 17 a=[ d i f f ( g , y ) , d i f f ( g , z ) ] ; 18 c=subs ( inv ( a ) , { y z } , { nghiem ( , ) nghiem ( , ) }) ; 19 nghiem=nghiem−c ∗b ; 20 dcx=sqrt ( ( nghiem ( , )−diemdau ( , ) ) ^2+( nghiem ( , )−diemdau ( , ) ) ^2) ; 21 end 22 %−−−−−−−−−−−−−−−−−−−−−− 23 y01 ( i +1)=nghiem ( , ) ; (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Phụ lục 24 z01 ( i +1)=nghiem ( , ) ; end 26 % kq so sanh 27 s s y = [ ] ; s s z = [ ] ; 28 for i =1:11 29 s s y= [ s s y , y01 ( i )−exp ( ( i −1)∗h )−sin ( ( i −1)∗h ) ] ; 30 s s z =[ s s z , z01 ( i )−sin ( ( i −1)∗h ) ^ ] ; 31 end 32 t = : : ; 33 nghcxy=exp ( t )+sin ( t ) ; 34 nghcxz=sin ( t ) ∗ sin ( t ) ; 35 % l a y k e t qua t h e o t h u t u nghiem tim dc y , z nghiem c h i n h xac y , z so sanh nghiem y va so sanh nghiem z 36 [ y01 ’ , z01 ’ , nghcxy ’ , nghcxz ’ , ssy ’ , s s z ’ ] 25 Ví dụ 2.2 Phương pháp BDF bước, bước cho phương trình vi phân đại số số Phương pháp BDF bước c l e ar a l l syms y t z y01 =[1 ] ; z01 =[0 ] ; h=0.1; f=y−y∗ z+cos ( t )−sin ( t )+exp ( t ) ∗( sin ( t ) ) ^2+( sin ( t ) ) ^ ; nghiem = [ ; ] ; for i =2:10 T=h∗ i ; 10 % phuong phap b d f buoc 11 g=[y−4/3∗y01 ( i ) +1/3∗ y01 ( i −1)−2/3∗h∗ subs ( f , t ,T) ;− z+y^2−2∗exp (T) ∗ sin (T)−exp (2∗T) ] ; 12 dcx =1; 13 % Phuong phap newton 14 while dcx >0.00001 15 diemdau=nghiem ; 16 b=subs ( g , { y z } , { nghiem ( , ) nghiem ( , ) } ) ; 17 a=[ d i f f ( g , y ) , d i f f ( g , z ) ] ; 18 c=subs ( inv ( a ) , { y z } , { nghiem ( , ) nghiem ( , ) 49 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Phụ lục 19 20 }) ; nghiem=nghiem−c ∗b ; dcx=sqrt ( ( nghiem ( , )−diemdau ( , ) ) ^2+( nghiem ( , )−diemdau ( , ) ) ^2) ; end 22 %−−−−−−−−−−−−−−−−−−−−−− 23 y01 ( i +1)=nghiem ( , ) ; 24 z01 ( i +1)=nghiem ( , ) ; 25 end 26 % kq so sanh 27 s s y = [ ] ; s s z = [ ] ; 28 for i =1:11 29 s s y= [ s s y , y01 ( i )−exp ( ( i −1)∗h )−sin ( ( i −1)∗h ) ] ; 30 s s z =[ s s z , z01 ( i )−sin ( ( i −1)∗h ) ^ ] ; 31 end 32 t = : : ; 33 nghcxy=exp ( t )+sin ( t ) ; 34 nghcxz=sin ( t ) ∗ sin ( t ) ; 35 % l a y k e t qua t h e o t h u t u nghiem tim dc y , z nghiem c h i n h xac y , z so sanh nghiem y va so sanh nghiem z 36 [ y01 ’ , z01 ’ , nghcxy ’ , nghcxz ’ , ssy ’ , s s z ’ ] 21 Phương pháp BDF bước c l e ar a l l syms y t z y01 =[1 ] ; z01 =[0 ] ; h=0.1; f=y−y∗ z+cos ( t )−sin ( t )+exp ( t ) ∗( sin ( t ) ) ^2+( sin ( t ) ) ^ ; nghiem = [ ; ] ; for i =3:10 T=h∗ i ; 10 % phuong phap b d f buoc 11 g=[y −(18/11) ∗ y01 ( i ) +(9/11) ∗ y01 ( i −1) −(2/11) ∗ y01 ( i −2) −(6/11) ∗h∗ subs ( f , t ,T) ;− z+y^2−2∗exp (T) ∗ sin ( T)−exp (2∗T) ] ; 12 dcx =1; 50 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Phụ lục % Phuong phap newton 14 while dcx >0.00001 15 diemdau=nghiem ; 16 b=subs ( g , { y z } , { nghiem ( , ) nghiem ( , ) } ) ; 17 a=[ d i f f ( g , y ) , d i f f ( g , z ) ] ; 18 c=subs ( inv ( a ) , { y z } , { nghiem ( , ) nghiem ( , ) }) ; 19 nghiem=nghiem−c ∗b ; 20 dcx=sqrt ( ( nghiem ( , )−diemdau ( , ) ) ^2+( nghiem ( , )−diemdau ( , ) ) ^2) ; 21 end 22 %−−−−−−−−−−−−−−−−−−−−−− 23 y01 ( i +1)=nghiem ( , ) ; 24 z01 ( i +1)=nghiem ( , ) ; 25 end 26 % kq so sanh 27 s s y = [ ] ; s s z = [ ] ; 28 for i =1:11 29 s s y= [ s s y , y01 ( i )−exp ( ( i −1)∗h )−sin ( ( i −1)∗h ) ] ; 30 s s z =[ s s z , z01 ( i )−sin ( ( i −1)∗h ) ^ ] ; 31 end 32 t = : : ; 33 nghcxy=exp ( t )+sin ( t ) ; 34 nghcxz=sin ( t ) ∗ sin ( t ) ; 35 % l a y k e t qua t h e o t h u t u nghiem tim dc y , z nghiem c h i n h xac y , z so sanh nghiem y va so sanh nghiem z 36 [ y01 ’ , z01 ’ , nghcxy ’ , nghcxz ’ , ssy ’ , s s z ’ ] 13 Ví dụ 3.1 Phương pháp Euler ẩn cho phương trình vi phân đại số số c l e ar a l l syms x1 x2 z x y1 t r e a l ld =0.5; x01 =2; h=0.1; t =0; x2=(x1 −1) ^ ; 51 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Phụ lục x02=(x01 −1) ^ ; z0 = ; 10 nghiemcantim =[ x01 , x02 , z0 ] ; 11 f =[ l d ∗x1−z ; ( ∗ ld −( sin ( t ) ) ^2) ∗x2+( sin ( t ) ) ^2∗( x1 −1) ^2]; 12 nghiem = [ ; ] ; 13 for i =1:10 14 t=t 0+h∗ i ; 15 % phuong t r i n h e u l e r an 16 g=[x1−x01 ; x2−x02]−h∗ f ; 17 dcx =1; 18 % Phuong phap newton 19 while dcx >0.0001 20 diemdau=nghiem ; 21 b=subs ( g , { x1 z } , { nghiem ( , ) nghiem ( , ) } ) ; 22 a=[ d i f f ( g , x1 ) d i f f ( g , z ) ] ; 23 c=subs ( inv ( a ) , { x1 z } , { nghiem ( , ) nghiem ( , ) }) ; 24 nghiem=nghiem−c ∗b ; 25 dcx=sqrt ( ( nghiem ( , )−diemdau ( , ) ) ^2+( nghiem ( , )−diemdau ( , ) ) ^2) ; 26 dcx=eval ( dcx ) ; 27 end 28 %−−−−−−−−−−−−−−−−−−−−−− 29 nghiem=eval ( nghiem ) ; 30 x01=nghiem ( , ) ; 31 x02=(x01 −1) ^2 ; 32 z1=nghiem ( , ) ; 33 nghiemcantim =[ nghiemcantim ; x01 , x02 , z1 ] ; 34 end 35 nghiemcantim 36 hold on 37 plot ( nghiemcantim ( : , ) , ’ g ’ ) 38 plot ( nghiemcantim ( : , ) ) 39 plot ( nghiemcantim ( : , ) , ’ r ’ ) Ví dụ 3.2 Phương pháp BDF bước, bước, bước cho phương trình vi phân đại số số 52 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc Phụ lục Phương pháp BDF bước % phuong phap buoc cho b a i c l e ar a l l syms x1 x2 z x y1 t r e a l ld =0.1; x01 =[2 1 ] ; for k =1:2 x02 ( k ) =(x01 ( k ) −1) ^ 2; end h=0.1; 10 t =0; 11 z0 = [ ] ; 12 %phuong t r i n h 13 f =[ l d ∗x1−z ; ( ∗ ld −( sin ( t ) ) ^2) ∗x2+( sin ( t ) ) ^2∗( x1 −1) ^2 ; x2−(x1 −1) ^ ] ; 14 nghiem = [ ; ; ] ; 16 for i =2:10 17 T=t 0+h∗ i ; 18 % phuong phap BDF buoc 19 g=[ x1 ; x2 ; ] − ( / ) ∗ [ x01 ( i ) ; x02 ( i ) ; ] + ( / ) ∗ [ x01 ( i −1) ; x02 ( i −1) ; ] − ( / ) ∗h∗ subs ( f , t ,T) ; 20 dcx =1; 21 % Phuong phap newton 22 while dcx >0.0001 23 diemdau=nghiem ; 24 b=subs ( g , { x1 x2 z } , { nghiem ( , ) nghiem ( , ) nghiem ( , ) } ) ; 25 a=[ d i f f ( g , x1 ) d i f f ( g , x2 ) d i f f ( g , z ) ] ; 26 c=subs ( inv ( a ) , { x1 x2 z } , { nghiem ( , ) nghiem ( , ) nghiem ( , ) } ) ; 27 nghiem=nghiem−c ∗b ; 28 dcx=sqrt ( ( nghiem ( , )−diemdau ( , ) ) ^2+( nghiem ( , )−diemdau ( , ) ) ^2+(nghiem ( , )− diemdau ( , ) ) ^2) ; 15 29 30 end 53 (LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc(LUAN.van.THAC.si).giai.so.phuong.trinh.vi.phan.dai.so.bang.phuong.phap.da.buoc TIEU LUAN MOI download : skknchat@gmail.com