一次函数最小二乘法的c语言代码,最小二乘法C语言和Matlab实现

最小二乘法的两种方法实现

包括Matlab和C语言的实现

最小二乘法的Matlab实现

① 一次函数使用polyfit(x,y,1) 　　②多项式函数使用 polyfit(x,y,n)，n为次数 　　拟合曲线 　　x=[0.5,1.0,1.5,2.0,2.5,3.0]， 　　y=[1.75,2.45,3.81,4.80,7.00,8.60]。 　　解：MATLAB程序如下： 　　x=[0.5,1.0,1.5,2.0,2.5,3.0]; 　　y=[1.75,2.45,3.81,4.80,7.00,8.60]; 　　p=polyfit(x,y,2) 　　x1=0.5:0.05:3.0; 　　y1=polyval(p,x1); 　　plot(x,y,'*r',x1,y1,'-b') 　　计算结果为： 　　p =0.5614 0.8287 1.1560 　　即所得多项式为y=0.5614x^2+0.08287x+1.15560 　　③非线性函数使用 lsqcurvefit(fun,x0,x,y)

///

下面是实现Gauss-Jordan法实矩阵求逆。

#include

#include

#include

int brinv(double a[], int n)

{ int *is,*js,i,j,k,l,u,v;

double d,p;

is=malloc(n*sizeof(int));

js=malloc(n*sizeof(int));

for (k=0; k<=n-1; k++)

{ d=0.0;

for (i=k; i<=n-1; i