Matlab随机种子问题：仅生成三类目标中的小型陨石坑

核心问题 #

• 数组索引覆盖：三类陨石坑共用x_crater、y_crater、t_crater数组，但大型/中型陨石坑生成后，小型陨石坑的循环从smallcrater=0开始重新填充数组，直接覆盖了前面的结果，最终只保留小型陨石坑的数据。
• 循环范围可能无效：如果nstep是奇数，nstep/2会是小数，MATLAB的for循环会自动取整，但会导致前后半段的划分不符合预期（比如nstep=101时，nstep/2=50.5，循环50.5:101实际从51开始，和预期的后半段范围有偏差）。

修正方案 #

1. 使用全局计数器避免数组覆盖 #

不要为每类陨石坑单独设置计数器，改用一个全局计数器统一管理数组索引，确保三类陨石坑的数据依次存入数组，不会互相覆盖。

2. 修正循环范围 #

用floor(nstep/2)明确划分前后半段，确保循环范围是整数区间，避免因nstep奇偶性导致的范围错误。

修正后的代码 #

% 初始化全局计数器和数组
crater_count = 0;
x_crater = [];
y_crater = [];
t_crater = [];

% medium craters - 后半段生成
half_step = floor(nstep/2);
for k = half_step+1:nstep
    Time_in_Myr = k*dt/1e6;
    med_c = rand;
    if med_c > 0.98
        crater_count = crater_count + 1;
        x_crater(crater_count) = dx*(rand*nx);
        y_crater(crater_count) = dy*(rand*ny);
        t_crater(crater_count) = k;
    end
end

% big craters - 前半段生成
for k = 1:half_step
    Time_in_Myr = k*dt/1e6;
    big_c = rand;
    if big_c > 0.991
        crater_count = crater_count + 1;
        x_crater(crater_count) = dx*(rand*nx);
        y_crater(crater_count) = dy*(rand*ny);
        t_crater(crater_count) = k;
    end
end

% small craters - 全程生成
for k = 1:nstep
    Time_in_Myr = k*dt/1e6;
    small_c = rand;
    if small_c > 0.97
        crater_count = crater_count + 1;
        x_crater(crater_count) = dx*(rand*nx);
        y_crater(crater_count) = dy*(rand*ny);
        t_crater(crater_count) = k;
    end
end

额外优化建议 #

• 预先估算陨石坑的最大可能数量，提前初始化数组（比如x_crater = zeros(1, estimated_max)），避免动态扩容带来的性能损耗。
• 新增crater_type数组，存入对应类型标识（如1=大型，2=中型，3=小型），方便后续区分处理各类陨石坑。

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