最小二乘法C语言的实现.doc

传感器技术题目要求：使用最小二乘法求线性度的拟合直线方程程序：#include #include #include #include float average(int n,float *x) int i; float av; av=0; for(i=0;in;i+) av+=*(x+i); av=av/n; return(av); /平方和 float spfh(int n,float *x) int i; float a,b; a=0; for(i=0;in;i+) a+=(*(x+i)*(*(x+i); return(a); /和平方 float shpf(int n,float *x) int i; float a,b; a=0; for(i=0;in;i+) a=a+*(x+i); b=a*a/n; return(b); /两数先相乘，再相加 float dcj(int n,float *x,float *y) int i; float a; a=0; for(i=0;in;i+) a+=(*(x+i)*(*(y+i); return(a); /两数先相加，再相乘 float djc(int n,float *x,float *y) int i; float a=0,b=0; for(i=0;in;i+) a=a+*(x+i); b=b+*(y+i); a=a*b/n; return(a); /系数a float xsa(int n,float *x,float *y) float a,b,c,d,e; a=spfh(n,x); b=shpf(n,x); c=dcj(n,x,y); d=djc(n,x,y); e=(c-d)/(a-b); /printf(%f %f %f %f,a,b,c,d); return(e); float he(int n,float *y) int i; float a; a=0; for(i=0;in;i+) a=a+*(y+i); return(a); float xsb(int n,float *x,float *y,float a) float b,c,d; b=he(n,y); c=he(n,x); d=(b-a*c)/n; return(d); void main() int n,i; float *x,*y,a,b; printf(请输入将要输入的有效数值组数n的值:); scanf(%d,&n); x=(float*)calloc(n,sizeof(float); if(x=NULL) printf(内存分配失败); exit(1); y=(float*)calloc(n,sizeof(float); if(y=NULL) printf(内存分配失败); exit(1); printf(请输入x的值n); for(i=0;in;i+)scanf(%f,x+i); printf(请输入y的值,请注意与x的值一一对应:n); for(i=0;in;i+)scanf(%f,y+i); for(i=0;in;i+) printf(x%d=%3.2f ,i,*(x+i); printf(y%d=%3.2fn,i,*(y+i); a=xsa(n,x,y); b=xsb(n,x,y,a); printf(经最小二乘法拟合得到的一元线性方程为:n); printf(f(x)=%3.2fx+%3.2f