如何在MATLAB中求解包含符号的方程组？

我试图求解“p，x，y，z”的四个方程，其中包含两个参数“b”，“a”，我使用了solve命令，如代码所示：

clc,clear;

syms  p x y z b a
 
eqn1 = (1)*p  + (exp(0.5*b^(1/3)))*x + (exp(-0.25*b^(1/3)) * cos(0.25*sqrt(3)*b^(1/3)))*y + (exp(-0.25*b^(1/3)) * sin(0.25*sqrt(3)*b^(1/3)))*z == 0;

eqn2 = (1)*p  + (exp(-0.5*b^(1/3)))*x + (exp(0.25*b^(1/3)) * cos(0.25*sqrt(3)*b^(1/3)))*y - (exp(0.25*b^(1/3)) * sin(0.25*sqrt(3)*b^(1/3)))*z == 0;

eqn3 = (b^(1/3)*exp(0.5*b^(1/3))+a*b^(2/3)*exp(0.5*b^(1/3)))*x + (-.866*b^.333*sin(.433*b^.333)*exp(-.25*b^.333)-.5*cos(.433*b^.333)*b^.333*exp(-.25*b^.333) + a(-.5*b^.666*cos(.433*b^.333)*exp(-.25*b^.333)+.866*b^.666*sin(.433*b^.333)*exp(-.25*b^.333)))*y + (.866*cos(.433*b^.333)*b^.333*exp(-.25*b^.333)-.5*b^.333*sin(.433*b^.333)*exp(-.25*b^.333)+a(-.5*b^.666*sin(.433*b^.333)*exp(-.25*b^.333)-.866*b^.666*cos(.433*b^.333)*exp(-.25*b^.333)))*z == 0;
                                                                 
eqn4 = (b^(1/3)*exp(-0.5*b^(1/3))-a*b^(2/3)*exp(-0.5*b^(1/3)))*x + (.866*b^.333*sin(.433*b^.333)*exp(.25*b^.333)-.5*cos(.433*b^.333)*b^.333*exp(.25*b^.333) - a(-.5*b^.666*cos(.433*b^.333)*exp(.25*b^.333)-.866*b^.666*sin(.433*b^.333)*exp(.25*b^.333)))*y + (.866*cos(.433*b^.333)*b^.333*exp(.25*b^.333)+.5*b^.333*sin(.433*b^.333)*exp(.25*b^.333)-a(.5*b^.666*sin(.433*b^.333)*exp(.25*b^.333)-.866*b^.666*cos(.433*b^.333)*exp(.25*b^.333) ))*z == 0;
                                                                       
sol = solve([eqn1, eqn2, eqn3, eqn4 ], [p, x, y, z]);
pSol = sol.p
xSol = sol.x
ySol = sol.y
zSol = sol.z

，但它不断显示错误，请您对此问题提出任何解决方案。

让我们看一下：

eqn3 = (b^(1/3)*exp(0.5*b^(1/3))+a*b^(2/3)*exp(0.5*b^(1/3)))*x + (-.866*b^.333*sin(.433*b^.333)*exp(-.25*b^.333)-.5*cos(.433*b^.333)*b^.333*exp(-.25*b^.333) + a(-.5*b^.666*cos(.433*b^.333)*exp(-.25*b^.333)+.866*b^.666*sin(.433*b^.333)*exp(-.25*b^.333)))*y + (.866*cos(.433*b^.333)*b^.333*exp(-.25*b^.333)-.5*b^.333*sin(.433*b^.333)*exp(-.25*b^.333)+a(-.5*b^.666*sin(.433*b^.333)*exp(-.25*b^.333)-.866*b^.666*cos(.433*b^.333)*exp(-.25*b^.333)))*z == 0;

我选错了一个零件，一共4个

a(-.5*b^.666

你看，如果

a（index）

存在，

index

必须是正整数。但是，在代码中，

-.5*b^.666*cos（.433*b^.333）*exp（-.25*b^.333）+.866*b^.666*sin（.433*b^.333）*exp（-.25*b^.333）

的值可能不是正整数。A甚至可能不是矩阵，更不用说索引了

---

另外，在你发帖之前，你最好检查一下基本的语法问题以节省彼此的时间。

你会遇到什么错误？@ThomasSablik，这个错误“使用sym/subsref的错误（第49行）使用maplemex的错误错误错误，（在MTM:-subsref中）索引必须是正整数，got-.500000000000000000*exp（.25000000000000000*b^.333000000000000018）*cos（.439999999999996*b^.333000000018）*b^.66600000000000036+.86599999999999992*exp（-.25000000000000000*b^.333000000018）*sin（.43999999999999996*b^.333000000000000018）*b^.66600000000000036“索引必须是正整数”