Matlab使用poly9拟合时出现条件数不佳警告的问题求助

问题描述 #

我有数据集FR，绘制置信水平时弹出警告：

Warning: Equation is badly conditioned. Remove repeated data points or try centering and scaling.

生成的图像显示拟合曲线波动剧烈，置信区间异常宽大。使用的Matlab代码如下：

%创建拟合模型与置信区间
 FitVec = fit(freq_tl,T121,'poly9')
 pRvecConf = predint(FitVec,freq_tl,0.95,'observation','off');
 
 %获取拟合值
 fitY=feval(FitVec,freq_tl);
 %将置信区间与sqrt(xVec(i)).^2 .*xVec'相乘
 ci=(fitY-pRvecConf(:,1));
 %获取加权置信区间
 Conf_new=[fitY-ci,fitY+ci];
 %绘图
 figure
 plot(FitVec,freq_tl,T121) 
 hold on
 plot(freq_tl,Conf_new,'m--')
 legend('Data','Fitted curve','Confidence','Location','se')
 xlabel('Length')
 ylabel('Strength')

警告原因 #

• 高阶多项式的数值不稳定性：你使用了9阶多项式（poly9），x的各次幂（x¹到x⁹）之间高度线性相关，导致拟合的系数矩阵条件数极高，求逆时数值误差被大幅放大。
• 数据未做归一化：如果freq_tl的数值范围较大，高阶项的数值会远超过低阶项，进一步加剧计算误差。
• 重复数据点：若freq_tl存在大量重复值，会让拟合矩阵出现线性相关的行，恶化数值稳定性。

解决方法 #

1. 降低多项式阶数（优先尝试） #

9阶多项式大概率存在过拟合问题，且数值稳定性极差。先换成低阶多项式（比如3阶或5阶）：

FitVec = fit(freq_tl,T121,'poly5'); % 替换poly9为更低阶数

观察拟合曲线和置信区间是否合理，警告通常会随之消失。

2. 启用中心化/标准化 #

直接在拟合时开启Matlab的归一化选项，这是解决条件数问题的标准方案：

FitVec = fit(freq_tl,T121,'poly9','Normalize','on');

该选项会自动将自变量缩放到均值为0、标准差为1的范围，消除不同阶次项的数值差异，大幅改善矩阵条件数。

3. 移除重复数据点 #

检查并清理重复的自变量值：

% 去重：保留唯一freq_tl，对应T121取均值
[unique_freq, idx] = unique(freq_tl);
unique_T121 = accumarray(idx, T121, [], @mean);
% 用去重后的数据拟合
FitVec = fit(unique_freq, unique_T121,'poly9');

4. 简化置信区间绘图代码 #

你当前的Conf_new计算是冗余的，predint返回的pRvecConf已经是置信区间的上下界，直接绘图即可：

plot(freq_tl, pRvecConf(:,1), 'm--', freq_tl, pRvecConf(:,2), 'm--');

避免额外计算引入不必要的误差。

内容的提问来源于stack exchange，提问作者Blob