无法从Matlabs一个答案解决功能

Question

当我尝试使用带有我收到以下错误消息，我的预定义功能和载体解决：

警告：显式的解决方案不能被发现。
在169
在cardanAngles解决在33

我有三个未知变量和三个方程式。我应该能够找到这个答案。什么可能是错的？

V0 = transpose([-0.7969, 0.1778, 0.5774]); 
V = [-0.7976 , 0.2282 , 0.5583]; 

syms a1 a2 a3 


F1 = [ cos(a3 * pi/180).*cos(a2 * pi/180), cos(a3 * pi/180).*sin(a2 * pi/180).*sin(a1 * pi/180)+sin(a3 * pi/180).*cos(a1 * pi/180), -cos(a3 * pi/180).*sin(a2 * pi/180).*cos(a1 * pi/180)+sin(a3 * pi/180).*sin(a1 * pi/180)]; 
F2 = [-sin(a3 * pi/180).*cos(a2 * pi/180), -sin(a3 * pi/180).*sin(a2 * pi/180).*sin(a1 * pi/180)+cos(a3 * pi/180).*cos(a1 * pi/180), sin(a3 * pi/180).*sin(a2 * pi/180).*cos(a1 * pi/180)+cos(a3 * pi/180).*sin(a1 * pi/180)]; 
F3 = [ sin(a2 * pi/180)     , -cos(a2 * pi/180).*sin(a1 * pi/180)               , cos(a2 * pi/180).*cos(a1 * pi/180)               ]; 


[a1, a2,a3] = solve(V(1)==F1*V0,V(2)==F2*V0, V(3)==F3*V0, a1,a1,a3); 

Answer 1

如果有解决方案，您可以从解决最后的eq。即F3*V0-V(3)=0。 情商写着：

F3 = [sin(a2), -cos(a2).*sin(a1), cos(a2).*cos(a1)]; 
V0 = [-0.7969, 0.1778, 0.5774]'; 

F=V0(1)*sin(a2) - V0(2)*cos(a2)*sin(a1) + V0(3)*cos(a1)*cos(a2) - V(3); 

现在我们需要寻找一个解决方案家族的一些地方F=0，快速分析会产生这样的：

现在你可以使用这些作为如注释中所建议的fsolve程序的初始猜测值。可以按摩了一下原来的形式，写了以下功能：

function F = qfun(x) 

V0 = transpose([-0.7969, 0.1778, 0.5774]); 
V = [-0.7976 , 0.2282 , 0.5583]; 
F1 = [ cos(x(3)).*cos(x(2)), cos(x(3)).*sin(x(2)).*sin(x(1))+sin(x(3)).*cos(x(1)), -cos(x(3)).*sin(x(2)).*cos(x(1))+sin(x(3)).*sin(x(1))]; 
F2 = [-sin(x(3)).*cos(x(2)), -sin(x(3)).*sin(x(2)).*sin(x(1))+cos(x(3)).*cos(x(1)), sin(x(3)).*sin(x(2)).*cos(x(1))+cos(x(3)).*sin(x(1))]; 
F3 = [ sin(x(2))     , -cos(x(2)).*sin(x(1))               , cos(x(2)).*cos(x(1))               ]; 
F=(F1*V0-V(1))+abs(F2*V0-V(2))+abs(F3*V0-V(3)); 

，然后在一个单独的文件运行类似：

options = optimset('Display','off'); 
for ... 
x0 = [a1_found a2_found 0]; % guess values 
[x,Fval,exitflag] = fsolve(@qfun,x0,options); 
end... 

每次猜测值向量将输出最近的解决方案，因此循环对角度值进行很好的采样以得到a3。