[求助]曲线拟合的问题

[心际依然]

已知函数x=exp(q/tan(B))*cos(q-qa)+pa*sin(q-qa+B);
y=exp(q/tan(B))*sin(q-qa)-pa*cos(q-qa+B);
其中只有q是变量，且q=[1.3,1.5];B,qa,pa都是常量;
现在要把x,y拟合为以下的函数形式:
x=exp(q/tan(B+tB))*cos(q+tq);
y=exp(q/tan(B+tB))*sin(q+tq);
怎么求出tB和tq呢?
如果用 p = curvefit（‘Fun’p0，xdata，ydata）这个,具体要怎么写.
我的x,y都是隐函数.

[watcher]

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)%线性拟合斜率

不知道可不可以，呵呵

[心际依然]

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

[watcher]

right
:smile:

[心际依然]

这样求出来误差很大，2条曲线差很多呀。好像是累积了很多误差

[zhanhuichong]

使用lsqcurvefit进行最小二乘拟合