给我三个自变量的连续局部枚举法MATLAB代码

以下是三个自变量的连续局部枚举法MATLAB代码：

1. 三个自变量的连续局部枚举法（循环嵌套版本）

% 定义目标函数
fun = @(x, y, z) x^2 + y^2 + z^2;

% 定义搜索范围
x_range = [-5, 5];
y_range = [-5, 5];
z_range = [-5, 5];

% 定义搜索步长
dx = 0.1;
dy = 0.1;
dz = 0.1;

% 初始化最小值
min_val = Inf;
min_x = NaN;
min_y = NaN;
min_z = NaN;

% 开始搜索
for x = x_range(1):dx:x_range(2)
    for y = y_range(1):dy:y_range(2)
        for z = z_range(1):dz:z_range(2)
            % 计算目标函数值
            val = fun(x, y, z);
            % 更新最小值和最小值对应的自变量值
            if val < min_val
                min_val = val;
                min_x = x;
                min_y = y;
                min_z = z;
            end
        end
    end
end

% 输出结果
fprintf('最小值为 %.4f，最小值对应的自变量为 (%.4f, %.4f, %.4f)\n', min_val, min_x, min_y, min_z);

2. 三个自变量的连续局部枚举法（矢量化版本）

% 定义目标函数
fun = @(x, y, z) x.^2 + y.^2 + z.^2;

% 定义搜索范围
x_range = [-5, 5];
y_range = [-5, 5];
z_range = [-5, 5];

% 定义搜索步长
dx = 0.1;
dy = 0.1;
dz = 0.1;

% 生成自变量网格点
[x_grid, y_grid, z_grid] = meshgrid(x_range(1):dx:x_range(2), y_range(1):dy:y_range(2), z_range(1):dz:z_range(2));

% 计算目标函数值
val_grid = fun(x_grid, y_grid, z_grid);

% 找到最小值和最小值对应的自变量值
[min_val, idx] = min(val_grid(:));
[min_x_idx, min_y_idx, min_z_idx] = ind2sub(size(val_grid), idx);
min_x = x_grid(min_x_idx, min_y_idx, min_z_idx);
min_y = y_grid(min_x_idx, min_y_idx, min_z_idx);
min_z = z_grid(min_x_idx, min_y_idx, min_z_idx);

% 输出结果
fprintf('最小值为 %.4f，最小值对应的自变量为 (%.4f, %.4f, %.4f)\n', min_val, min_x, min_y, min_z);

3. 三个自变量的连续局部枚举法（向量化版本）

% 定义目标函数
fun = @(x, y, z) x.^2 + y.^2 + z.^2;

% 定义搜索范围
x_range = [-5, 5];
y_range = [-5, 5];
z_range = [-5, 5];

% 定义搜索步长
dx = 0.1;
dy = 0.1;
dz = 0.1;

% 生成自变量向量
x_vec = x_range(1):dx:x_range(2);
y_vec = y_range(1):dy:y_range(2);
z_vec = z_range(1):dz:z_range(2);

% 生成自变量矩阵
[x_mat, y_mat, z_mat] = meshgrid(x_vec, y_vec, z_vec);

% 计算目标函数值矩阵
val_mat = fun(x_mat, y_mat, z_mat);

% 找到最小值和最小值对应的自变量值
[min_val, idx] = min(val_mat(:));
[min_x_idx, min_y_idx, min_z_idx] = ind2sub(size(val_mat), idx);
min_x = x_vec(min_x_idx);
min_y = y_vec(min_y_idx);
min_z = z_vec(min_z_idx);

% 输出结果
fprintf('最小值为 %.4f，最小值对应的自变量为 (%.4f, %.4f, %.4f)\n', min_val, min_x, min_y, min_z);