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

最新推荐文章于 2021-11-06 20:07:51 发布

大菜市场

最新推荐文章于 2021-11-06 20:07:51 发布

阅读量1.3k

点赞数

CC 4.0 BY-SA版权

分类专栏：数值分析 Python

本文链接：https://blog.youkuaiyun.com/qq_36440163/article/details/89052509

这篇博客介绍了数值分析中的牛顿差商法，并提供了Python3.6环境下的实现代码，利用matplotlib和numpy库。欢迎大家进行交流和指正。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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

仅供参考

配置环境

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]):
            #对每行
            for i in range(j-1,W.shape[0]):
                W[i,j] = (W[i,j-1] - W[i-1,j-1]) / (W[