Chương 5 NỘI SUY VÀ XẤP XỈ HÀM - TaiLieu.VN

Trang chủ » Hàm Nội Suy Trong Matlab » Chương 5 NỘI SUY VÀ XẤP XỈ HÀM - TaiLieu.VN

Có thể bạn quan tâm

OPTADS360 intTypePromotion=1 zunia.vn Tuyển sinh 2024 dành cho Gen-Z zunia.vn zunia.vn tailieu.vn NÂNG CẤP Đăng Nhập | Đăng Ký Chủ đề »
  • Đề thi toán cao cấp 2
  • Đại số tuyến tính
  • Toán rời rạc
  • Xác suất thống kê
  • Phương trình vi phân
    • Toán cao cấp
    • Toán kinh tế
  • HOT
    • FORM.04: Bộ 240+ Biểu Mẫu Chứng Từ Kế...
    • CEO.24: Bộ 240+ Tài Liệu Quản Trị Rủi...
    • CEO.27: Bộ Tài Liệu Dành Cho StartUp...
    • FORM.07: Bộ 125+ Biểu Mẫu Báo Cáo...
    • CEO.29: Bộ Tài Liệu Hệ Thống Quản Trị...
    • FORM.08: Bộ 130+ Biểu Mẫu Thống Kê...
    • TL.01: Bộ Tiểu Luận Triết Học
    • LV.26: Bộ 320 Luận Văn Thạc Sĩ Y...
    • LV.11: Bộ Luận Văn Tốt Nghiệp Chuyên...
    CMO.03: Bộ Tài Liệu Hệ Thống Quản Trị Marketing...
TUYỂN SINH YOMEDIA ADSENSE Trang Chủ » Khoa Học Tự Nhiên » Toán học Chương 5 NỘI SUY VÀ XẤP XỈ HÀM

Chia sẻ: Tuanlong Long | Ngày: | Loại File: PDF | Số trang:7

Thêm vào BST Báo xấu 745 lượt xem 54 download Download Vui lòng tải xuống để xem tài liệu đầy đủ

Tham khảo tài liệu 'chương 5 nội suy và xấp xỉ hàm', tài liệu phổ thông, toán học phục vụ nhu cầu học tập, nghiên cứu và làm việc hiệu quả

AMBIENT/ Chủ đề:
  • giải tích số
  • phương pháp dạy học toán
  • phương pháp giải tích
  • hàm số
  • nội suy
  • xấp xỉ hàm số

Bình luận(0) Đăng nhập để gửi bình luận!

Đăng nhập để gửi bình luận! Lưu

