| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 3540 人关注过本帖
标题:编程.用牛顿方法计算浮点数的平方根.请纠错
只看楼主 加入收藏
天使的恶魔
Rank: 1
来 自:广东
等 级:新手上路
帖 子:26
专家分:0
注 册:2010-9-14
结帖率:83.33%
收藏
 问题点数:0 回复次数:7 
编程.用牛顿方法计算浮点数的平方根.请纠错
# include <stdio.h>
# include <stdlib.h>
# include <math.h>

int main ()
{
    double x,y,z,t;
    printf("Enter a positive number:");
    scanf("%lf",&x);
    y=1;
   
    do {
        z=x/y;               
        t=(y+(x/y))/2;                    //平均數
        if ( fabs(y-t)>0.0001*y)
        y=t;
    } while ( fabs(y-t)>0.0001*y);
   
    printf("Square root:%f",y);
                                 
    system("pause");
    return 0;
}        

运行不出来,哪里错了


搜索更多相关主题的帖子: 平方根 牛顿 点数 纠错 
2010-10-24 15:37
天使的恶魔
Rank: 1
来 自:广东
等 级:新手上路
帖 子:26
专家分:0
注 册:2010-9-14
收藏
得分:0 
这道题的题目:
设x是用户输入的数。牛顿方法需要先给出x平方根的猜测值y(我们使用1)。后续的猜测值通过计算y和x/y的平均值得到。
x       y              x/y              y和x/y的平均数
3       1               3                        2
3       2               1.5                      1.75
3       1.75            1.71429                  1.73214
3       1.73214         1.73196                  1.73205
3       1.73205         1.73205                  1.73205

当y的新旧值之差的绝对值小于0.0001和y的乘积时程序终止。(用fabs函数求double类型数值的绝对值)

记住我最好的方式就是忘了我
2010-10-24 16:12
beyondyf
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
等 级:贵宾
威 望:103
帖 子:3282
专家分:12654
注 册:2008-1-21
收藏
得分:0 
你代码中公式用的不对,看来楼主还不了解牛顿迭代法。简单推导一下吧。
函数y=f(x)在x0点的切线方程为

y - f(x0) = f'(x0)(x - x0)

在y = 0时,有

x1 = x0 - f(x0)/f'(x0)

要求平方根的话对应的函数 f(x) = x * x - C
C是代求其根的平方数。

f'(x) = 2 * x;

那么迭代公式为

x1 = x0 - (x * x - C) / (2 * x)

整理一下

x1 = (x0 + C / x0) / 2

这里的C就是你代码里的 x

下面是完整的示例:

#include <stdio.h>
#define AC  0.0000001
double MySqrt(double a)
{
    double x = 1;
    double t = x * x - a;
    while(t > AC || t < -AC)
    {
        x = (x + a / x) / 2;
        t = x * x - a;
    }
    return x;
}

int main()
{
    double a, sa;
    printf("Enter a positive number:");
    scanf("%lf", &a);
    sa = MySqrt(a);
    printf("Square root:%f",sa);

}

重剑无锋,大巧不工
2010-10-24 17:26
m21wo
Rank: 10Rank: 10Rank: 10
等 级:青峰侠
威 望:4
帖 子:440
专家分:1905
注 册:2010-9-23
收藏
得分:0 
想过没开根号 精度没!譬如精确到500位!!

If You Want Something, Go Get It, Period.
2010-10-24 18:16
天使的恶魔
Rank: 1
来 自:广东
等 级:新手上路
帖 子:26
专家分:0
注 册:2010-9-14
收藏
得分:0 
回复 3楼 beyondyf
.......

不明白

记住我最好的方式就是忘了我
2010-10-24 23:03
beyondyf
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
等 级:贵宾
威 望:103
帖 子:3282
专家分:12654
注 册:2008-1-21
收藏
得分:0 
回复 5楼 天使的恶魔
哪里不明白?需要解释哪一方面?

重剑无锋,大巧不工
2010-10-24 23:08
xx342508809
Rank: 2
等 级:论坛游民
帖 子:89
专家分:51
注 册:2010-7-28
收藏
得分:0 
LZ你能不能把赋值y=1放在定义后面
# include <stdio.h>
# include <stdlib.h>
# include <math.h>

int main ()
{
    double x,y,z,t;
    y=1;
    printf("Enter a positive number:");
    scanf("%lf",&x);
    do {
        z=x/y;               
        t=(y+(x/y))/2;                    //平均數
        if ( fabs(y-t)>0.0001*y)
        y=t;
    } while ( fabs(y-t)>0.0001*y);
   
    printf("Square root:%f",y);
                                 
    system("pause");
    return 0;
}        
2010-10-25 11:29
lovesunshine
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2013-10-11
收藏
得分:0 
牛顿法求浮点数的平方根,完全根据牛顿法来的
#include <stdio.h>
#include<math.h>
#define AC  0.0000001
double MySqrt(double a)
{
     double x = 1;
     double t = x * x - a;
     while(fabs(t)>AC)
     {
         x=x-t/(2*x);
         t = x * x - a;
     }
     return x;
}
int main()
{
     double a, sa;
     printf("Enter a positive number:");
     scanf("%lf", &a);
     sa = MySqrt(a);
     printf("Square root:%f\n",sa);
     return 0;
 }
2013-10-11 16:50
快速回复:编程.用牛顿方法计算浮点数的平方根.请纠错
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.038973 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved