tan计算

clc;clear;


%%调整精度
r=50;
table=zeros(r,r);

for y=0:r
   for x=0:r
       if(x==0) 
           table(y+1,x+1)=90;
       else
           table(y+1,x+1)=atan2(y,x)*180/pi;
       end 
   end
end


DX=35;
DY=-55;


%%查表计算
dy=abs(DY);
dx=abs(DX);
    
while(1) 
    if dx>=(r+1)||dy>=(r+1)
        dx=floor(dx/2);
        dy=floor(dy/2);
    else
        break;
    end
end
angle1=table(dy+1,dx+1);

%%通过dx和dy的正负判断象限
if(DX>0&&DY<0)
    angle1=360-angle1;
end

if(DX<0&&DY<0)
    angle1=angle1+180;
end

if(DX<0&&DY>0)
    angle1=180-angle1;
end


%调用函数计算，和查表结果比较近似相等
angle2=atan2(DY,DX)*180/pi;
if(angle2<0)
    angle2=angle2+360;
end