Nội dung Text: Chương 5 NỘI SUY VÀ XẤP XỈ HÀM

  1. Chương 5 NỘI SUY VÀ XẤP XỈ HÀM 5.1 NỘI SUY BẰNG ĐA THỨC Nội suy là cơ sở của nhiều khái niệm trong giải tích số. Đó là công cụ để khôi phục các đặc trưng liên tục của một hàm số y=f(x) từ các tập hợp dữ liệu rời rạc do đo đạc hay quan sát được. Khi f(x) là một hàm phức tạp, khó tính toán và khảo sát thì cũng cần được xấp xỉ bởi một đa thức. Nội suy đơn giản nhất là nội suy bằng đa thức. Lý do đa thức là một hàm đơn giản: dễ tính đạo đạo hàm và nguyên hàm… Nội suy bằng đa thức là tìm một đa thức P(xi) bậc n-1 qua n mốc nội suy xi với i  1, n thỏa mãn P(xi )= f(xi ). Nói cách khác, có thể mô tả tập các điểm dữ liệu rời rạc của hàm y = f(x) dưới dạng bảng: x x1 x2 ... xn y y1 y2 ... yn Sau đó tính các hệ số của đa thức P(x) bậc n-1 thỏa mãn: yi=P (xi), i  1, n (5.1) Bây giời ta cần xây dựng công thức tính các hệ số của đa thức P(x). Giả sử đa thức P(x) được viết dưới dạng tường minh: P(x) = p1xn-1 + p2xn-2+...+pn-1x+ pn Từ điều kiện (5.1), để tìm các hệ số pi của đa thức nội suy P(x) ta có thể giải hệ phương trình sau đây:  p x n 1  p2 x1  2  ... n pn 1 x1  pn  y1  11  p1 x2 1  n p2 x2  2  ... n pn 1x1  pn  y2   ... ... ... ...  n 1  p1 xn  p2 xn  2  ... n pn 1xn  pn  yn 111
  2.  x n 1... x 2 x1 1   p1   y1   1 1      x2 1... x2 n 2 x2 1   p2   y2  hay    .  ... ... ... ...   ...   ...   n 1      xn ... xn 2 xn 1   pn   yn  Ma trận hệ số của hệ phương trình chính là ma trận Vandermonde của vector x=( x1,x2,...xn). Do đó, để giải hệ phương trình trên có thể sử dụng một câu lệnh đơn giản trong Matlab: >> p =vander(x)\y (5.2) Đa thức nội suy được tính theo công thức trên đơn giản, nhưng khá cồng kềnh. Tính hệ số của đa thức nội suy dạng tường minh bằng giải hệ phương trình như trên sẽ mất nhiều công sức tính toán khi số nút nội suy lớn. Do đó ta cần phải nghiên cứu một số phương pháp tìm đa thức nội suy khác, đơn giản hơn. Định lý 5.1 (Tính duy nhất của đa thức nội suy). Đa thức nội suy bậc n-1 thoả mãn (5.1) là duy nhất. Chứng minh. Thật vậy. Giả sử có 2 đa thức nội suy bậc n-1 là P(x) và Q(x) cùng thoả mãn điều kiện (5.1). Nghĩa là yi=P(xi)=Q(x) với i  1, n . Xét đa thức R(x)=P(x)-Q(x). Rõ ràng là: R(xi)=P(xi)-Q(xi) =0 , i  1, n . Như vậy R(x) là đa thức bậc không quá n-1 nhưng có tới n nghiệm khác nhau x1,x2,...xn. Do đó R(x) =0 với mọi x, hay P(x)  Q(x). Điều đó chứng tỏ rằng đa thức nội suy bậc n-1 thỏa mãn (5.1) là duy nhất (đ.p.c.m). Khi số nút nội suy lớn, thì việc giải hệ phương trình như trên tốn rất nhiều công sức. Sau đây chúng ta sẽ nghiên cứu một số phương pháp khác để tìm đa thức nội suy mà không cần giải hệ phương trình. 5.1.1 Đa thức nội suy Lagrange Trước hết ta xây dựng các đa thức cơ bản như sau: Li  x    x  x1  x  x2  ...  x  xi 1  x  xi 1  ...  x  xN  với i  1, n .  xi  x1  xi  x2  ...  xi  xi 1  xi  xi 1  ...  xi  xN  112
  3. Dễ thấy các đa thức cơ bản có tính chất:  0 khi i  j   Li x j   1 khi i  j n Do đó nếu đặt P (x) =  yi Li ( x) (5.3) i 1 thì P(x) là đa thức bậc không quá n thoả mãn P(xi)=yi, với i  1, n . Do đó P(x) chính là đa thức nội suy bậc n-1 của hàm số đã cho. Đa thức dạng (5.3) còn gọi là đa thức nội suy Lagrange. Nó có dạng tổng của n đa thức bậc n-1. Thí dụ 1. Xây dựng đa thức nội suy Lagrange bậc 2 từ bảng dữ liệu có dạng sau đây: x x1 x2 x3 y y1 y2 y3 Giải: Ta có các đa thức nội suy cơ bản: L1 ( x)   x  x2   x  x3  , L ( x)   x  x1   x  x3   x1  x2   x1  x3  2  x2  x1   x2  x3  và L3 ( x)   x  x1  x  x2  .  x3  x1  x3  x2  Do đó P( x)   x  x2   x  x3  y   x  x1   x  x3  y   x  x1  x  x2  y .  x1  x2   x1  x3  1  x2  x1   x2  x3  2  x3  x1  x3  x2  3 5.1.2 Đa thức nội suy Newton Nội suy bằng đa thức Lagrange là một phương pháp khá đơn giản, sử dụng rất ít các kiến thức về đại số, nên dễ nhớ. Tuy nhiên đây lại là một phương pháp kém hiệu quả. Bây giờ xét trường hợp giữ bảng dữ liệu cũ và bổ sung thêm một nút nội suy mới (để hàm số được nội suy chính xác hơn) thì tất cả các đa thức nội suy cơ bản lại phải tính toán lại từ đầu. Thay cho công thức nội suy dạng Lagrange ta viết đa thức nội suy P(x) dưới dạng: P(x)= a 1+ a2(x-x1) + a3(x-x1)(x-x2)+…+ aN(x-x1)(x-x2)(x-x3)...(x-xn-1) (5.4) Các hệ số a i của đa thức có thể được tính trong bảng tỉ hiệu (tỉ sai phân) theo công thức qui nạp như sau: 113
  4. y y f1  xi , xi 1   i 1 i : Tỉ hiệu cấp 1 tại xi ; xi 1  xi f [ x , x ]  f1[ xi , xi 1] f 2  xi , xi 1, xi  2   1 i 1 i  2 : Tỉ hiệu cấp 2 tại xi ;… xi  2  xi f [ x , x ,..., xi  k ]  fk 1[ xi , xi 1,..., xi  k 1] fk  xi ,..., xi  k   k 1 i 1 i  2 :Tỉ hiệu cấp k tại xi; xi  k  xi fn  2 [ x2 , x3 ,..., xn ]  f n  2[ x1, x2 ,..., xn 1] … f n 1  x1, x2 ,..., xn   :Tỉ hiệu cấp n-1 tại x1. xn  x1 Tại nút xi chỉ phải tính các tỉ hiệu cấp 1 đến cấp n-i. Ta lập một bảng để thuận tiện khi tính toán các tỉ hiệu. Bảng 5-1 Bảng tỉ hiệu với 6 núy nội suy x y f1 f2 f3 f4 f5 x1 y1 f1[x1,x2] x2 y2 f2[x1,x2,x 3] f1[x2,x3] f3[x 1,x 2 ,x3,x4] x3 y3 f2[x2,x3,x4 ] f4[x1,x2,x 3,x 4,x5] f1[x3,x4] f3[x2,x3,x4,x5 ] f5[x 1,x2,x3,x4,x5,x6] x4 y4 f2[x3,x4,x5 ] f4[x2,x3,x4,x5,x6] f1[x4,x5] f3[x3,x4,x5,x6 ] x5 y5 f2[x4,x5,x6 ] f1[x5,x6] x6 y6 Khi đó các hệ số của đa thức nội suy (5.4) được xác định như sau: a 1=y1 , a2= f1[x1,x2], a3= f2[x1,x2,x3] ... an= fn-1 [x1,x2,x3,...,xn]. Công thức (5.4) với cách tính các hệ số ai như trên gọi là công thức nội suy Newton tiến xuất phát từ x1. Với công thức này, khi thêm một nút nội suy mới xn+1 thì ta chỉ cần tính thêm một hệ số mới an+1. Khi đó bảng tỉ hiệu chỉ phải thêm một dòng. Mặt khác, nếu chú ý ta thấy công thức Newton không đòi hỏi sự sắp xếp thứ tự về giá trị của dữ liệu, nên khi đảo ngược thứ tự của dữ liệu thì dạng mới của đa thức nội suy là: P(x)=b1+ b2(x-xn) + b 3(x-xn)(x-xn-1)+…+bn(x-xn)(x-xn-1)(x-xn-2)...(x-x2) ( 5.5) Khi đó các hệ số của đa thức nội suy dạng (5.5) được xác định như sau: b1=yn , b 2= f1[xn-1,xn], b3= f2[xn-2,xn-1 ,xn] ... bn= fn-1[x1,x2,x3,...,xn]. 114
  5. Công thức với cách tính các hệ số bi như trên gọi là công thức nội suy Newton lùi xuất phát từ xn. Do định lý 5.1 có thể thấy công thức Lagrange và các công thức Newton tiến hay lùi đều xác định cùng một đa thức, chỉ có hình thức thể hiện là khác nhau và thuận tiện áp dụng cho các trường hợp khác nhau. Thí dụ 2. Cho hàm y =f(x) dưới dạng bảng số sau: x 1 2 3 5 6 8 y 5,230 2,092 1,406 -1,202 -1,321 0,015 Hãy lập đa thức nội suy cho hàm F(x) đã cho dưới dạng: a. Tường minh; b. Lagrange; c. Newton tiến xuất phát từ x1 =1. d. Newton lùi xuất phát từ x6 =8. Giải. a. Để tìm hệ số tường minh của đa thức P(x) ta sử dụng công thức (5.2): >> x=[1 2 3 5 6 8]; >> y=[ 5.230 ; 2.092; 1.406; -1.202 ; -1.321; 0.015]; >> p=vander(x)\y p= -0.0187 0.4201 -3.4803 13.2919 -24.3719 19.3890 Do đó : P  x   0, 0187 x5  0, 4201x 4  3, 4803 x3 + 13.2919 x 2  24,3719 x  19,3890 b. Lập đa thức nội suy Lagrange: ( x  2)( x  3)( x  5)( x  6)( x  8) ( x  1)( x  3)( x  5)( x  6)( x  8) P( x)  5, 230 2, 092  (1  2)(1  3)(1  5)(1  6)(1  8) (2  1)(2  3)(2  5)(2  6)(2  8) 115
  6. BÀI TẬP A. Cài đặt chương trình và lập hàm. 1. Cho tập dữ liệu {(xi,yi)}, i  1, n . Cài đặt hàm tính vector hệ số: a=(a1,a2, a3,…, an) của đa thức nội suy Newton tiến bậc n-1, xuất phát từ x1: P(x) = a 1 +a 2(x-x1)+ a3(x-x1)(x-x2) +...+ a n(x-x1)(x-x2)…(x-xn-1) theo bảng tỉ hiệu. Lệnh gọi hàm có dạng: a = NewtonFit(x,y). 2. Cho tập dữ liệu (xi,yi)}, i  1, n . Cài đặt hàm tính xấp xỉ giá trị của đa thức nội suy P(x) bậc n-1 tại xx=(xx1, xx2,xx3,…,xxm) bằng phương pháp Lagrange. Lệnh gọi hàm có dạng: yy = Lagrange(x,y,xx) trong đó yy là vector giá trị hàm tại xx tính được. 3. Cho tập dữ liệu {(xi,yi)}, i  1, n . Cài đặt hàm tính hệ số a, b và c của quan hệ hàm y = ax2 + bx + c theo phương pháp bình phương tối thiểu. Lệnh gọi hàm có dạng: [a,b,c] =ParaboFit(x,y). 4. Cho tập dữ liệu {(xi,yi)}, i  1, n . Cài đặt hàm tính hệ số a, b và c của quan hệ hàm y = a +b cosx+ c sinx theo phương pháp bình phương tối thiểu. Lệnh gọi hàm có dạng: [a,b,c] =CosSinFit(x,y). 5. Cho tập dữ liệu {(xi,yi)}, i  1, n . Xây dựng công thức tính và cài đặt hàm tính hệ số a và b của quan hệ hàm y = aebx (a>0) theo phương pháp bình phương tối thiểu. Lệnh gọi hàm có dạng: [a,b] =ExpFit(x,y). 133
  7. B. Sử dụng các hàm nội trú của Matlab. 1. Tìm xấp xỉ các hệ số của đa thức P(x) bậc 3 được cho trong bảng bằng phương pháp bình phương tối thiểu và tính P(x) tại x= 1; 1,2; 1.4; ...; 2,8; 3. x 1 1,5 2 2,5 3 3,5 y 1,2341 3,9242 2,4563 - 0,2224 -1,3215 0,5506 2. Đa thức P(x)= ax3 + bx2 +cx +d đo được dưới dạng bảng: x 1 1,5 2 2,5 3 3,5 y 5,23401 2,09242 1,40563 - 1,20224 -1,32105 0,01501 - Tìm các hệ số a, b, c, d bằng phương pháp bình phương tối thiểu. - Tính P(x) tại x = 1,5; 1,6; 1,7; ... 3,2. 3. Hàm số y=f(x) được cho dưới dạng bảng: x 1,0 1,5 2,0 2,5 3,0 3,5 y -2,95001 -2,90250 -1,00020 2,75754 8,37002 15,83752 Tính xấp xỉ hàm f(x) tại x= 1; 1.1; 1.2; ...; 3 bằng phương pháp spline bậc 3. 4. Hàm số y=f(x) được cho bởi bảng số. Tính gần đúng f(x) tại x=1,25 bằng phương pháp spline bậc 3: x 1,0 1,1 1,2 1,3 1,4 1,5 y 3,85011 -2,30650 -1,12320 4,85754 8,36602 10,99154 5. Dùng ma trận pascal, đưa ra màn hình vector hệ số của khai triển nhị thức Newton bậc 10 . 6. Hàm số y=f(x) được cho dưới dạng bảng, Tính vector hệ số của đa thức nội suy bậc 5 bằng cách sử dụng ma trận Vandermonde x 1,0 1,5 2,0 2,5 3,0 3,5 y -2,30650 -1,12320 -1,00020 2,75754 8,37002 15,83752 7. Cho 2 vector: 5x10 + 12x9 -4x7 +3x6 -5 x2 -10 = 0 -3x8 + 6x5 +8x7 +3x3 +2x4 -11x -2 = 0 Hãy tính vector p, q và s tương ứng là hệ số của tổng tích, thương và phần dư phép chia của 2 đa thức đã cho. 134
