In[]:=
vars={X,Y,Z,M,A,B};(*定义参数-假设所有参数为正实数*)params={ωa>0,ωc>0,κ1>0,λx>0,λy>0};(*方程组*)equations={2ωaX-4λxZA==0,(*E1*)2ωcM-κ1A-16λyYB==0,(*E2*)2ωcA+κ1M+16λxXB==0,(*E3*)-κ1B-4λxXM-4λyYA==0,(*E4*)-2ωaY-4λyZM==0,(*E5*)4λxYA+4λyXM==0,(*E6*)X^2+Y^2+Z^2==1};(*方法1：尝试用Solve直接求解解析解*)solutions=Solve[equations,vars,Reals];Print["解析解："];solutions//Simplify