采用Matlab软件，计算当梁受一10N简谐力、作用频率1-1000Hz、作用点在距离一端0.15m处时的强迫振动，获取梁中心点处的振动加速度响应

下面是使用MATLAB计算梁受一10N简谐力、作用频率1-1000Hz、作用点在距离一端0.15m处时的强迫振动，获取梁中心点处的振动加速度响应的程序：

clear all; close all; clc;

% 梁的基本参数
L = 1; % 梁长
E = 2.1e11; % 弹性模量
I = 1.33e-5; % 截面惯性矩
rho = 7800; % 材料密度
w = 1; % 梁的宽度
h = 0.1; % 梁的厚度

% 计算固有频率
f = sqrt(E*I/rho/w/h^3)/(2*pi); % Hz

% 计算作用点的位置
x = 0.15; % m

% 计算简谐力的振动加速度响应
omega = 2*pi*(1:1000); % 强迫振动的角频率
F = 10; % 简谐力的振幅
a = F/E/I/(omega.^2-f^2).*sin(omega*x/sqrt(E*I/rho/w/h^3)); % 振动加速度响应

% 绘制图像
figure;
plot(omega/(2*pi), abs(a)); % 振动加速度响应的幅值
xlabel('频率 (Hz)');
ylabel('振动加速度 (m/s^2)');
title('梁中心点处的振动加速度响应');

程序中，我们首先根据梁的基本参数计算出固有频率。然后，利用简谐力的振动加速度公式，计算在不同频率下梁中心点处的振动加速度响应。最后，使用 plot 函数绘制振动加速度响应的幅值随频率变化的曲线图。

需要注意的是，在实际应用中，应根据具体情况调整参数和计算方法。例如，如果梁不是简支梁，计算方法会有所不同。