数组索引必须为正整数或逻辑值。 出错 untitled4 (第 31 行) a(i) = (1/(m + theta*Delta_t*c + theta^2*Delta_t^2*k))*(F - c*(1 - theta)*v(i-1) - k*(1 - theta)*u(i-1) - theta*c*v(i-2) - theta^2*k*u(i-2));

这个错误提示通常表示在 MATLAB 中使用了负数的索引值，或者使用了非逻辑的索引值。在这个具体的代码中，可能是因为在第三个时间步长（即 i=3）时，会使用到 v(i-2) 和 u(i-2) 的值，但是这些索引值为负数，因此会出现错误提示。

为了解决这个问题，可以在循环开始之前初始化前两个时间步长的加速度、速度和位移向量，以确保在计算 a(i) 时不会出现负数的索引值。

下面是修改后的 MATLAB 代码：

% 定义参数
Delta_t = 0.01; % 时间步长
T = 10; % 总计算时间
N = T/Delta_t; % 计算时间步数
m = 10; % 质量
c = 3; % 阻尼
k = 10; % 刚度
theta = 0.25; % Wilson-θ 法参数

% 初始化加速度、位移和速度向量
a = zeros(N, 1); % 加速度向量
v = zeros(N, 1); % 速度向量
u = zeros(N, 1); % 位移向量

% 初始位移和速度
u(1) = 0; % 初始位移
v(1) = 0; % 初始速度

% 计算初始加速度
a(1) = (1/m)*(-c*v(1) - k*u(1));

% 计算第二个时间步长的加速度
t = Delta_t;
F = 5*t;
a(2) = (1/(m + theta*Delta_t*c + theta^2*Delta_t^2*k))*(F - c*(1 - theta)*v(1) - k*(1 - theta)*u(1) - theta*c*v(1) - theta^2*k*u(1));
v(2) = v(1) + Delta_t*((1 - theta)*a(1) + theta*a(2));
u(2) = u(1) + Delta_t*v(1) + ((1 - 2*theta)/2)*Delta_t^2*a(1);

% 计算第三个时间步长开始的加速度、速度和位移向量
for i = 3:N
    t = i*Delta_t;
    F = 0;
    if t <= 1
        F = 5*t;
    elseif t <= 2
        F = 10 - 5*t;
    end
    a(i) = (1/(m + theta*Delta_t*c + theta^2*Delta_t^2*k))*(F - c*(1 - theta)*v(i-1) - k*(1 - theta)*u(i-1) - theta*c*v(i-2) - theta^2*k*u(i-2));
    v(i) = v(i-1) + Delta_t*((1 - theta)*a(i-1) + theta*a(i));
    u(i) = u(i-1) + Delta_t*v(i-1) + ((1 - 2*theta)/2)*Delta_t^2*a(i-1);
end

% 输出动力学响应结果
plot(0:Delta_t:T-Delta_t, u);
xlabel('Time (s)');
ylabel('Displacement (m)');
title('Wilson-θ Method');

这样，就可以避免使用负数的索引值，从而消除了这个错误提示。