已知质量矩阵和刚度矩阵数据的两个excel文件，动能法matlab代码

以下是使用MATLAB进行动能法计算的示例代码：

% 从Excel文件中读取质量矩阵和刚度矩阵数据
massFile = '质量矩阵.xlsx';  % 替换为质量矩阵的Excel文件路径
stiffnessFile = '刚度矩阵.xlsx';  % 替换为刚度矩阵的Excel文件路径
massMatrix = xlsread(massFile);
stiffnessMatrix = xlsread(stiffnessFile);

% 计算质量矩阵的逆矩阵
invMassMatrix = inv(massMatrix);

% 计算刚度矩阵的特征值和特征向量
[eigVec, eigVal] = eig(stiffnessMatrix);

% 计算模态质量
modalMass = eigVec' * massMatrix * eigVec;

% 计算模态动能
modalKineticEnergy = diag(eigVec' * invMassMatrix * eigVec);

% 打印每个模态的模态质量和模态动能
for i = 1:size(eigVal, 1)
    fprintf('Mode %d: Modal Mass = %.4f, Modal Kinetic Energy = %.4f\n', i, modalMass(i, i), modalKineticEnergy(i));
end

在上述代码中，我们假设质量矩阵数据存储在名为"质量矩阵.xlsx"的Excel文件中，刚度矩阵数据存储在名为"刚度矩阵.xlsx"的Excel文件中。您需要根据实际情况替换这两个文件的路径。

首先，我们使用`inv`函数计算质量矩阵的逆矩阵。然后，使用`eig`函数计算刚度矩阵的特征值和特征向量。接下来，通过矩阵乘法和转置操作，计算模态质量和模态动能。最后，使用循环遍历打印每个模态的模态质量和模态动能。

请确保您的Excel文件中的数据与代码中指定的文件名和格式相匹配，并且包含正确的质量矩阵和刚度矩阵数据。

如果您有任何其他问题，请随时提问。