Matlab双运动圆模拟：循环中初始条件加载问题排查

核心原因分析 #

你遇到的“索引超出数组元素数目”错误，本质是数组未预分配足够长度或者循环起始索引错误——比如当i=1时，x(i-1)指向x(0)，但Matlab数组索引从1开始，直接触发越界。

具体修复步骤 #

1. 提前预分配所有数组长度 #

先根据仿真总时长和时间步长计算总步数，给位置、速度、加速度数组分配对应长度的空间，同时设置初始值：

dt = 0.01; % 时间步长
T_total = 10; % 总仿真时间
N = T_total / dt; % 总步数

% 预分配两个圆的位置、速度、加速度数组
x1 = zeros(1, N); y1 = zeros(1, N);
v1x = zeros(1, N); v1y = zeros(1, N);
a1x = zeros(1, N); a1y = zeros(1, N);

x2 = zeros(1, N); y2 = zeros(1, N);
v2x = zeros(1, N); v2y = zeros(1, N);
a2x = zeros(1, N); a2y = zeros(1, N);

% 设置初始位置（i=1对应t=0时刻）
x1(1) = 1; y1(1) = 2; % 圆1初始坐标
x2(1) = 5; y2(1) = 3; % 圆2初始坐标
% 设置初始速度
v1x(1) = 0.5; v1y(1) = 0;
v2x(1) = -0.3; v2y(1) = 0.2;

2. 从i=2开始执行循环计算 #

第1个元素已经是初始状态，后续每一步都基于前一时刻（i-1）的数据计算，循环起始索引设为2：

for i = 2:N
    % 计算当前时刻的加速度（替换为你的受力/运动模型）
    a1x(i) = get_acc_x1(x1(i-1), y1(i-1), x2(i-1), y2(i-1));
    a1y(i) = get_acc_y1(x1(i-1), y1(i-1), x2(i-1), y2(i-1));
    a2x(i) = get_acc_x2(x1(i-1), y1(i-1), x2(i-1), y2(i-1));
    a2y(i) = get_acc_y2(x1(i-1), y1(i-1), x2(i-1), y2(i-1));
    
    % 更新速度：基于前一时刻的速度和加速度
    v1x(i) = v1x(i-1) + a1x(i-1)*dt;
    v1y(i) = v1y(i-1) + a1y(i-1)*dt;
    v2x(i) = v2x(i-1) + a2x(i-1)*dt;
    v2y(i) = v2y(i-1) + a2y(i-1)*dt;
    
    % 更新位置：基于前一时刻的位置和速度
    x1(i) = x1(i-1) + v1x(i-1)*dt;
    y1(i) = y1(i-1) + v1y(i-1)*dt;
    x2(i) = x2(i-1) + v2x(i-1)*dt;
    y2(i) = y2(i-1) + v2y(i-1)*dt;
end

3. 避坑注意事项 #

• 禁止在循环中动态扩展数组（比如x1 = [x1, new_val]），这种操作不仅效率低，还容易导致索引逻辑混乱。
• 确保所有数组长度一致（都是N），避免某类数组只初始化了单个元素，循环到i>1时触发越界。
• 如果加速度基于当前位置计算，速度更新可改为v1x(i) = v1x(i-1) + a1x(i)*dt，但要保证a1x(i)已提前计算完成。

测试不同初始位置 #

直接修改x1(1)、y1(1)、x2(1)、y2(1)的数值即可，比如：

% 测试另一组初始位置
x1(1) = 0; y1(1) = 0;
x2(1) = 10; y2(1) = 10;

重新运行仿真就能观察不同初始位置下的运动行为。

内容的提问来源于stack exchange，提问作者Ytwo Monk