Matlab最小二乘法曲线拟合

1、最小二乘法在曲线拟合中比较普遍。拟合的模型主要有1 .直线型2 .多项式型3 .分数函数型4 .指数函数型5 .对数线性型6 .高斯函数型一般对于LS问题，通常利用反斜杠运算、fminsearch或优化工具箱提供的极小化函数求解。在Matlab中，曲线拟合工具箱也提供了曲线拟合的图形界面操作。在命令提示符后键入：cftool ,即可根据数据，选择适当的拟合模型。命令1 .假设要拟合的多项式是：y=a+b*x+c*xA2.首先建立设计矩阵 X：X=ones(size(x) x xA2;执行：para=Xypara 中包含了三个参数：para(1)=a;para(2)=b;para(3)=c;这种

2、方法对于系数是线性的模型也适应。2 .假设要才以合:y=a+b*exp(x)+cx*exp(xA2)设计矩阵X为X=ones(size(x) exp(x) x.*exp(x.A2);para=Xy3 .多重回归(乘积回归)设要拟合：y=a+b*x+c*t ,其中x和t是预测变量，y是响应变量。设计矩阵为X=ones(size(x) x t % 注意x,t大小相等！para=Xypolyfit 函数polyfit函数不需要输入设计矩阵，在参数估计中，polyfit会根据输入的数据生成设计矩阵。1 .假设要拟合的多项式是：y=a+b*x+c*xA2p=polyfit(x,y,2)然后可以使用pol

3、yval在t处预测：y_hat=polyval(p,t)polyfit函数可以给出置信区间。p S=polyfit(x,y,2) %S中包含了 标准差y_fit,delta = polyval(p,t,S) %按照拟合模型在 t处预测在每个 t 处的 95%CI 为：*delta, y_fit+*delta)2 .指数模型也适应假设要才以合:y = a+b*exp(x)+c*exp(x.2)p=polyfit(x,log(y),2)fminsearch 函数fminsearch是优化工具箱的极小化函数。LS问题的基本思想就是残差的平方和(一种范数，由此，LS产生了许多应用)最小，因此可以利用f

4、minsearch函数进行曲线拟合。假设要才以合:y = a+b*exp(x)+c*exp(x.2)首先建立函数，可以通过m文件或函数句柄建立：x=.1；y=;f=(p,x) p(1)+p(2)*exp(x)+p(3)*exp(x.2) % 注意向量化:p(1)=a;p(2)=b;p(3)=c;%可以根据需要选择是否优化参数%opt=options() p0=ones(3,1);% 初值para=fminsearch(p) (y-f(p,x),A2,p0) %可以输出Hessian矩阵res=y-f(para,x)%拟合残差曲线拟合工具箱提供了很多拟合函数，对大样本场合比较有效！非线性拟合nl

5、infit 函数clear all;x1=;x2=;x=x1 x2;y=;f=(p,x)*p(1)*(1-1/p(2)*(1-(1-x(:,1)A(1/p(2)Ap(2)A2.*(x(:,1)A(-l/p(2)-1)A(-p(2).*x(:,1)A(-1/p(2).*x(:,2);p0=8 ;opt=optimset(TolFun,1e-3,TolX,1e-3);%p R=nlinfit(x,y,f,p0,opt)例子例子例子例子例子例子例子例子例子例子例子例子例子例子例子例子直线型例子2.多项式型的一个例子1900-2000年的总人口情况的曲线拟合 clear all;close all;%

6、cftool提供了可视化的曲线拟合！t=1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 2000;y=;%t太大，以t的嘉作为基函数会导致设计矩阵尺度太差，列变量几乎线性相依。变换为 -1 1上 s=(t-1950)/50;%plot(s,y,ro);刎归线：y=a+bxmx=mean(s);my=mean(y);sx=std(s);sy=std(y);r=corr(s,y);b=r*sy/sx;a=my-b*mx;rline=a+b.*s;figure;subplot(3,2,1 2)plot(s,y,ro,s,rline,k);%titl

7、e( 多项式拟合);set(gca,XTick,s,XTickLabel,sprintf(%d|,t);%hold on;n=4;PreYear=2010 2015 2020;% 预测年份tPreYear=(PreYear-1950)/50;Y=zeros(length(t),n);res=zeros(size(Y);delta=zeros(size(Y);PrePo=zeros(length(PreYear),n);Predelta=zeros(size(PrePo);for i=1:nP S(i)=polyfit(s,y,i);Y(:,i) delta(:,i)=polyval(p,s,S

8、(i);%拟合的 YPrePo(:,i) Predelta(:,i)=polyval(p,tPreYear,S(i);%预测res(:,i)=y-Y(:,i);% 残差end% plot(s,Y);%2009a自动添加不同颜色% legend(data,regression line,1st poly,2nd poly,3rd poly,4th poly,2)% plot(tPreYear,PrePo,);% hold off% plot(Y,res,o);% 残差图r=corr(s,Y),A2 %RA2%以合误差估计CIYearAdd=t;PreYear;tYearAdd=s;tPreYea

9、r;CFtit=一阶拟合,二阶才合,三阶拟合,四阶拟合);for col=1:nsubplot(3,2,col+2);plot(s,y,ro,s,Y(:,col),g-);%原始数据和拟合数据legend(Original,Fitted,2);hold on;plot(s,Y(:,col)+2*delta(:,col),r:);%95% CIplot(s,Y(:,col)-2*delta(:,col),r:);plot(tPreYear,PrePo(:,col),);%预测值plot(tPreYear,PrePo(:,col)+2*Predelta(:,col);%预测 95% CIplot(tPreYear,PrePo(:,col)-2*Predelta(:,col);axis( 0 400);set(gca,XTick,tYearAdd,XTickLabel,sprintf(%d|,YearAdd);title(CFtitcol);hold off;endfigure;%残差图for col=1:nsubplot(2,2,col);plot(Y(:,i),res(:,i),o);end一个非线性的 应用例子（多元情况）在百度知道 中，要拟合 y=a*x1An1+b*x2An2+c*x3An3%4:只是作为应用，模型不一定正确！%x2=x3!2,;y=;x1=1x2