最小二乘法Matlab自编函数实现及示例(共4页).doc

一、 最小二乘拟合原理x= x1 x2 xny= y1y2 yn求m次拟合a0 a1 amT=(ATA)-1ATy所以m次拟合曲线为 y=a0+a1x+a2x2+amxm二、 Matlab实现程序function p=funLSM(x,y,m)%x,y为序列长度相等的数据向量，m为拟合多项式次数format short;A=zeros(m+1,m+1);for i=0:m for j=0:m A(i+1,j+1)=sum(x.(i+j);endb(i+1)=sum(x.i.*y);enda=Ab;p=fliplr(a);三、作业题1：给出如下数据，使用最小二乘法球一次和二次拟合多项式（取小数点后3位）X 1.36 1.49 1.73 1.81 1.95 2.16 2.28 2.48Y 14.094 15.069 16.844 17.378 18.435 19.949 20.963 22.495