给我一个matlab代码，要求可以通过孔隙率，纤维直径来控制的随机纤维模型

在MATLAB中，我们可以创建一个简单的随机纤维模型，它基于给定的孔隙率（porosity）和纤维直径（fiber diameter）。以下是一个基本的示例代码，它使用了`random`函数生成随机坐标点，并模拟纤维填充空间的过程：

% 定义常量
domain_size = [100 100 100]; % 空间域大小
fiber_diameter = 5; % 单位长度下的纤维直径
max_fibers = 100; % 最大纤维数
porosity = 0.4; % 孔隙率

% 计算纤维间隔
fiber_spacing = fiber_diameter / sqrt(3); % 三维空间中，理想情况下的纤维间距为直径的根号三倍

% 创建一个空矩阵用于存储纤维位置
fiber_positions = zeros(max_fibers, 3);

% 随机生成纤维位置
num_fibers = floor(porosity * domain_size(1) * domain_size(2) * domain_size(3) / (pi * fiber_diameter^3)); % 约束在给定孔隙率内的纤维数量
for i = 1:num_fibers
    while true
        % 随机选择一个点
        random_point = rand([1, 3]) * domain_size;
        
        % 检查是否与已有纤维冲突
        if isin(random_point, fiber_positions)
            continue;
        end
        
        % 更新纤维位置矩阵
        fiber_positions(i, :) = random_point;
        break;
    end
end

% 打印纤维数量和孔隙率
disp(['Generated ' num2str(num_fibers) ' fibers for a porosity of ' num2str(porosity*100) '%']);

% 可视化纤维模型（假设你有绘图库如Image Processing Toolbox）
% figure;
% hold on;
% scatter3(fiber_positions(:,1), fiber_positions(:,2), fiber_positions(:,3), 'r', 'filled');
% axis equal;
% xlabel('X');
% ylabel('Y');
% zlabel('Z');
% title('Random Fiber Model');