用matlab求解振型，质量矩阵M是8行8列的矩阵，刚度矩阵K是8行8列的矩阵，零件模型是一个梁，求解前三阶振型

在MATLAB中，计算结构（如梁）的前三阶振型通常涉及到动力学分析中的模态分析。这个过程一般包括以下几个步骤：

1. **定义矩阵**：
- 质量矩阵 \( M \) 和刚度矩阵 \( K \) 分别表示系统的质量和弹性性质。

2. **构建系统矩阵**：
- 结合 \( M \) 和 \( K \)，形成系统矩阵 \( A = K - M \cdot \omega^2 \)，其中 \( \omega^2 \) 是角频率的平方，对于简谐振动，它可以用0代替，因为初始条件通常是静止状态。

3. **特征值和特征向量计算**：
- 使用 `eig(A)` 函数找到矩阵 \( A \) 的特征值（即固有频率的平方）和对应的特征向量（代表振型），`[V,D] = eig(A)`，其中 \( V \) 是振型矩阵，\( D \) 是对角矩阵，对角线元素是特征值。

4. **前几阶振型**：
- 挑选前三大的特征值（对应的是低阶振型），然后取对应的正交特征向量作为前三个振型。可以使用 `sort(D, 'descend')` 来排序并选取前几个。

5. **振型表达式**：
- 每一阶振型 \( X_n \) 可以通过 \( X_n = V(:,n) \) 计算得到，这里的 `:` 表示所有行，`n` 是阶数。

[V,D] = eig(K-M*0); % 假设静止初始条件
[~,idx] = sort(diag(D),'descend'); % 排序
前三阶振型 = V(:, idx(1:3)); % 提取前三阶振型