请各位高手帮我看一下程序(定点DSP中求浮点数的平方根--使用牛顿迭代法)

2019-03-26 16:36发布

以下main.c 程序是小弟写的,不使用中断。编译,链接都没问题,就是输入数值后,输出一直为0,感觉是浮点转定点后在Sqrt函数里面运算时出了问题了,但不知道问题出在哪,小弟是菜鸟,请各位大侠帮帮忙,多谢多谢!

#include <stdio.h>
#include <stdlib.h>
#include <std.h>
#include <math.h>

long Sqrt(long);
//float fabs(float);

float Num;   //输入浮点值
long temp; //输入值转化为定点值
double result; //求得的平方根


void main()
{

   printf("请输入一个数:");
   scanf("%f",&Num);
  
   temp=(long)(Num*(1<<15)); //浮点数转为定点数,Q15

   //检查输入的数是否能开方
   if(temp<0)
   {printf("输入错误,负数不能开方,请输入一个非负数 ");
    return ;
   }
   
   result=(double)(Sqrt(temp)/(1<<15)); //计算得到平方根

   printf("输入的数%f的平方根为%f ",Num,result);       

   printf("平方根计算完毕。");
}

/*牛顿迭代法求非负实数x的平方根*/
long Sqrt(long x)
{

long x0;
long x1=x/2;

float error=(float)((1e-5)*(1<<15));  //误差阈值:0.00001

do
{
x0=x1;
x1=(x0+x/x0)/2;
}
while(fabs(x1-x0)>=error);

return(x1);
}

//*求绝对值*/
//float Abs(float a)
//{
// if( a < 0 )
// {
//  a = -a;
// }
// return(a);
//} 此帖出自小平头技术问答
友情提示: 此问题已得到解决,问题已经关闭,关闭后问题禁止继续编辑,回答。