```matlab
% 定义质量、阻尼比和刚度参数
m = 1; % 质量
c = 0.2; % 阻尼比
k = 100; % 刚度
% 计算自振频率
omega_n = sqrt(k / m); % 无阻尼自振角频率
% 输出自振频率
fprintf('系统的自振频率为：%f rad/s\n', omega_n);
% 如果需要计算阻尼自振频率，可以增加阻尼项
% 阻尼自振角频率 omega_d = omega_n * sqrt(1 - zeta^2)
zeta = c / (2 * sqrt(m * k)); % 阻尼比
omega_d = omega_n * sqrt(1 - zeta^2);
fprintf('考虑阻尼后的自振频率为：%f rad/s\n', omega_d);
```
【考研刷题通】——您的考研刷题好帮手！无论是政治、英语还是数学，我们的考研刷题小程序都能为您提供全面的练习题库和模拟测试，助您高效备考，轻松应对考研挑战！立即加入我们，开启高效刷题之旅！【考研刷题通】小程序，助力考研成功！