最小二乘法拟合C语言实现

a=0; for(i=0;i<n;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;i<n;i++)
a+=(*(x+i))*(*(y+i)); return(a); }
}
�
av=av/n; return(av); } //平方和 float spfh(int n,float *x)
{int i; float a; a=0; for(i=0;i<n;i++)
a+=(*(x+i))*(*(x+i)); return(a); } //和平方
float shpf(int n,float *x) {int i; float a,b;
{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)
//两数先相加，再相乘
float djc(int n,float *x,float *y)
{int i; float a=0,b=0; for(i=0;i<n;i++)
{a=a+*(x+i); b=b+*(y+i); } a=a*b/n; return(a); } //系数a
float xsa(int n,float *x,float *y)
} void main()
{ int n,i; float *x,*y,a,b;
printf("请输入将要输入的有效数值组数n的值");
scanf("%d",&n);
printf("n=%d\n",n);
x=(float*)calloc(n,sizeof(float));
if(x==NULL) {printf("内存分配失败"); exit(1); }
for(i=0;i<n;i++)scanf("%f",y+i);
//x[0]=0.1;x[1]=0.3;x[2]=0.4;x[3]=0.55;x[4]=0.7;x[5]=0.8;x[6]=0.95;
//y[0]=15;y[1]=18;y[2]=19;y[3]=21;y[4]=22.6;y[5]=23.8;y[6]=26; for(i=0;i<n;i+&#gt;
#include <stdio.h>
#include <stdlib.h>
#include<malloc.h>
float average(int n,float *x)
{int i; float av; av=0; for(i=0;i<n;i++) av+=*(x+i);
{ printf("x[%d]=%5.4f ",i,*(x+i));
printf("y[%d]=%5.4f\n",i,*(y+i));
} a=xsa(n,x,y); b=xsb(n,x,y,a);
printf("经最小二乘法拟合得到的一元线性方程为:\n");
printf("f(x)=%5.4fx+%5.4f\n",a,b);
y=(float*)calloc(n,sizeof(float));
if(y==NULL) {printf("内存分配失败"); exit(1);
} printf("请输入x的值");
for(i=0;i<n;i++)scanf("%f",x+i);
printf("请输入y的值,请注意与x的值一一对应");
{int i; float a;
a=0; for(i=0;i<n;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);