Lagrange插值，自写py代码

yinda_peng

1. import numpy as np
   
2. 
   
3. class Interpolation:
   
4.     def Lagrange(self, xx, x, y):
   
5.         m = len(x)
   
6.         n = len(y)
   
7.         if m != n:
   
8.             print("x与y维数不同")
   
9. 
   
10.         s = np.ones(len(x))
    
11.         for i in range(0, n):
    
12.             d = 1
    
13.             for j in range(0, n):
    
14.                 if i != j:
    
15.                     d = d * (xx - x[j]) / (x[i] - x[j])
    
16.             s[i] = d
    
17.         f = y * s
    
18.         return sum(f)
    
19. 
    
20. class Main:
    
21.     def run(self):
    
22.         input_str = input("请输入一个x列表（用空格分隔各元素）: ")
    
23.         x = [int(i) for i in input_str.split()]
    
24.         
    
25.         input_str = input("请输入一个y列表（用空格分隔各元素）: ")
    
26.         y = [int(i) for i in input_str.split()]
    
27. 
    
28.         print("输入的列表为:\n", x, '\n', y)
    
29. 
    
30.         interpolator = Interpolation()
    
31. 
    
32.         while True:
    
33.             a = input("是否输入插值点（Y/N）\n")
    
34.             if a.lower() == 'n':
    
35.                 break
    
36.             xx = int(input("请输入插值点: "))
    
37.             print(interpolator.Lagrange(xx, x, y))
    
38. 
    
39. #如果以文件形式，文件内以元组((x1,y1),(x2,y2)...)形式传入
    
40. '''
    
41. class Main:
    
42.     def read_data_from_file(self, file_path):
    
43.         data = []
    
44.         with open(file_path, 'r') as file:
    
45.             lines = file.readlines()
    
46.             for line in lines:
    
47.                 line = line.strip()
    
48.                 if line:
    
49.                     pair = tuple(map(int, line.split(',')))
    
50.                     data.append(pair)
    
51.         return data
    
52. 
    
53.     def run(self):
    
54.         file_path = input("请输入包含数据的文件路径：")
    
55.         data = self.read_data_from_file(file_path)
    
56. 
    
57.         x = []
    
58.         y = []
    
59.         for pair in data:
    
60.             x.append(pair[0])
    
61.             y.append(pair[1])
    
62. 
    
63.         print("输入的列表为:\n", x, '\n', y)
    
64. 
    
65.         interpolator = Interpolation()
    
66. 
    
67.         while True:
    
68.             a = input("是否输入插值点（Y/N）\n")
    
69.             if a.lower() == 'n':
    
70.                 break
    
71.             xx = int(input("请输入插值点: "))
    
72.             print(interpolator.Lagrange(xx, x, y))
    
73. '''
    
74. 
    
75. if __name__ == "__main__":
    
76.     main = Main()
    
77.     main.run()

复制代码

FishC

你好！下面是使用Lagrange插值方法的Python代码示例：

1. import numpy as np
   
2. 
   
3. class Interpolation:
   
4.     def Lagrange(self, xx, x, y):
   
5.         m = len(x)
   
6.         n = len(y)
   
7.         if m != n:
   
8.             print("x与y维数不同")
   
9. 
   
10.         s = np.ones(len(x))
    
11.         for i in range(0, n):
    
12.             d = 1
    
13.             for j in range(0, n):
    
14.                 if i != j:
    
15.                     d = d * (xx - x[j]) / (x[i] - x[j])
    
16.             s = s * d
    
17.         f = y * s
    
18.         return sum(f)
    
19. 
    
20. class Main:
    
21.     def run(self):
    
22.         input_str = input("请输入一个x列表（用空格分隔各元素）: ")
    
23.         x = list(map(float, input_str.split()))
    
24. 
    
25.         input_str = input("请输入一个y列表（用空格分隔各元素）: ")
    
26.         y = list(map(float, input_str.split()))
    
27. 
    
28.         print("输入的列表为:\n", x, '\n', y)
    
29. 
    
30.         interpolator = Interpolation()
    
31. 
    
32.         while True:
    
33.             a = input("是否输入插值点（Y/N）\n")
    
34.             if a.lower() == 'n':
    
35.                 break
    
36.             xx = float(input("请输入插值点: "))
    
37.             print(interpolator.Lagrange(xx, x, y))
    
38. 
    
39. if __name__ == "__main__":
    
40.     main = Main()
    
41.     main.run()

复制代码

这段代码实现了Lagrange插值方法。在运行代码时，首先需要输入一个x列表和一个y列表，其中各元素之间用空格分隔。然后，程序会不断询问是否输入插值点，如果输入'Y'则继续，如果输入'N'则结束程序。如果需要输入插值点，则需要输入一个插值点的值。程序会根据输入的x列表、y列表和插值点的值，使用Lagrange插值方法计算出插值点对应的y值，并将结果输出。

以上回复来自 FishC AI 小助理，如未能正确解答您的问题，请继续追问。