[mathjiang]

表达式
D=solve('S^3-I1*S^2+I2*S-I3')
中，原I1，I2，I3的值并没有传入方程
S^3-I1*S^2+I2*S-I3=0
因为solve是符号求解器。
破解之法是求数值解：
clear all
clc
s=eye(3);
%计算应力不变量
I1=trace(s);%计算第一应力不变量I1
I2=det(s(2:3,2:3))+det(s(1:2:3,1:2:3))+det(s(1:2,1:2));%计算I2
I3=det(s);%计算I3
%解三次方程，求主应力。
c = [1 -I1 I2 -I3];%这就是三次多项式 S^3-I1*S^2+I2*S-I3
D=roots(c);