MATLAB：具有阻尼项的齐次二阶方程组的特征值分析

Question

我有一个动态方程系统，最终可以用着名的弹簧质量阻尼器＃34;形式：

[M]{q''}+[C]{q'}+[K]{q}={0}

[M], [C], [K]: n-by-n Coefficient Matrices
{q}: n-by-1 Vector of the Degrees of Freedom
(the ' mark represents a time derivative)

我想找到这个系统的特征值和特征向量。显然由于术语[C] {q＆＃39;}，标准的MATLAB函数eig（）将没有用。

有没有人知道一个简单的MATLAB例程来确定该系统的特征值，特征向量？系统是同质的，因此有效的特征值分析应该是非常可行的，但我有点挣扎。

显然，我可以使用蛮力和符号计算软件来找到巨大的特征多项式。但这对我来说似乎效率低下，尤其是因为我通过代码的其他部分将其循环以确定频率作为其他变化参数的函数。