double f(double x) {
return x * x * x - 5 * x * x + 10 * x - 80;
}
const double N = 1e-6;
int main() {
int num = 1;
double l = 0, r = 100;
double root = l + (r - l) / 2;
double y = f(root);
while (fabs(y)>= N) {
if (y>0) {//根的值大于0
r = root;
}
else l = root;
root = l + (r - l) / 2;
y = f(root);
++num;
}
cout <<"二分查找根和次数"<< root << " " << num;
//上面为二分查找,下面为牛顿迭代法。
cout << endl;
root = 6;
num = 0;
while (f(fabs(root)) >N) {
root = root - ((root * root * root - 5 * root * root + 10 * root - 80) / (6*root*root-10*root+10));
++num;
}
cout << "牛顿迭代法的根"<<root<<" "<<num;
return 0;
}
二分查找不用多讲了,牛顿迭代法记住一个公式就行:
记住这个公式,就完全莫得问题啦!简化一下
root=root-(f(root)/f’(root));这样进行迭代就OK啦!此方程根结果一样,二分查找更快一些。
接下来说开平方根,一想到开根除了用sqrt()函数就真的不知道还能怎么做,昨天做了一道关于开根的题目
题目链接:https://leetcode-cn.com/problems/sqrtx/
用二分查找法怎么做都做出不来,后来看别人的题解,整了一个牛顿迭代法的公式,上来我就蒙逼了,公式如下
r=(r+x/r)/2;
求平方根的函数可以看成f(x)=x^2-n,为什么要-n?为了确保函数得到一个根,n相当于题目里的x。
我看的第一眼心里想的就是这啥玩意啊,咋所有人都知道这是啥,我咋啥也不知道。又翻了许多牛顿迭代法的资料,看也看不懂,后来自己发现上面的这个公式 r=(r+x/r)/2;就是由一开始第一个公式推导而出。第一个公式又好记也好写,个人觉得看起来比第二个公式舒心一些。r=(r+x/r)/2;
这种写法应该是为了防止数溢出。