编写matlab程序，计算简支梁的参数为：长0.6m，矩形截面尺寸为0.01×0.01m，材料为 钢计算当梁受一10N简谐力、作用频率1-1000Hz、作用点在距离一端0.15m处时的强迫振动，获取梁中心点处的振动加速度响应。

好的，我会为您编写一个 MATLAB 程序，用于计算简支梁的参数并计算其在给定条件下的振动加速度响应。请确保您已经安装了 MATLAB 软件。

程序如下：

% 定义常量
L = 0.6; % 梁的长度，单位：m
b = 0.01; % 梁的宽度，单位：m
h = 0.01; % 梁的高度，单位：m
E = 2.1e11; % 梁的弹性模量，单位：Pa
rho = 7800; % 梁的密度，单位：kg/m^3
f = 1:1000; % 振动频率范围，单位：Hz
F = 10; % 外力，单位：N
x = 0.15; % 外力作用点距离梁左端的长度，单位：m

% 计算惯性矩和刚度
I = b*h^3/12; % 惯性矩，单位：m^4
k = 3*E*I/L^3; % 刚度，单位：N/m

% 计算阻尼比和阻尼系数
zeta = 0.01; % 阻尼比
c = 2*zeta*sqrt(k*rho*I); % 阻尼系数，单位：N*s/m

% 计算振动加速度响应
omega = 2*pi*f; % 角频率
X = F./(k-i*c*omega-(omega.^2)*rho*I); % 位移响应函数
A = -omega.^2.*X; % 加速度响应，单位：m/s^2

% 绘制图形
plot(f,A);
xlabel('频率（Hz）');
ylabel('振动加速度（m/s^2）');
title('简支梁在外力作用下的振动加速度响应');

程序说明：

程序首先定义了简支梁的各种参数，包括长度、宽度、高度、弹性模量、密度等。然后定义了振动频率范围和外力的大小，以及外力作用点距离梁左端的长度。

接下来，程序计算了简支梁的惯性矩和刚度，以及阻尼比和阻尼系数。

最后，程序根据振动频率范围计算了梁中心点处的振动加速度响应，并绘制了图形。

您可以将以上代码复制到 MATLAB 软件中执行，即可得到简支梁在给定条件下的振动加速度响应图形。