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的實部值 實在不清楚虛部從哪來的
請問是哪裡有問題嗎?