ADSENSE

CÓ THỂ BẠN MUỐN DOWNLOAD

  • Giáo trình Giải tích số - Lê Minh Lưu

    pdf 77 p | 513 | 185

  • Giáo trình Giải tích số: Phần 1

    pdf 96 p | 396 | 89

  • Giáo trình tính toán khoa học - Chương 5

    pdf 24 p | 304 | 86

  • Bài giảng Chương 5 - Nội suy và xấp xỉ hàm

    pdf 24 p | 180 | 31

Thêm tài liệu vào bộ sưu tập có sẵn: Đồng ý Thêm vào bộ sưu tập mới: *Tên bộ sưu tập Mô Tả: *Từ Khóa: Tạo mới Báo xấu
  • Hãy cho chúng tôi biết lý do bạn muốn thông báo. Chúng tôi sẽ khắc phục vấn đề này trong thời gian ngắn nhất.
  • Không hoạt động
  • Có nội dung khiêu dâm
  • Có nội dung chính trị, phản động.
  • Spam
  • Vi phạm bản quyền.
  • Nội dung không đúng tiêu đề.
Hoặc bạn có thể nhập những lý do khác vào ô bên dưới (100 ký tự): Vui lòng nhập mã xác nhận vào ô bên dưới. Nếu bạn không đọc được, hãy Chọn mã xác nhận khác.. Đồng ý LAVA AANETWORK THÔNG TIN
  • Về chúng tôi
  • Quy định bảo mật
  • Thỏa thuận sử dụng
  • Quy chế hoạt động
TRỢ GIÚP
  • Hướng dẫn sử dụng
  • Upload tài liệu
  • Hỏi và đáp
HỖ TRỢ KHÁCH HÀNG
  • Liên hệ
  • Hỗ trợ trực tuyến
  • Liên hệ quảng cáo
Theo dõi chúng tôi

Chịu trách nhiệm nội dung:

Nguyễn Công Hà - Giám đốc Công ty TNHH TÀI LIỆU TRỰC TUYẾN VI NA

LIÊN HỆ

Địa chỉ: P402, 54A Nơ Trang Long, Phường 14, Q.Bình Thạnh, TP.HCM

Hotline: 093 303 0098

Email: support@tailieu.vn

Giấy phép Mạng Xã Hội số: 670/GP-BTTTT cấp ngày 30/11/2015 Copyright © 2022-2032 TaiLieu.VN. All rights reserved.

Đang xử lý... Đồng bộ tài khoản Login thành công! AMBIENT

Từ khóa » Hàm Nội Suy Trong Matlab