就是一件小小的源码

				% 最小二乘法多项式拟合 
				function p = nihe(x,y,m)    
				% 用途：多项式拟合
				% 格式：p=nihe(x,y,m),x,y是数据向量，m为拟合次数，p为返回多项式系数降幂排列
				A=zeros(m+1,m+1);
				for i =0 : m
				    for j =0 : m
				        A(i+1,j+1)=sum(x.^(i+j));
				    end
				    b(i+1)=sum(x.^i.*y);
				end
				a=A\b';
				p=fliplr(a');