matlab计算悬臂梁前两阶频率的代码

以下是使用MATLAB计算悬臂梁前两阶频率的代码：

%% 定义参数
L = 1; % 梁的长度
b = 0.1; % 梁的宽度
h = 0.2; % 梁的高度
E = 2e11; % 梁的弹性模量
rho = 7800; % 梁的密度
I = b*h^3/12; % 梁的截面惯性矩
A = b*h; % 梁的截面面积

%% 计算前两阶频率
syms x
f(x) = E*I*x^4/L^3 - rho*A*x^2/L;
f1 = diff(f,x);
f2 = diff(f1,x);
w = sqrt(f2/(rho*A));
w1 = double(subs(w,x,1));
w2 = double(subs(w,x,2));

disp(['第一阶频率：',num2str(w1),' Hz']);
disp(['第二阶频率：',num2str(w2),' Hz']);

运行以上代码即可得到悬臂梁的前两阶频率，结果会在命令窗口中显示出来。