关于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