cordic.m

来自「MATLAB6P5下的GUI设计－实时参数二阶时域曲线的自动绘制.rar」· M 代码 · 共 34 行

a=5,b=8
%z=a+bi
over=0
for n=0:8
    k(n+1)=1/2^n
end
if(b>0)
    sign=-1
else
    sign=1
end
ph=sign*pi/2
temp=b*(-1*sign)
b=a*sign
a=temp          %旋转90度
for n=0:7
   if(b>0)
    sign=-1
   else
    sign=1
   end
 %k=sign*1/(2^n)
 temp=a-sign*b*k(n+1)
 b=sign*a*k(n+1)+b
 a=temp
 ph=ph+sign*atan(k(n+1))
    if(abs(temp)>=32768 | abs(b)>=32768)        %检查溢出
        over=1
    end
end
absz=a/1.6467
ph=(-1)*ph