%x,y表示直角坐标,a,b表示极坐标
function [a,b]=trans(x,y)
a = sqrt(x^2+y^2);
if x>0
b = arsin(y/a);
if b