开方运算的DSP实现

(2014-04-07 10:59:36)
sg_trans.gif转载▼
标签：

c

dsp

分类：DSP

  1. //===============================================

  3. //函数名：VSqrt3

  4. //功能： 实现对32位定点数的开方

  5. //性能： 60M主频28015硬件下运行时间小于10us

  6. //转自：http://read.pudn.com/downloads180/sourcecode/mpi/840129/sqrt.c__.htm

  8. unsigned longVSqrt3(unsigned longx )

  9. {

  10. unsigned longx1;

  11. ints=1;

  12. unsigned longg0,g1;

  14. if(x<=1) returnx;

  16. x1=x-1;

  17. if(x1>65535)

  18. {

  19. s+=8;

  20. x1>>=16;

  21. }

  22. if(x1>255)

  23. {

  24. s+=4;

  25. x1>>=8;

  26. }

  27. if(x1>15)

  28. {

  29. s+=2;

  30. x1>>=4;

  31. }

  32. if(x1>3)

  33. {

  34. s+=1;

  35. }

  36. g0=1;

  37. g0=g0<<s;

  38. g1 =(g0 +(x>>s))>>1;

  40. while(g1

  41. {

  42. g0=g1;

  43. g1=(g0+x/g0)>>1;

  44. }

  45. returng0;

  46. }

  48. //===============================================