牛顿法推导平方根迭代公式

平方根迭代公式

假设√a = x
则平方根迭代公式为：
x_n+1 = 1/2 (x_n + a/x_n)

推导过程

假设有如上图的曲线，x_n点的直线为曲线的切线，容易得到x_n的方程为 f(x_n) + f^’(x_n)(x - x_n)

上面的方程是根据泰勒级数展开后，去除高阶项后得到。

这样让切线方程等于0的点的x设为x_n+1，方程变为了f(x_n) + f^’(x_n)(x_n+1 - x_n)(粗体是导数) = 0

解得x_n+1 = x_n - f(x_n)/f^’(x_n)(粗体是导数)。

进入今天的主题，已知√a = x，所以x² - a = 0，曲线如下图，


我们任意取x_n点，假设是上面的两种情况。
上面的推导公式同样适用于当前曲线情况。通过上面的迭代公式x_n+1 = x_n - f(x_n)/f^’(x_n)(粗体是导数)和图像不难看出，x_n+1总是向y = 0的x点趋近，这样只要我们选择了x_n后，不断地迭代求x_n+1，它总是能够到达y = 0的点，此时的x便是平方根的值。

我们将x² - a = 0带入到上面的迭代公式，得到了x_n+1 = 1/2 (x_n + a/x_n)。
需要注意的是我们的根号后得到的x为正值，而二次函数显然有两个值，一正一负，这就需要我们选择的初始x_n点应该为正值。

实现代码

代码就是简单的实现了一下，还有漏洞，希望读者自己查缺补漏吧。

#include <iostream>
#include <cmath>

using namespace std;

double f(double x, double a)
{
	x = (x + a/x)/2;
	
	if (abs(a-x*x) < 0.000001)
		return x; 
		
	return f(x, a);
}

int main()
{
	double n;
	
	cin >> n;
	
	cout << f(5, n) << endl;
	
	return 0;
}