采用中心差分法计算其位移响应，并作出前10s内每个质点的时间位移曲线。m1=m2=1，m3=2，c1=c2=c3=1，k1=1,k2=5,k3=8,f1=0,f2=sin(w2*t),f3=0,w1=w2=w3=1,x1=x1'=x2=x2'=x3=x3'=0,matlab代码

好的，根据您的要求，我来给您提供一个MATLAB代码，该代码可以用于计算所述系统的位移响应并绘制前10s内每个质点的时间位移曲线。以下是代码：

% 定义系统参数
m = [1; 1; 2];  % 质量向量
c = [1; 1; 1];  % 阻尼向量
k = [1, -1, 0; -1, 6, -5; 0, -5, 13];  % 刚度矩阵
f = [0; sin(t); 0];  % 外力向量

% 定义初始条件
x0 = [0; 0; 0; 0; 0; 0];  % 初位移和初速度向量

% 定义时间步长和时间向量
dt = 0.01;  % 时间步长
t = 0:dt:10;  % 时间向量

% 初始化位移和速度向量
x = zeros(6, length(t));

% 迭代计算位移和速度
for i = 2:length(t)
    % 计算加速度向量
    a = (f - c.*(x(4:6, i-1)) - k*x(1:3, i-1))./m;
    
    % 计算位移和速度
    x(1:3, i) = x(1:3, i-1) + x(4:6, i-1)*dt + 0.5*a*dt^2;
    x(4:6, i) = x(4:6, i-1) + 0.5*(a + (f - c.*(x(4:6, i-1)) - k*x(1:3, i))./m)*dt;
end

% 绘制时间位移曲线
figure;
plot(t, x(1,:), t, x(2,:), t, x(3,:));
xlabel('时间（s）'); ylabel('位移（m）');
legend('质点1', '质点2', '质点3');

希望这个代码可以帮助到您！如果您还有其他问题，请随时问我。