求平方根的算法 牛顿迭代法和二分法

牛顿法:

public double sqr(double n){
        double x=n,y=0.0;
        while (Math.abs(x-y)>0.0001){
            y=x;
            x=(x+n/x)/2;
        }
        return x;
    }

设r是f(x) = 0的根,选取x0作为r初始近似值,过点(x0,f(x0))做曲线y = f(x)的切线L,L的方程为y = f(x0)+f‘(x0)(x-x0),求出L与x轴交点的横坐标 x1 = x0-f(x0)/f‘(x0),称x1为r的一次近似值。

过点(x1,f(x1))做曲线y = f(x)的切线,并求该切线与x轴交点的横坐标 x2 = x1-f(x1)/f‘(x1),称x2为r的二次近似值。重复以上过程,得r的近似值序列,其中x(n+1)=x(n)-f(x(n))/f‘(x(n)),称为r的n+1次近似值,上式称为牛顿迭代公式。

二分法:

#define eps 0.00001
float SqrtByDichotomy(float n)
{
    if (n < 0)
    {
        return -1.0;
    }
    else
    {
        float low, up, mid, last;
        low = 0, up = (n>=1?n:1);
        mid = (low + up) / 2;
        do {
            if (mid*mid>n)
                up = mid;
            else
                low = mid;
            mid = (up+low)/2;
        } while (fabsf(mid - last) > eps);
        return mid;

    }

原文地址:https://www.cnblogs.com/10zhang/p/10557649.html

时间: 2024-11-08 11:59:20

求平方根的算法 牛顿迭代法和二分法的相关文章

求方程解,牛顿迭代和二分

牛顿迭代 #include<iostream> #include<string.h> #include<math.h> using namespace std; float f(float x){ return (pow(x,3)-5*pow(x,2)+16*x+80); } float f1(float x){ return (3*pow(x,2)-5*x+16); } int main(){ // x*x*x-5*x*x+16*x+80; float x=1,x1,

求平方根关键算法

求n以内(不包括n)同时能被3和7整除的所有自然数之和的平方根s,然后将结果s输出.例如若n为1000时,则s=153.909064. 注意:使用循环语句结构实现. ②n由键盘输入,且100 ≤ n ≤10000. package seven; import java.util.Scanner; public class two { public static void main(String[] args) { Scanner sc=new Scanner(System.in); System

利用牛顿迭代法求平方根

求n的平方根,先如果一推測值X0 = 1,然后依据下面公式求出X1,再将X1代入公式右边,继续求出X2…通过有效次迭代后就可以求出n的平方根,Xk+1 先让我们来验证下这个巧妙的方法准确性,来算下2的平方根 (Computed by Mathomatic) 1-> x_new = ( x_old + y/x_old )/2 y (x_old + -----) x_old #1: x_new = --------------- 2 1-> calculate x_old 1 Enter y: 2

24.用牛顿迭代法求平方根

假设a.欲求a的平方根,首先猜测一个值X1=a/2,然后根据迭代公式X(n+1)=(Xn+a/Xn)/2,算出X2,再将X2代公式的右边算出X3等等,直到连续两次算出的Xn和X(n+1)的差的绝对值小于某个值,即认为找到了精确的平方根.例算步骤如下. 1.假设求6的平方根,当Xn和X(n+1)的差值小于0.001时,可以认为已经找到了精确值. 2.根据牛顿迭代法的步骤,首先猜测一个值X1,猜测X1=6/2=3. 求6的平方根:6/2=3 3.将X1=3代入公式X(n+1)=(Xn+a/Xn)/2

求平方根算法 Heron’s algorithm

求平方根问题 概述:本文介绍一个古老但是高效的求平方根的算法及其python实现,分析它为什么可以快速求解,并说明它为何就是牛顿迭代法的特例. 问题:求一个正实数的平方根. 给定正实数 \(m\),如何求其平方根\(\sqrt{m}\)? 你可能记住了一些完全平方数的平方根,比如\(4, 9, 16, 144, 256\)等.那其它数字(比如\(105.6\))的平方根怎么求呢? 实际上,正实数的平方根有很多算法可以求.这里介绍一个最早可能是巴比伦人发现的算法,叫做Heron's algorit

速求平方根倒数

在游戏3D建模方面很多时候要用到求平方根的倒数,而本文章打算介绍的算法会比正常算法快上4倍左右.这对于产品性能将是一个大幅度的提高. 那我们要从哪里开始呢?首先不得不提一提 idsoftware.这是一个创建之初只有13个人的小公司,但它推出的毁灭战士(DOOM)系列游戏可以说改变了游戏世界,极大地推动了游戏产业的发展,因为在当时贫瘠的电脑性能的支撑下,一个开发者能够在游戏中加入一段流畅的动画都会让人惊叹不已,而我们所说的idsoftware,在那个年代就已经做出了画面远超同代其余作品的游戏,像

Leetcode 69. Sqrt(x) 解题报告【C库函数sqrt(x)模拟-求平方根】

69. Sqrt(x) Total Accepted: 93296 Total Submissions: 368340 Difficulty: Medium 提交网址: https://leetcode.com/problems/sqrtx/ Implement int sqrt(int x). Compute and return the square root of x. 分析: 解法1:牛顿迭代法(牛顿切线法) Newton's Method(牛顿切线法)是由艾萨克·牛顿在<流数法>(M

求平方根——分治法

题目描述: 不用sqrt(x)库函数,实现求平方根. 解题思路: 采用二分法 假定要求数num的平方根,那么首先取1~num之间的中点mid. 若 mid * mid > num,那么 根在 1~mid-1之间: 若 mid * mid < num,那么根在 mid+1~num 之间: 若 mid * mid == num,直接输出 mid: 由于整数int求平方根是向下取整,所以,若mid * mid < x情况下,根可能是mid. 依据上面假设根在mid+1~num之间,那么mid+

经典算法:牛顿迭代法求平方根

//牛顿迭代法求平方根 1 double mysqrt(double num) 2 { 3 double x = num/2; 4 double y = 0; 5 do{ 6 x = x/2+num/(2*x); 7 y = x*x-num; 8 if(y<0) y = -y; 9 } while(y>0.0001); 10 return x; 11 } 12 int main(int argc, char* argv[]) 13 { 14 printf("%.3f",my