matlab传递矩阵法求临界转速

MATLAB传递矩阵法可以用于求解临界转速。具体步骤如下：

1. 定义转子和轴承系统的动力学方程。
2. 将动力学方程转换为矩阵形式，包括质量矩阵、刚度矩阵和阻尼矩阵。
3. 将矩阵形式的方程组简化为特征值问题，即求解特征值和特征向量。
4. 根据特征值和特征向量计算出临界转速。

下面是一个简单的示例代码：

% 定义系统参数
m = 1; % 转子质量
k = 10; % 转子刚度
c = 0.1; % 转子阻尼
J = 0.5; % 转子惯量
omega = 0:0.1:50; % 转速范围

% 定义矩阵
M = [m 0; 0 J]; % 质量矩阵
K = [k -k; -k k]; % 刚度矩阵
C = [c -c; -c c]; % 阻尼矩阵

% 求解特征值和特征向量
[V,D] = eig(inv(M)*K); % D为特征值矩阵，V为特征向量矩阵

% 计算临界转速
wc = sqrt(D(1,1)); % 第一个特征值对应的临界转速

% 绘制临界转速曲线
plot(omega,sqrt(omega.^2-wc^2),'r');
xlabel('转速');
ylabel('非稳定运动幅值');
title('临界转速曲线');

此代码可以绘制临界转速曲线，其中x轴为转速，y轴为非稳定运动幅值。曲线上方表示系统处于稳定状态，曲线下方表示系统处于非稳定状态。