+ Trả lời chủ đề
Hiện kết quả từ 1 tới 2 của 2

Chủ đề: Đa thức nội suy - Phương pháp Newton tiến

  1. #1
    Status : Administrator đang ẩn
    Tham gia ngày : Aug 2011
    Đến từ : Ninh Kiều - Cần Thơ
    Bài gửi : 1.450

    Đa thức nội suy - Phương pháp Newton tiến

    Mã:
    TyHieu[X_, Y_, i_, j_, Cap_] :=
        		If[Cap == 1,
          				Return[(Y[[i]] - Y[[j]])/(X[[i]] - X[[j]])];
          			,
          				Return[(TyHieu[X, Y, i, j - 1, Cap - 1] - TyHieu[X, Y, i + 1, j, Cap - 1])/(X[[i]] -X[[j]])];
          		];
    NoiSuyNewtonTien[X_, Y_, n_] := Module[{P, i, j, temp},
          		P = Y[[1]];
          		temp = 1;
          		For[i = 2, i <= n, i++,
            			temp = temp*(x - X[[i - 1]]);
            			P = P + temp *TyHieu[X, Y, 1, i, i - 1];
            		];
          		Return[Expand[P]];
          ];
    
    Clear[X, Y, n];
    X = Input["Nhap danh sach toa do X"];
    Y = Input["Nhap danh sach toa do Y"];
    n = Length[X];
    Print["X = ", X];
    Print["Y = ", Y];
    Print["Da thuc noi suy Newton:"];
    Print["P(x) = ", NoiSuyNewtonTien[X, Y, n] ];

  2. #2
    Status : Administrator đang ẩn
    Tham gia ngày : Aug 2011
    Đến từ : Ninh Kiều - Cần Thơ
    Bài gửi : 1.450
    Phương pháp Newton lùi:
    Mã:
    NoiSuyNewtonLui[X_, Y_, n_] := Module[{P, i, j, temp},
          		P = Y[[n]];
          		temp = 1;
          		For[i = n - 1, i >= 1, i--,
            			temp = temp*(x - X[[i + 1]]);
            			P = P + temp *TyHieu[X, Y, i, n, n - i];
            		];
          		Return[Expand[P]];
          ];

+ Trả lời chủ đề

Từ khóa (Tag) của chủ đề này

Quyền viết bài

  • Bạn không thể gửi chủ đề mới
  • Bạn không thể gửi trả lời
  • Bạn không thể gửi file đính kèm
  • Bạn không thể sửa bài viết của mình
Trang Chủ Việc Làm Gia Sư Gia sư