Question

怎样用MATLAB批量求解一元三次方程x^3+a1*x^2+a2*x+a3=0，三个系数已知，但是是一系列值，请问怎样求解

如题，相当于很多个方程，谢谢！

Answer 1

>> a=magic(3);
X=solve('x^3+a1*x^2+a2*x+a3=0','x'宴缺者晌薯扮简);
X=inline(X);
a1=a(1,:);a2=a(2,:);a3=a(3,:);
root=feval(X,a1,a2,a3)

root =

-0.1615 - 0.7035i -1.5421 -0.4384 - 0.0000i
-0.1615 + 0.7035i 0.2711 - 2.4005i -1.0000 + 0.0000i
-7.6771 - 0.0000i 0.2711 + 2.4005i -4.5616

追问

能解释一下吗，我QQ：1402266560。我的悬赏加倍，多谢！

追答

a=magic(3);%3阶幻方 生成系数矩阵
X=solve('x^3+a1*x^2+a2*x+a3=0','x');%解方程
X=inline(X);%
a1=a(1,:);a2=a(2,:);a3=a(3,:);%系数赋值
root=feval(X,a1,a2,a3)%带入求根公式

Answer 2

新版本的matlab已经不建议使用inline函数了，因为其对于较复杂的表差晌达式效率较虚滚锋低，建议备腊用匿名函数代替：
a=magic(3);
f=@(a1,a2,a3)eval(solve('x^3+a1*x^2+a2*x+a3'));
a1=a(1,:);a2=a(2,:);a3=a(3,:);
root=f(a(1,:),a(2,:),a(3,:))

root =

-0.1615 - 0.7035i -1.5421 -0.4384 - 0.0000i
-0.1615 + 0.7035i 0.2711 - 2.4005i -1.0000 + 0.0000i
-7.6771 - 0.0000i 0.2711 + 2.4005i -4.5616

Answer 3

可用偏微分方法化成线性矩阵,然后在按线性矩阵求解就行了