matlab计算悬臂梁前两阶频率的代码并备注各物理量单位且不使用syms函数

以下是MATLAB计算悬臂梁前两阶频率的代码，并且已注明各物理量单位，不使用syms函数：

clear all;

% 物理量输入
L = 1; % 梁长，单位：m
b = 0.02; % 梁宽，单位：m
h = 0.004; % 梁高，单位：m
E = 2.1e11; % 弹性模量，单位：Pa
rho = 7850; % 材料密度，单位：kg/m^3

% 计算相关物理量
I = b * h^3 / 12; % 惯性矩，单位：m^4
A = b * h; % 梁截面积，单位：m^2
m = rho * A * L; % 梁质量，单位：kg

% 求解悬臂梁前两阶固有频率
f1 = 1 / (2 * pi) * sqrt(E * I / (m * L^3)); % 第一阶固有频率，单位：Hz
f2 = 1 / (2 * pi) * sqrt(3 * E * I / (m * L^3)); % 第二阶固有频率，单位：Hz

% 显示结果
fprintf('第一阶固有频率为：%.2f Hz\n', f1);
fprintf('第二阶固有频率为：%.2f Hz\n', f2);

其中，E为弹性模量，I为惯性矩，A为梁截面积，m为梁质量，f1和f2分别为第一阶和第二阶固有频率。单位分别为：

- L：m
- b：m
- h：m
- E：Pa
- rho：kg/m^3
- I：m^4
- A：m^2
- m：kg
- f1：Hz
- f2：Hz