已知单自由度系统固有频率wn的公式，编写matlab程序绘制wn随刚度k和质量m的变化规律。 m=1,k取0-1000；k=1，m取0-10

单自由度系统的固有频率 \( w_n \) 可以通过经典物理公式计算，对于弹簧-质量系统，它与质量和弹性系数 (刚度) 成反比，并与质量的平方根成正比，即：

\[ w_n = \sqrt{\frac{k}{m}} \]

在 MATLAB 中，你可以使用 for 循环来创建一个二维图形，展示 \( w_n \) 随着 \( k \) 和 \( m \) 的变化情况。这里是一个简单的示例代码：

% 定义变量范围
m_range = 0:0.1:10; % 质量范围，从0到10，步长0.1
k_range = 0:100:1000; % 刚度范围，从0到1000，步长100

% 初始化固有频率矩阵
wn_matrix = zeros(length(m_range), length(k_range));

% 计算每个 (m, k) 对应的 wn
for i = 1:length(m_range)
    for j = 1:length(k_range)
        wn_matrix(i,j) = sqrt(k_range(j) / m_range(i));
    end
end

% 创建散点图
figure;
surf(k_range, m_range, wn_matrix);
xlabel('刚度 k');
ylabel('质量 m');
zlabel('\( w_n \)');
title('固有频率 \( w_n \) 随 k 和 m 变化');

% 添加网格线和颜色标尺
grid on;
colorbar;

% 提问部分