晶粒生长的matlab代码

您好！对于晶粒生长的 Matlab 代码，我可以给您提供一个简单的例子，供参考：

% 设定初始条件
gridSize = 100; % 网格大小
numGrains = 20; % 晶粒数量
numIterations = 100; % 迭代次数

% 创建初始晶粒分布
grainMap = zeros(gridSize);
for i = 1:numGrains
    x = randi([1, gridSize]);
    y = randi([1, gridSize]);
    grainMap(x, y) = i;
end

% 迭代进行晶粒生长
for iter = 1:numIterations
    for x = 1:gridSize
        for y = 1:gridSize
            if grainMap(x, y) == 0 % 如果当前格点为空
                neighbors = [];
                % 寻找当前格点的邻居晶粒
                if x > 1 && grainMap(x-1, y) ~= 0
                    neighbors = [neighbors, grainMap(x-1, y)];
                end
                if x < gridSize && grainMap(x+1, y) ~= 0
                    neighbors = [neighbors, grainMap(x+1, y)];
                end
                if y > 1 && grainMap(x, y-1) ~= 0
                    neighbors = [neighbors, grainMap(x, y-1)];
                end
                if y < gridSize && grainMap(x, y+1) ~= 0
                    neighbors = [neighbors, grainMap(x, y+1)];
                end
                
                % 如果有邻居晶粒，则在当前格点随机选择一个邻居晶粒进行生长
                if ~isempty(neighbors)
                    grainMap(x, y) = neighbors(randi([1, length(neighbors)]));
                end
            end
        end
    end
end

% 可视化晶粒分布
figure;
imagesc(grainMap);
colormap(jet(numGrains));
colorbar;

这是一个简单的二维晶粒生长模拟代码，首先随机生成一些初始晶粒，然后通过迭代生长的方式，晶粒会扩展到空白区域。请根据实际需求进行修改和扩展。