如何在没有eig的情况下在matlab中计算特征向量

我有一个矩阵，我需要得到特征向量。我已经计算了特征值，假设我们有特征值，我写了这个

  for i = 1:length(c)
  syms y
  cal_vec = (c-eig_Val(i)*I)*y == 0;
  eigVec(:,i) = double(solve(cal_vec,y));
end

---

现在我得到了0作为y，但是我需要得到y1和y2，你要找的方法是，计算c-eig_Val（I）*I的零空间

---

请注意，如果您的特征值具有多个特征向量，或者如果您的特征值（足够）不准确，那么它将无法开箱即用。

请创建一个“如何计算零空间？”