程序：
A = sym('[s, -1, 0, 0;-0.8832*k1, s+0.0883-0.8832*k2, -0.8832*k3-0.6293, -0.8832*k4;0, 0, s, -1;-2.3566*k1, 0.2357-2.3566*k2, -27.8285-2.3566*k3, s-2.3566*k4]')
d = det(A)
说明：第一句是定义符号矩阵，第二句是计算符号矩阵的行列式。
结果：行列式为d =s^4-2.3566*s^3*k4-27.8285*s^2-2.3566*s^2*k3+.883e-1*s^3+.8246e-4*s^2*k4-2.30893054*s+.8246e-4*s*k3-.8832*k2*s^3+23.09512282*s*k2-.8832*k1*s^2+23.09512282*k1。
