Newton法、弦截法

一、简介

牛顿迭代法（Newton’s method）又称为牛顿-拉夫逊（拉弗森）方法（Newton-Raphson method），它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。多数方程不存在求根公式，因此求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要。方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根，此时线性收敛，但是可通过一些方法变成超线性收敛。另外该方法广泛用于计算机编程中。

弦截法是一种求方程根的基该方法，在计算机编程中常用。
他的思路是这样的：任取两个数x1、x2，求得对应的函数值f(x1）、f(x2）。如果两函数值同号，则重新取数，直到这两个函数值异号为止。连接（x1,f(x1））与（x2,f(x2））这两点形成的直线与x轴相交于一点x，求得对应的f(x），判断其与f(x1）、f(x2）中的哪个值同号。如f(x）与f(x1）同号，则f(x）为新的f(x1)。将新的f(x1）与f(x2）连接，如此循环。体现的是极限的思想。

牛顿法详解

二、实现

# -*- coding: utf-8 -*-
"""
Created on Fri Dec 16 12:41:37 2016

    Newton法、弦截法

@author: Administrator
"""

from numpy import *
from matplotlib import pyplot as plt 

def f(x):
    return x**3 + 2 * x**2 + 10 *x - 100

def f_derivative(x):
    return 3 * x**2 + 4 *x + 10

#牛顿法
def newton():
    x0 = 0
    x1 = 1.1
    count = 0
    while True:
        x0 = x1
        f_derivative_num = f_derivative(x0)
        if f_derivative_num != 0:
            x1 = x0 - f(x0) / f_derivative_num
        else:
            return
        if abs(x1 - x0) < 1e-10:
            break
        else:
            print 'Newton:',count
            print abs(x1 - x0)
            count += 1
    return x1

#弦截法
def secant():
    x0 = 0
    x1 = 1.1
    count = 0
    while True:       
        x2 = x1 - f(x1) * (x1 - x0) / (f(x1) - f(x0))         
        if abs(x2-x1) < 1e-10:             
            break
        else:
            print 'Secant:',count
            print abs(x2 - x1)
            x0 = x1
            x1 = x2
            count += 1
    return x2

if __name__ == '__main__':
    ans1 = newton()
    ans2 = secant()

    x = linspace(-5,5,1000)
    y = f(x)
    fig = plt.figure(figsize=(8,4))
    ax = fig.add_subplot(111)
    ax.plot(x,y,color='r',linestyle='-',label='f(x)')
    ax.scatter(ans1,f(ans1),label='Newton',color='y')
    ax.scatter(ans2,f(ans2),label='Secant',color='k')
    ax.legend(loc='lower right')
    fig.show()
    fig.savefig('a.png')