加上上述修正的代码呈现一个完整的代码

以下是加上修正后的完整代码，它解决了将符号值转换为数值进行计算的问题：

clc;
clear all;

syms x1 x2;

% 定义多项式函数
f = 2*x1^3 + 4*x1*x2^3 - 10*x1*x2 + x2^2;

% 寻找极值点
df_x1 = diff(f, x1);
df_x2 = diff(f, x2);

system_of_equations = [df_x1 == 0, df_x2 == 0];
extrema_points = solve(system_of_equations, [x1, x2], 'Real', true);

if isempty(extrema_points)
    disp('没有找到实数解。');
else
    disp('extrema_points的结构：');
    struct_info = structfun(@class, extrema_points, 'UniformOutput', false);
    disp(struct_info);

    % 初始化数值结果矩阵
    f_values = zeros(size(extrema_points, 1), 1);

    for i = 1:size(extrema_points, 1)
        if isfield(extrema_points(i), 'x1') && isfield(extrema_points(i), 'x2')
            x1_val = double(extrema_points(i).x1);
            x2_val = double(extrema_points(i).x2);
            
            % 将x1_val和x2_val转换为符号形式，然后替换到f中
            symbolic_f = subs(f, [x1, x2], [x1_val, x2_val]);
            f_values(i) = double(subs(symbolic_f)); % 计算并存储数值结果
        else
            disp(['在第 ', num2str(i), ' 个元素处无法正确提取x1和x2的值，跳过。']);
            continue;
        end
    end

    % 找到最小值及其索引
    [min_value, min_index] = min(f_values);
    
    % 提取最小值对应的点
    min_point_x1 = double(extrema_points(min_index).x1);
    min_point_x2 = double(extrema_points(min_index).x2);
    min_point = [min_point_x1, min_point_x2];

    disp(['函数在(0,0)附近的最小值为: ', num2str(min_value)]);
    disp(['最小值对应的点为: ', num2str(min_point)]);
end

这个版本的代码应该能成功计算并输出函数的极小值及其位置。如果你在运行此代码时遇到其他问题，请随时提问。