Bài giảng Phương Pháp Tính - Chương II

Bài giảng Phương Pháp Tính - Chương II

Chương 2NỘI SUY VÀ LẤY XẤP XỈ HÀM SỐI. Nội suy. 1. Đặt vấn đề.- không biết biểu thức giải tích của hàm;y = f ( x );- biết giá trị của hàm tại một số hữu hạn điểm trên đoạn [ a, b] ( bằng đo đạc hoặc thực nghiệm):xy=f( x )xo x1 . . . xi . . . xn-1 xnyo y1 . . . yi . . . yn-1 yn - Tìm giá trị của hàm số tại một số điểm trung gian khác.-Xây dựng một hàm φ( x ) có biểu thức đơn giản, có giá trị trùng với giá trị của f ( x ) tại các điểm , còn tại các điểm khác trên đoạn [a, b] thì φ( x ) khá gần f( x ) [phản ảnh gần đúng quy luật f( x ) ] có thể suy ra giá trị gần đúng của f( x ) tại các giá trị x bất kỳ thoả mãn xo < x < xn.xo x1 . . . xn Bài toán nội suy:- φ( x ) – hàm nội suy của f( x ) trên đoạn [a, b]. -Ý nghĩa hình học: xây dựng đường cong y = φ( x ) đi qua các điểm cho trước (xi, yi), I = 0, 1, . . . , n. 2. Đa thức nội suy.Thường chọn đa thức làm hàm nội suy vì:- Đa thức là loại hàm đơn giản;- Luôn có đạo hàm và nguyên hàm;- Việc tính giá trị của chúng đơn giản.- Trên đoạn a ≤ x ≤ b cho một lưới các điểm chia (điểm nút)xi; i = 0, 1, 2, . . ., n; vớia ≤ xo, x1, x2, . . . .xn ≤ b.- Cho giá trị tương ứng của hàm y = f( x ) tại các nút:xy=f( x )xo x1 . . . xi . . . xn-1 xnyo y1 . . . yi . . . yn-1 yn- Cần xây dựng một đa thức bậc n:Pn( x ) = aoxn + a1xn-1 + . . . + an-1x + an.sao cho Pn(x) trùng với f(x) tại các nút xi,Pn(xi) = yi ; i= 0, 1, 2, . . . , n.Bài toán: )1( Sự duy nhất của đa thức nội suy: đa thức nội suy Pn( x ) của hàm f( x ) định nghĩa ở trên nếu có thì chỉ có một mà thôi.Đa thức nội suy có thể xây dựng theo nhiều cách nhưng do tính duy nhất nên các dạng của nó đều có thể quy về nhau.3. Đa thức nội suy Lagrăng.trong đó;)) .()() .(()) .()() .(()(1111niiiiioiniioixxxxxxxxxxxxxxxxxl−−−−−−−−=+−+−li(x) – đa thức bậc n Pn(x) – đa thức bậc n =)(jixl1 khi j = i;0 j ≠ I ; Pn(xi)= yi; i = 0, 1, 2, . . ., n.);()(0xlyxPiniin∑==( 2 )( 2 ) – đa thức nội suy Lagrăng. */ Nội suy tuyến tính. ( n = 1 )xyxox1yoy1( 2 ));()()(111xlyxlyxPoo+=( 3 )ooooxxxxxlxxxxxl−−=−−=1111)(;)(;)(1111oooonxxxxyxxxxyxP−−+−−=Có dạng P1(x) = Ax + B - bậc nhất đối với x.*/ Nội suy bậc 2. ( n = 2 )xyxox1 x2yoy1 y2);()()()(22112xlyxlyxlyxPoo++=( 4 )))(())(()(;))(())(()(21212121xxxxxxxxxlxxxxxxxxxloooooo−−−−=−−−−=;))(())(()(12212xxxxxxxxxloo−−−−=P2(x) có dạng : P2(x) = Ax2 + Bx + C - bậc 2 đối với x. */ Sai số nội suy.Định lý. Nếu hàm f(x) liên tục trên [a, b] và có trong [a, b] đạohàm liên tục đến cấp n+1 thì sai số nội suy rn(x) =f(x) – Pn(x) có biểu thức:[ ]bacnxcfxrnn,;)!1()()()()1(∈+=+π( 5 ))) .()(()(1 noxxxxxxx −−−=π([a, b] - khoảng chứa các nút xi)Gọi[ ]baxxfMn,;)(max)1(∈=+thì;)()!1()( xnMxrnπ+≤-Ưu điểm của đa thức nội suy Lagrăng : đơn giản;- Nhược điểm : thêm một nút thì phải tính lại toàn bộ. 4. Đa thức nội suy Niutơn.a/ Sai phân hữu hạn.y = f(x) có giá trị yi = f(xi) tại các nút xi cách đều nhau vớixi+1 – xi = h = const; i = 0, 1, 2, . . ., nx0 ; x1 = xo + h; x2 = x0 + 2h; . . .xi = xo + ih . . .Định nghĩa sai phân hữu hạn của hàm y = f(x):Sai phân cấp 1(hạng 1): Δyi = yi+1 – yi;Sai phân cấp hai:Δ2yi = Δyi+1 –Δyi;Δnyi = Δ(Δn-1yi) = Δn-1yi+1 –Δn-1yi;Sai phân cấp n là sai phân của sai phân cấp n-1:Sai phân cấp 3:Δ3yi = Δ2yi+1 –Δ2yi;;1 ooyyy−=∆;2)()(121122oooyyyyyyyy +−=−−−=∆;331233ooyyyyy−+−=∆;)1()1(!3)2)(1(!2)1(!111321 onnnnnnonyyynnnynnynyy−+−+−−−−+−=∆−−−−. . . . . . . . . . . . . . . . . . . . . b/ Đa thức nội suy Niutơn tiến (nội suy về phía phải).Trường hợp các nút cách đều, Niutơn đa thức có dạng:)()())(()()(1121−−⋅⋅⋅−+⋅⋅⋅+−−+−+=nonooonxxxxaxxxxaxxaaxP- x = xo ao = Pn(xo) = yo;- x = x1;)(1111hyhyyxxaxPaoooon∆=−=−−=- x = x2;!2222)(2.22))(()()(2221221221222122hyhyyyhyyyyhhhhyyyxxxxxxaaxPaooooooooon∆=+−=−−−=⋅∆−−=−−−−−=;!ioiihiya∆=;)(iinyxP=noaaaa ,,,,21⋅⋅⋅Xác địnhtừ điều kiện );()(!))((!2)(!1)(1122−−⋅⋅⋅−∆+⋅⋅⋅+−−∆+−∆+=nononooooonxxxxhnyxxxxhyxxhyyxPĐổi biến, đặt hxxto−=x = xo + t.hxi = xo + i.hx - xi = x - xo - i.h = t.h - i.h = (t-i)h;;!)1) .(1(!2)1()(2onooonynntttyttytyxP∆+−−+⋅⋅⋅+∆−+∆+=Thường dùng để tính các giá trị của hàm ở gần xo đầu bảng.3/ Đa thức nội suy Niutơn lùi (nội suy về phía phải)Với cách làm tương tự nhưng bắt đầu với x = xn);()(!))((!2)(!1)(112221xxxxhnyxxxxhyxxhyyxPnnonnnnnnnn−⋅⋅⋅−∆+⋅⋅⋅+−−∆+−∆+=−−−;!)1) .(1(!2)1(!1)(221 onnnnnynntttyttytyxP∆−+++⋅⋅⋅+∆++∆+=−−Ưu điểm của công thức nội suy Niutơn: thêm nút chỉ cần thêm số hạng, không cần phải tính lại.Để thuận tiện tính toán thường lập bảng sai phân đường chéo. x y ΔyΔ2yΔ3y Δ4y Δ5y Δ6yxox1x2x3x4x5x6yoy1y2y3y4y5y6ΔyoΔy1Δy2Δy3Δy4Δy5Δ2yoΔ2y1Δ2y2Δ2y3Δ2y4Δ3yoΔ3y1Δ3y2Δ3y3Δ4yoΔ4y1Δ4y2Δ5yoΔ5y1Δ6yoBảng sai phân đường chéo của công thức nội suy tiến x y ΔyΔ2yΔ3y Δ4y Δ5y Δ6yxn-6xn-5xn-3xn-2xn-1xnΔyn-6Δyn-5Δyn-4Δyn-3Δyn-2Δyn-1Δ2yn-6Δ2yn-5Δ2yn-4Δ2yn-3Δ2yn-2Δ3yn-6Δ3yn-5Δ3yn-4Δ3yn-3Δ4yn-6Δ4yn-5Δ4yn-4Δ5yn-6Δ5yn-5Δ6yn-6Bảng sai phân đường chéo của công thức nội suy lùixn-4yn-6yn-5yn-3yn-2yn-1ynyn-4 [...]... dạng : P 2 (x) = Ax 2 + Bx + C - bậc 2 đối với x. Sự duy nhất của đa thức nội suy: đa thức nội suy P n ( x ) của hàm f( x ) định nghĩa ở trên nếu có thì chỉ có một mà thơi. Đa thức nội suy có thể xây dựng theo nhiều cách nhưng do tính duy nhất nên các dạng của nó đều có thể quy về nhau. 3. Đa thức nội suy Lagrăng. trong đó ; )) ()() (( )) ()() (( )( 11 11 niiiiioi niio i xxxxxxxx xxxxxxxx xl −−−− −−−− = +− +− l i (x)... VÀ LẤY XẤP XỈ HÀM SỐ I. Nội suy. 1. Đặt vấn đề. - không biết biểu thức giải tích của hàm; y = f ( x ); - biết giá trị của hàm tại một số hữu hạn điểm trên đoạn [ a, b] ( bằng đo đạc hoặc thực nghiệm): x y=f( x ) x o x 1 . . . x i . . . x n-1 x n y o y 1 . . . y i . . . y n-1 y n - Tìm giá trị của hàm số tại một số điểm trung gian khác. - Xây dựng một hàm φ( x ) có biểu thức đơn giản,... nội suy Lagrăng. - Nội dung của phương pháp bình phương nhỏ nhất là tìm cực tiểu của hàm [ ] ;)()( 0 2 ∑ = −= n i ii xQxfM ( c ) f(x i ), i = 0, 1, . . ., n là những giá trị đã biết của f(x) trên tập điểm x o , x 1 , . . ., x n . Coi các hệ số C k là các ẩn phải tìm, để M cực tiểu phải có: ), ,1,0(;0 mk C M k == ∂ ∂ ( d ) Từ ( d ) tìm được các hệ số C k . f(x i ) – Q(x i ) - Sai số tại x i khi... Δy i+1 –Δy i ; Δ n y i = Δ(Δ n-1 y i ) = Δ n-1 y i+1 –Δ n-1 y i ; Sai phân cấp n là sai phân của sai phân cấp n-1: Sai phân cấp 3: Δ 3 y i = Δ 2 y i+1 –Δ 2 y i ; ; 1 oo yyy −=∆ ;2)()( 12112 2 ooo yyyyyyyy +−=−−−=∆ ;33 123 3 oo yyyyy −+−=∆ ;)1()1( !3 )2)(1( !2 )1( !1 1 1 321 o nn nnnno n yyy nnn y nn y n yy −+−+ −− − − +−=∆ − −−− . . . . . . . . . . . . . . . . . . . . . Chương 2 NỘI SUY VÀ LẤY XẤP... dạng: )()())(()()( 1121 − −⋅⋅⋅−+⋅⋅⋅+−−+−+= nonooon xxxxaxxxxaxxaaxP - x = x o a o = P n (x o ) = y o ; - x = x 1 ; )( 1 1 1 1 h y h yy xx axP a oo o on ∆ = − = − − = - x = x 2 ; !22 2 2 )(2 .2 2 ))(( )()( 2 2 2 12 2 12 2 122 212 2 h y h yyy h yyyy hh h h y yy xxxx xxaaxP a oooo o o o oon ∆ = +− = −−− = ⋅ ∆ −− = −− −−− = ; ! i o i i hi y a ∆ = ;)( iin yxP = no aaaa ,,,, 21 ⋅⋅⋅ Xác định từ điều kiện 4. Đa... có biểu thức: [ ] bac n x cfxr n n ,; )!1( )( )()( )1( ∈ + = + π ( 5 ) )) ()(()( 1 no xxxxxxx −−−= π ([a, b] - khoảng chứa các nút x i ) Gọi [ ] baxxfM n ,;)(max )1( ∈= + thì ;)( )!1( )( x n M xr n π + ≤ - Ưu điểm của đa thức nội suy Lagrăng : đơn giản; - Nhược điểm : thêm một nút thì phải tính lại tồn bộ. ... x ) [phản ảnh gần đúng quy luật f( x ) ] có thể suy ra giá trị gần đúng của f( x ) tại các giá trị x bất kỳ thoả mãn x o < x < x n . x o x 1 . . . x n Bài toán nội suy: - φ( x ) – hàm nội suy của f( x ) trên đoạn [a, b]. - Ý nghĩa hình học: xây dựng đường cong y = φ( x ) đi qua các điểm cho trước (x i , y i ), I = 0, 1, . . . , n. x y Δy Δ 2 y Δ 3 y Δ 4 y Δ 5 y Δ 6 y x o x 1 x 2 x 3 x 4 x 5 x 6 y o y 1 y 2 y 3 y 4 y 5 y 6 Δy o Δy 1 Δy 2 Δy 3 Δy 4 Δy 5 Δ 2 y o Δ 2 y 1 Δ 2 y 2 Δ 2 y 3 Δ 2 y 4 Δ 3 y o Δ 3 y 1 Δ 3 y 2 Δ 3 y 3 Δ 4 y o Δ 4 y 1 Δ 4 y 2 Δ 5 y o Δ 5 y 1 Δ 6 y o Bảng... Δ 6 y x o x 1 x 2 x 3 x 4 x 5 x 6 y o y 1 y 2 y 3 y 4 y 5 y 6 Δy o Δy 1 Δy 2 Δy 3 Δy 4 Δy 5 Δ 2 y o Δ 2 y 1 Δ 2 y 2 Δ 2 y 3 Δ 2 y 4 Δ 3 y o Δ 3 y 1 Δ 3 y 2 Δ 3 y 3 Δ 4 y o Δ 4 y 1 Δ 4 y 2 Δ 5 y o Δ 5 y 1 Δ 6 y o Bảng sai phân đường chéo của công thức nội suy tiến */ Nội suy tuyến tính. ( n = 1 ) x y x o x 1 y o y 1 ( 2 ) );()()( 111 xlyxlyxP oo += ( 3 ) o o o o xx xx xl xx xx xl − − = − − = 1 1 1 1 )(;)( ;)( 1 1 1 1 o o o on xx xx y xx xx yxP − − + − − = Có dạng P 1 (x) = Ax + B - bậc nhất đối với x. */ Nội suy bậc 2. ( n = 2 ) x y x o x 1 x 2 y o y 1 y 2 );()()()( 22112 xlyxlyxlyxP oo ++= (... b/ Đa thức xấp xỉ là những hàm phi tuyến. Trong nhiều trường hợp có thể tuyến tính hố để việc tính tốn thuận lợi: */ ; 1 bax y + = Đặt Y = 1/y Y = ax + b; */ ; bax x y + = Đặt Y = 1/y và X = 1/x Y = a + bX; */ ; bx aey = với a > 0; Lấy logarit thập phân cả 2 vế: logy = loga + xbloge; Đặt . Δ5y Δ6yxn-6xn-5xn-3xn-2xn-1xnΔyn-6Δyn-5Δyn-4Δyn-3Δyn-2Δyn-1Δ2yn-6Δ2yn-5Δ2yn-4Δ2yn-3Δ2yn-2Δ3yn-6Δ3yn-5Δ3yn-4Δ3yn-3Δ4yn-6Δ4yn-5Δ4yn-4Δ5yn-6Δ5yn-5Δ6yn-6Bảng. a, b. */ Q(x) = ax2 + bx + c.( );022∑=−−−=niiiicbxaxyM;0;0;0 =∂∂=∂∂=∂∂cMbMaM;;;2232234iiiiiiiiiiiiiycnxbxaxyxcxbxaxyxcxbxa∑=+∑+∑∑=∑+∑+∑∑=∑+∑+∑xiyixi2xi2yiΣ.