Chương trình C ++ để thực hiện phép nội suy ngược sử dụng công thức Lagrange

Trong hướng dẫn này, chúng ta sẽ thảo luận về một chương trình thực hiện phép nội suy ngược sử dụng công thức Lagrange.

Nội suy ngược được định nghĩa là phương pháp tìm giá trị của một biến độc lập từ giá trị đã cho của giá trị phụ thuộc nằm giữa hai tập giá trị được lập bảng cho một hàm chưa biết.

Ví dụ

#include <bits/stdc++.h>
using namespace std;
//structuring the values of x and y
struct Data {
   double x, y;
};
//calculating inverse interpolation
double calc_invinter(Data d[], int n, double y){
   double x = 0;
   int i, j;
   for (i = 0; i < n; i++) {
      double xi = d[i].x;
      for (j = 0; j < n; j++) {
         if (j != i) {
            xi = xi * (y - d[j].y) / (d[i].y - d[j].y);
         }
      }
      x += xi;
   }
   return x;
}
int main(){
   Data d[] = {
      { 1.27, 2.3 },
      { 2.25, 2.95 },
      { 2.5, 3.5 },
      { 3.6, 5.1 }
   };
   int n = 6;
   double y = 4.5;
   cout << "Value of x (y = 4.5) : " << calc_invinter(d, n, y) << endl;
   return 0;
}

Đầu ra

Value of x (y = 4.5) : 2.51602