Matlab数据拟合：曲线拟合、线性/非线性拟合、贝塞尔曲线、B样条曲线及迭代法、梯度法、牛顿...
Matlab数据拟合曲线拟合线性/非线性拟合贝塞尔曲线B样条曲线迭代法梯度法牛顿插值拉格朗日插值三次样条等。
某个深夜调试温度传感器时，我盯着采集到的13组离散数据点突然意识到——数据拟合本质上是在混乱中创造秩序的行为。就像用乐高积木搭建曲线桥，不同的拼接手法会呈现出截然不同的建筑美学。今天咱们来聊聊Matlab工具箱里的这些"积木块"，看看它们如何把离散点变成优雅的曲线。
先来点基本功。线性拟合就像用直尺画线，polyfit函数是得力助手：
x = 1:0.5:7; y = 2.3*x + randn(size(x))*0.8; % 带噪声的线性数据 p = polyfit(x,y,1); fit_line = polyval(p,x); plot(x,y,'o',x,fit_line,'-r')那个神秘数字'1'指定多项式阶数，改成2就变成抛物线拟合。但真实世界往往更复杂——某次处理卫星轨道数据时，指数型衰减模型让我卡壳三天：
modelfun = @(b,x) b(1)*exp(-b(2)*x) + b(3); beta0 = [1, 0.1, 0.2]; % 初始猜想值 beta = lsqcurvefit(modelfun,beta0,x,y);这里beta0的选择充满玄学，有次误把衰减系数初始值设为10，结果迭代直接飞向无穷大。后来发现用半衰期倒推初始值更靠谱，这大概就是理论与实战的差距。
说到插值，拉格朗日像是用橡皮泥捏曲线，每个点都精确穿过：
xi = linspace(min(x),max(x),100); yi_lagrange = interp1(x,y,xi,'linear'); % 其实这是分段线性 yi_spline = interp1(x,y,xi,'spline');但遇到陡峭变化时，高次多项式会像过山车般震荡。有次处理ECG信号时，三次样条插值成功避免了这种"龙格现象"，秘诀在于分段低次多项式+连续导数的约束。
工业设计领域更偏爱贝塞尔曲线这种矢量玩家。虽然Matlab没有原生支持，但可以魔改绘图函数：
P = [0 0; 1 3; 4 2; 5 0]; % 控制点 t = linspace(0,1); Bezier = (1-t).^3*P(1,:) + 3*(1-t).^2.*t*P(2,:) + ... 3*(1-t).*t.^2*P(3,:) + t.^3*P(4,:);这种基于伯恩斯坦基函数的构造法，控制点就像磁铁般牵引曲线走向。而B样条更聪明地引入节点向量概念，局部修改特性让汽车曲面设计师们爱不释手：
knots = [0 0 0 0 1 2 2 2 2]; % 非均匀节点 t = linspace(0,2,50); coeffs = [1 3 5 2; 0 2 4 1]; curve = spcol(knots,3,t) * coeffs';注意那个魔法数字3表示三次样条，调整节点向量的稠密度能控制曲线刚度。有次模拟飞机机翼时，非均匀节点分布让前缘更尖锐而后缘平滑，这比均匀节点方案节省了60%的控制点。
当遇到超大规模数据拟合时，梯度下降法如同愚公移山。有次处理百万级气象数据，我这样构建迭代过程：
learning_rate = 0.001; for iter = 1:1000 grad = compute_gradient(x,y,theta); % 自定义梯度计算 theta = theta - learning_rate*grad; if norm(grad) < 1e-6 break; end end这里的学习率设定充满血泪史——设大了参数在山谷间跳跃，设小了迭代三天三夜。后来加入动量项后，收敛速度提升就像给老自行车装上电动马达。
从线性到非线性，从插值到逼近，每种工具都暗含着设计者的哲学。就像木匠选择不同刨刀，数据工匠也要懂得在过拟合与欠拟合的钢丝上保持平衡。下次当你的拟合曲线开始跳舞时，不妨换个基函数试试，说不定就能找到那个让数据安睡的摇篮曲。