关于matlab一次求解很多个形式相同的方程组

代码不长,问题可不少:

1、第二行:

定义符号变量应该用syms t,而不是sym t(如果用sym函数,则应该返回参数t=sym(t));

2、第三行:

eq=p1(:).t^3+p2(:).t^2+p3(:).t-x

首先,相乘应该用“*”,而不是“.”;

其次,如果你把方程组写到一起,则求解时会被认为联立求解。我猜测你的1000个方程应该是独立求解的,这与你的希望不符。

3、如果想求数值解,可考虑用roots,速度比起求符号解相差无数倍。只需把循环体里面的内容改成下面的代码即可:

T(i,:)=roots([p1(i) p2(i) p3(i) -x]);

以下是改后的参考代码(每个方程的根保存为一行):

p1=rand(1000,1); p2=rand(1000,1); p3=rand(1000,1); x=0.1;

T = sym(zeros(1000,3));

for i=1:length(p1)

syms t

eq=p1(i)*t^3+p2(i)*t^2+p3(i)*t-x;

T(i,:)=solve(eq).';

end