matlab中eig函数求解特征值与特征向量有五种标准调用方式：一、e=eig(a)得特征值列向量；二、[v,d]=eig(a)满足av=vd；三、[v,d]=eig(a,'nobalance')禁用平衡预处理；四、[v,d]=eig(a,b)解广义特征值问题；五、支持'vector'或'matrix'输出格式指定。
如果您在MATLAB中需要求解方阵的特征值与特征向量，但结果不符合预期或报错，则可能是由于输入矩阵非方阵、调用语法不匹配或未处理复数/广义情形。以下是eig函数多种标准调用方式的具体解析：
一、基础调用：仅获取特征值列向量
该方式适用于只需数值结果、无需特征向量的场景，返回一个包含全部特征值的列向量，顺序不保证单调性，且允许含复数。
1、确保待计算矩阵A为n×n方阵，例如定义 A = [4, 3; 2, 1];
2、执行命令 E = eig(A)；
3、输出变量E为n×1列向量，其元素即为A的所有特征值。
二、标准双输出：同时获取特征向量与对角特征值矩阵
该调用返回正交（或可逆）特征向量矩阵V与对角矩阵D，满足关系式 A*V = V*D，是工程分析中最常用的格式。
1、定义方阵A，如 A = [1, 1, 0.5; 1, 1, 0.25; 0.5, 0.25, 2];
2、执行命令 [V, D] = eig(A)；
3、检查V的每一列是否为单位向量：norm(V(:,1)) 应近似等于1；
4、验证计算正确性：max(abs(A*V - V*D)) 应接近机器精度（如1e-15量级）。
三、禁用平衡预处理：提升病态矩阵精度
默认情况下eig会对输入矩阵先做相似变换以提升数值稳定性，但对某些严重缩放失衡或已平衡的矩阵，该步骤反而引入误差，此时应显式禁用。
1、当A存在极端数量级差异（如含1e-6与1e6元素），且观察到特征值虚部异常偏大时启用本方法；
2、执行命令 [V, D] = eig(A, 'nobalance')；
3、注意：若A为对称矩阵，此选项将被自动忽略，因对称矩阵本身已具备良好平衡性。
四、广义特征值问题：求解AV = BVD形式
该模式用于质量-刚度特征值问题、微分方程离散系统等物理建模场景，其中B通常为正定矩阵（如质量矩阵），A为刚度或系数矩阵。
1、定义同阶方阵A与B，要求B可逆且数值性质良好，例如 K = [2,-1,0; -1,2,-1; 0,-1,2]; M = eye(3);
2、执行命令 [V, D] = eig(K, M)；
3、验证关系式：max(abs(K*V - M*V*D)) 应处于浮点误差范围内；
4、若B为对称正定矩阵，MATLAB默认采用Cholesky分解算法；如需强制使用QZ算法，可追加flag参数：[V, D] = eig(K, M, 'qz')。
五、控制输出形式：向量或对角矩阵格式指定
从R2021a起，eig支持outputForm参数，便于后续排序或索引操作，避免手动diag(D)转换。
1、若需直接获得按列排列的特征值向量（而非对角矩阵），使用：[V, D] = eig(A, 'vector')；
2、若明确需要对角矩阵形式（即使已用[V,D]=eig(A)），可显式指定：[V, D] = eig(A, 'matrix')；
3、该参数必须置于最后位置，不可与' nobalance'或'qz'等其他选项混用。
