imax=find(max(y)==y);%% 这一步求出来的imax有两个值，是个1X2的矩阵[1 50]theta=atan2(fy(imax),fx(imax));%%所以这一步求出来的theta也是个1X2的矩阵yy(i)=theta;%%这一步等号两侧维度不匹配，所以出错%%换成atan也不行啊，atan输入变量只能有一个呀？
