Question

matlab解方程组应该有4组解，但只解出一组。求解！！！

syms f1 f2 a t;
s=solve('f1 =5196.4+9081.3*cos(a)+10919*cos(t)*cos(a)-6062*sin(t)*sin(a)','f2 =9081.3*sin(a)+10919*cos(t)*sin(a)+6062*sin(t)*cos(a)','f1^2/11282881-f2^2/15720000=1','(-10919*sin(t)*sin(a)+6062*cos(t)*cos(a))/(-10919*sin(t)*cos(a)-6062*cos(t)*sin(a))=(15720000*f1)/(11282881*f2)');
s=[s.f1 s.f2 s.a s.t]
这是求椭圆与双曲线相切的问题，明显有四组解，但解出来的只有一组，求高人指点。

Answer 1

可以得到四组解。
>> syms f1 f2 a t;s=solve('f1 =5196.4+9081.3*cos(a)+10919*cos(t)*cos(a)-6062*sin(t)*sin(a)','f2 =9081.3*sin(a)+10919*cos(t)*sin(a)+6062*sin(t)*cos(a)','f1^2/11282881-f2^2/15720000=1','(-10919*sin(t)*sin(a)+6062*cos(t)*cos(a))/(-10919*sin(t)*cos(a)-6062*cos(t)*sin(a))=(15720000*f1)/(11282881*f2)');s=[s.f1 s.f2 s.a s.t]
s =
[ -5099.4606699720349415169380782717, -4528.9136572364426638675729120767, 475.36365690426645418837220559876, 29.88139092416406081904881270519]

追问

不是吧 分别是 f1 f2 a t的值吧 是一组的吧。

追答

噢，是一组解。问题可能出在a、t变量上?

追问

什么问题 可否赐教

追答

因为 椭圆与双曲线相切点位置不仅仅只限于一个象限，它有四个点分别 在四个象限位置上。所以a、t有四个，那么f1 f2 也有四个。

追问

我擦，大哥 我不是让你来给我解释为啥有四个切点 我是让你解释为啥解不出四个解。我是初学matlab，数学功底自认为还是可以的。。。
额 算了 还是谢谢你回答了好几次 但真心说 你的回答对我没一点用额。