[討論] solve後電腦會當掉
各位神人好
想請大家幫我看一下程式
最後建出來的方程式沒有錯誤
但求解會一直當在那裏
我猜想可能是方程式比較複雜
不知大家有沒有建議的方法
程式如下
------------------------------------------------
Acz=167.7598;
Acy=189.1289;
Sz=-7.6482;
Sy=210.4124;
R=[166.5252 165.2098 168.6954];
d=1;
syms x;
length_SE=abs(Acy-Sy-(Sz-Acz)*x/(1-x^2)^0.5);
length_AE=((Acz-Sz+R(d)*(1-x^2)^0.5)^2+((Sz-Acz)*x/(1-x^2)^0.5-R(d)*x)^2)^0.5;
length_SA=(Acz(d)-Sz+R(d)*(1-x^2)^0.5)^2+(Acy-Sy-R(d)*x^2)^0.5;
alpha=acos((length_AE^2+length_SA^2-length_SE^2)/(2*length_AE*length_SA));
length_AB=((Acz-Bz+R(d)*(1-x^2)^0.5)^2+(Acz-By-R(d)*x)^2)^0.5;
length_AF=((Acz-Bz+R(d)*(1-x^2)^0.5)^2+((Bz-Acz)*x/(1-x^2)^0.5-R(d)*x)^2)^0.5;
length_BF=abs(Acy-By-(Bz-Acz)*x/(1-x^2)^0.5);
beta=acos((length_AF^2+length_AB^2-length_BF^2)/(2*length_AF*length_AB));
eq=subs('sin(alpha)-n*sin(beta)');
sol=solve(eq);
------------------------------------------------
我有試過使用最佳化的方法
如下
opts=optimest('Display','iter');
fzero(eq,3,opts);
但解最佳化時它會跑出
??? Undefined function or method 'optimest' for input arguments of type
'char'.
似乎是建方程式時要先給定未知數前的參數數值
而不能使用讀取的方式
希望大家能給我點方向
感謝!!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.114.57.216
→
04/05 17:52, , 1F
04/05 17:52, 1F
→
04/06 23:08, , 2F
04/06 23:08, 2F
討論串 (同標題文章)