作者ax80097 (零度空间)
看板MATLAB
标题[讨论] solve後电脑会当掉
时间Thu Apr 5 14:34:08 2012
各位神人好
想请大家帮我看一下程式
最後建出来的方程式没有错误
但求解会一直当在那里
我猜想可能是方程式比较复杂
不知大家有没有建议的方法
程式如下
------------------------------------------------
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
1F：→ mp19990920:参数设定function名称打错了(optimset) 04/05 17:52
2F：→ maplefog:n是什麽? 04/06 23:08