需要運行下面的代碼，得到矩陣E，然后才能進行仿真，更具體思路的可以參考 https://blog.csdn.net/XSTX1996/article/details/83117733 A = [1 1; 0 0]; C = [1 0]; syms e1 e2 lambda; E = [e1; e2]; T=0.01; eigPloy = det(A-E*C - lambda*eye(2)); lambdaVal = roots([2*T^2, 2*T, 1]); eigPloy = subs(eigPloy, lambda, lambdaVal); [e1, e2] = solve(eigPloy, e1, e2); E = double([e1 e2])';