x_int=[50:1:250];
y_int=spline(xc,yc,x_int);
.
.
.
.
thedam(1,i)=(y_int(1,i)-ya)./x_int(1,i);
thedam1(1,i)=(L1.^2+(x_int(1,i)-xa).^2+(y_int(1,i)-ya).^2
-L2.^2)/(2*L1*((x_int(1,i)-xa).^2+(y_int(1,i)-ya).^2).^(1/2));
theda1(1,i)=atan(thedam(1,i))+acos(thedam1(1,i));
在做机构设计，要找出杆件的角度
找到一组参数 thedam=2 thedam1=1.00000
可是代入theda1 求角度的时候会有I值，导致polt不出图形
theda1=1.10714871779409 + 2.10734242554470e-08i
因为知道acos里输入的值 介于 1~-1
也试过 atan(2)+acos(1) 出来为theda1的实部值 实在不清楚虚部从哪来的
请问是哪里有问题吗?