如何补全Matlab代码求解polyfit拟合二次函数的极值点

嘿，我来帮你搞定这个二次拟合函数的极值求解问题！首先咱们先理清楚逻辑：你用polyfit(t,y,2)得到的p是二次多项式的系数数组，对应函数形式是 y = p(1)*t² + p(2)*t + p(3)。二次函数的极值点就在导数为0的位置，咱们一步步来实现：

核心原理 #

1. 二次函数的导数为：dy/dt = 2*p(1)*t + p(2)
2. 令导数等于0，解出极值点的t值：t_extreme = -p(2)/(2*p(1))
3. 将t_extreme代入拟合函数，得到对应的极值y_extreme
4. 通过p(1)的正负判断极值类型：
   • 若p(1) > 0：该点是极小值点
   • 若p(1) < 0：该点是极大值点

补全后的完整代码 #

path3='D:\SelectedCycle2\001\nm-01\';
list3 = dir(path3);
fName3 = {list3.name};
[~,y3]=size(fName3);
sumvec=double([]);
t=[];
for fff_no=3:y3
    image=double(imread(char(strcat(path3,fName3(fff_no)))));
    sumvec=[sumvec,sum(image(:))];
    t=[t,fff_no-2];
end
y=sumvec;
p = polyfit(t,y,2);

% ------------------- 新增的极值求解代码 -------------------
% 提取二次多项式系数
a = p(1);
b = p(2);
c = p(3);

% 处理特殊情况：如果a接近0，说明拟合结果近似线性，无极值
if abs(a) < 1e-8
    disp('拟合结果近似线性，不存在极值点');
else
    % 计算极值点的t值
    t_extreme = -b/(2*a);
    % 计算对应的极值
    y_extreme = polyval(p, t_extreme);
    
    % 判断极值类型并输出结果
    if a > 0
        fprintf('极小值点：t=%.4f，对应的y值=%.4f\n', t_extreme, y_extreme);
    else
        fprintf('极大值点：t=%.4f，对应的y值=%.4f\n', t_extreme, y_extreme);
    end
end

补充说明 #

• 我加了一个小判断abs(a) < 1e-8，用来避免当拟合结果接近线性时（a几乎为0）出现除以0的错误，这是实际编程里很实用的小细节~
• 用polyval(p, t_extreme)来计算极值比手动代入a*t_extreme^2 + b*t_extreme + c更稳妥，尤其是当你后续调整拟合阶数时，这个函数能自动适配。

内容的提问来源于stack exchange，提问作者Sanjay Kumar Gupta Res Scholar