【数值分析】牛顿差商法 -Python实现

本篇为数值分析课程代码实现-牛顿差商法 的实现

仅供参考

配置环境

• Python3.6
• matplotlib
• numpy

• # -*- coding: utf-8 -*-
  """
  Created on Tue Dec  4 21:04:35 2018
  @author: hhuaf
  """
  import matplotlib.pyplot as plt
  import numpy as np
  
  #可以显示中文
  plt.rcParams["font.sans-serif"] = ["SimHei"]
  plt.rcParams['axes.unicode_minus'] = False
  
  # 设置风格
  plt.style.use('ggplot')
  
  class Newton():
  
      def __init__(self):
          
          self.W = None
          self.X = None
          
      def fit(self, X, y):#0-n-1
          
          if type(X) != np.ndarray:
              X=np.array(X)
          if type(y) != np.ndarray:
              y=np.array(y)
          W = np.zeros([len(X),len(X)+1])
          W_list=[]
          W[:,0] = X
          W[:,1] = y
          #针对每列
          W_list.append(W[0,1])
          for j in range(2,W.shape[1]):
              #对每行