1，由q=[1.3,1.5]获得x,y的离散序列
2，由x=exp(q/tan(B+tB))*cos(q+tq);y=exp(q/tan(B+tB))*sin(q+tq);得到
atan(y/x)=q+tq%线性拟合截距
ln(x^2+y^2)=q/tan(B+tB)%线性拟合斜率
你的意思是这样吗？
令x1=q；
y1=atan(y/x)+ln(x^2+y^2)=q+tq+2q/tan(B+tB）=（1+2/tan(B+tB））*x+tq
p=polyfit(x1y11);
求出来的p（1）=1+2/tan(B+tB）
p（2）=tq
|