帮我编写一个外径25cm内径10cm，磨粒为1cm的砂轮表面磨粒分布的matlab代码

好的，以下是一份简单的 Matlab 代码，可以生成一个外径为 25cm，内径为 10cm，磨粒为 1cm 的砂轮表面磨粒分布。其中，我们将砂轮表面分成若干个方格，每个方格内随机分布一定数量的磨粒，以模拟真实的磨粒分布情况。代码如下：

% 砂轮参数定义
outer_r = 0.25; % 外径
inner_r = 0.1; % 内径
grain_size = 0.01; % 磨粒大小
grid_size = 0.005; % 方格大小
num_grains = 100; % 每个方格内的磨粒数

% 计算砂轮表面上的方格数
grid_num = ceil((outer_r^2 - inner_r^2) / grid_size^2);

% 计算每个方格内的磨粒数
num_grains_per_grid = poissrnd(num_grains, [grid_num, 1]);

% 计算每个方格内磨粒的位置
offset = (outer_r + inner_r) / 2;
x = repmat(-offset:grid_size:offset, grid_num, 1);
y = repmat((-offset:grid_size:offset)', 1, grid_num);
x = x(:) + grid_size * rand(grid_num^2, num_grains);
y = y(:) + grid_size * rand(grid_num^2, num_grains);

% 计算每个磨粒的坐标
r = sqrt(x.^2 + y.^2);
theta = atan2(y, x);
theta(theta < 0) = theta(theta < 0) + 2*pi;
x = r .* cos(theta);
y = r .* sin(theta);

% 绘制砂轮表面磨粒分布图像
scatter(x, y, 1, 'k', 'filled');
axis equal;
xlim([-outer_r, outer_r]);
ylim([-outer_r, outer_r]);

运行上述代码后，将生成一个随机分布的砂轮表面磨粒分布图像。您可以根据需要调整参数，如方格大小、每个方格内的磨粒数等，以获得不同的砂轮表面磨粒分布图像。