求解三元一次方程时，​得出的解是一个代数式​而不是数。

2 次查看（过去 30 天）
xebad
xebad 2022-11-18
回答： lefow 2022-11-18
deltax=200;
deltay=-300;
deltaz=50;
taoxy=500;
taoyz=0;
taozx=300;%x,y,z方向的正应力和切应力
I1=deltax+deltay+deltaz;
I2=deltax*deltay+deltay*deltaz+deltaz*deltax-taoxy^2-taoyz^2-taozx^2;
I3=deltax*deltay*deltaz-deltaz*taoxy^2-deltax*taoyz^2-deltaz*taozx^2+2*taoxy*taoyz*taozx;%I1，I2，I3为有正应力和切应力组成的系数
syms x;
x = solve(x^3-I1*x^2+I2*x-I3 == 0,'x'),
x1 = 50;
x2 = 50*161^(1/2) - 50;
x3 = - 50*161^(1/2) - 50;%x1,x2,x3为上述方程的解
a1= deltax-x1;a2= deltay-x1;a3=deltaz-x1;%a1,a2,a3上述方程的解组成的系数
[l1,m1,n1]=solve('a1*l1+taoxy*m1+taozx*n1=0','taoxy*l1+a2*m1+taoyz*n1=0','l1^2+m1^2+n1^2=1','l1','m1','n1')%求解l1,m1,n1

请先登录，再回答此问题。

采纳的回答

lefow
lefow 2022-11-18
最后加上：
l1=subs(l1)
m1=subs(m1)
n1=subs(n1)

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!