GiẢI PHƯƠNG TRÌNH PHI TUYẾN

CHƯƠNG 4 GiẢI PHƯƠNG TRÌNH PHI TUYẾN Nội dung Đặt vấn đề 1. Phương pháp chia đôi 2. Phương pháp dây cung 3. Phương pháp Newton 4. Phương pháp cát tuyến 5. Phương pháp lặp 6. Phương pháp Bairstow Đặt vấn đề • Phương trình phi tuyến (PTPT) – VD1: x 2 = 0 – VD2: 1 + 2x + x 2 - 3x 3 + 7x 4 = 0 – VD3: ln(x+1) = 0 – VD4: tg(x) – artg(2x) = 0 – Tổng quát: f(x) = 0 • Giải phương trình phi tuyến (root finding) – Tìm x để f(x) = 0 – X được gọi là nghiệm của PT, cũng được gọi là không điểm của hàm f • Tìm nghiệm dưới dạng công thức hiện: Khó, một số không tồn tại ( VD PT đa thức bậc lớn hơn 4) => sử dụng PP số dựa trên thủ tục lặp Giải PTPT: Một số khái niệm (1) • Sự tồn tại nghiệm – Định lý: Cho hàm f:R->R; [a,b] là đoạn phân ly nghiệm nếu f(a) và f(b) trái dấu. Nếu thêm điều kiện f liên tục trên [a,b] thì tồn tại nghiệm x * ϵ [a,b] sao cho f(x * )=0. – VD: e x + 1 = 0 vô nghiệm 2x + 3 = 0 có một nghiệm x 2 + 3x + 1 = 0 có hai nghiệm sin(x) = 0 có vô số nghiệm • Độ nhạy và điều kiện của bài toán giải PTPT – Số điều kiện của bài toán tìm nghiệm x * : )( 1 *' xf Giải PTPT: Một số khái niệm (2) • Giải PTPT bằng phương pháp lặp – Điều kiện dừng • • ɛ là độ chính xác cho trước – Tốc độ hội tụ: • Gọi sai số ở bước lặp k là e k = x k - x * ; x k là lời giải xấp xỉ tại bước k, x * là nghiệm chính xác. • Dãy {e k } hội tụ với tốc độ r nếu: C ≠ 0 – r = 1: hội tụ tuyến tính – r > 1: hội tụ trên tuyến tính – r = 2: hội tụ bình phương  )(xf   xx * ; 1 lim C e e r k k k    Giải PTPT: Phương pháp chia đôi (1) • Ý tưởng: nếu [a,c] chỉ chứa một nghiệm của PT f(x)=0 thì f(a)*f(c)≤0; [a,c]-khoảng phân ly nghiệm • Phương pháp chia đôi: Chia đôi khoảng phân ly nghiệm liên tục cho đến khi đủ nhỏ, như sau: – Chia đôi: b = (a+c)/2 – Kiểm tra: • Nếu f(b) = 0, => b là nghiệm • Nếu f(a)*f(b)≤0 thì đặt [a,b] là khoảng phân ly nghiệm mới • Nếu f(c)*f(b)≤0 thì đặt [b,c] là khoảng phân ly nghiệm mới – Lặp cho đến khi khoảng phân ly nghiệm nhỏ hơn độ chính xác ɛ cho trước Giải PTPT: Phương pháp chia đôi (2) • Độ dài khoảng phân ly nghiệm sau mỗi bước lặp: – Bước 1: (c-a)/2 1 – Bước 2: (c-a)/2 2 – Bước n: (c-a)/2 n • Cho trước độ chính xác ɛ, thì số bước lặp cần thiết là số nguyên n thỏa mãn: • Vậy số bước lặp cần thiết là: 7    n ac 2  ac n   2 log          ac n 2 log Giải PTPT: Phương pháp chia đôi (3) • VD: PT e x – 2 = 0 có nghiệm nằm trong khoảng [0,2]. Tìm nghiệm với sai số cho phép 0.01 – Đặt a = 0, c = 2, => f(a)*f(c) = -1*5.389 < 0 – Bước lặp 1: • Đặt b = (2+0)/2 = 1; f(b) = 0.718 • Kiểm tra: f(a)*f(b) < 0, => [0,1] là khoảng phân ly nghiệm mới – Bước lặp 2: • Đặt b = (1+0)/2 = 0.5; f(b) = - 0.351 • Kiểm tra: f(b)*f(c) < 0, => [0.5,1] là khoảng phân ly nghiệm mới – ………. Giải PTPT: Phương pháp chia đôi (4) 24-Nov-13 9 Lần lặp a b c f(a) f(b) f(c) Sai số (độ dài khoảng PLN) 1 0 1 2 -1 0.718 5.3890 1 2 0 0.5 1 -1 -0.351 0.718 0.5 3 0.5 0.75 1 -0.351 0.117 0.718 0.25 4 0.5 0.625 0.75 -0.351 -0.132 0.117 0.125 5 0.625 0.688 0.75 -0.132 -0.011 0.117 0.0625 6 0.688 0.719 0.75 -0.011 0.058 0.117 0.03125 7 0.688 0.703 0.719 -0.011 0.020 0.052 0.015625 8 0.688 0.695 0.703 -0.011 0.004 0.020 0.0078125 • Ghi chú: số bước lặp:   8200log 01.0 02 log 22          n • Yêu cầu và tính năng: – Yêu cầu phải biết trước khoảng phân ly nghiệm – Không đòi hỏi tính liên tục của đạo hàm bậc nhất – Có thể giải kiểu PTPT bất kỳ – Có thể áp dụng cho hàm không biểu diễn dưới dạng giải tích 24-Nov-13 10 Giải PTPT: Phương pháp chia đôi (5) [...].. .Giải PTPT: Phương pháp chia đôi (6) • Bài tập: Viết chương trình Matlab giải phương trình phi tuyến bằng phương pháp chia đôi 24-Nov-13 11 Giải PTPT: Phương pháp dây cung (1) • Thay vì chia đôi khoảng phân ly nghiệm, phương pháp dây cung sử dụng đoạn thẳng đi qua hai đầu mút của khoảng phân ly nghiệm để tìm... phương Giải PTPT: Phương pháp Newton (6) • VD: Giải PT sau: f(x) = x2 – 4 sin(x) = 0: – Ta có: f ’(x) = 2x – 4 cos(x) – Suy ra công thức lặp Newton: xn  4 sin( xn )  xn  2 xn  4 cos( xn ) 2 xn 1 – Lấy x0 = 3, ta có kết quả như bảng sau: Bước lặp 0 1 2 x 3 2.153 1.954 f(x) 8.346 1.295 0.108 Giải PTPT: Phương pháp Newton (7) • Bài tập: Viết chương trình Matlab giải PTPT bằng phương pháp Newton Giải. .. = 0 khi đó ta sẽ có: R(x) = 0 Giải PTPT: Phương pháp Bairstow (4) – Thay (3) và (4) vào (2) và viết phương trình thu được theo dạng chuỗi lũy thừa Bởi vì phương trình này phải bằng phương trình (1), nên bằng cách cân bằng hệ số ta có: aN = bN aN-1 = bN-1 + p*bN aN-2 = bN-2 + p*bN-1 + q*bN ………………… a2 = b2 + p*b3 + q*b4 a1 = b1 + p*b2 + q*b3 a0 = b0 + q*b2 (6) Giải PTPT: Phương pháp Bairstow (5) – Viết... f (a ) Giải PTPT: Phương pháp dây cung (2) y A(a,f(a)) b2 b1 c a B(c,f(c)) x Giải PTPT: Phương pháp dây cung (3) • Khác so với phương pháp chia đôi: – Không đặt b=(c+a)/2 – Đặt: af (c)  cf (a ) b f (c)  f (a ) Giải PTPT: Phương pháp dây cung (4) • Yêu cầu và tính năng: – Yêu cầu phải biết trước khoảng phân ly nghiệm – Có thể giải kiểu PTPT bất kỳ – Hội tụ nhanh nếu hàm có dạng phép nội suy tuyến. .. PTPT: Phương pháp cát tuyến • Ý tưởng: Thay việc tính đạo hàm trong phương pháp Newton bằng việc tính sai phân xấp xỉ dựa trên hai bước lặp liên tiếp • Phương pháp Newton f ( xk 1 ) xk  xk 1  , k  1,2, f ' ( xk 1 ) • Phương pháp cát tuyến: f ( xk 1 )  f ( xk 2 ) x  x  f ( xk 1 ) , k  1,2, Sk 1  f ( xk 1 )  ; k k 1 Sk 1 xk 1  xk 2 ' – Cần hai điểm xuất phát: x0 và x1 Giải PTPT: Phương. .. điểm bất động Giải PTPT: Phương pháp lặp (2) • Các ví dụ: – Phương pháp Newton, vì f ( xk 1 ) xk  xk 1  , k  1,2, f ' ( xk 1 ) f ( x) , ta được phương pháp lặp Nên có thể đặt g ( x )  x  f ' ( x) – Tìm nghiệm của PT: f(x) = x - ex , => g(x) = ex – Tìm nghiệm của PT: f(x) = x2 – x – 2, => g(x) = x2 – 2 • Công thức giải PTPT bằng phương pháp lặp xk = g(xk-1); k = 1, 2, … Giải PTPT: Phương pháp... ( x0 )  ( x  x0 ) f ' ( x0 ) • Vậy giải PT f(x)=0  giải PT f ( x0 ) f ( x0 )  ( x  x0 ) f ' ( x0 )  0  x  x0  f ' ( x0 ) Giải PTPT: Phương pháp Newton (4) • Thủ tục lặp để giải PTPT bẳng phương pháp Newton: – Chọn nghiệm xấp xỉ x0 – Tìm nghiệm theo công thức lặp f ( xk 1 ) xk  xk 1  , k  1,2, f ' ( xk 1 ) – Kết thúc khi: f ( xk )   Giải PTPT: Phương pháp Newton (5) • Nhận xét: – Đòi... Hội tụ nhanh nếu hàm có dạng phép nội suy tuyến tính; hội tụ chậm nếu khoảng phân ly nghiệm lớn 24-Nov-13 15 Giải PTPT: Phương pháp Newton (1) • Ý tưởng: – Thay PTPT f(x) = 0 bằng một phương trình tuyến tính với x – Yêu cầu biết nghiệm xấp xỉ ban đầu – Dựa trên khai triển Taylor Giải PTPT: Phương pháp Newton (2) • Khai triển Taylor: Giả sử f, f ’,…,f(n) liên tục trên [a,b]; f(n+1)(x) tồn tại với mọi... (11) Giải PTPT: Phương pháp Bairstow (8) • Tóm tắt phương pháp Bairstow: – (1) Khởi tạo giá trị p, q; tính b0, b1 theo (7); – (2) Tính (b0)p, (b1)p, (b0)q, (b1)q theo (10) và (11); – (3) Giải (9) để tìm Δp và Δq; – (4) Tìm p* và q* theo công thức: p*=p+ Δp;q*=q+ Δq Các bước trên được lặp lài nhờ sử dụng p*, q* của bước trước như là giá trí khởi tạo p, q của bước sau Một số hàm trên MatLab để giải phương. .. bằng tìm nghiệm của các đa thức bậc 2 Giải PTPT: Phương pháp Bairstow (2) • Mô tả phương pháp Bairstow: – Xét đa thức bậc N: y = a0 + a1x + a2x2 + …+ aNxN Ta có thể viết dưới dạng: y = (x2 + px + q)*G(x) + R(x) • • • p, q chọn tùy ý G(x) là đa thức bậc N-2: G(x) = b2 + b3x + b4x2 + …+ bNxN-2 R(x) là phần dư, thường là bậc 1: R(x) = b0 + b1x (1) (2) (3) (4) Giải PTPT: Phương pháp Bairstow (3) – Nếu chọn . CHƯƠNG 4 GiẢI PHƯƠNG TRÌNH PHI TUYẾN Nội dung Đặt vấn đề 1. Phương pháp chia đôi 2. Phương pháp dây cung 3. Phương pháp Newton 4. Phương pháp cát tuyến 5. Phương pháp lặp 6. Phương pháp. Phương trình phi tuyến (PTPT) – VD1: x 2 = 0 – VD2: 1 + 2x + x 2 - 3x 3 + 7x 4 = 0 – VD3: ln(x+1) = 0 – VD4: tg(x) – artg(2x) = 0 – Tổng quát: f(x) = 0 • Giải phương trình phi tuyến. thể giải kiểu PTPT bất kỳ – Có thể áp dụng cho hàm không biểu diễn dưới dạng giải tích 24-Nov-13 10 Giải PTPT: Phương pháp chia đôi (5) • Bài tập: Viết chương trình Matlab giải phương trình